Hi Zhang,<br><br>since you are new to GrADS I recommend you to take a look at grib2ctl.pl, it's a Perl program that try to catch informations from files in GRIB format and build a control file (.ctl).<br>( <a href="http://www.cpc.ncep.noaa.gov/products/wesley/grib2ctl.html">
http://www.cpc.ncep.noaa.gov/products/wesley/grib2ctl.html</a> )<br>Line <br> APCP 0 99 total precipitation<br>is configured to access a flat binary file, and not a grib file. If grib2ctl does not work you can try to change this, or at least use the wgrib to convert it to a binary file.
<br><br>regards,<br><br>Ricardo.<br><br><div class="gmail_quote">On Nov 15, 2007 11:02 AM, Baolin Zhang <<a href="mailto:baolin@alrc.tottori-u.ac.jp">baolin@alrc.tottori-u.ac.jp</a>> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Heiner,<br><br>Thank you for your reply.<br>Grads is new to me. I made the .ctl by myself, according to wgrib.<br><div class="Ih2E3d">I used wgrib a.grb, the last line is:<br>1812:10894976:d=00121612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=0:sfc:anl:NAve=0
<br></div>Using wgrin -V, I can find the last record as follow:<br>rec 1812:10894976:date 2100121612 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)<br>grid=255 sfc anl:<br>ACPC=total precipitation<br>timerange 10 P1 0 P2 0 TimeU0 nx 57 nny 52 GDS grid 4 num_in_ave 0 missing 0
<br>center 0 subcenter 0 process 0 Table 0<br>gaussian: lat -56.890000 to 38.238000<br> long -33.750000 to 71.250000 by 1.875000, (57X52) scan 0 mode<br>128 bdsgrid 1<br>min/max data 0 8.88966e-05 num bits 16 BDS_Ref 0 DecScale 0 BinScale -29
<br><br>So, I made a .ctl file like this:<br><div class="Ih2E3d">dset ^a.grb<br>dtype grib<br> index ^a.grb.gmp<br> undef -9.99E+33<br> title a.grb<br> xdef 57 linear -56.890000 1.875000<br> ydef 52 linear -33.750000 1.875000
<br> zdef 1 levels 1 1<br> tdef 1812 linear jan1950 1mon<br> vars 1<br> APCP 0 99 total precipitation<br> endvars<br></div>But it does not work.<br>Any suggestion?<br><br>Thank you very much<br><br><br>Baolin Zhang<br><br>
Arid Land Research Center<br>Tottori University<br><div><div></div><div class="Wj3C7c"><br><br><br>----- Original Message -----<br>From: "Heiner Körnich" <<a href="mailto:heiner@MISU.SU.SE">heiner@MISU.SU.SE</a>
><br>To: <<a href="mailto:GRADSUSR@LIST.CINECA.IT">GRADSUSR@LIST.CINECA.IT</a>><br>Sent: Thursday, November 15, 2007 9:48 PM<br>Subject: Re: help me the IPCC AR4 dada<br><br><br>> Hi,<br>><br>> have you run "gribmap" over the ctl-file?
<br>> The entry for index is a bit strange. It should be ^ instead of ~.<br>><br>> Regards,<br>> Heiner<br>><br>> On Thu, 2007-11-15 at 20:11 +0900, Baolin Zhang wrote:<br>>> Dear all<br>>> I want to open IPCC AR4 precipitation data over Africa.
<br>>><br>>> I try the follow:<br>>> dset ^a.grb<br>>> dtype grib<br>>> index ~a.grb.gmp<br>>> undef -9.99E+33<br>>> title a.grb<br>>> xdef 57 linear -56.890000 1.875000<br>
>> ydef 52 linear -33.750000 1.875000<br>>> zdef 1 levels 1 1<br>>> tdef 1812 linear jan1950 1mon<br>>> vars 1<br>>> APCP 0 99 total precipitation<br>>> endvars<br>>><br>>> But the result is: all grid undefined.
<br>>><br>>> Anyone knows the format? I am not sure about the bold lines.<br>>> I used wgrib a.grb, the last line is:<br>>> 1812:10894976:d=00121612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=0:sfc:anl:NAve=0
<br>>> Any suggestions?<br>>><br>>> Thank you<br>>><br>>> Best regards<br>>><br>>><br>>> Baolin Zhang<br>> --<br>> Heiner Körnich<br>> Dept. of Meteorology Tel: +46 8 164333
<br>> Stockholms University, SE-106 91 Stockholm, Sweden<br>> Email: <a href="mailto:heiner@misu.su.se">heiner@misu.su.se</a> <a href="http://www.misu.su.se/%7Eheiner/" target="_blank">www.misu.su.se/~heiner/</a><br>
><br>><br></div></div></blockquote></div><br><br clear="all"><br>-- <br>============================================================<br>Ricardo Marcelo da Silva<br>LAMMA - Laboratório de Modelagem de Processos Marinhos e Atmosféricos
<br>Depto de Meteorologia - Instituto de Geociências <br>Centro de Ciências Matemáticas e da Natureza<br>Universidade Federal do Rio de Janeiro<br>e-mail: <a href="mailto:ricardom@acd.ufrj.br">ricardom@acd.ufrj.br</a> <br>
phone: 55 21 2598-9470 r 24<br>============================================================