[gradsusr] Change of month

Mon Sep 28 19:23:15 EDT 2015

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 here <http://www.iges.org/grads/gadoc/templates.html>. Then
adjust the TDEF line of the control file accordingly. The last thing to do
is to run gribmap on that control file.

Jeff

On Mon, Sep 28, 2015 at 5:48 PM, Miquel Bernis <bernisfeina at gmail.com>
wrote:

> 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}
> <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>
> "
>
> # 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
>>
>>
>
> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150928/5f23033a/attachment.html 