From 56918bd0ff10def10d6063aa9a66955e009f7380 Mon Sep 17 00:00:00 2001 From: allen Date: Mon, 4 Dec 2000 23:38:09 +0000 Subject: Added back output of imaginary parts which had been accidently removed git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@25 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129 --- src/Extract.F | 30 ++++++++++++++++++++++++++++-- 1 file 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 -- cgit v1.2.3