An EinsteinExact version of qlm-ks.par. Generates NaN's inside theHEADmaster

horizon, but those get's masked by the carpet mask and so does
not affect the norms. Passes the test when compared with the
data produced by qlm-ks.par (at least on my laptop).


git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/QuasiLocalMeasures/trunk@40 ef6f4158-a8ee-47d1-ba14-cb73256398e0

1 files changed, 243 insertions, 0 deletions

diff --git a/test/qlm-ks-EE.par b/test/qlm-ks-EE.par
new file mode 100644
index 0000000..2605ac5
--- /dev/null
+++ b/test/qlm-ks-EE.par
@@ -0,0 +1,243 @@
+Cactus::cctk_run_title = "Quasi-Local Measures for Kerr-Schild data"
+
+Cactus::cctk_full_warnings         = yes
+Cactus::highlight_warning_messages = no
+
+Cactus::cctk_itlast = 3
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir       = $parfile
+IO::out_fileinfo  = "axis labels"
+IO::parfile_write = "no"
+
+
+
+ActiveThorns = "AEILocalInterp"
+
+ActiveThorns = "Fortran"
+
+ActiveThorns = "GSL"
+
+ActiveThorns = "GenericFD"
+
+ActiveThorns = "LoopControl"
+
+ActiveThorns = "Slab"
+
+
+
+ActiveThorns = "InitBase"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce"
+
+Carpet::domain_from_coordbase = yes
+driver::ghost_size            = 3
+Carpet::init_fill_timelevels  = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 1
+NaNChecker::action_if_found = "just warn"
+NaNChecker::check_vars      = "
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+"
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry90 SymBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmin =  0.00
+CoordBase::ymin =  0.00
+CoordBase::zmin =  0.00
+CoordBase::xmax = +4.00
+CoordBase::ymax = +4.00
+CoordBase::zmax = +4.00
+CoordBase::dx   =  0.125 
+CoordBase::dy   =  0.125
+CoordBase::dz   =  0.125
+
+CoordBase::boundary_size_x_lower     = 3
+CoordBase::boundary_size_y_lower     = 3
+CoordBase::boundary_size_z_lower     = 3
+CoordBase::boundary_size_x_upper     = 3
+CoordBase::boundary_size_y_upper     = 3
+CoordBase::boundary_size_z_upper     = 3
+
+CoordBase::boundary_shiftout_x_lower = 1
+CoordBase::boundary_shiftout_y_lower = 1
+CoordBase::boundary_shiftout_z_lower = 1
+
+CartGrid3D::type = "coordbase"
+
+ReflectionSymmetry::reflection_z   = yes
+ReflectionSymmetry::avoid_origin_z = no
+
+
+
+ActiveThorns = "SphericalSurface"
+
+SphericalSurface::nsurfaces = 1
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+
+
+
+ActiveThorns = "CarpetMask"
+
+CarpetMask::excluded_surface       [0] = 0
+CarpetMask::excluded_surface_factor[0] = 1.0
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.40
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"
+
+
+
+ActiveThorns = "KerrSchild" 
+
+ADMBase::initial_data  = "KerrSchild"
+ADMBase::initial_lapse = "KerrSchild"
+ADMBase::initial_shift = "KerrSchild"
+ADMBase::initial_dtlapse = "KerrSchild"
+ADMBase::initial_dtshift = "KerrSchild"
+
+KerrSchild::M      = 1.0
+KerrSchild::a      = 0.8
+
+ADMBase::evolution_method       = "KerrSchild"
+#ADMBase::lapse_evolution_method = "static"
+#ADMBase::shift_evolution_method = "static"
+
+
+
+ActiveThorns = "AHFinderDirect"
+
+AHFinderDirect::find_every = 1
+
+AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation"
+AHFinderDirect::geometry_interpolator_pars = "order=4"
+AHFinderDirect::surface_interpolator_name  = "Lagrange polynomial interpolation"
+AHFinderDirect::surface_interpolator_pars  = "order=4"
+
+AHFinderDirect::output_h_every = 0
+
+AHFinderDirect::N_horizons = 1
+
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0
+AHFinderDirect::which_surface_to_store_info        [1] = 0
+AHFinderDirect::reset_horizon_after_not_finding    [1] = no
+
+
+
+ActiveThorns = "QuasiLocalMeasures"
+
+QuasiLocalMeasures::verbose              = yes
+QuasiLocalMeasures::interpolator         = "Lagrange polynomial interpolation"
+QuasiLocalMeasures::interpolator_options = "order=4"
+QuasiLocalMeasures::spatial_order        = 4
+QuasiLocalMeasures::surface_index[0]     = 0
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        QuasiLocalMeasures::qlm_spin[0]
+        QuasiLocalMeasures::qlm_radius[0]
+        QuasiLocalMeasures::qlm_mass[0]
+        QuasiLocalMeasures::qlm_3det[0]
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::one_file_per_group = yes
+
+IOScalar::outScalar_every      = 1
+IOScalar::outScalar_reductions = "average count maximum minimum norm1 norm2 norm_inf"
+IOScalar::outScalar_vars       = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        SphericalSurface::sf_radius
+        QuasiLocalMeasures::qlm_shapes
+        QuasiLocalMeasures::qlm_coordinates
+        QuasiLocalMeasures::qlm_tetrad_l
+        QuasiLocalMeasures::qlm_tetrad_n
+        QuasiLocalMeasures::qlm_ricci_scalars
+        QuasiLocalMeasures::qlm_twometric
+        QuasiLocalMeasures::qlm_killing_vector
+        QuasiLocalMeasures::qlm_killed_twometric
+        QuasiLocalMeasures::qlm_invariant_coordinates
+        QuasiLocalMeasures::qlm_3determinant
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::out0D_every = 1
+IOASCII::out0D_vars  = "
+        CarpetReduce::weight
+        SphericalSurface::sf_valid
+        SphericalSurface::sf_info
+        SphericalSurface::sf_origin
+        SphericalSurface::sf_coordinate_descriptors
+        QuasiLocalMeasures::qlm_state
+        QuasiLocalMeasures::qlm_grid_int
+        QuasiLocalMeasures::qlm_grid_real
+        QuasiLocalMeasures::qlm_scalars
+        QuasiLocalMeasures::qlm_multipole_moments
+"
+
+IOASCII::out1D_every = 1
+IOASCII::out1D_vars  = "
+        CarpetReduce::weight
+        SphericalSurface::sf_radius
+        QuasiLocalMeasures::qlm_shapes
+        QuasiLocalMeasures::qlm_coordinates
+        QuasiLocalMeasures::qlm_tetrad_l
+        QuasiLocalMeasures::qlm_tetrad_n
+        QuasiLocalMeasures::qlm_tetrad_m
+        QuasiLocalMeasures::qlm_newman_penrose
+        QuasiLocalMeasures::qlm_weyl_scalars
+        QuasiLocalMeasures::qlm_ricci_scalars
+        QuasiLocalMeasures::qlm_twometric
+        QuasiLocalMeasures::qlm_killing_vector
+        QuasiLocalMeasures::qlm_killed_twometric
+        QuasiLocalMeasures::qlm_invariant_coordinates
+        QuasiLocalMeasures::qlm_3determinant
+"