Problem in Rebind 2D disks
I am working on 2D simulations, I have tried to rebind the disk ( probably for the first time) I have the following error:
/scratch/elena/Test2DAdiabatic/AdiaStep4M10alpha5HR$
/home/elena/fargOCA2020/build/utils/rebind_disk -i disk90000.h5 -g grid1000x760x1.h5 --interpolation linear -o disk90000.h5reb
terminate called after throwing an instance of 'std::domain_error'
what(): Approximation order must be < data length.
I guess the reason is that in the vertical direction we have 1 layer only, so the interpolation fails, may be we only have to avoid interpolating along the vertical direction, @alainm I do not know if we can add something like:
if(grid.dim() == 3){
before
myPhi.rebind(nextGrid, interpolators);
in vectorField.cpp
pobably is more comlicated than this,