diff options
| -rw-r--r-- | src/Extract.F | 30 |
1 files changed, 28 insertions, 2 deletions
diff --git a/src/Extract.F b/src/Extract.F index eddd0dc..b1c0f23 100644 --- a/src/Extract.F +++ b/src/Extract.F @@ -650,11 +650,11 @@ c Output gauge invariant variables & CHAR(im+48)//timestring c Write even parity waveforms - call Extract_write(cctkGH,openfile,out1,time,Qeven(:,il,im)) + call Extract_write2(cctkGH,openfile,out1,time,Qeven(:,il,im)) c Only write odd parity waveforms if full grid IF (igrid == 0) THEN - call Extract_write(cctkGH,openfile,out2,time,Qodd(:,il,im)) + call Extract_write2(cctkGH,openfile,out2,time,Qodd(:,il,im)) END IF END DO loop_m5 @@ -723,3 +723,29 @@ c Only write odd parity waveforms if full grid close(444) END SUBROUTINE Extract_write + + + SUBROUTINE Extract_write2(cctkGH, openfile,filename,value1,value2) + + implicit none + + DECLARE_CCTK_PARAMETERS + + CCTK_POINTER cctkGH + CCTK_REAL value1 + CCTK_REAL value2(2) + integer openfile + character*(*) filename + + if (openfile==1) then + OPEN(UNIT= 444,FILE=filename,STATUS="unknown") + call Extract_Advertise(cctkGH,filename) + else + OPEN(UNIT=444,FILE=filename,STATUS="old",POSITION="append") + end if + + write(444,*) value1,value2(1),value2(2) + + close(444) + + END SUBROUTINE Extract_write2 |