Hi all.
I'm trying something very simple -- running the "qc0-mclachlan" parameter file (specifically, the version found in Cactus/par after checking out EinsteinToolkit -- the version in Cactus/repos/McLachlan/par seems almost identical).
After getting rid of one thorn I can't check out anyway (LSUThorns/SummationByParts), the parfile seems to run OK. I've since modified it by switching to a full x-domain, removing the Pi-symmetry, and moving the punctures out a bit. That all seems to work.
But now I'm getting stuck trying to modify it in a simple way -- enlarging the grid while adding refinement levels "on the outside". Diffing my new "bad" parfile with the old "good" one yields the differences below. I'm attaching the bad parfile here also.
Can anyone tell me where I've gone wrong, please? Are there assumed parameter values never set in the test parfile?
Thanks,
Bernard
P.S. The errors all seem to be from NaNChecker. Here's the first printed sample:
WARNING level 1 in thorn NaNChecker processor 29 host borgh098.prv.cube (line 875 of /gpfsm/dnb31/bjkelly1/Cactus/configs/mcl/build/NaNChecker/NaNCheck.c): -> There were 35 NaN/Inf value(s) found in variable 'ML_BSSN::phi' (NaNmask bitfield 0)
---------- results of diffing "good" (<) and "bad" (>) parfiles -----------------
47c47 < Carpet::max_refinement_levels = 10 ---
Carpet::max_refinement_levels = 11
99,100c99,100 < CoordBase::xmin = -120.00 < CoordBase::ymin = -120.00 ---
CoordBase::xmin = -240.00 CoordBase::ymin = -240.00
102,107c102,107 < CoordBase::xmax = +120.00 < CoordBase::ymax = +120.00 < CoordBase::zmax = +120.00 < CoordBase::dx = 2.00 < CoordBase::dy = 2.00 < CoordBase::dz = 2.00 ---
CoordBase::xmax = +240.00 CoordBase::ymax = +240.00 CoordBase::zmax = +240.00 CoordBase::dx = 4.00 CoordBase::dy = 4.00 CoordBase::dz = 4.00
170c170 < CarpetRegrid2::num_levels_1 = 7 ---
CarpetRegrid2::num_levels_1 = 8
172,177c172,178 < CarpetRegrid2::radius_1[ 1] = 64.0 < CarpetRegrid2::radius_1[ 2] = 16.0 < CarpetRegrid2::radius_1[ 3] = 8.0 < CarpetRegrid2::radius_1[ 4] = 4.0 < CarpetRegrid2::radius_1[ 5] = 2.0 < CarpetRegrid2::radius_1[ 6] = 1.0 ---
CarpetRegrid2::radius_1[ 1] = 128.0 CarpetRegrid2::radius_1[ 2] = 64.0 CarpetRegrid2::radius_1[ 3] = 16.0 CarpetRegrid2::radius_1[ 4] = 8.0 CarpetRegrid2::radius_1[ 5] = 4.0 CarpetRegrid2::radius_1[ 6] = 2.0 CarpetRegrid2::radius_1[ 7] = 1.0
180c181 < CarpetRegrid2::num_levels_2 = 7 ---
CarpetRegrid2::num_levels_2 = 8
182,187c183,189 < CarpetRegrid2::radius_2[ 1] = 64.0 < CarpetRegrid2::radius_2[ 2] = 16.0 < CarpetRegrid2::radius_2[ 3] = 8.0 < CarpetRegrid2::radius_2[ 4] = 4.0 < CarpetRegrid2::radius_2[ 5] = 2.0 < CarpetRegrid2::radius_2[ 6] = 1.0 ---
CarpetRegrid2::radius_2[ 1] = 128.0 CarpetRegrid2::radius_2[ 2] = 64.0 CarpetRegrid2::radius_2[ 3] = 16.0 CarpetRegrid2::radius_2[ 4] = 8.0 CarpetRegrid2::radius_2[ 5] = 4.0 CarpetRegrid2::radius_2[ 6] = 2.0 CarpetRegrid2::radius_2[ 7] = 1.0
190c192 < CarpetRegrid2::num_levels_3 = 7 ---
CarpetRegrid2::num_levels_3 = 8
192,197c194,200 < CarpetRegrid2::radius_3[ 1] = 64.0 < CarpetRegrid2::radius_3[ 2] = 16.0 < CarpetRegrid2::radius_3[ 3] = 8.0 < CarpetRegrid2::radius_3[ 4] = 4.0 < CarpetRegrid2::radius_3[ 5] = 2.0 < CarpetRegrid2::radius_3[ 6] = 1.0 ---
CarpetRegrid2::radius_3[ 1] = 128.0 CarpetRegrid2::radius_3[ 2] = 64.0 CarpetRegrid2::radius_3[ 3] = 16.0 CarpetRegrid2::radius_3[ 4] = 8.0 CarpetRegrid2::radius_3[ 5] = 4.0 CarpetRegrid2::radius_3[ 6] = 2.0 CarpetRegrid2::radius_3[ 7] = 1.0
206c209 < Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256]" ---
Carpet::time_refinement_factors = "[1,1,1,2,4,8,16,32,64,128,256]"
----------------------------------------------------------------