[Carpet] interpolator problem

[Yosef Zlochower]

Hi,


   I am having some difficulty interpolating analysis variables (
specifically psi4).
I give the variables global storage of 1 time level. I calculate
the analysis quantities with the schedule option
OPTIONS: GLOBAL loop-local
I then determine which level is the finest level containing
a given sphere of points (for a later spherical harmonic
decomposition).
I then call a routine which will interpolate these points with the
the options
OPTIONS: Global loop-level.
Inside the routine I check if the refinement level is the one
I found previously, and if it is, call CCTK_InterpGridArrays.
To complicate the issue a little, the gfs are actually calculated
and interpolated twice during the analysis step (because I use
two different tetrads).

This works for most spheres, but in the case of a sphere
with radius 5, the global interpolator seems to pass the interpolation
points to the wrong processor. The interpolation then fails at the
local level with errors like:

WARNING[L1,P0] (AEILocalInterp):
    CCTK_InterpLocalUniform():
         interpolation point is either outside the grid,
         or inside but too close to the grid boundary!
         (this may be caused by a global interpolation with
          driver::ghost_size too small)
         0-origin interpolation point number pt=4 of N_interp_points=16137
         interpolation point (x,y,z)=(0.294297,-0.00577926,4.99133)
         grid x_min(delta_x)x_max = -5.87305(0.125)1.12695
         grid y_min(delta_y)y_max = -5.87305(0.125)1.12695
         grid z_min(delta_z)z_max = -5.99805(0.125)1.00195

(note that the interpolation point is in the +z bitant while
this processor owns the -z bitant)

I have attached the relevant sections of the parfile and
schedule file.

Yosef
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: parfile
Url: /archives/developers/attachments/20060927/c3f93df5/attachment.txt 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: schedule.ccl
Url: /archives/developers/attachments/20060927/c3f93df5/attachment-0001.txt