diff options
Diffstat (limited to 'src/SourceData.F77')
| -rw-r--r-- | src/SourceData.F77 | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/SourceData.F77 b/src/SourceData.F77 index 17b3efb..1cacee5 100644 --- a/src/SourceData.F77 +++ b/src/SourceData.F77 @@ -71,7 +71,7 @@ c Set up all coefficients and initial guess for solution Mcoeff(i,j,k) = 0.0d0 - if (r(i,j,k) <= radius) then + if (r(i,j,k) .le. radius) then Ncoeff(i,j,k) = charge_factor else Ncoeff(i,j,k) = 0.0d0 @@ -144,12 +144,12 @@ c Set up last timestep ... assume (first order) time symmetry c Output exact solution if required - if (output_tmp==1) then + if (output_tmp .eq. 1) then do k=1, cctk_lsh(3) do j=1, cctk_lsh(2) do i=1, cctk_lsh(1) - if (r(i,j,k) >= radius) then + if (r(i,j,k) .ge. radius) then temp(i,j,k) = charge/r(i,j,k) else temp(i,j,k) = charge/(2.0d0*radius**3)* |