<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body style='font-size: 10pt; font-family: Tahoma,Arial,Helvetica,sans-serif'>
<p> </p>
<p>Hi Guys,</p>
<p> </p>
<p>How would this work around apply to windows users who use perl!!!! Normally a ctl file would be produced like such:</p>
<p>@ECHO Please wait, concetenating grb files...<br />cat *.grb > gfs_grib.grib2<br />@ECHO Please wait, creating ctl file...<br />perl g2ctl.pl -verf gfs_grib.grib2 > gfs_grib.ctl<br />@ECHO Please wait, creating idx file...<br />gribmap -i gfs_grib.ctl -e</p>
<p>regards,</p>
<p>Stuart</p>
<p> </p>
<p> </p>
<p>On 17-12-2015 22:17, Brian Gaze wrote:</p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px"><!-- html ignored --><!-- head ignored --><!-- meta ignored -->
<p>Wesley,</p>
<p>Many thanks for the explanation and solution using alt_g2ctl and alt_gmp. </p>
<p>Brian</p>
<p> </p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px">On 17 December 2015 at 13:14 Wesley Ebisuzaki - NOAA Federal <wesley.ebisuzaki@noaa.gov> wrote: <br /> <br />
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>Brian, <br /> </div>
I am being pedantic because I hope that the ARPEGE and AROME </div>
people will get word and fix their grib2 files. <br /> </div>
In order for GrADS to plot the ARPEGE file, GrADS needs to know the </div>
time of the fields. For Product Definition Template (PDT) 4.8 which is <br />used to describe averages and accumulations, there are 3 different times. <br /> </div>
Reference time: start of the forecast <br />
<div>Section 1, octets 13-19 <br /> </div>
<div>Start of the forecast time interval: </div>
<div>reference time + forecast time (Section 3, octets 18-22) <br /> </div>
<div>End of the forecast time interval: </div>
<div>Section 3, octets 35-41 <br /> </div>
<div>For your ARPEGE file, the only unique time should be the end of the forecast </div>
<div>time interval. However when you examine the ARPEGE file. <br /> <br />ebis@linux-landing:~/Downloads$ wgrib2 ARPEGE_0.1_SP1_37H48H_201512150000.grib2 -s -start_ft -end_ft -match PRATE <br />133:32079786:d=2015121500:TPRATE:surface:0-37 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />134:32391856:d=2015121500:TPRATE:surface:0-38 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />135:32706662:d=2015121500:TPRATE:surface:0-39 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />136:33024062:d=2015121500:TPRATE:surface:0-40 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />137:33344053:d=2015121500:TPRATE:surface:0-41 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />138:33666613:d=2015121500:TPRATE:surface:0-42 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />139:33991684:d=2015121500:TPRATE:surface:0-43 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />140:34319034:d=2015121500:TPRATE:surface:0-44 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />141:34648323:d=2015121500:TPRATE:surface:0-45 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />142:34979540:d=2015121500:TPRATE:surface:0-46 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />143:35312524:d=2015121500:TPRATE:surface:0-47 hour acc fcst::start_ft=2015121500:end_ft=2015121600 <br />144:35647622:d=2015121500:TPRATE:surface:0-2 day acc fcst::start_ft=2015121500:end_ft=2015121700 <br /> </div>
<div>You see that the TPRATE (totat precipitation rate) has the end of </div>
<div>the forecast time interval of 2015121600 which is wrong. This ARPEGE file </div>
<div>needs to be fixed. A few weeks ago, someone sent me a AROME </div>
<div>forecast file which had the same problem. <br /> </div>
<div>So g2ctl and gribmap will not work until the grib file is fixed. <br /> </div>
<div>You can display the fields using <br /> <br /> alt_g2ctl -0t ARPEGE_0.1_SP1_37H48H_201512150000.grib2 >a.ctl <br /> alt_gmp -i a.ctl <br /> </div>
<div>This solution gives each field a unique name and doesn't use the </div>
<div>end of the forecast time. <br /> </div>
<div>Maybe someone will inform the ARPEGE and AROME people and they <br />will fix their files. </div>
<div> </div>
<div>Wesley </div>
<div><br />
<div>
<div>
<div>
<div>
<div>
<div class="gmail_extra"><br />
<div class="gmail_quote">On Thu, Dec 17, 2015 at 6:13 AM, Brian Gaze < <a href="mailto:brian.gaze@ntlworld.com">brian.gaze@ntlworld.com</a>> wrote: <br />
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px"><span style="text-decoration: underline;"></span>
<div>
<p>Wesley,</p>
<p>You can ftp down the grib file here using anonymous:</p>
<p>137.135.248.188</p>
<p>ARPEGE_0.1_SP1_37H48H_201512150000.grib2</p>
<p>PS: It's about 40 meg and I tested the ftp download successfully from a command line rather than gui client</p>
<p>Cheers</p>
<p>Brian</p>
<p> </p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px">On 16 December 2015 at 15:03 Wesley Ebisuzaki - NOAA Federal < <a href="mailto:wesley.ebisuzaki@noaa.gov">wesley.ebisuzaki@noaa.gov</a>> wrote: <br /> <br />
<div dir="ltr">
<div>Brian, <br /> </div>
<div>It would be much easier if you made your grib file available </div>
<div>by FTP. If the file is huge or is a proprietary forecast, you can </div>
<div>zero out the data by <br /> </div>
<div>wgrib2 file.grib -rpn 0 -grib_out file00.grb </div>
<div> </div>
<div>Wesley </div>
<br /> <br />
<div>
<div>
<div class="gmail_extra"><br />
<div class="gmail_quote">On Tue, Dec 15, 2015 at 4:32 PM, Brian Gaze < <a href="mailto:brian.gaze@ntlworld.com">brian.gaze@ntlworld.com</a>> wrote: <br />
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px">
<div>
<p>Hi,</p>
<p>I'm working with Arpège grib2 files. These should I think each contain up to 13 hourly forecast timesteps. If I run...</p>
<p><a href="http://g2ctl.pl">g2ctl.pl</a> -verf -b file.grib2 > file.ctl<br />gribmap -i file.ctl</p>
<p>...I'm able to open the .ctl file in GrADS but it only has 2 timesteps. T1 seems to be the time the model was initialized and T2 the time of the first forecast step.</p>
<p> </p>
<p>If I try:</p>
<p><a href="http://g2ctl.pl">g2ctl.pl</a> -b -ts1hr file.grib2 > file.ctl<br />gribmap -i file.ctl</p>
<p>It opens ok in GrADS and the correct number of timesteps are present but they begin from the model initialization hour rather than incrementing from the time of the first forecast step. For example:</p>
<p>Model initialization: 0z 15/12<br />Forecast time range contained in the file is +37 hours to +48 hours after initialization</p>
<p>t2 = 1z 15/12<br />t3 = 2z 15/12<br />.. etc<br /><br />As a consequence, when I try to plot a variable the grid is undefined.</p>
<p> </p>
<p>What should I think be present is:<br />t2 = 13z 16/12<br />t3 = 14z 16/12<br />..etc</p>
<p>Any suggestions / ideas on what I'm doing wrong?</p>
<p>Thanks</p>
<p>BWG</p>
</div>
<br />_______________________________________________ <br /> gradsusr mailing list <br /> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <br /> <a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a> <br /> </blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<p> </p>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<p><br /> </p>
<!-- html ignored --><br />
<pre>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<p> </p>
<div>-- <br />
<p>Stuart Markham</p>
<p><img src="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG" alt="" /></p>
</div>
</body></html>