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

15 Nov 2013

User: szilagyi
Date: 2013/11/15 02:09 PM
Modified:
 /branches/tapir/src/
  NullSHRE_WTSphHarm.F90
Log:
 add more diagnostic output (time-deriv of ell)
File Changes:
Directory: /branches/tapir/src/
===============================
File [modified]: NullSHRE_WTSphHarm.F90
Delta lines: +136 -103
===================================================================

--- branches/tapir/src/NullSHRE_WTSphHarm.F90	2013-11-15 20:03:03 UTC (rev 31)
+++ branches/tapir/src/NullSHRE_WTSphHarm.F90	2013-11-15 20:09:38 UTC (rev 32)
@@ -72,12 +72,12 @@
   CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: fx_s,fy_s,fz_s
CCTK_REAL, dimension(null_lsh(1), null_lsh(2)) :: &
-    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
+       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
+       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 :: &
@@ -88,7 +88,9 @@
   logical, parameter :: output_lapse = .false.
   logical, parameter :: output_shift = .false.
   logical, parameter :: output_dgij = .false.
-  logical, parameter :: output_ell = .false.
+  logical, parameter :: output_ell = .true.
+  logical, parameter :: output_w0parts = .false.
+  logical, parameter :: output_u1parts = .false.
CCTK_INT :: a, b
   angular_tensor_1_n(1)%d => tmp1(:,:,1)
@@ -420,19 +422,7 @@
      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)
