[RegCNET] SST interpolation crossing date line
Raymond Arritt
rwarritt at gmail.com
Thu Oct 10 05:47:53 CEST 2013
Dear RegCNET friends,
I'm attempting to run RegCM-4.4rc8 for a large North America domain
using HadGEM2 initial and boundary conditions. When running the 'sst'
preprocessing step the code fails with the error message
NEST DOMAIN TOO NEAR TO PARENT.
I've traced this back to the part of mod_interp.f90 that searches for
the global data grid points that surround the RegCM4 grid points.
Through reading the code and some trial-and-error I'm almost certain
that this arises because the northwest corner of the domain crosses the
date line, and the code does not account for the wraparound of
longitudes at +/-180. That is, the code does not recognize that (for
example) RegCM4 grid longitude 179.888 is near the global data longitude
-180.
I have already done a large suite of current climate runs for the same
domain that were driven by ERA-Interim. These all worked fine, because
'sst' performs the interpolation for ERA-Interim without having to
search for surrounding points.
Before rewriting mod_interp.f90 to account for such cases I would like
to ask if anyone else already has done so, or if there is another fix
you could recommend.
Thanks for whatever help you can suggest.
Ray Arritt
More information about the RegCNET
mailing list