[Commits] [svn:einsteintoolkit] NullNews/branches/tapir/src/ (Rev. 23)

[bela at caltech.edu]

User: szilagyi
Date: 2014/01/07 04:36 PM

Modified:
 /branches/tapir/src/
  NullNews_SphHarmDecomp.F90

Log:
 work on diagnostic code

File Changes:

Directory: /branches/tapir/src/
===============================

File [modified]: NullNews_SphHarmDecomp.F90
Delta lines: +13 -82
===================================================================
--- branches/tapir/src/NullNews_SphHarmDecomp.F90	2014-01-07 22:36:05 UTC (rev 22)
+++ branches/tapir/src/NullNews_SphHarmDecomp.F90	2014-01-07 22:36:50 UTC (rev 23)
@@ -17,24 +17,10 @@
   implicit none
 
   CCTK_COMPLEX, dimension (:,:,:), allocatable:: cmplxtmp
-  ! CCTK_COMPLEX, dimension (:,:,:), allocatable:: NewsError
   CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: NewsSphCoeff
   CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: NewsBSphCoeff
   CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: Psi4SphCoeff
-  CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: JnSphCoeff
-  CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: Jn_lSphCoeff
-  CCTK_COMPLEX, dimension(1:lmax, -lmax:lmax) :: QnSphCoeff
-  CCTK_COMPLEX, dimension(1:lmax, -lmax:lmax) :: Qn_lSphCoeff
-  CCTK_COMPLEX, dimension(1:lmax, -lmax:lmax) :: UnSphCoeff
-  CCTK_COMPLEX, dimension(1:lmax, -lmax:lmax) :: Un_lSphCoeff
-  CCTK_COMPLEX, dimension(0:lmax, -lmax:lmax) :: omeganSphCoeff
-  CCTK_COMPLEX, dimension(0:lmax, -lmax:lmax) :: betanSphCoeff
-  CCTK_COMPLEX, dimension(0:lmax, -lmax:lmax) :: WnSphCoeff
-  CCTK_COMPLEX, dimension(0:lmax, -lmax:lmax) :: uBondiSphCoeff
   CCTK_COMPLEX, dimension(2:lmax, -lmax:lmax) :: linStrainSphCoeff
-  CCTK_REAL :: NewsNorm
-  CCTK_INT :: spin, lmode, nchar
-  character(len=256) :: io_out_dir, filestring_re, filestring_im, filestring, tmp
   CCTK_REAL :: time
 
   DECLARE_CCTK_FUNCTIONS
@@ -46,30 +32,6 @@
 
   truncate = (IO_TruncateOutputFiles(cctkGH) .ne. 0) .and. first_time
 
-!  write (*,*) "computing ylm coefs: time=", cctk_time, ", J=", maxval(abs(Jn))
-! fill the waveforms with Ylm at Scri before decomposition
-! if (Ylm_at_Scri == 1) then
-
-!    News = YlmScri_2
-!    NewsB = YlmScri_2
-!    Psi4 = YlmScri_2
-!    News_uBondi = YlmScri_2
-!    NewsB_uBondi = YlmScri_2
-!    Psi4_uBondi = YlmScri_2
-!    Jn = YlmScri_2
-!    Jn_l = YlmScri_2
-!    Qn = YlmScri_1
-!    Qn_l = YlmScri_1
-!    Un = YlmScri_1
-!    Un_l = YlmScri_1
-!    omegan = YlmScri_0
-!    deltan = YlmScri_0
-!    betan = YlmScri_0
-!    uBondi = YlmScri_0
-
-! end if
-
-
   ! allocate(NewsError(null_lsh(1), null_lsh(2), 2) )
 
   ! SpinDecompCoefs(GH, nq, np, s, z, f, sphcoef)
@@ -106,37 +68,6 @@
       endif
   endif
 
-  if (debug_output == 1) then
-      allocate(cmplxtmp(null_lsh(1), null_lsh(2), 2) )
-  
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 2, &
-         zeta, Jn, JnSphCoeff)
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 2, &
-         zeta, Jn_l, Jn_lSphCoeff)
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 1, &
-         zeta, Qn, QnSphCoeff)
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 1, &
-         zeta, Qn_l, Qn_lSphCoeff)
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 1, &
-         zeta, Un, UnSphCoeff)
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 1, &
-         zeta, Un_l, Un_lSphCoeff)
-      cmplxtmp = omegan
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 0, &
-         zeta, cmplxtmp, omeganSphCoeff)
-      cmplxtmp = betan
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 0, &
-         zeta, cmplxtmp, betanSphCoeff)
-      cmplxtmp = Wn
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 0, &
-         zeta, cmplxtmp, WnSphCoeff)
-      cmplxtmp = uBondi
-      call SpinDecompCoefs(cctkGH, null_lsh(1), null_lsh(2), 0, &
-         zeta, cmplxtmp, uBondiSphCoeff)
-      
-      deallocate(cmplxtmp)
-  endif
-
   ! deallocate(NewsError)
 
   ! Only output if this is processor 0
@@ -153,26 +84,26 @@
      call NullDecomp_WriteCoefFile('linStrain_scri', truncate, 2, lmax, linStrainSphCoeff, time)
   endif
   
-  if (debug_output == 1) then
+! if (debug_output == 1) then
 
-     time = cctk_time - 0.5 * cctk_delta_time
+!    time = cctk_time - 0.5 * cctk_delta_time
 
-     call NullDecomp_WriteCoefFile('J_scri', truncate, 2, lmax, JnSphCoeff, time)
-     call NullDecomp_WriteCoefFile('Jl_scri', truncate, 2, lmax, Jn_lSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('J_scri', truncate, 2, lmax, JnSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('Jl_scri', truncate, 2, lmax, Jn_lSphCoeff, time)
 
-     call NullDecomp_WriteCoefFile('U_scri', truncate, 1, lmax, UnSphCoeff, time)
-     call NullDecomp_WriteCoefFile('Ul_scri', truncate, 1, lmax, Un_lSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('U_scri', truncate, 1, lmax, UnSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('Ul_scri', truncate, 1, lmax, Un_lSphCoeff, time)
 
-     call NullDecomp_WriteCoefFile('Q_scri', truncate, 1, lmax, QnSphCoeff, time)
-     call NullDecomp_WriteCoefFile('Ql_scri', truncate, 1, lmax, Qn_lSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('Q_scri', truncate, 1, lmax, QnSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('Ql_scri', truncate, 1, lmax, Qn_lSphCoeff, time)
 
-     call NullDecomp_WriteCoefFile('B_scri', truncate, 0, lmax, betanSphCoeff, time)
-     call NullDecomp_WriteCoefFile('W_scri', truncate, 0, lmax, WnSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('B_scri', truncate, 0, lmax, betanSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('W_scri', truncate, 0, lmax, WnSphCoeff, time)
 
-     call NullDecomp_WriteCoefFile('omega_scri', truncate, 0, lmax, omeganSphCoeff, time)
-     call NullDecomp_WriteCoefFile('uBondi_scri', truncate, 0, lmax, uBondiSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('omega_scri', truncate, 0, lmax, omeganSphCoeff, time)
+!    call NullDecomp_WriteCoefFile('uBondi_scri', truncate, 0, lmax, uBondiSphCoeff, time)
 
-  endif
+! endif
 
 
   first_time = .FALSE.