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

[bela at caltech.edu]

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