diff options
| -rw-r--r-- | src/SourceData.F77 | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/SourceData.F77 b/src/SourceData.F77 index d5fc8da..17b3efb 100644 --- a/src/SourceData.F77 +++ b/src/SourceData.F77 @@ -60,7 +60,7 @@ c Get variable indices for all grid functions c Set up all coefficients and initial guess for solution - pi = 4.0*atan(1.0) + pi = 4.0d0*atan(1.0d0) charge_factor = 4.0d0*pi*charge*3.0d0/(4.0d0*pi*radius**3) do k=1, cctk_lsh(3) @@ -74,7 +74,7 @@ c Set up all coefficients and initial guess for solution if (r(i,j,k) <= radius) then Ncoeff(i,j,k) = charge_factor else - Ncoeff(i,j,k) = 0d0 + Ncoeff(i,j,k) = 0.0d0 end if end do @@ -95,9 +95,9 @@ c Set tolerance for stopping elliptic solve c Set any options needed for the elliptic solve - call Ell_SetStrKey(ierr, "yes","EllLinFlat::Bnd::Robin") - call Ell_SetRealKey(ierr, 0d0, "EllLinFlat::Bnd::Robin::inf") - call Ell_SetIntKey (ierr, 1, "EllLinFlat::Bnd::Robin::falloff") + call Ell_SetStrKey (ierr, "yes", "EllLinFlat::Bnd::Robin") + call Ell_SetRealKey(ierr, 0.0d0, "EllLinFlat::Bnd::Robin::inf") + call Ell_SetIntKey (ierr, 1, "EllLinFlat::Bnd::Robin::falloff") c Set parameters for specific solvers if (CCTK_EQUALS(solver,"sor")) then |