<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000099">
<font color="#000000">Hi, Lynette.<br>
<br>
The default output for fwrite is "stream". Your descriptor file for
your new binary file includes "options sequential", so, grads is
skipping what it thinks are fortran control words at the beginning and
end of each record. Instead it's skipping valid data points resulting
in your shift.<br>
<br>
If you want your new file to be sequential, include the option "-sq" in
your call to "set fwrite". Otherwise, remove the sequential option
from your descriptor file.<br>
<br>
For what it's worth, you can also convert grib to binary using wgrib.
(<a class="moz-txt-link-freetext" href="http://www.cpc.ncep.noaa.gov/products/wesley/wgrib.html">http://www.cpc.ncep.noaa.gov/products/wesley/wgrib.html</a>).<br>
<br>
Good luck.<br>
Diane<br>
</font><br>
Lynette Bettio wrote:
<blockquote
cite="mid656E711DAB9AA74FBB22DE320668E99201385EA3@officeho2.bom.gov.au"
type="cite">
<meta http-equiv="Content-Type" content="text/html; ">
<meta content="MSHTML 6.00.2900.2912" name="GENERATOR">
<div dir="ltr" align="left"><font face="Arial" size="2"><span
class="357344306-19062006">Hi,</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">hopefully
I have made a simple mistake and this is an easy fix...</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><span class="357344306-19062006"><font face="Arial"><font
size="2">I have a grib and ctl file (see attached). with a 2.5X2.5
grid, 145 lons, 73 lats, time=312<span class="455272207-19062006"><font
color="#0000ff"> </font></span></font></font></span></div>
<div><span class="357344306-19062006"><font face="Arial"><font
size="2"><span class="455272207-19062006"></span></font></font></span> </div>
<div><span class="357344306-19062006"><font face="Arial"><font
color="#0000ff" size="2"><span class="455272207-19062006">couldnt
attach grib file as too big but this is the output from a wgrib -V call
on the last record: </span></font></font></span></div>
<div><span class="357344306-19062006"><font face="Arial"><font
color="#0000ff" size="2"><span class="455272207-19062006"> </span></font></font></span></div>
<div><span class="357344306-19062006"><font face="Arial"><font
size="2"><font color="#0000ff"><span class="455272207-19062006"> </span>rec
312:6613726:date 2004120100 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=255 sfc anl:<br>
APCP=Total precipitation [kg/m^2]<br>
timerange 10 P1 0 P2 0 TimeU 1 nx 145 ny 73 GDS grid 0 num_in_ave 0
missing 0<br>
center 54 subcenter 0 process 45 Table 2<br>
latlon: lat -90.000000 to 90.000000 by 2.500000 nxny 10585<br>
long 0.000000 to 360.000000 by 2.500000, (145 x 73) scan 64
mode 136 bdsgrid 1<br>
min/max data -1 1557.09 num bits 16 BDS_Ref -1e+06 DecScale 6
BinScale 15<br>
</font></font></font></span></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">I
want the output in binary format for reading into fortran. I have a
simple script to do this:</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">---------------------------------------------------------------------------------</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">'set
gxout fwrite'<br>
'set fwrite pcpn_xie_arkin_1979_2005_2.5x2.5_ms.dat'</span></font></div>
<div> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">'set
y 1 73'<br>
'set x 1 145'</span></font></div>
<div> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">tvar=1<br>
while (tvar<=312)<br>
'set t 'tvar<br>
'd apcpsfc'<br>
tvar=tvar+1<br>
endwhile</span></font></div>
<div> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">'disable
fwrite'</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">-----------------------------------------------------------------------------------</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">and
my ctl file for the binary data is this:</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"><span
class="455272207-19062006">-----------------------------------------------------------------------------------</span></span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">DSET
pcpn_xie_arkin_1979_2005_2.5x2.5_ms.dat<br>
TITLE xie arkin surf total precip kg/m^2 3mon avg<br>
OPTIONS sequential<br>
UNDEF 99999.9<br>
XDEF 145 LINEAR 0.0 2.5<br>
YDEF 73 LINEAR -90.000 2.5<br>
ZDEF 1 LINEAR 1 1<br>
TDEF 312 LINEAR JAN1979 1mo<br>
VARS 1<br>
prp 1 99 DMP output variable<br>
ENDVARS</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">
<div><font face="Arial" size="2"><span class="357344306-19062006"><span
class="455272207-19062006">-----------------------------------------------------------------------------------</span></span></font></div>
</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">When
I run the script it all looks fine but unfortunately when I examine the
resultant binary file, at each time step the grid shifts slightly in
the x direction. So at a later time step e.g. set t 200, when compared
with the original grib data, though the binary displays the same
pattern it is quite obviously shifted in the x direction. I was
thinking it may have something to do with the 145 lons - as the zero is
repeated twice (0/360)??</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">I
would really appreciate any suggestions as I have also tried decoding
with wgrib and again it looks like it does everything correctly but the
outputted binary pattern at the various time steps, though on the
correct grid this time, displays a different pattern to the original
grib data. So again I would welcome any suggestions from anyone more
familiar with grib. </span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">Thankyou
for your time.</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
<div><font face="Arial" size="2"><span class="357344306-19062006">Regards,</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006">Lyn
Bettio</span></font></div>
<div><font face="Arial" size="2"><span class="357344306-19062006"></span></font> </div>
</blockquote>
</body>
</html>