[Commits] [svn:einsteintoolkit] NullSHRExtract/branches/tapir/src/ (Rev. 19)
bela at caltech.edu
bela at caltech.edu
Wed Nov 6 19:02:36 CST 2013
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)
More information about the Commits
mailing list