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

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