diff options
| author | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2005-09-23 22:45:02 +0000 |
|---|---|---|
| committer | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2005-09-23 22:45:02 +0000 |
| commit | e30b0d281e291c38396305d1c1eb8b374f1114d2 (patch) | |
| tree | 7b6376e506a87278d4af3cba9f0d92e90230fec6 | |
| parent | 09a50cc3030484bb052b8fe1301ed3dcd518237c (diff) | |
Replace pause and stop statements with CCTK_WARN (0, ...)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@236 e296648e-0e4f-0410-bd07-d597d9acff87
| -rw-r--r-- | src/include/slice_normal.inc | 2 | ||||
| -rw-r--r-- | src/metrics/boost_rotation_symmetric.F77 | 12 | ||||
| -rw-r--r-- | src/xyz_blended_boundary.F77 | 2 |
3 files changed, 8 insertions, 8 deletions
diff --git a/src/include/slice_normal.inc b/src/include/slice_normal.inc index b8bb3d6..e56b884 100644 --- a/src/include/slice_normal.inc +++ b/src/include/slice_normal.inc @@ -72,7 +72,7 @@ C Debugging test. write(6,*) m, n, gu(m,n) end do end do - STOP + call CCTK_WARN (0, "aborting") end if C And now n^A = - 1/N g^AB n_B. (Sign so that it is future-pointing.) do m=1,4 diff --git a/src/metrics/boost_rotation_symmetric.F77 b/src/metrics/boost_rotation_symmetric.F77 index 8ab6fa6..1346c3b 100644 --- a/src/metrics/boost_rotation_symmetric.F77 +++ b/src/metrics/boost_rotation_symmetric.F77 @@ -69,7 +69,9 @@ C Intermediate quantities. C Make sure we are not sitting on one of the two source wordlines, C given by x = y = 0, z = +/- sqrt(h^2 + t^2) - if (num / h**4 .le. numlim) stop 'too close to source wordline' + if (num / h**4 .le. numlim) then + call CCTK_WARN (0, "too close to source wordline") + end if div = 1.d0 / sqrt(num**3) f = d**2 * ((0.25d0 * (a + b)**2 - h**2)**2 @@ -114,7 +116,7 @@ C we only need to check that gdzz is positive. write(*,*) 'WARNING 3-metric not spacelike in boostrot at' write(*,*) 't =', t, 'z =', z write(*,*) 'x =', x, 'y =', y - pause + call CCTK_WARN (0, "aborting") end if C Calculate inverse metric. That is not too difficult as it is @@ -123,8 +125,7 @@ c in block-diagonal form. tmp = gdtt * gdzz - gdtz**2 if (tmp .eq. 0.d0) then - write(*,*) 'boostrot metric inverse failed in tz plane' - STOP + call CCTK_WARN (0, "boostrot metric inverse failed in tz plane") end if gutt = gdzz / tmp @@ -134,8 +135,7 @@ c in block-diagonal form. tmp = gdxx * gdyy - gdxy**2 if (tmp .eq. 0.d0) then - write(*,*) 'boostrot metric inverse failed in xy plane' - STOP + call CCTK_WARN (0, "boostrot metric inverse failed in xy plane") end if guxx = gdyy / tmp diff --git a/src/xyz_blended_boundary.F77 b/src/xyz_blended_boundary.F77 index 66bc924..d8357bb 100644 --- a/src/xyz_blended_boundary.F77 +++ b/src/xyz_blended_boundary.F77 @@ -190,7 +190,7 @@ c OK so 6 blending cases. If frac = 1 we get all exact if (ninterps .eq. 0 .or. ninterps .gt. 3) then print *,"NINTERPS error", ninterps - STOP + call CCTK_WARN (0, "aborting") endif sfrac = (xfrac + yfrac + zfrac) * oonints |