diff options
| author | eschnett <eschnett@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2013-10-30 15:20:38 +0000 |
|---|---|---|
| committer | eschnett <eschnett@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2013-10-30 15:20:38 +0000 |
| commit | 445e2850879ea9da712bda3cd5aee7a9de4b294a (patch) | |
| tree | 75a8a0e7e0e21e32e77479a87074f989f0813f83 /src/setupbrilldata3D.F | |
| parent | 21f2f0f02ce20a45b0033aa6c3ec424dc74b9217 (diff) | |
Make code compile when CCTK_REAL is real*4svn
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@129 a678b1cf-93e1-4b43-a69d-d43939e66649
Diffstat (limited to 'src/setupbrilldata3D.F')
| -rw-r--r-- | src/setupbrilldata3D.F | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/setupbrilldata3D.F b/src/setupbrilldata3D.F index b37c28c..8101cce 100644 --- a/src/setupbrilldata3D.F +++ b/src/setupbrilldata3D.F @@ -100,16 +100,16 @@ c is not regular on the axis. gzz(i,j,k) = e2q gxy(i,j,k) = - (one - e2q)*x1*y1/rho2 - brillMlinear(i,j,k) = 0.25D0/e2q + brillMlinear(i,j,k) = 0.25/e2q . *(brillq(rho1,z1+eps,phi) . + brillq(rho1,z1-eps,phi) . + brillq(rho1+eps,z1,phi) . + brillq(rho1-eps,z1,phi) - . - 4.0D0*brillq(rho1,z1,phi)) + . - 4.0*brillq(rho1,z1,phi)) . / eps**2 - brillMlinear(i,j,k) = brillMlinear(i,j,k) + 0.25D0/rho2 - . *(three*0.25D0*(brillq(rho1,z1,phi+eps) + brillMlinear(i,j,k) = brillMlinear(i,j,k) + 0.25/rho2 + . *(three*0.25*(brillq(rho1,z1,phi+eps) . - brillq(rho1,z1,phi-eps))**2 . + two*(brillq(rho1,z1,phi+eps) . - two*brillq(rho1,z1,phi) @@ -122,7 +122,7 @@ c is not regular on the axis. gzz(i,j,k) = one gxy(i,j,k) = zero - brillMlinear(i,j,k) = 0.25D0/e2q + brillMlinear(i,j,k) = 0.25/e2q . *(brillq(rho1,z1+eps,phi) . + brillq(rho1,z1-eps,phi) . + two*brillq(eps,z1,phi) |