diff options
| author | jthorn <jthorn@0a4070d5-58f5-498f-b6c0-2693e757fa0f> | 2005-05-26 11:30:08 +0000 |
|---|---|---|
| committer | jthorn <jthorn@0a4070d5-58f5-498f-b6c0-2693e757fa0f> | 2005-05-26 11:30:08 +0000 |
| commit | d70d7eb0fcf42ee8ec7135bf54492d00021f6a40 (patch) | |
| tree | d110e7180748296974c9a970fc1afcd8ca1666dd | |
| parent | 5c4a2df0af4927748a8fa062a53f6dfb9bcf2fdf (diff) | |
s/1.5/1.5d0/ and similarly for a few other single-precision constants
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@71 0a4070d5-58f5-498f-b6c0-2693e757fa0f
| -rw-r--r-- | src/IDAxiBrillBH.F | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/src/IDAxiBrillBH.F b/src/IDAxiBrillBH.F index 9e15dab..b4784de 100644 --- a/src/IDAxiBrillBH.F +++ b/src/IDAxiBrillBH.F @@ -192,7 +192,7 @@ c Initialize some arrays deta = etamax/(neb-3) do j=1,nqb - qgrd(j) = (j-1.5)*dq + qgrd(j) = (j-1.5d0)*dq do i=1,neb etagrd(i) = (i-2)*deta #include "CactusEinstein/IDAxiBrillBH/src/bhbrill.x" @@ -201,18 +201,18 @@ c Initialize some arrays c Boundary conditions do j=1,nqb ce(2,j)=ce(2,j)+cw(2,j) - cw(2,j)=0.0 + cw(2,j)=0.0d0 cw(neb-1,j)=cw(neb-1,j)+ce(neb-1,j) cc(neb-1,j)=cc(neb-1,j)-deta*ce(neb-1,j) - ce(neb-1,j)=0.0 + ce(neb-1,j)=0.0d0 enddo do i=1,neb cc(i,2)=cc(i,2)+cs(i,2) - cs(i,2)=0.0 + cs(i,2)=0.0d0 cc(i,nqb-1)=cc(i,nqb-1)+cn(i,nqb-1) - cn(i,nqb-1)=0.0 + cn(i,nqb-1)=0.0d0 enddo c Do the solve @@ -522,12 +522,12 @@ c dxgxx = 1/2 \partial gxx / \partial x enddo c Extrinsic Curvature is identically zero - kxx = 0.0D0 - kxy = 0.0D0 - kxz = 0.0D0 - kyy = 0.0D0 - kyz = 0.0D0 - kzz = 0.0D0 + kxx = 0.0d0 + kxy = 0.0d0 + kxz = 0.0d0 + kyy = 0.0d0 + kyz = 0.0d0 + kzz = 0.0d0 if (debug .ge. 6) then print *, '### final results (again at this 3-D grid point) ###' |