19990804: gdbiint

Brett,

I'm trying to duplicate your problem here by interpolating the
AVN thinned grids to the GAREA and PROJ you list below and assuming
a 1x1 degree interpolation grid.

In gdcfil I create the file using:
 PROJ     = ced
 GRDAREA  = 29;-161;58;-90
 KXKY     = 30;72

First point I noticed when trying to look at the TMPK@800 grid
you mention below- I noticed that the AVN doesn't contain that level.

So I tried 850mb and get:
    TIME1             TIME2         LEVL1 LEVL2   VCORD PARM
990804/0000F000                       850          PRES TMPK        

 GAREA:    25;-165;61;-85                                  
 SCALE FACTOR : 10** 0
 OUTPUT:    TERM/              

 MINIMUM AND MAXIMUM VALUES   275.02   303.02


So In this case, it looks like things are fine.


So, the one piece of information I didn't have is what your kxky
value is for this interpolation domain...presumably you don't have
more than the maximum number of grid points allowed in the interpolation
array (should be ~97,000).

What OS are you running on (solaris, Dec Unix, etc.).

I didn't see any division by zero errors that you had previously mentioned,
but I was wondering if this may be specific to the OS/compiler settings used.

Any more information would help.

Thanks,

Steve


>From: Brett Newkirk <address@hidden>
>Organization: .
>Keywords: 199907192220.QAA02818

>Hi Steve:
>       I have been working with the GDBIINT program but I have a
>question. I create a file (i.e. my own file with PROJ=CED and a
>particular region of the Pacific/West Coast: 29;-161;58;-90). I am trying
>to interpolate the AVN thinned grids to this new grid. The interior points
>look fine, but the boundary points have RMISSD values along with some
>other large negative values in TMPK at 800mb values for example. Did you
>notice problems in the testing of this routine or is there something that
>I am missing in using the program? Can one use the program to go from one
>projection to another with varying resolutions as long as the navigation
>blocks are defined?
>       I am noticing this problem in converting the ETA grids to this
>grid as well. However, the NOGAPS grid (CED, global coverage) does not
>have this boundary problem. Any insight that you could give on this issue
>I would greatly appreciate.
>       Also, I spoke to David Ovens and he had a bicubic spline routine
>written (usable by GR_INTP, for example). Thanks again for your help---
>Brett
>
>Brett Newkirk    E-MAIL: address@hidden
>Office: ATG 424  Atmospheric Science/Geophysics Building
>Mailing Address: Department of Atmospheric Sciences, University of Washington
>                 Box 351640
>                 Seattle, WA 98195-1640
>Office Phone: (206) 685-2183
>