[gradsusr] Unable To Open Concatanated Regional GEM GRIB2s With GrADS

[Perry, Aaron @ LSC]

Hi Greg,

Thanks again for taking the time and effort to help me out with this issue. I'm actually using Wesley's get_inv/get_grib perl scripts to pull specific variables from the CMC server.

Here's the bash script that I'm using to try to pull the 10km RGEM grib2s with: http://paste2.org/DjhcB8UI.

I also cat all of the grib2s together into one big grib2. I just went back to check and I also get the same exact grid description in the ctl file as you do.

Do you see anything obviously wrong that I'm doing in my bash script?



Onward and Upward,

Aaron D. Perry
Atmospheric Sciences
Lyndon State College
Class of 2016
Boston, MA
Twitter: @AaronPerryLSC<https://twitter.com/AaronPerryLSC>
Mobile: 617-780-4312
________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on behalf of Greg West [gregory.l.west at gmail.com]
Sent: Tuesday, March 25, 2014 12:27 AM
To: gradsusr at gradsusr.org
Subject: Re: [gradsusr] Unable To Open Concatanated Regional GEM GRIB2s With GrADS

Aaron,

I also download it from that source. When I use g2ctl.pl<http://g2ctl.pl> to create ctl files, if I remove CAPE and 4LFTX from the grib2 files, this is the grid description in the ctl file:

pdef 935 824 nps 456.142068646443 732.42845661651 249.000000 10
xdef 4000 linear 0.000000 0.0899962568142035
ydef 808 linear 17.342726 0.0900337967781908

If I separate out CAPE and 4LFTX into their own grib2 file, this is what I get:

pdef 858 665 nps 361.120321747126 683.553026459468 249.000000 10
xdef 1945 linear 161.949691 0.0899750853909465
ydef 767 linear 19.241056 0.0900398276762402

You can see that the convective indices are on a different grid for some reason.

These are the commands in my perl script I use to work around it...
Datamart, as you know, has every variable and fcst hr in its own file. First take the files containing the 2 special vars and put them in their own grib2 file...
`cat \`ls $data_root/$modcap/$date$hh/*4LFTX*_P$fhrname.grib2 $data_root/$modcap/$date$hh/*CAPE*_P$fhrname.grib2\`
> $data_root/$modcap/$date$hh/$mod.$date$hh.conv.f$fhrname`
Then take all the other files (except the 2 special vars) and put them into a grib2 file.
`cat \`ls $data_root/$modcap/$date$hh/*_P$fhrname.grib2 | grep -v '.*4LFTX.*' | grep -v '.*CAPE.*'\` > $data_root/$
modcap/$date$hh/$mod.$date$hh.master.f$fhrname`
Then I create different template ctl/idx files for each set of files. Your plotting scripts will have to open the 2nd ctl file if you're plotting the special vars.

Greg

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140324/48e84100/attachment-0001.html

------------------------------

Message: 3
Date: Mon, 24 Mar 2014 22:49:45 +0000
From: "Perry, Aaron @ LSC" <Aaron.Perry at lsc.vsc.edu<mailto:Aaron.Perry at lsc.vsc.edu>>
Subject: Re: [gradsusr] Unable To Open Concatanated Regional GEM
        GRIB2s  With    GrADS
To: GrADS Users Forum <gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>>
Message-ID:
        <DEE427CC71C1184E9F7C903BF1EEE7A80B65D4BC at LSC-EX10MBX03.lsc.vsc.edu<mailto:DEE427CC71C1184E9F7C903BF1EEE7A80B65D4BC at LSC-EX10MBX03.lsc.vsc.edu>>
Content-Type: text/plain; charset="iso-8859-1"

Hi Greg,

Thank you very much for chiming in! I've been using the usual http://dd.weatheroffice.gc.ca/model_gem_regional/10km/grib2/ datasource to pull the GRIB2s and then run them through g2ctl/gribmap.

Truly stumped right now...



Onward and Upward,

Aaron D. Perry
Atmospheric Sciences
Lyndon State College
Class of 2016
Boston, MA
Twitter: @AaronPerryLSC<https://twitter.com/AaronPerryLSC>
Mobile: 617-780-4312<tel:617-780-4312>
________________________________



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140325/0e869f54/attachment.html