Thanks Bary and Diane. When I am running the following gs file for regridding, it is<br> 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. I tried with options sequential, and the blank region (no data) is moving (spatially) with the time. I give the details here...<br><br>1. this is the regrid script....<br> open rain_1996-2000.ctl'<br> 'set gxout fwrite'
<br> 'set fwrite rain_2.5x2.5_2.dat'<br> 'set z 1'<br> 'set lat -28.75 43.75'<br> 'set lon 1.25 358.75'<br> 'set t 1 1826'<br> 'define p55=regrid2(rain,2.5,2.5,ba_p1,1.25,- 28.75)'<br> 'd p55'<br> 'disable fwrite'<br><br>
2. the rain_1996-2000.ctl file is the following...<br> dset ^rain_1996-2000.dat<br>undef 1.e35<br>xdef 1111 levels <br>0.1620178223<br> 0.4860534668<br> .........<br> 359.8379822<br>ydef 254 levels<br>-29.97691345<br>-29.69582748
<br>.........<br>45.02470398<br>zdef 1 linear 1 1<br>tdef 1826 linear 0z01jan1996 1dy<br>vars 1<br>rain 0 0 ACCUMULATED TOTAL PRECIPITATION<br>endvars<br><br>3. The rain_2.5x2.5.ctl file is theis.....
<br>set ^rain_2.5x2.5.dat<br>undef 1.e35<br>xdef 144 linear 1.25 2.5<br>ydef 30 linear -28.75 2.5<br>zdef 1 linear 1 1<br>tdef 1826 linear 0z01jan1996 1dy<br>vars 1<br>rain 1 99 ACCUMULATED TOTAL PRECIPITATION<br>
endvars<br><br>4. File sizes are following....<br>The rain_1996-2000.dat file size is 2061144976 (1111 x 254 x 1 x 1826 x 4)<br>The rain_2.5x2.5.dat file size is 31553280 (144 x 30 x 1 x 1826 x 4 ) <br>The file sizes are what I expect.
<br><br> Any help will b e great.<br><br> thanks.<br><br><br><br><br><br><br><br><div><span class="gmail_quote">On 10/17/06, <b class="gmail_sendername">Barry</b> <<a href="mailto:Barry.E.Schwartz@noaa.gov">Barry.E.Schwartz@noaa.gov
</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">John<br><br>For what it might be worth- I had similar problems<br>with a missing point and it turned out that since I wrote
<br>my binary data using FORTRAN, it creates an extra record<br>that C does not automatically understand. I added<br><br>OPTIONS SEQUENTIAL<br><br>to my control (ctl) file to account for it-<br><br>Barry<br><br><br>John Guhin wrote:
<br>> Hi,<br>> When I am running the following gs file for regridding, it is<br>> working fine, except that I am getting missing values at the four<br>> corners of the domain. I do not have any missing values in the
<br>> rain_1996-2000.dat file.<br>><br>> 'open rain_1996-2000.ctl'<br>> 'set gxout fwrite'<br>> 'set fwrite rain_2.5x2.5_2.dat'<br>> 'set z 1'<br>> 'set lat -28.75 43.75'<br>> 'set lon 1.25 358.75'
<br>> 'set t 1 1826'<br>> 'define p55=regrid2(rain,2.5,2.5,ba_p1,1.25,- 28.75)'<br>> 'd p55'<br>> 'disable fwrite'<br>><br>> John Guhin<br></blockquote></div><br><br clear="all"><br>-- <br>John Guhin