[gradsusr] meteosat grib2 (grb) files

Thu Apr 28 03:58:11 EDT 2016

Dear Wesley

Thanx for your help.

This time it worked! I really don't know what happened in my previous try, but your solution seems to be working perfectly.

Kind Regards

Athanasios F. Karagiannidis
Physicist, PhD Meteorology and Climatology
E-mail: thankar at live.com

Date: Wed, 27 Apr 2016 15:11:18 -0400
From: wesley.ebisuzaki at noaa.gov
To: gradsusr at gradsusr.org
Subject: Re: [gradsusr] meteosat grib2 (grb) files

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.

WesleyPS 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

_______________________________________________
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/20160428/c622cd35/attachment.html 