User: szilagyi Date: 2013/11/15 02:03 PM
Modified: /branches/tapir/src/ NullSHRE_WTSphHarm.F90
Log: add more diagnostic output (bits of u1, w)
File Changes:
Directory: /branches/tapir/src/ ===============================
File [modified]: NullSHRE_WTSphHarm.F90 Delta lines: +71 -1 =================================================================== --- branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-14 20:11:15 UTC (rev 30) +++ branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-15 20:03:03 UTC (rev 31) @@ -75,6 +75,9 @@ w_1_n, w_2_n, w_3_n, w_4_n, w_5_n, & w_1_s, w_2_s, w_3_s, w_4_s, w_5_s
+ 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
CCTK_REAL, dimension(null_lsh(1),null_lsh(2),4),target :: tmp1 CCTK_REAL, dimension(null_lsh(1),null_lsh(2),6),target :: & @@ -87,6 +90,7 @@ logical, parameter :: output_dgij = .false. logical, parameter :: output_ell = .false.
+ CCTK_INT :: a, b angular_tensor_1_n(1)%d => tmp1(:,:,1) angular_tensor_1_n(2)%d => tmp1(:,:,2) angular_tensor_1_s(1)%d => tmp1(:,:,3) @@ -416,7 +420,19 @@ call NullDecomp_WriteCoefsForRealFunc('ell_t_wt', cctkGH, truncate, & null_lsh, zeta, ell_n(4)%d, ell_s(4)%d, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc('dell_x_dt_wt', cctkGH, truncate, & + null_lsh, zeta, j1_n(1,4)%d, j1_s(1,4)%d, cctk_time)
+ call NullDecomp_WriteCoefsForRealFunc('dell_y_dt_wt', cctkGH, truncate, & + null_lsh, zeta, j1_n(2,4)%d, j1_s(2,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dell_z_dt_wt', cctkGH, truncate, & + null_lsh, zeta, j1_n(3,4)%d, j1_s(3,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dell_t_dt_wt', cctkGH, truncate, & + null_lsh, zeta, j1_n(4,4)%d, j1_s(4,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('s_x_wt', cctkGH, truncate, & null_lsh, zeta, sa_n(1)%d, sa_s(1)%d, cctk_time)
@@ -492,7 +508,7 @@ null_lsh, 1, zeta, u_wt_n%d, u_wt_s%d, cctk_time)
call NullDecomp_WriteCoefsForComplexFunc('Ul_wt', cctkGH, truncate, & - null_lsh, 2, zeta, u_l_n%d, u_l_s%d, cctk_time) + null_lsh, 1, zeta, u_l_n%d, u_l_s%d, cctk_time)
call NullDecomp_WriteCoefsForRealFunc('B_wt', cctkGH, truncate, & null_lsh, zeta, beta_wt_n%d, beta_wt_s%d, cctk_time) @@ -561,6 +577,60 @@ ! maxval(abs( -w_wt_s%d + w_1_s + w_2_s + w_3_s + w_4_s + w_5_s))
+ ! bits of U_{,\lambda} + + + ul_1_s = - etaup1_s(1,2)%d * qa(2)%d - etaup1_s(1,3)%d * qa(3)%d + ul_1_n = - etaup1_n(1,2)%d * qa(2)%d - etaup1_n(1,3)%d * qa(3)%d + + call NullDecomp_WriteCoefsForComplexFunc(& + 'Ul_part_negative_detaup_dlambda_lA_qA', cctkGH, truncate, & + null_lsh, 1, zeta, ul_1_n, ul_1_s, cctk_time) + + + ul_2_s = 0 + ul_2_n = 0 + ul_3_s = 0 + ul_3_n = 0 + do a = 2, 3 + do b = 2, 3 + ul_2_s = ul_2_s - dr1_s(b)%d * etaup0_s(a,b)%d * qa(a)%d / dr0_s(1)%d + ul_2_n = ul_2_n - dr1_n(b)%d * etaup0_n(a,b)%d * qa(a)%d / dr0_n(1)%d + ul_3_s = ul_3_s - dr0_s(b)%d * etaup1_s(a,b)%d * qa(a)%d / dr0_s(1)%d + ul_3_n = ul_3_n - dr0_n(b)%d * etaup1_n(a,b)%d * qa(a)%d / dr0_n(1)%d + end do + end do + + call NullDecomp_WriteCoefsForComplexFunc(& + 'Ul_part_negative_etaup_AB_r1_B_qA_over_r1', cctkGH, truncate, & + null_lsh, 1, zeta, ul_2_n, ul_2_s, cctk_time) + + call NullDecomp_WriteCoefsForComplexFunc(& + 'Ul_part_negative_detaup_dlambda_AB_r0_B_qA_over_r1', cctkGH, truncate, & + null_lsh, 1, zeta, ul_3_n, ul_3_s, cctk_time) + + ul_4_s = 2.d0 * beta_l_s%d * u_wt_s%d + ul_4_n = 2.d0 * beta_l_n%d * u_wt_n%d + + call NullDecomp_WriteCoefsForComplexFunc(& + 'Ul_part_twice_beta_l_U', cctkGH, truncate, & + null_lsh, 1, zeta, ul_4_n, ul_4_s, cctk_time) + + ul_5_s = 2.d0 * beta_l_s%d * (etaup0_s(1,2)%d * qa(2)%d +& + & etaup0_s(1,3)%d * qa(3)%d) + ul_5_n = 2.d0 * beta_l_n%d * (etaup0_n(1,2)%d * qa(2)%d +& + & etaup0_n(1,3)%d * qa(3)%d) + + call NullDecomp_WriteCoefsForComplexFunc(& + 'Ul_part_twice_beta_l_etaup_lA_qA', cctkGH, truncate, & + null_lsh, 1, zeta, ul_5_n, ul_5_s, cctk_time) + + +! write (*,*) 'diff_ul = ' , & +! maxval(abs( -u_l_n%d + ul_1_n + ul_2_n + ul_3_n + ul_4_n + ul_5_n)), & +! maxval(abs( -u_l_s%d + ul_1_s + ul_2_s + ul_3_s + ul_4_s + ul_5_s)) + + end subroutine NullSHRE_WTSphHarm