#2601: Add four more parameters for carpet

Reporter: Liwei Ji Status: new Milestone: Version: development version Type: enhancement Priority: major Component: Carpet

Comment (by Liwei Ji):

Hi @{557058:59e031ba-9bb5-4298-a472-7b99d0ae6f22} , here is a test case

``` ActiveThorns = " Boundary Carpet CarpetIOASCII CarpetLib # CarpetRegrid2 CartGrid3D CoordBase IOUtil InitBase MoL SymBase "

$ghost_size = 1 $boundary_size = 1 $num_levels = 1

IOASCII::out1D_vars = " grid::x "

IO::out_dir = $parfile IO::out_fileinfo = "none" IO::parfile_write = no

IOASCII::out1D_every = 1 IOASCII::out1D_d = no IOASCII::out1D_y = no IOASCII::out1D_z = no

Carpet::max_refinement_levels = 1 Carpet::use_buffer_zones = yes Carpet::prolongation_order_space = 1 Carpet::prolongation_order_time = 1

Carpet::init_fill_timelevels = yes

CartGrid3D::type = "CoordBase" CartGrid3D::domain = "full" CartGrid3D::avoid_origin = no

CoordBase::xmin = -1. CoordBase::xmax = +1 - 0.125

CoordBase::ymin = -0.125 CoordBase::ymax = +0.125 CoordBase::zmin = -0.125 CoordBase::zmax = +0.125 CoordBase::dx = 0.125 CoordBase::dy = 0.125 CoordBase::dz = 0.125

driver::ghost_size = $ghost_size Cactus::cctk_itlast = 0

Carpet::domain_from_CoordBase = yes #Carpet::split_interior_points_x = yes Carpet::processor_topology = "manual" Carpet::processor_topology_3d_x = 4 Carpet::processor_topology_3d_y = 1 Carpet::processor_topology_3d_z = 1

CoordBase::boundary_size_x_lower = $boundary_size CoordBase::boundary_size_y_lower = $boundary_size CoordBase::boundary_size_z_lower = $boundary_size CoordBase::boundary_size_x_upper = $boundary_size CoordBase::boundary_size_y_upper = $boundary_size CoordBase::boundary_size_z_upper = $boundary_size CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_shiftout_x_upper = 1 #CoordBase::boundary_shiftout_y_lower = 1 #CoordBase::boundary_shiftout_y_upper = 1 #CoordBase::boundary_shiftout_z_lower = 1 #CoordBase::boundary_shiftout_z_upper = 1

#CarpetRegrid2::regrid_every = 0 #CarpetRegrid2::verbose = yes

#CarpetRegrid2::num_centres = 1 #CarpetRegrid2::num_levels_1 = $num_levels #CarpetRegrid2::position_x_1 = 0.0 #CarpetRegrid2::radius_1[1] = 0.2 ```

I made the following changes:

1. remove `CarpetRegrid2`, it seems that the new pars are not consistent with it, 2. set `Carpet::max_refinement_levels = 1` the new pars only can work with no refinement case for now, 3. add `CoordBase::boundary_shiftout_x_lower = 1` and `CoordBase::boundary_shiftout_x_upper = 1`, for some reason only when I turn on these. The boundary ghost points are counted as global points. And these option are on in `SphericalNR`.

After these changes, when `Carpet::split_interior_points_x = yes`, I got

``` ipos = [0,0,0], nprocs_dir = [4,1,1], cub = [5,3,3], rub0 = [18,3,3] ipos = [1,0,0], nprocs_dir = [4,1,1], cub = [9,3,3], rub0 = [18,3,3] ipos = [2,0,0], nprocs_dir = [4,1,1], cub = [13,3,3], rub0 = [18,3,3] ipos = [3,0,0], nprocs_dir = [4,1,1], cub = [18,3,3], rub0 = [18,3,3] ```

when `Carpet::split_interior_points_x = no`, I got

``` ipos = [0,0,0], nprocs_dir = [4,1,1], cub = [5,3,3], rub0 = [18,3,3] ipos = [1,0,0], nprocs_dir = [4,1,1], cub = [10,3,3], rub0 = [18,3,3] ipos = [2,0,0], nprocs_dir = [4,1,1], cub = [14,3,3], rub0 = [18,3,3] ipos = [3,0,0], nprocs_dir = [4,1,1], cub = [18,3,3], rub0 = [18,3,3] ```

basally, the distribution of `(clb, cub)` change in `x` change from `[(0,5), (5,9), (9,13), (13,18)]` to `[(0,5), (5,10), (10,14), (14,18)]`

-- Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2601/add-four-more-para...