diff options
Diffstat (limited to 'src/cartesian_to_spherical.F')
-rw-r--r-- | src/cartesian_to_spherical.F | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/cartesian_to_spherical.F b/src/cartesian_to_spherical.F index d55f84f..f75aea4 100644 --- a/src/cartesian_to_spherical.F +++ b/src/cartesian_to_spherical.F @@ -56,7 +56,7 @@ c ------------------------------------------------------------------ DO it = 1, Nt - ct = COS(theta(it)) ; ct2 = ct*ct + ct = COS(theta(it)) ; ct2 = ct*ct st = SIN(theta(it)) ; st2 = st*st DO ip = 1, Np @@ -64,9 +64,9 @@ c ------------------------------------------------------------------ cp = COS(phi(ip)) ; cp2 = cp*cp sp = SIN(phi(ip)) ; sp2 = sp*sp - gxx = gxxs(it,ip) ; dgxx = dgxxs(it,ip) - gxy = gxys(it,ip) ; dgxy = dgxys(it,ip) - gxz = gxzs(it,ip) ; dgxz = dgxzs(it,ip) + gxx = gxxs(it,ip) ; dgxx = dgxxs(it,ip) + gxy = gxys(it,ip) ; dgxy = dgxys(it,ip) + gxz = gxzs(it,ip) ; dgxz = dgxzs(it,ip) gyy = gyys(it,ip) ; dgyy = dgyys(it,ip) gyz = gyzs(it,ip) ; dgyz = dgyzs(it,ip) gzz = gzzs(it,ip) ; dgzz = dgzzs(it,ip) @@ -99,7 +99,7 @@ c ------------------------------------------------------------------ & -two*cp*sp*dgxy+cp2*dgyy) END DO - + END DO |