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!!
_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users
