Added a VTK writer to simplify visualization with VisIt.

git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/QuasiLocalMeasures/trunk@7 ef6f4158-a8ee-47d1-ba14-cb73256398e0

1 files changed, 12 insertions, 0 deletions

diff --git a/src/qlm_calculate.F90 b/src/qlm_calculate.F90
index 19b0ffb..17b4177 100644
--- a/src/qlm_calculate.F90
+++ b/src/qlm_calculate.F90
@@ -18,6 +18,9 @@ subroutine qlm_calculate (CCTK_ARGUMENTS)
   integer   :: h0, hn
   
   character :: msg*1000
+  character :: slabel*2, ilabel*8
+  character(len=200) :: odir
+  integer   :: nchars
  
   logical   :: did_allocate
   
@@ -109,6 +112,15 @@ subroutine qlm_calculate (CCTK_ARGUMENTS)
         if (qlm_have_killing_vector(hn) /= 0) then
            call qlm_multipoles_normalise (CCTK_PASS_FTOF, hn)
         end if
+
+        if (output_vtk /= 0) then
+          if (mod(cctk_iteration,output_vtk_every) == 0) then
+            write(slabel,'(I2.2)'), hn
+            write(ilabel,'(I8.8)'), cctk_iteration
+            call CCTK_ParameterValString (nchars, "out_dir", "IOUtil", odir) 
+            call qlm_outputvtk (CCTK_PASS_FTOF, hn, odir(1:nchars)//'/surface'//slabel//'_'//ilabel//'.vtk', 1)
+          end if
+        end if
         
 9999    continue