[Commits] [svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 39)

[bela at caltech.edu]

User: szilagyi
Date: 2013/11/18 06:36 PM

Modified:
 /branches/tapir/src/
  NullSHRE_WTSphHarm.F90, NullSHRE_modAnaCoord.F90

Log:
 add output for Jt_wt

File Changes:

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

File [modified]: NullSHRE_WTSphHarm.F90
Delta lines: +9 -3
===================================================================
--- branches/tapir/src/NullSHRE_WTSphHarm.F90	2013-11-18 22:58:35 UTC (rev 38)
+++ branches/tapir/src/NullSHRE_WTSphHarm.F90	2013-11-19 00:36:47 UTC (rev 39)
@@ -77,11 +77,9 @@
 
   CCTK_COMPLEX, dimension(null_lsh(1), null_lsh(2)) :: &
        ul_1_n, ul_2_n, ul_3_n, ul_4_n, ul_5_n, &
-       ul_1_s, ul_2_s, ul_3_s, ul_4_s, ul_5_s
+       ul_1_s, ul_2_s, ul_3_s, ul_4_s, ul_5_s, jt_n, jt_s
 
   CCTK_REAL, dimension(null_lsh(1),null_lsh(2),4),target :: tmp1
-  CCTK_REAL, dimension(null_lsh(1),null_lsh(2),6),target :: &
-       tmp2, tmp3, tmp4, tmp5
   type (gf2d), dimension(2) :: angular_tensor_1_n, angular_tensor_1_s
 
   logical, parameter :: output_gij = .false.
@@ -546,6 +544,14 @@
      call NullDecomp_WriteCoefsForComplexFunc('J_wt', cctkGH, truncate, &
           null_lsh, 2, zeta, j_wt_n%d, j_wt_s%d, cctk_time)
 
+     if(cctk_iteration.gt.0) then
+       jt_n = (j_wt(:,:,1)-j_wt_p(:,:,1))/(null_time-null_time_p)
+       jt_s = (j_wt(:,:,2)-j_wt_p(:,:,2))/(null_time-null_time_p)
+
+       call NullDecomp_WriteCoefsForComplexFunc('Jt_wt', cctkGH, truncate, &
+            null_lsh, 2, zeta, jt_n, jt_s, 0.5*(null_time+null_time_p))
+     end if
+
      call NullDecomp_WriteCoefsForComplexFunc('Jl_wt', cctkGH, truncate, &
           null_lsh, 2, zeta, j_l_n%d, j_l_s%d, cctk_time)
 

File [modified]: NullSHRE_modAnaCoord.F90
Delta lines: +5 -0
===================================================================
--- branches/tapir/src/NullSHRE_modAnaCoord.F90	2013-11-18 22:58:35 UTC (rev 38)
+++ branches/tapir/src/NullSHRE_modAnaCoord.F90	2013-11-19 00:36:47 UTC (rev 39)
@@ -30,6 +30,11 @@
      dspherical_dstereo(2,2)%d =   dstereo_dspherical(1,1)%d / det
      dspherical_dstereo(2,1)%d = - dstereo_dspherical(2,1)%d / det
 
+!   InvF(1,1) =  F(2,2)/DetF;
+!   InvF(1,2) = -F(1,2)/DetF;
+!   InvF(2,2) =  F(1,1)/DetF;
+!   InvF(2,1) = -F(2,1)/DetF;
+
    end subroutine wt_dspherical_dstereo
 
    subroutine wt_dstereo_dspherical(ip, dstereo_dspherical)