[Commits] [svn:einsteintoolkit] NullNews/trunk/ (Rev. 16)

[reisswig at tapir.caltech.edu]

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
 
+