[Commits] [svn:einsteintoolkit] NullEvolve/branches/tapir/src/ (Rev. 16)

[bela at caltech.edu]

User: szilagyi
Date: 2013/11/20 12:57 PM

Modified:
 /branches/tapir/src/
  NullEvol_Diag.F90

Log:
 add output of Jt

File Changes:

Directory: /branches/tapir/src/
===============================

File [modified]: NullEvol_Diag.F90
Delta lines: +13 -0
===================================================================
--- branches/tapir/src/NullEvol_Diag.F90	2013-10-22 18:26:19 UTC (rev 15)
+++ branches/tapir/src/NullEvol_Diag.F90	2013-11-20 18:57:18 UTC (rev 16)
@@ -27,6 +27,19 @@
 
   call NullEvol_DiagImArray(cctkGH, 'dxJ_ev', truncate, 2, dxjcn, dxjcs, diagtmp, null_xb, cctk_time)
 
+  if(.not. first_time) then
+
+    dxjcn = (jcn-jcn_p)/cctk_delta_time
+    dxjcs = (jcs-jcs_p)/cctk_delta_time
+
+    call NullEvol_DiagImArray(cctkGH, 'dtJ_ev', truncate, 2, &
+      dxjcn, dxjcs, diagtmp, null_xb, cctk_time-0.5*cctk_delta_time)
+
+    dxjcn = 1.e16
+    dxjcs = 1.e16
+
+  end if
+
   call NullEvol_DiagImArray(cctkGH, 'Q_ev', truncate, 1, qcn, qcs, diagtmp, null_xb, cctk_time)
 
   call NullEvol_DiagImArray(cctkGH, 'U_ev', truncate, 1, ucn, ucs, diagtmp, null_xbh, cctk_time)