[Users] adding levels to Carpet parameter file

Bernard Kelly bernard.j.kelly at nasa.gov
Mon Oct 11 12:52:07 CDT 2010


  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]"

----------------------------------------------------------------

-- 
-----------------------------------------------
Bernard J. Kelly

NASA Goddard Space Flight Center, Code 663
8800 Greenbelt Road
Greenbelt, MD 20771, U.S.A.

phone: +1 (301) 286-7243
-----------------------------------------------

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: qc0-mclachlan.par.6.FAILED
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20101011/64c1b3a0/attachment.pl 


More information about the Users mailing list