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.