[gradsusr] meteosat grib2 (grb) files

Wed Apr 27 15:11:18 EDT 2016

Athanasios,

I recreated the plot using the directions given.  See the attachment.
I am using wgrib2 v2.0.4 on a Redhat Linux system.

Wesley

On Wed, Apr 27, 2016 at 1:20 PM, Athanasios Karagiannidis <thankar at live.com>
wrote:

> Dear Wesley
>
> Firstly I tried exactly the same thing that you typed (-90:1801:0.1) and I
> got the wrong visualisation of "test1.png" I sent. Then I experimented with
> other values (e.g. 0:1801:01) and did not succeed either.
> As you can see in Test1.png the image seems to be shifted to the north by
> 90 degrees.
> Could that be corrected somehow?
>
>
> Kind regards
>
> Athanasios F. Karagiannidis
> Physicist, PhD Meteorology and Climatology
> E-mail: thankar at live.com
>
>
>
> ------------------------------
> Date: Tue, 26 Apr 2016 13:23:23 -0400
> From: wesley.ebisuzaki at noaa.gov
> To: gradsusr at gradsusr.org
> Subject: Re: [gradsusr] meteosat grib2 (grb) files
>
> Dear Athanasios,
>
> This is exactly what I typed
>
> wgrib2 201604250700-fir.grb -set_pdt +0 -grib junk.grb
> wgrib2 junk.grb  -set_grib_type c2 -lola 0:3600:0.1 -90:1801:0.1 junk2.grb  grib
> g2ctl junk2.grb >junk2.ctl
> gribmap -i junk2.ctl
>
> The plot should be of a region that is centered over 0N OE.
>
> In the second line, latitude range is     -90:1801:0.1
>
> I.e, the grid latitude starts at 90S, and goes northward by 0.1 degrees for 1801 points.
>
> In your version, you had 0:1801:0.1
>
> So the latitude goes from the equator to 180N. I forgot the case
> when the latitude is larger than 90N.
>
>
> Wesley
>
> PS to get rid of missing fcst statement, modify the line to
>
> wgrib2 201604250700-fir.grb -set_pdt +0 -set_ftime anl -grib new_fir.grb2
>
> On Tue, Apr 26, 2016 at 12:00 PM, Athanasios Karagiannidis <
> thankar at live.com> wrote:
>
> Dear Wesley
>
> When I try to convert to PTD 4.0 I get the following
>
>
> *1:0:d=2016042507:FIREDI:no_level:-2147483647 missing fcst:*
> but i guess that is not a problem.
>
> However, after converting to the new grid something goes wrong because the
> image seems shifted by 90 degrees to the north (please see test1.png). I
> tried some tests but nothing worked. For example when I asked for the
> regriding to start from 0 degrees north (wgrib2 new-fir.grb -set_grib_type
> c2 -lola 0:3600:0.1 *0*:1801:0.1 lola-fir.grb grib), I got the image of
> test2.png. The image there seems cut and not shifted to the north!
>
> What can I do for that?
>
> Kind Regards
>
> Athanasios F. Karagiannidis
> Physicist, PhD Meteorology and Climatology
> E-mail: thankar at live.com
>
>
>
> ------------------------------
> Date: Tue, 26 Apr 2016 09:34:58 -0400
> From: wesley.ebisuzaki at noaa.gov
> To: gradsusr at gradsusr.org
> Subject: Re: [gradsusr] meteosat grib2 (grb) files
>
> Dear Athanasios,
>
> There are two problems that need to overcome before displaying with GrADS.
>
> 1)  The file is using grib2 Product Definition Template (PDT) 4.30.
> GrADS doesn't
>       support PDT 4.30.   So the first step to convert to a PDT that GrADS
> can handle
>       such as PDT 4.0.   (You should consider making the grib file with
> PDT 4.0 unless
>       you want to document the spectral bands that are being used.)
>
>       to convert to PDT 4.0:    wgrib2 201604250700-fir.grb -set_pdt +0
> -grib new-fir.grb
>       requires a newer wgrib2
>
> If you want to see the raw grid (space view perspective)
>
>      g2ctl -raw new-fir.grb >new-fir.ctl
>      gribmap -i new-fir.ctl
>
> 2)  The other problem is that GrADS doesn't have the code to handle
> space-view perspective
>        (orographic).   Wgrib2 to the rescue again.  Here we convert the
> field to a
>        longitude-latitude 0.1 degree global grid.
>
>         wgrib2 new-fir.grb -set_grib_type c2 -lola 0:3600:0.1 -90:1801:0.1
> lola-fir.grb grib
>
>  GrADS/g2ctl can handle lola-fir.grb
>
>
> Yours,
> Wesley
>
>
> On Tue, Apr 26, 2016 at 5:10 AM, Athanasios Karagiannidis <
> thankar at live.com> wrote:
>
> Hello to all.
>
> I'm trying to open and visualize the contents of the "fire detection
> indicator" parameter that is included in meteosat grib2 file like the one I
> attach.
>
> When trying to create the ctl file i get the following message.
>
> [...]$ g2ctl.pl -0 201604250700-fir.grb >201604250700-fir.ctl
> *** script needs to be modified ***
> unknown user-defined grid
>
> As a consequence the ctl file that is created can not be opened by gribmap
>
> [...]$ gribmap -0 -i 201604250700-fir.ctl
> Open Error:  missing XDEF record
> Open Error:  missing YDEF record
>   The data file was not opened
>
> and no idx file is created.
>
> Could someone provide any help?
>
> Athanasios F. Karagiannidis
> Physicist, PhD Meteorology and Climatology
> E-mail: thankar at live.com
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> _______________________________________________ gradsusr mailing list
> gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> _______________________________________________ gradsusr mailing list
> gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160427/acc9a4b5/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fir.png
Type: image/png
Size: 14064 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20160427/acc9a4b5/attachment-0001.png 