<div dir="ltr"><div><div><div>Thanks Jennifer, gràcies Oriol!<br><br></div>I finally solved the problem. Now I can generate correctly the template and it seems the best way to work. About that, my only problem was that if you don't download also the first file "f000", the template don't get correctly assembled. Once I realized that it's running ok.<br><br></div>Thanks a lot,<br></div>Miquel <br></div><div class="gmail_extra"><br><div class="gmail_quote">2015-10-05 14:53 GMT+02:00 Oriol Ciurana <span dir="ltr"><<a href="mailto:uri@uriens.cat" target="_blank">uri@uriens.cat</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hola Miquel,<div><br></div><div>Sorry, i haven't seen your mail until now. For GFS0.25 I use templates too, is the best way, and I never haven't had any problem.</div><div>I generate the CTL and IDX with these commands, like yours:</div><div><br></div><div>${config_g2ctl} gfs.t${hini}z.pgrb2.0p25.f%f3 > gfs.t${hini}z.pgrb2.0p25.f000.ctl<br></div><div> gribmap -i gfs.t${hini}z.pgrb2.0p25.f000.ctl<br></div><div><br></div><div>And this is the generated CTL header:</div><div><br></div><div><div>dset ^gfs.t00z.pgrb2.0p25.f%f3</div><div>index ^gfs.t00z.pgrb2.0p25.f000.idx</div><div>undef 9.999E+20</div><div>title gfs.t00z.pgrb2.0p25.f000</div><div>* produced by g2ctl v0.0.9</div><div>* command line options: gfs.t00z.pgrb2.0p25.f%f3</div><div>* griddef=1:0:(1440 x 721):grid_template=0:winds(N/S): lat-lon grid:(1440 x 721) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.250000 lon 0.000000 to 359.750000 by 0.250000 #points=1038240:winds(N/S)</div><div><br></div><div>dtype grib2<br></div><div>options template</div><div>ydef 721 linear -90.000000 0.25</div><div>xdef 1440 linear 0.000000 0.250000</div><div>tdef 41 linear 00Z05oct2015 3hr</div><div>* PROFILE hPa</div><div>zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000</div><div>options pascals</div><div>vars 165</div><div>[...]</div><div>ENDVARS</div></div><div><br></div><div>[CATALAN]</div><div>Com baixes els fitxers de GFS? M'ha semblat entendre que tu només baixes una zona concreta i cada 6h, m'equivoco? Suposo que per guanyar velocitat en la descàrrega.<br></div><div>Prova, ni que només sigui un cop, de baixar els fitxers de tot el món i cada 3h i fes servir les comandes que et poso a veure si et funciona bé i digue'ns. Potser el problema pot venir per aquí...</div><div><br></div><div>Salut!</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Oriol</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">2015-09-29 16:06 GMT+02:00 Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">When you have <b><i>monthly</i></b> data, GrADS is happiest when the hour and day in the TDEF statement are "00z01". Even if your data variable is an accumulation over the entire month, and those variables usually take the time stamp of the end of the accumulation period, everything in GrADS works better (for monthly data) if you use the first of the month in the date stamp instead. If your data is in grib format, use -b with gribmap. <div><br></div><div><div><span><div><br></div><blockquote type="cite">
<div bgcolor="#FFFFFF" text="#000000">Here's the problem. Your first file TDEF is:<br>
<br>
tdef 1 linear 18Z30sep2015 1mo<br></div></blockquote><div><br></div></span><div>Yes, this is a problem. It should be “TDEF 1 linear 00z01sep2015 1mo"</div><span><br><blockquote type="cite"><div bgcolor="#FFFFFF" text="#000000">
<br>
and the second file TDEF is:<br>
<br>
tdef 1 linear 00Z01oct2015 1mo<br></div></blockquote></span>This is okay. </div><div><span><br><blockquote type="cite"><div bgcolor="#FFFFFF" text="#000000">
<br>
So, when you open the first file, the first time value will be set
to 1800 on Sept. 30, 2015. The increment is define as one month
(1mo), so t=2 would presumably be 1800 Oct. 30, 2015, t=3 same for
Nov., and so on. The combined file has:<br>
<br>
tdef 30 linear 06Z29sep2015 6hr<br></div></blockquote><div><br></div></span><div>There are two problems with this entry: the start time and the increment. This is monthly data, correct? </div><div>Try “TDEF 30 linear 00z01sep2015 1mo” instead for the case where you are templating files together. But be careful about your file naming convention and your DSET entry — make sure you are templating over month, not forecast hour. If you need help figuring out templating, please have a look at <a href="http://iges.org/grads/gadoc/templates.html" target="_blank">http://iges.org/grads/gadoc/templates.html</a></div><div><br></div><div><br></div><div>Also, when diagnosing problems, please be sure to run gribmap in the verbose mode (-v) and make sure that all records in your data file are matching. If not, there will be some informative messages about why the record didn’t match. If you’re still having problems, please send the gribmap output along with your descriptor files. </div><div><br></div><div>—Jennifer</div><div><div><div><br></div><br><blockquote type="cite"><div bgcolor="#FFFFFF" text="#000000">
<br>
Here, the initial time is 0600 on Sept. 29, 2015 and the increment
is 6 hours, so the next time step should be 1200 on Sept. 29. I
suspect that your files either have too many, or too few, time
steps. Maybe the files contain forecast information that overlaps?
It looks like the names of the files are f066 then f072, so maybe 6
time-steps in each, which means (complete speculation here) that
they have 6 time-steps in each for a total of 36 hours? If true,
the first file starts at 0600 Sept 29, the next file should start
1200 Sept 30?<br>
<br>
Jim<br>
<br>
<div>On 9/29/15 2:32 AM, Miquel Bernis
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">Thanks for your instructions Jeff and James, I
think I understood what you said, and I tried that way. It seems
to be everything right, but when I open de .ctl and try to
display a var, I get the error:<br>
<br>
Cannot contour grid - all undefined values <br>
<div>
<div><br>
</div>
<div>These are commands that I execute to create de template,
and after that de .ctl header that results. I don't know
what might the mistake.<br>
<br>
</div>
<div>Sorry for being so hopeless on that.<br>
</div>
<div>Miquel<br>
</div>
<div><br>
>>><br>
<br>
#el convertim a .ctl<br>
perl <a href="http://g2ctl.pl/" target="_blank">g2ctl.pl</a>
GFS025_f%f3.grb >GFS025.ctl<br>
<br>
#creem arxiu .idf<br>
gribmap -i GFS025.ctl<br>
<br>
>>><br>
<br>
dset ^GFS025_f%f3.grb<br>
index ^GFS025_f006.grb.idx<br>
undef 9.999E+20<br>
title GFS025_f006.grb<br>
* produced by g2ctl v0.1.1<br>
* command line options: GFS025_f%f3.grb<br>
* griddef=1:0:(53 x 41):grid_template=0:winds(N/S): lat-lon
grid:(53 x 41) units 1e-06 input WE:SN output WE:SN res 48
lat 35.000000 to 45.000000 by 0.250000 lon 355.000000 to
8.000000 by 0.250000 #points=2173:winds(N/S)<br>
<br>
dtype grib2<br>
options template<br>
ydef 41 linear 35.000000 0.25<br>
xdef 53 linear -5 0.250000<br>
tdef 30 linear 06Z29sep2015 6hr<br>
* PROFILE hPa<br>
zdef 3 levels 100000 85000 50000<br>
options pascals<br>
vars 27<br>
<br>
<br>
</div>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-09-29 1:23 GMT+02:00 Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>As James mentioned and I suspected, doing this one
file at a time is harder than using templates. What I
would do is download all the GFS data files, then use
one single control file to open all of them. Use
templating, as described <a href="http://www.iges.org/grads/gadoc/templates.html" target="_blank">here</a>. Then adjust the TDEF line of
the control file accordingly. The last thing to do is to
run gribmap on that control file.<span><font color="#888888"><br>
<br>
</font></span></div>
<span><font color="#888888">Jeff<br>
</font></span></div>
<div>
<div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Mon, Sep 28, 2015 at 5:48
PM, Miquel Bernis <span dir="ltr"><<a href="mailto:bernisfeina@gmail.com" target="_blank"></a><a href="mailto:bernisfeina@gmail.com" target="_blank">bernisfeina@gmail.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>Thanks again Jeff and James. <br>
<br>
I've pasted here a pair of my .ctl file
example headers. After that, you have the
commands I use for downloading the GFS 0.25
data, and for the conversion to .ctl. As you
can see I use the perl g2ctl. <br>
<br>
</div>
<div>I really appreciate your help. <br>
</div>
<div>Miquel <br>
</div>
<div><br>
>>>><br>
<br>
dset ^GFS025_f066.grb<br>
index ^GFS025_f066.grb.idx<br>
undef 9.999E+20<br>
title GFS025_f066.grb<br>
* produced by g2ctl v0.1.1<br>
* command line options: GFS025_f066.grb<br>
* griddef=1:0:(53 x
41):grid_template=0:winds(N/S): lat-lon
grid:(53 x 41) units 1e-06 input WE:SN output
WE:SN res 48 lat 35.000000 to 45.000000 by
0.250000 lon 355.000000 to 8.000000 by
0.250000 #points=2173:winds(N/S)<br>
<br>
dtype grib2<br>
ydef 41 linear 35.000000 0.25<br>
xdef 53 linear -5 0.250000<br>
tdef 1 linear 18Z30sep2015 1mo<br>
* PROFILE hPa<br>
zdef 3 levels 100000 85000 50000<br>
options pascals<br>
</div>
vars 27<br>
....<br>
<br>
<div>>>>><br>
dset ^GFS025_f072.grb<br>
index ^GFS025_f072.grb.idx<br>
undef 9.999E+20<br>
title GFS025_f072.grb<br>
* produced by g2ctl v0.1.1<br>
* command line options: GFS025_f072.grb<br>
* griddef=1:0:(53 x
41):grid_template=0:winds(N/S): lat-lon
grid:(53 x 41) units 1e-06 input WE:SN output
WE:SN res 48 lat 35.000000 to 45.000000 by
0.250000 lon 355.000000 to 8.000000 by
0.250000 #points=2173:winds(N/S)<br>
<br>
dtype grib2<br>
ydef 41 linear 35.000000 0.25<br>
xdef 53 linear -5 0.250000<br>
tdef 1 linear 00Z01oct2015 1mo<br>
* PROFILE hPa<br>
zdef 3 levels 100000 85000 50000<br>
options pascals<br>
</div>
<div>vars 27<br>
....<br>
</div>
<div><br>
>>>><br>
<br>
URL="<a href="http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_0p25.pl?file=gfs.t$%7Bfhr%7Dz.pgrb2.0p25.f$%7Bhr%7D&lev_10_m_above_ground=on&lev_2_m_above_ground=on&lev_500_mb=on&lev_850_mb=on&lev_1000_mb=on&lev_entire_atmosphere=on&lev_high_cloud_layer=on&lev_low_cloud_layer=on&lev_middle_cloud_layer=on&lev_surface=on&var_APCP=on&var_CAPE=on&var_CFRZR=on&var_CICEP=on&var_CRAIN=on&var_CSNOW=on&var_DPT=on&var_GUST=on&var_HGT=on&var_LFTX=on&var_PRATE=on&var_RH=on&var_TCDC=on&var_TMAX=on&var_TMIN=on&var_TMP=on&var_UGRD=on&var_VGRD=on&subregion=&leftlon=-5&rightlon=8&toplat=45&bottomlat=35&dir=%2Fgfs.$%7Bdate%7D$%7Bfhr%7D" target="_blank">http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_0p25.pl?file=gfs.t${fhr}z.pgrb2.0p25.f${hr}&lev_10_m_above_ground=on&lev_2_m_above_ground=on&lev_500_mb=on&lev_850_mb=on&lev_1000_mb=on&lev_entire_atmosphere=on&lev_high_cloud_layer=on&lev_low_cloud_layer=on&lev_middle_cloud_layer=on&lev_surface=on&var_APCP=on&var_CAPE=on&var_CFRZR=on&var_CICEP=on&var_CRAIN=on&var_CSNOW=on&var_DPT=on&var_GUST=on&var_HGT=on&var_LFTX=on&var_PRATE=on&var_RH=on&var_TCDC=on&var_TMAX=on&var_TMIN=on&var_TMP=on&var_UGRD=on&var_VGRD=on&subregion=&leftlon=-5&rightlon=8&toplat=45&bottomlat=35&dir=%2Fgfs.${date}${fhr}</a>"<br>
<br>
# download file<br>
curl "$URL" -o GFS025_f"$hr".grb<br>
sleep 1<br>
<br>
#convert to .ctl<br>
perl <a href="http://g2ctl.pl/" target="_blank">g2ctl.pl</a>
GFS025_f"$hr".grb > GFS025_f"$hr".ctl<br>
<br>
#create .idf file<br>
gribmap -i GFS025_f"$hr".ctl<br>
<br>
</div>
</div>
<div>
<div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-09-28 21:22
GMT+02:00 James T. Potemra <span dir="ltr"><<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Miquel,<br>
<br>
We'd need to see you control files,
but I suspect the TDEF line is causing
problems. I believe this gets set
with your first open file. For
example, if you open a file and "set t
1", time will get set to some date as
defined in that control file. If the
subsequent files don't have this date
defined you'll get the grid coordinate
error. You could issue a "set dfile"
command to re-set the coordinates.<br>
<br>
Jim
<div>
<div><br>
<br>
<div>On 9/28/15 3:12 PM, Miquel
Bernis wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>Thanks a lot for your
answer Jeff. I'm not sure
how can I give more clues
about the error. There's no
problem when I work with the
files separately, only when
I operate between files in a
different month. I realized
that there's something wrong
about the time vars. When I
open consecutive files 24z,
30z, 36z, 42z... the first
one has the right time
value, but the next ones
shows a change of month,
instead of just a 6h
advance:<br>
<br>
ga-> set t 1<br>
Time values set: 2015:9:30:6
2015:9:30:6 <br>
ga-> set t 2<br>
Time values set:
2015:10:30:6 2015:10:30:6 <br>
ga-> set t 3<br>
Time values set:
2015:11:30:6 2015:11:30:6 <br>
ga-> set t 4<br>
Time values set:
2015:12:30:6 2015:12:30:6 <br>
ga-> </div>
<div><br>
</div>
<div>If open the files
separately there's no time
problem. Kind of weird.</div>
<div><br>
</div>
<div>Anyway, you suggested:</div>
<div><br>
</div>
<div>"<span style="font-size:12.8px">I
would just use templates
and do the sum over
different time index
values"</span></div>
<div><br>
</div>
<div>I'm sorry, I don't
understand that. Maybe you
can tell the same in another
way or with more details. It
would be helpful.</div>
<div><br>
</div>
<div>Thanks a lot again.</div>
<div>Miquel</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<br>
<br>
</div>
<br>
<div><br>
<br>
, <br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-09-28
17:51 GMT+02:00 Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank"></a><a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>Perhaps some more
information on the
specifics of the
problem would help us
to troubleshoot.
However, if you're
summing precip totals
from a single GFS run,
I would just use
templates and do the
sum over different
time index values.<br>
<br>
</div>
Jeff Duda<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">
<div>
<div>On Mon, Sep 28,
2015 at 4:19 AM,
Miquel Bernis <span dir="ltr"><<a href="mailto:bernisfeina@gmail.com" target="_blank"></a><a href="mailto:bernisfeina@gmail.com" target="_blank">bernisfeina@gmail.com</a>></span>
wrote:<br>
</div>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div>
<div dir="ltr">
<div style="font-size:12.8px">
<div>
<div>Hi
everyone,<br>
<br>
</div>
I'm stucked in
a problem that
might be
simple to fix,
maybe any of
you could help
me. The point
is that I
can't operate
with files
when there's
change of
month. I'm
working with
GFS 0.25, and
when I try to
sum the rain
of the last
file of month
with the rain
of the first
file of the
moth ahead, it
returns and
error like
this:<br>
<br>
Invalid grid
coordinates<br>
World
coordinates
convert to
non-integer
grid
coordinates<br>
Variable =
apcpsfc
Dimension = 3 <br>
<br>
</div>
Has any of you
a solution for
that? <br>
Thank you very
much.<br>
</div>
<div style="font-size:12.8px">I'll
really
appreciate
your help.</div>
<span><font color="#888888">
<div style="font-size:12.8px"><br>
</div>
<div style="font-size:12.8px">Miquel</div>
</font></span></div>
<br>
</div>
</div>
_______________________________________________<br>
gradsusr mailing
list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank"></a><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"></a><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
</blockquote>
</div>
<span><font color="#888888"><br>
<br clear="all">
<br>
-- <br>
<div>Jeff Duda<br>
Graduate research
assistant<br>
University of
Oklahoma School of
Meteorology<br>
Center for
Analysis and
Prediction of
Storms<br>
</div>
</font></span></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>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
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>
</pre>
</blockquote>
<br>
</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>
</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>
<br clear="all">
<br>
-- <br>
<div>Jeff Duda<br>
Graduate research assistant<br>
University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms<br>
</div>
</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>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
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>
</pre>
</blockquote>
<br>
</div>
_______________________________________________<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" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br></blockquote></div></div></div><br><div>
<span style="border-collapse:separate;border-spacing:0px"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div>--</div><div>Jennifer M. Adams<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>111 Research Hall, Mail Stop 2B3<br>George Mason University<br>4400 University Drive<br>Fairfax, VA 22030 <br><br></div><div><br></div><br></span></span><br>
</div>
<br></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><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" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>