diff options
| author | tradke <tradke@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 2006-09-27 16:19:13 +0000 |
|---|---|---|
| committer | tradke <tradke@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 2006-09-27 16:19:13 +0000 |
| commit | ba1f67d962fefa4ded627da879a841b94fbe1713 (patch) | |
| tree | 6590ce291011cbd7d69b706d0d8e5b9773e28111 | |
| parent | 0cb3a8838bbb1a9986866ebff44591a7ead621f1 (diff) | |
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
| -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)* |