@@ -458,6 +448,46 @@
      call NullDecomp_WriteCoefsForRealFunc('elldenom_wt', cctkGH, truncate, &
           null_lsh, zeta, elld_n%d, elld_s%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('ds_x_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dsa_n(1,4)%d, dsa_s(1,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('ds_y_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dsa_n(2,4)%d, dsa_s(2,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('ds_z_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dsa_n(3,4)%d, dsa_s(3,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('dn_x_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dna_n(1,4)%d, dna_s(1,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('dn_y_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dna_n(2,4)%d, dna_s(2,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('dn_z_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dna_n(3,4)%d, dna_s(3,4)%d, cctk_time)
+
+     call NullDecomp_WriteCoefsForRealFunc('dn_t_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, dna_n(4,4)%d, dna_s(4,4)%d, cctk_time)
+
+
+     call NullDecomp_WriteCoefsForRealFunc('delldenom_dt_wt', cctkGH, truncate, &
+          null_lsh, zeta, delld_n(4)%d, delld_s(4)%d, cctk_time)
+
+
   end if
! xb, rb, r_l, r_t
@@ -488,10 +518,10 @@
        null_lsh, zeta, etaup0_n(1,4)%d, etaup0_s(1,4)%d, cctk_time)
call NullDecomp_WriteCoefsForComplexFunc(&
-    'AffineInvMetric_rA_qA', cctkGH, truncate, null_lsh, 1, zeta, &
-     etaup0_n(1,2)%d*qa(2)%d + etaup0_n(1,3)%d*qa(3)%d,&
-     etaup0_s(1,2)%d*qa(2)%d + etaup0_s(1,3)%d*qa(3)%d,&
-     cctk_time)
+       'AffineInvMetric_rA_qA', cctkGH, truncate, null_lsh, 1, zeta, &
+       etaup0_n(1,2)%d*qa(2)%d + etaup0_n(1,3)%d*qa(3)%d,&
+       etaup0_s(1,2)%d*qa(2)%d + etaup0_s(1,3)%d*qa(3)%d,&
+       cctk_time)
! J, beta, U, Q, W
@@ -524,113 +554,116 @@
! bits of W
+  if(output_w0parts) then
+     w_1_n = dr0_n(1)%d * etaup0_n(1,1)%d / r0_n%d
+     w_2_n = - 2.d0 * dr0_n(4)%d / r0_n%d
+     w_3_n =  &
+          + 2.d0 * ( dr0_n(2)%d * etaup0_n(1,2)%d &
+          + dr0_n(3)%d * etaup0_n(1,3)%d  ) / r0_n%d
+     w_4_n = &
+          + ( dr0_n(2)%d ** 2 * etaup0_n(2,2)%d &
+          + 2.d0 * dr0_n(2)%d * dr0_n(3)%d * etaup0_n(2,3)%d &
+          + dr0_n(3)%d ** 2 * etaup0_n(3,3)%d &
+          ) / (r0_n%d*dr0_n(1)%d)
+     w_5_n = - 1.d0 / r0_n%d
-  w_1_n = dr0_n(1)%d * etaup0_n(1,1)%d / r0_n%d
-  w_2_n = - 2.d0 * dr0_n(4)%d / r0_n%d
-  w_3_n =  &
-    + 2.d0 * ( dr0_n(2)%d * etaup0_n(1,2)%d &
-             + dr0_n(3)%d * etaup0_n(1,3)%d  ) / r0_n%d
-  w_4_n = &
-    + ( dr0_n(2)%d ** 2 * etaup0_n(2,2)%d &
-    + 2.d0 * dr0_n(2)%d * dr0_n(3)%d * etaup0_n(2,3)%d &
-    + dr0_n(3)%d ** 2 * etaup0_n(3,3)%d &
-      ) / (r0_n%d*dr0_n(1)%d)
-  w_5_n = - 1.d0 / r0_n%d
+     w_1_s = dr0_s(1)%d * etaup0_s(1,1)%d / r0_s%d
+     w_2_s = - 2.d0 * dr0_s(4)%d / r0_s%d
+     w_3_s =  &
+          + 2.d0 * ( dr0_s(2)%d * etaup0_s(1,2)%d &
+          + dr0_s(3)%d * etaup0_s(1,3)%d  ) / r0_s%d
+     w_4_s = &
+          + ( dr0_s(2)%d ** 2 * etaup0_s(2,2)%d &
+          + 2.d0 * dr0_s(2)%d * dr0_s(3)%d * etaup0_s(2,3)%d &
+          + dr0_s(3)%d ** 2 * etaup0_s(3,3)%d &
+          ) / (r0_s%d*dr0_s(1)%d)
+     w_5_s = - 1.d0 / r0_s%d
-  w_1_s = dr0_s(1)%d * etaup0_s(1,1)%d / r0_s%d
-  w_2_s = - 2.d0 * dr0_s(4)%d / r0_s%d
-  w_3_s =  &
-    + 2.d0 * ( dr0_s(2)%d * etaup0_s(1,2)%d &
-             + dr0_s(3)%d * etaup0_s(1,3)%d  ) / r0_s%d
-  w_4_s = &
-    + ( dr0_s(2)%d ** 2 * etaup0_s(2,2)%d &
-    + 2.d0 * dr0_s(2)%d * dr0_s(3)%d * etaup0_s(2,3)%d &
-    + dr0_s(3)%d ** 2 * etaup0_s(3,3)%d &
-      ) / (r0_s%d*dr0_s(1)%d)
-  w_5_s = - 1.d0 / r0_s%d
+     call NullDecomp_WriteCoefsForRealFunc(&
+          'W_part_rl_etaup_ll_over_r', cctkGH, truncate, &
+          null_lsh, zeta, w_1_n, w_1_s, cctk_time)
-  call NullDecomp_WriteCoefsForRealFunc(&
-    'W_part_rl_etaup_ll_over_r', cctkGH, truncate, &
-       null_lsh, zeta, w_1_n, w_1_s, cctk_time)
+     call NullDecomp_WriteCoefsForRealFunc(&
+          'W_part_negative_2ru_over_r', cctkGH, truncate, &
+          null_lsh, zeta, w_2_n, w_2_s, cctk_time)
-  call NullDecomp_WriteCoefsForRealFunc(&
-    'W_part_negative_2ru_over_r', cctkGH, truncate, &
-       null_lsh, zeta, w_2_n, w_2_s, cctk_time)
+     call NullDecomp_WriteCoefsForRealFunc(&
+          'W_part_2rA_etaup_lA_over_r', cctkGH, truncate, &
+          null_lsh, zeta, w_3_n, w_3_s, cctk_time)
-  call NullDecomp_WriteCoefsForRealFunc(&
-    'W_part_2rA_etaup_lA_over_r', cctkGH, truncate, &
-       null_lsh, zeta, w_3_n, w_3_s, cctk_time)
+     call NullDecomp_WriteCoefsForRealFunc(&
+          'W_part_rA_rB_etaup_AB_over_rl_over_r', cctkGH, truncate, &
+          null_lsh, zeta, w_4_n, w_4_s, cctk_time)
-  call NullDecomp_WriteCoefsForRealFunc(&
-    'W_part_rA_rB_etaup_AB_over_rl_over_r', cctkGH, truncate, &
-       null_lsh, zeta, w_4_n, w_4_s, cctk_time)
+     call NullDecomp_WriteCoefsForRealFunc(&
+          'W_part_negative_one_over_r', cctkGH, truncate, &
+          null_lsh, zeta, w_5_n, w_5_s, cctk_time)
-  call NullDecomp_WriteCoefsForRealFunc(&
-    'W_part_negative_one_over_r', cctkGH, truncate, &
-       null_lsh, zeta, w_5_n, w_5_s, cctk_time)
+  end if
+  ! write (*,*) 'diff_n = ' , &
+  !   maxval(abs( -w_wt_n%d + w_1_n + w_2_n + w_3_n + w_4_n + w_5_n))
-! write (*,*) 'diff_n = ' , &
-!   maxval(abs( -w_wt_n%d + w_1_n + w_2_n + w_3_n + w_4_n + w_5_n))
-  
-! write (*,*) 'diff_s = ' , &
-!   maxval(abs( -w_wt_s%d + w_1_s + w_2_s + w_3_s + w_4_s + w_5_s))
-  
+  ! write (*,*) 'diff_s = ' , &
+  !   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}
+  if(output_w0parts) then
-    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
+     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)
+     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
+     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_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)
+     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
+     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)
+     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)
+     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)
+     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)
+  end if
-! 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))
-  
+  ! 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

    