[gradsusr] Change of month

Mon Sep 28 18:48:42 EDT 2015

Thanks again Jeff and James.

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.

I really appreciate your help.
Miquel

>>>>

dset ^GFS025_f066.grb
index ^GFS025_f066.grb.idx
undef 9.999E+20
title GFS025_f066.grb
* produced by g2ctl v0.1.1
* command line options: GFS025_f066.grb
* 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)

dtype grib2
ydef 41 linear 35.000000 0.25
xdef 53 linear -5 0.250000
tdef 1 linear 18Z30sep2015 1mo
* PROFILE hPa
zdef 3 levels 100000 85000 50000
options pascals
vars 27
....

>>>>
dset ^GFS025_f072.grb
index ^GFS025_f072.grb.idx
undef 9.999E+20
title GFS025_f072.grb
* produced by g2ctl v0.1.1
* command line options: GFS025_f072.grb
* 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)

dtype grib2
ydef 41 linear 35.000000 0.25
xdef 53 linear -5 0.250000
tdef 1 linear 00Z01oct2015 1mo
* PROFILE hPa
zdef 3 levels 100000 85000 50000
options pascals
vars 27
....

>>>>

URL="
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}
"

# download file
curl "$URL" -o GFS025_f"$hr".grb
sleep 1

#convert to .ctl
perl g2ctl.pl GFS025_f"$hr".grb > GFS025_f"$hr".ctl

#create .idf file
gribmap -i GFS025_f"$hr".ctl

2015-09-28 21:22 GMT+02:00 James T. Potemra <jimp at hawaii.edu>:

> Miquel,
>
> 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.
>
> Jim
>
>
> On 9/28/15 3:12 PM, Miquel Bernis wrote:
>
> 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:
>
> ga-> set t 1
> Time values set: 2015:9:30:6 2015:9:30:6
> ga-> set t 2
> Time values set: 2015:10:30:6 2015:10:30:6
> ga-> set t 3
> Time values set: 2015:11:30:6 2015:11:30:6
> ga-> set t 4
> Time values set: 2015:12:30:6 2015:12:30:6
> ga->
>
> If open the files separately there's no time problem. Kind of weird.
>
> Anyway, you suggested:
>
> "I would just use templates and do the sum over different time index
> values"
>
> 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.
>
> Thanks a lot again.
> Miquel
>
>
>
>
>
>
>
>
> ,
>
> 2015-09-28 17:51 GMT+02:00 Jeff Duda <jeffduda319 at gmail.com>:
>
>> 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.
>>
>> Jeff Duda
>>
>> On Mon, Sep 28, 2015 at 4:19 AM, Miquel Bernis < <bernisfeina at gmail.com>
>> bernisfeina at gmail.com> wrote:
>>
>>> Hi everyone,
>>>
>>> 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:
>>>
>>> Invalid grid coordinates
>>> World coordinates convert to non-integer  grid coordinates
>>> Variable = apcpsfc  Dimension = 3
>>>
>>> Has any of you a solution for that?
>>> Thank you very much.
>>> I'll really appreciate your help.
>>>
>>> Miquel
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
>>
>> --
>> Jeff Duda
>> Graduate research assistant
>> University of Oklahoma School of Meteorology
>> Center for Analysis and Prediction of Storms
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> _______________________________________________
> gradsusr mailing listgradsusr at gradsusr.orghttp://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/20150929/54dae96e/attachment-0001.html 