From ba1f67d962fefa4ded627da879a841b94fbe1713 Mon Sep 17 00:00:00 2001 From: tradke Date: Wed, 27 Sep 2006 16:19:13 +0000 Subject: Replaced '==' by '.eq.' and '<='/'>=' by '.le.'/'.ge.'. Amazing how long this syntax error went annoticed. git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveElliptic/trunk@72 41e88fdd-2190-4c69-9c84-4659c8cf322e --- src/SourceData.F77 | 6 +++--- 1 file 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)* -- cgit v1.2.3