help me the IPCC AR4 dada

[Ricardo Marcelo]

Hi Zhang,

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).
( http://www.cpc.ncep.noaa.gov/products/wesley/grib2ctl.html )
Line
 APCP 0 99 total precipitation
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.

regards,

Ricardo.

On Nov 15, 2007 11:02 AM, Baolin Zhang <baolin at alrc.tottori-u.ac.jp> wrote:

> Dear Heiner,
>
> Thank you for your reply.
> Grads is new to me. I made the .ctl by myself, according to wgrib.
> I used wgrib a.grb, the last line is:
>
> 1812:10894976:d=00121612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=0:sfc:anl:NAve=0
> Using wgrin -V, I can find the last record as follow:
> rec 1812:10894976:date 2100121612 APCP kpds5=61 kpds6=1 kpds7=0
> levels=(0,0)
> grid=255 sfc anl:
> ACPC=total precipitation
> timerange 10 P1 0 P2 0 TimeU0 nx 57 nny 52 GDS grid 4 num_in_ave 0 missing
> 0
> center 0 subcenter 0 process 0 Table 0
> gaussian: lat -56.890000 to 38.238000
>             long -33.750000 to 71.250000 by 1.875000, (57X52) scan 0 mode
> 128 bdsgrid 1
> min/max data 0 8.88966e-05 num bits 16 BDS_Ref 0 DecScale 0 BinScale -29
>
> So, I made a .ctl file like this:
> dset ^a.grb
> dtype grib
>  index ^a.grb.gmp
>  undef -9.99E+33
>  title a.grb
>  xdef 57 linear -56.890000 1.875000
>  ydef 52 linear -33.750000 1.875000
>  zdef 1 levels 1 1
>  tdef 1812 linear jan1950 1mon
>  vars 1
>  APCP 0 99 total precipitation
>  endvars
> But it does not work.
> Any suggestion?
>
> Thank you very much
>
>
> Baolin Zhang
>
> Arid Land Research Center
> Tottori University
>
>
>
> ----- Original Message -----
> From: "Heiner Körnich" <heiner at MISU.SU.SE>
> To: <GRADSUSR at LIST.CINECA.IT>
> Sent: Thursday, November 15, 2007 9:48 PM
> Subject: Re: help me the IPCC AR4 dada
>
>
> > Hi,
> >
> > have you run "gribmap" over the ctl-file?
> > The entry for index is a bit strange. It should be ^ instead of ~.
> >
> > Regards,
> > Heiner
> >
> > On Thu, 2007-11-15 at 20:11 +0900, Baolin Zhang wrote:
> >> Dear all
> >> I want to open IPCC AR4 precipitation data over Africa.
> >>
> >> I try the follow:
> >> dset ^a.grb
> >> dtype grib
> >> index ~a.grb.gmp
> >> undef -9.99E+33
> >> title a.grb
> >> xdef 57 linear -56.890000 1.875000
> >> ydef 52 linear -33.750000 1.875000
> >> zdef 1 levels 1 1
> >> tdef 1812 linear jan1950 1mon
> >> vars 1
> >> APCP 0 99 total precipitation
> >> endvars
> >>
> >> But the result is: all grid undefined.
> >>
> >> Anyone knows the format? I am not sure about the bold lines.
> >> I used wgrib a.grb, the last line is:
> >>
> 1812:10894976:d=00121612:APCP:kpds5=61:kpds6=1:kpds7=0:TR=10:P1=0:P2=0:TimeU=0:sfc:anl:NAve=0
> >> Any suggestions?
> >>
> >> Thank you
> >>
> >> Best regards
> >>
> >>
> >> Baolin Zhang
> > --
> > Heiner Körnich
> > Dept. of Meteorology            Tel:  +46 8 164333
> > Stockholms University, SE-106 91 Stockholm, Sweden
> > Email: heiner at misu.su.se   www.misu.su.se/~heiner/<http://www.misu.su.se/%7Eheiner/>
> >
> >
>



-- 
============================================================
Ricardo Marcelo da Silva
LAMMA - Laboratório de Modelagem de Processos Marinhos e Atmosféricos
Depto de Meteorologia - Instituto de Geociências
Centro de Ciências Matemáticas e da Natureza
Universidade Federal do Rio de Janeiro
e-mail: ricardom at acd.ufrj.br
phone: 55 21 2598-9470 r 24
============================================================
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20071115/b94bce29/attachment.html