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)