Regrid problem : missing values at the corner after regridding

[Diane Stokes]

John,

I'm wondering about the start/end points of your input file and what
regrid2 does with box averaging if not all points in the new box (2.5
deg box in your case) are available.  Perhaps it needs a certain
percentage of the new box to be defined and your corner points don't
meet that criteria.

It's good that you included the script.  You might likely get a quicker
and more definite response if you also include your initial descriptor
file (rain_1996-2000.ctl), your new ctl file (the one for
rain_2.5x2.5_2.dat) and the file size (in bytes) of the output file from
fwrite (rain_2.5x2.5_2.dat).  Those might help someone see the problem
right off the bat.

    Diane


John Guhin wrote:
> Hi,
>      When I am running the following gs file for regridding, it is
> working fine, except that I am getting missing values at the four
> corners of the domain. I do not have any missing values in the
> rain_1996-2000.dat file.
>
> 'open rain_1996-2000.ctl'
> 'set gxout fwrite'
> 'set fwrite rain_2.5x2.5_2.dat'
> 'set z 1'
> 'set lat -28.75 43.75'
> 'set lon 1.25 358.75'
> 'set t 1 1826'
> 'define p55=regrid2(rain,2.5,2.5,ba_p1,1.25,- 28.75)'
> 'd p55'
> 'disable fwrite'
>
> John Guhin