User: szilagyi Date: 2013/11/06 07:02 PM
Modified: /branches/tapir/src/ NullSHRE_WTSphHarm.F90
Log: added more diagnostic output
File Changes:
Directory: /branches/tapir/src/ ===============================
File [modified]: NullSHRE_WTSphHarm.F90 Delta lines: +271 -0 =================================================================== --- branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-10-31 21:25:34 UTC (rev 18) +++ branches/tapir/src/NullSHRE_WTSphHarm.F90 2013-11-07 01:02:36 UTC (rev 19) @@ -4,6 +4,21 @@ #include "cctk_Parameters.h" #include "cctk_Functions.h"
+subroutine NullSHRE_ConvertDerivs(null_lsh, fr, fq, fp, fx, fy, fz, j0inv) + use NullSHRE_modGFdef + implicit none + CCTK_INT, dimension(2), intent(in) :: null_lsh + CCTK_REAL, dimension(null_lsh(1), null_lsh(2)), intent(in) :: fr, fq, fp + CCTK_REAL, dimension(null_lsh(1), null_lsh(2)), intent(out) :: fx, fy, fz + + type (gf2d), dimension (4,4), intent(in) :: j0inv + + fx = j0inv(1,1)%d*fr + j0inv(2,1)%d*fq + j0inv(3,1)%d*fp + fy = j0inv(1,2)%d*fr + j0inv(2,2)%d*fq + j0inv(3,2)%d*fp + fz = j0inv(1,3)%d*fr + j0inv(2,3)%d*fq + j0inv(3,3)%d*fp + +end subroutine NullSHRE_ConvertDerivs + subroutine NullSHRE_WTSphHarm(CCTK_ARGUMENTS)
use NullInterp @@ -21,6 +36,9 @@ DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS
+ CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: fx_n,fy_n,fz_n + CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: fx_s,fy_s,fz_s + truncate = (IO_TruncateOutputFiles(cctkGH) .ne. 0) .and. first_time first_time = .FALSE.
@@ -65,7 +83,260 @@ call NullDecomp_WriteCoefsForRealFunc('alpha_wt', cctkGH, truncate, & null_lsh, zeta, alpha_n%d, alpha_s%d, cctk_time)
+ ! g_ab{,c}
+ call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(1,1,1)%d, dg_n(1,1,2)%d, dg_n(1,1,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(1,1,1)%d, dg_s(1,1,2)%d, dg_s(1,1,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_xx_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xx_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xx_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xx_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(1,1,4)%d, dg_s(1,1,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xx_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(1,1)%d, g1_s(1,1)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(1,2,1)%d, dg_n(1,2,2)%d, dg_n(1,2,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(1,2,1)%d, dg_s(1,2,2)%d, dg_s(1,2,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_xy_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xy_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xy_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xy_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(1,2,4)%d, dg_s(1,2,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xy_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(1,2)%d, g1_s(1,2)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(1,3,1)%d, dg_n(1,3,2)%d, dg_n(1,3,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(1,3,1)%d, dg_s(1,3,2)%d, dg_s(1,3,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_xz_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xz_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xz_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xz_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(1,3,4)%d, dg_s(1,3,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xz_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(1,3)%d, g1_s(1,3)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(1,4,1)%d, dg_n(1,4,2)%d, dg_n(1,4,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(1,4,1)%d, dg_s(1,4,2)%d, dg_s(1,4,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_xt_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xt_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xt_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xt_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(1,4,4)%d, dg_s(1,4,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_xt_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(1,4)%d, g1_s(1,4)%d, cctk_time) + + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(2,2,1)%d, dg_n(2,2,2)%d, dg_n(2,2,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(2,2,1)%d, dg_s(2,2,2)%d, dg_s(2,2,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_yy_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yy_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yy_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yy_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(2,2,4)%d, dg_s(2,2,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yy_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(2,2)%d, g1_s(2,2)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(2,3,1)%d, dg_n(2,3,2)%d, dg_n(2,3,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(2,3,1)%d, dg_s(2,3,2)%d, dg_s(2,3,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_yz_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yz_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yz_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yz_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(2,3,4)%d, dg_s(2,3,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yz_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(2,3)%d, g1_s(2,3)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(2,4,1)%d, dg_n(2,4,2)%d, dg_n(2,4,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(2,4,1)%d, dg_s(2,4,2)%d, dg_s(2,4,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_yt_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yt_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yt_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yt_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(2,4,4)%d, dg_s(2,4,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_yt_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(2,4)%d, g1_s(2,4)%d, cctk_time) + + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(3,3,1)%d, dg_n(3,3,2)%d, dg_n(3,3,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(3,3,1)%d, dg_s(3,3,2)%d, dg_s(3,3,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_zz_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zz_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zz_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zz_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(3,3,4)%d, dg_s(3,3,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zz_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(3,3)%d, g1_s(3,3)%d, cctk_time) + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(3,4,1)%d, dg_n(3,4,2)%d, dg_n(3,4,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(3,4,1)%d, dg_s(3,4,2)%d, dg_s(3,4,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_zt_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zt_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zt_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zt_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(3,4,4)%d, dg_s(3,4,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_zt_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(3,4)%d, g1_s(3,4)%d, cctk_time) + + + + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_n(4,4,1)%d, dg_n(4,4,2)%d, dg_n(4,4,3)%d, & + fx_n, fy_n, fz_n, j0inv_n) + + call NullSHRE_ConvertDerivs(null_lsh, & + dg_s(4,4,1)%d, dg_s(4,4,2)%d, dg_s(4,4,3)%d, & + fx_s, fy_s, fz_s, j0inv_s) + + call NullDecomp_WriteCoefsForRealFunc('dg_tt_x_wt', cctkGH, truncate, & + null_lsh, zeta, fx_n, fx_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_tt_y_wt', cctkGH, truncate, & + null_lsh, zeta, fy_n, fy_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_tt_z_wt', cctkGH, truncate, & + null_lsh, zeta, fz_n, fz_s, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_tt_t_wt', cctkGH, truncate, & + null_lsh, zeta, dg_n(4,4,4)%d, dg_s(4,4,4)%d, cctk_time) + + call NullDecomp_WriteCoefsForRealFunc('dg_tt_l_wt', cctkGH, truncate, & + null_lsh, zeta, g1_n(4,4)%d, g1_s(4,4)%d, cctk_time) + + ! null vector and parts of it call NullDecomp_WriteCoefsForRealFunc('ell_x_wt', cctkGH, truncate, & null_lsh, zeta, ell_n(1)%d, ell_s(1)%d, cctk_time)