undef_val = numpy.array([-1.0E10], dtype=numpy.float64)
# Create the expected results on the curvilinear grid
curv_expect_data = curv_data.copy()
curv_expect_data[curv_center_ignore] = undef_val
# Create the expected results on the rectilinear grid
rect_expect_data = rect_data.copy()
rect_expect_data[rect_center_ignore] = undef_val
# Initialize ESMP
if not ESMPControl().startCheckESMP():
raise RuntimeError("Unexpected failure to start ESMP")
# Create the regridder
regridder = CurvRect3DRegridder()
if direction in ('cw2r', 'r2cw'):
# Create the curvilinear grid with corner and center points
regridder.createCurvGrid(curv_center_lons, curv_center_lats,
curv_center_depths, curv_center_ignore, True,
curv_corner_lons, curv_corner_lats,
curv_corner_depths, curv_corner_ignore)
elif direction in ('co2r', 'r2co'):
# Create the curvilinear grid with only center points
regridder.createCurvGrid(curv_center_lons, curv_center_lats,
curv_center_depths, curv_center_ignore, True)
else:
raise ValueError("unexpected direction of %s" % direction)
# Create the rectilinear grid with corner and center points