From 46baef91c3fe3a7107ba168bf1e6e0a5ded007f2 Mon Sep 17 00:00:00 2001
From: allen <allen@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7>
Date: Wed, 10 Nov 1999 09:20:19 +0000
Subject: Adding CCTK_INFOs and CCTK_WARNs instead of print statements and
 stops

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@41 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7
---
 README           |  2 --
 src/planewaves.F | 21 +++++++++++++++------
 src/teukwaves.F  | 42 +++++++++++++++++++++---------------------
 3 files changed, 36 insertions(+), 29 deletions(-)

diff --git a/README b/README
index 71f9143..2a27ef7 100644
--- a/README
+++ b/README
@@ -3,8 +3,6 @@ Authors: Malcolm Tobias, Joan Masso
 CVS info   : $Header$
 --------------------------------------------------------------------------
 
-1. Purpose of the thorn
-
 Provides analytic data for linear gravitational waves.
 
 
diff --git a/src/planewaves.F b/src/planewaves.F
index 0c8d858..95b119a 100644
--- a/src/planewaves.F
+++ b/src/planewaves.F
@@ -53,7 +53,8 @@ c Using macro definitions from Einstein
       CCTK_REAL wave,wavep
       CCTK_REAL kx,ky,kz,w
       CCTK_REAL kxp,kyp,kzp,wp
-      
+      CHARACTER*200 :: infoline
+
 c     local arrays (scalars! Man, what a sweat!)
       CCTK_REAL plus,minus,plusp,minusp,ain,aout
       
@@ -94,11 +95,19 @@ c     and convert it from degrees to radians
       the = pi*wavetheta/180.d0
       phi = pi*wavephi/180.d0
       
-      write(*,*)'Plane waves'
-      write(*,*)'amplitude  = ',amplitude
-      write(*,*)'wavecenter = ',ra
-      write(*,*)'wavelength = ',wave
-      write(*,*)'wavepulse  = ',wavep
+      call CCTK_INFO('Plane waves')
+      write(infoline,'(A13,G12.7)')
+     &    ' amplitude = ',amplitude
+      call CCTK_INFO(infoline)
+      write(infoline,'(A14,G12.7)')
+     &    ' wavecenter = ',ra
+      call CCTK_INFO(infoline)
+      write(infoline,'(A14,G12.7)')
+     &    ' wavelength = ',wave
+      call CCTK_INFO(infoline)
+      write(infoline,'(A14,G12.7)')
+     &    ' wavepulse = ',wavep
+      call CCTK_INFO(infoline)
 
 c     precalc
       kx = 2*pi*sin(the)*cos(phi)/wave
diff --git a/src/teukwaves.F b/src/teukwaves.F
index e5df827..24b8ab7 100644
--- a/src/teukwaves.F
+++ b/src/teukwaves.F
@@ -37,7 +37,6 @@ c Using macro definitions from Einstein
       DECLARE_CCTK_FARGUMENTS
       DECLARE_CCTK_PARAMETERS
       
-      
       CCTK_REAL amp,m,ra,pi
       CCTK_REAL wave,wave2,wave3,wave4,wave5,wave6,wave7,wave8
       INTEGER ipacket,iparity
@@ -62,8 +61,8 @@ c     uffff....
 c     from old spheretocart
       CCTK_REAL drx,dry,drz,dtx,dty,dtz,dpx,dpy,dpz
       
+      CHARACTER*200 infoline
       INTEGER i,j,k
-      
       INTEGER CCTK_Equals
       
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -90,27 +89,31 @@ c     by kappa, because they use that in the packet and
 c     4/3 because they use I not Teuk. F  
       if (CCTK_Equals(packet,'eppley') == 1) then
          ipacket = 1
-         write(*,*)'Teuk.Packet= eppley' 
+         call CCTK_INFO('Teukolsky Packet = eppley') 
       elseif (CCTK_Equals(packet,'evans') == 1) then
          ipacket = 2
-         write(*,*)'Teuk.Packet= evans'
+         call CCTK_INFO('Teukolsky Packet = evans')
 c     to be consistent with Evans & Abrahams, rescale the amplitude
 c     by kappa, because they use that in the packet and
 c     4/3 because they use I not Teuk. F
          amp = amp*(4.0d0/3.0d0)*kappa
       elseif (CCTK_Equals(packet,'square') == 1) then
          ipacket = 3
-         write(*,*)'Teuk.Packet= square'
-      else
-         write(*,*)'Teukwaves: unsupported packet'
-         STOP
-      endif      
-
+         call CCTK_INFO('Teukolsky Packet = square')
+      end if
       
-      write(*,*)'amplitude  = ',amp
-      write(*,*)'     m     = ',m
-      write(*,*)'wavecenter = ',ra
-      write(*,*)'wavelength = ',wave
+      write(infoline,'(A13,G12.7)')
+     &    ' amplitude = ',amp
+      call CCTK_INFO(infoline)
+      write(infoline,'(A11,G12.7)')
+     &    ' m value = ',m
+      call CCTK_INFO(infoline)
+      write(infoline,'(A14,G12.7)')
+     &    ' wavecenter = ',ra
+      call CCTK_INFO(infoline)
+      write(infoline,'(A14,G12.7)')
+     &    ' wavelength = ',wave
+      call CCTK_INFO(infoline)
       
       if (CCTK_Equals(parity,'even') == 1) then    
          iparity = 0
@@ -133,9 +136,7 @@ c     4/3 because they use I not Teuk. F
       
       if (ipacket.eq.1) then
          if (ingoing.eq.0.or.outgoing.eq.0) then
-            write(*,*)'WARNING: Epply packet is only non-singular at'
-            write(*,*)'the origin for ingoing-outgoing combination'
-            write(*,*)'of waves.  Expect problems ;-) '
+            call CCTK_WARN(4,'Epply packet is only non-singular at the origin for ingoing-outgoing combination of waves')
          endif
       endif
       
@@ -303,8 +304,8 @@ c     evans package  --> w^4(1-x^2/w^2)^6
 c     square package --> (1-x^2/w^2)^2
                if(ipacket.eq.3) then 
                   
-                  write(*,*)'WARNING: need to calculate fpe,fne'
-                  write(*,*)'WARNING: need conditionals for in out waves'
+                  call CCTK_WARN(4,'Need to calculate fpe,fne')
+                  call CCTK_WARN(4,'Need conditionals for in out waves')
                   
                   if (abs(teuk_tp).lt.wave) then
                      fp = amp * (1-(tp2/wave2))**2
@@ -437,8 +438,7 @@ c     mvalue
                   dtp = (2-sint**2)*2*sinp*cosp
                   dpp = 2*cost*(1-2*sinp**2)
                else
-                  write(*,*) 'teukwaves: m should be one of -2,-1,0,1,2'
-                  STOP
+                  call CCTK_WARN(0,'m should be one of -2,-1,0,1,2')
                endif
                
 c     parity
-- 
cgit v1.2.3