<div dir="ltr"><pre>Dear Athanasios,<br><br>This is exactly what I typed<br><br>wgrib2 201604250700-fir.grb -set_pdt +0 -grib junk.grb<br>wgrib2 junk.grb -set_grib_type c2 -lola 0:3600:0.1 -90:1801:0.1 junk2.grb grib<br>g2ctl junk2.grb >junk2.ctl<br>gribmap -i junk2.ctl<br><br>The plot should be of a region that is centered over 0N OE.<br><br></pre><div class="gmail_extra">In the second line, latitude range is -90:1801:0.1<br></div><div class="gmail_extra"><pre>I.e, the grid latitude starts at 90S, and goes northward by 0.1 degrees for 1801 points.<br><br></pre><pre>In your version, you had 0:1801:0.1<br><br>So the latitude goes from the equator to 180N. I forgot the case <br>when the latitude is larger than 90N.<br></pre><pre><br></pre><pre>Wesley</pre><pre>PS to get rid of missing fcst statement, modify the line to <br> <br>wgrib2 201604250700-fir.grb -set_pdt +0 -set_ftime anl -grib new_fir.grb2<br><br></pre><div class="gmail_quote">On Tue, Apr 26, 2016 at 12:00 PM, Athanasios Karagiannidis <span dir="ltr"><<a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div><div dir="ltr">Dear Wesley<br><br>When I try to convert to PTD 4.0 I get the following<br><br><i>1:0:d=<a href="tel:2016042507" value="+12016042507" target="_blank">2016042507</a>:FIREDI:no_level:-<a href="tel:2147483647" value="+12147483647" target="_blank">2147483647</a> missing fcst:<br></i><br>but i guess that is not a problem.<br><br>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 <b><u>0</u></b>: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!<i><br></i><br>What can I do for that?<br><br>Kind Regards <br><pre>Athanasios F. Karagiannidis<br>Physicist, PhD Meteorology and Climatology<br>E-mail: <a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a></pre><br><br><div><hr>Date: Tue, 26 Apr 2016 09:34:58 -0400<br>From: <a href="mailto:wesley.ebisuzaki@noaa.gov" target="_blank">wesley.ebisuzaki@noaa.gov</a><br>To: <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>Subject: Re: [gradsusr] meteosat grib2 (grb) files<br><br><div dir="ltr"><div><div><div><div>Dear Athanasios,<br><br></div>There are two problems that need to overcome before displaying with GrADS.<br><br></div>1) The file is using grib2 Product Definition Template (PDT) 4.30. GrADS doesn't<br></div> support PDT 4.30. So the first step to convert to a PDT that GrADS can handle<br></div><div> such as PDT 4.0. (You should consider making the grib file with PDT 4.0 unless<br></div><div> you want to document the spectral bands that are being used.)<br><br></div><div> to convert to PDT 4.0: wgrib2 201604250700-fir.grb -set_pdt +0 -grib new-fir.grb<br></div><div> requires a newer wgrib2<br><br></div><div>If you want to see the raw grid (space view perspective)<br><br></div><div> g2ctl -raw new-fir.grb >new-fir.ctl<br></div><div> gribmap -i new-fir.ctl<br><br></div><div>2) The other problem is that GrADS doesn't have the code to handle space-view perspective<br></div><div> (orographic). Wgrib2 to the rescue again. Here we convert the field to a<br></div><div> longitude-latitude 0.1 degree global grid.<br></div><div><br></div><div> wgrib2 new-fir.grb -set_grib_type c2 -lola 0:3600:0.1 -90:1801:0.1 lola-fir.grb grib<br><br></div> GrADS/g2ctl can handle lola-fir.grb<br><div><div><div><br><br><div>Yours,<br></div><div>Wesley<br></div><div><br><div><div><div><div><br><div>On Tue, Apr 26, 2016 at 5:10 AM, Athanasios Karagiannidis <span dir="ltr"><<a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a>></span> wrote:<br><blockquote style="border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div><div dir="ltr">Hello to all.<br><br>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.<br><br>When trying to create the ctl file i get the following message.<br><br>[...]$ <a href="http://g2ctl.pl" target="_blank">g2ctl.pl</a> -0 201604250700-fir.grb >201604250700-fir.ctl<br>*** script needs to be modified ***<br>unknown user-defined grid<br><br>As a consequence the ctl file that is created can not be opened by gribmap<br><br>[...]$ gribmap -0 -i 201604250700-fir.ctl <br>Open Error: missing XDEF record <br>Open Error: missing YDEF record <br> The data file was not opened<br><br>and no idx file is created.<br><br>Could someone provide any help?<br><pre>Athanasios F. Karagiannidis<br>Physicist, PhD Meteorology and Climatology<br>E-mail: <a href="mailto:thankar@live.com" target="_blank">thankar@live.com</a></pre> </div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div></div></div></div></div></div></div></div></div>
<br>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a></div> </div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div></div>