User: reisswig Date: 2013/06/16 04:10 AM

Modified: /trunk/ param.ccl /trunk/src/ NullNews_Bondi.F90, NullNews_GetNews.F90

Log: Allow linearized inertial frame interpolation via a Cactus parameter. Fixed a bug in linearized strain computation.

File Changes:

Directory: /trunk/src/ ======================

File [modified]: NullNews_Bondi.F90 Delta lines: +1 -2 =================================================================== --- trunk/src/NullNews_Bondi.F90 2013-06-05 07:58:20 UTC (rev 15) +++ trunk/src/NullNews_Bondi.F90 2013-06-16 09:10:32 UTC (rev 16) @@ -67,7 +67,6 @@ CCTK_INT, dimension(5) :: out_array_type_codes CCTK_INT, dimension(10), save :: gf_indices

- logical, parameter :: linearized = .false.

CCTK_COMPLEX, dimension(:), allocatable, save ::& Iarr1_this, Iarr1_other,& @@ -198,7 +197,7 @@ ! write(message, *) "maxval Uyo = ", maxval(abs(Uyo)) ! CALL CCTK_INFO(trim(message))

- if (.NOT. linearized) then + if (linearized_inertial_frame .eq. 0) then Uyh = 0

do ip = 1, 2

File [modified]: NullNews_GetNews.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/NullNews_GetNews.F90 2013-06-05 07:58:20 UTC (rev 15) +++ trunk/src/NullNews_GetNews.F90 2013-06-16 09:10:32 UTC (rev 16) @@ -206,7 +206,7 @@ u0 = uBondi call NullInterp_d2 (etheth_u0(:,:,1), dcmplx(u0(:,:,1)), 0, 1, 1) call NullInterp_d2 (etheth_u0(:,:,2), dcmplx(u0(:,:,2)), 0, 1, 1) - call NullInterp_cinterp(cctkGH, tmp_cgfn, tmp_cgfs, etheth_u0(:,:,1), etheth_u0(:,:,2), 1) + call NullInterp_cinterp(cctkGH, tmp_cgfn, tmp_cgfs, etheth_u0(:,:,1), etheth_u0(:,:,2), 2) call NullNews_ResetInactive(null_lsh, etheth_u0) linStrain(:,:,1) = J_l(:,:,1) + tmp_cgfn linStrain(:,:,2) = J_l(:,:,2) + tmp_cgfs

Directory: /trunk/ ==================

File [modified]: param.ccl Delta lines: +4 -0 =================================================================== --- trunk/param.ccl 2013-06-05 07:58:20 UTC (rev 15) +++ trunk/param.ccl 2013-06-16 09:10:32 UTC (rev 16) @@ -109,4 +109,8 @@ { } no

+BOOLEAN linearized_inertial_frame "Interpolation to inertial angular coordinates is linearized" +{ +} no

+