[Users] Issue with ReadInterpolate

[Luciano Combi]

Hi people,

I have used ReadInterpolate successfully for different applications in the past. I'm now trying to use it to interpolate onto a new grid that has 5 ghost zones i.e.:

CoordBase::boundary_size_x_lower     = 5
CoordBase::boundary_size_y_lower     = 5
CoordBase::boundary_size_z_lower     = 5
CoordBase::boundary_size_x_upper     = 5
CoordBase::boundary_size_y_upper     = 5
CoordBase::boundary_size_z_upper     = 5
Driver::ghost_size                   = 5
GRHydro::GRHydro_stencil             = 5
But I'm getting this error:

CCTK_Proc569.err <==
  while executing schedule bin ReadInterpolate_ReadData, routine ReadInterpolate::ReadInterpolate_Read
  in thorn ReadInterpolate, file /home/mvxcombi/ET_Zelmani/Cactus/arrangements/EinsteinInitialData/ReadInterpolate/src/read.c:559:
  -> Skipping integer variable 'SPACEMASK::space_mask it=335925 tl=0 rl=0 c=569'. Do not know how to interpolate integers.
ERROR from host bcn1814.usr.hlrn.de process 569
  while executing schedule bin ReadInterpolate_ReadData, routine ReadInterpolate::ReadInterpolate_Read
  in thorn ReadInterpolate, file /home/mvxcombi/ET_Zelmani/Cactus/arrangements/EinsteinInitialData/ReadInterpolate/src/interpolate.cc:130:
  -> There were 70227 points in 27 variables that could not be set.
cactus_etk_devel: /home/mvxcombi/ET_Zelmani/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const cGH*, int): Assertion `0' failed.
Rank 569 with PID 1675935 received signal 6
Writing backtrace to AIC_interp/backtrace.569.txt
The original grid has 3 ghost zones. I'm using the following parfile:
# What to interpolate
# -------------------
ReadInterpolate::files                        = "checkpoint.chkpt.it_335925"
ReadInterpolate::max_number_of_read_variables = 30
ReadInterpolate::only_these_datasets          = ":Y_e, :rho, :temperature, :vel, :Ax, :Ay, :Az, :psi6phi, :kxx, :kxy, :kxz, :kyy, :kyz, :kzz, :gxx, :gxy, :gxz, :gyy, :gyz, :gzz, :betax, :betay, :betaz, :alp"
ReadInterpolate::minimum_reflevel             = 0     # Default
ReadInterpolate::maximum_reflevel             = 1000  # Default

# Interpolator settings
# ---------------------
ReadInterpolate::epsilon                 = 1.e-12 #2  # Default
ReadInterpolate::interpolator_half_width = -1      # Default, use cctk_nghostzones
ReadInterpolate::interpolator_name       = "Lagrange polynomial interpolation"
ReadInterpolate::interpolator_pars       = "
    order = 3
    boundary_off_centering_tolerance = {1. 1. 1. 1. 1. 1.}
    boundary_extrapolation_tolerance = {1.e-12 1.e-12 1.e-12 1.e-12 1.e-12 1.e-12}
"

# Fixing the read-in data
# -----------------------
ReadInterpolate::shift_read_datasets_by[0] = 0.
ReadInterpolate::shift_read_datasets_by[1] = 0.
ReadInterpolate::shift_read_datasets_by[2] = 0.

ReadInterpolate::enforce_symmetries_after_reading = no  # Default
ReadInterpolate::fix_cell_centered_origins        = no  # Default

Any suggestions on how to solve this issue?

Thanks a lot!!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20250725/87444213/attachment-0001.htm>