Regrid problem : missing values at the corner after regridding

[John Guhin]

Thanks Bary and Diane.  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.datfile. I tried with options sequential, and the blank
region (no data) is
moving (spatially) with the time.  I give the details here...

1. this is the regrid script....
 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'

2. the rain_1996-2000.ctl file is the following...
 dset ^rain_1996-2000.dat
undef 1.e35
xdef 1111 levels
0.1620178223
 0.4860534668
 .........
 359.8379822
ydef  254 levels
-29.97691345
-29.69582748
.........
45.02470398
zdef   1 linear 1 1
tdef    1826 linear  0z01jan1996  1dy
vars    1
rain      0  0 ACCUMULATED TOTAL PRECIPITATION
endvars

3.   The  rain_2.5x2.5.ctl   file is theis.....
set ^rain_2.5x2.5.dat
undef 1.e35
xdef 144 linear 1.25 2.5
ydef 30  linear -28.75 2.5
zdef 1   linear 1 1
tdef 1826 linear  0z01jan1996  1dy
vars    1
rain    1 99 ACCUMULATED TOTAL PRECIPITATION
endvars

4. File sizes are following....
The rain_1996-2000.dat file size is 2061144976 (1111 x 254 x 1 x 1826 x 4)
The rain_2.5x2.5.dat file size is 31553280 (144 x 30 x 1 x 1826 x 4 )
The file sizes are what I expect.

 Any help will b e great.

 thanks.







On 10/17/06, Barry <Barry.E.Schwartz at noaa.gov> wrote:
>
> John
>
> For what it might be worth- I had similar problems
> with a missing point and it turned out that since I wrote
> my binary data using FORTRAN, it creates an extra record
> that C does not automatically understand. I added
>
> OPTIONS SEQUENTIAL
>
> to my control (ctl) file to account for it-
>
> Barry
>
>
> 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
>



--
John Guhin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20061017/aca16a53/attachment.html