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)