[gradsusr] Issue with duplicate Varible in GRIB Data

Wesley Ebisuzaki - NOAA Federal wesley.ebisuzaki at noaa.gov
Wed Jan 18 16:02:02 EST 2023


Mike,

    The directions in README.ICON.DWD
<https://www.ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/_README.ICON.DWD> are for
getting the lat/lon
values associated with each grid point.  Once you do that, you can
use -lola to create a lat/lon grid using nearest neighbor interpolation.
Then you can use that lat/lon file with GrADS.  It wouldn't be that hard to
modify wgrib2 to generate a pdef file so that you could read the ICON
file directly.

Wesley





On Fri, Jan 13, 2023 at 3:24 AM <mike at weatherwatch.net.au> wrote:

> Hi Wesley,
>
>
>
> That worked a treat thanks with alt_g2ctl.  I ended up using that and it’s
> working great (had to use -short as the variable length exceeded 15
> characters):
>
> alt_g2ctl -short IDY25001.APS3.syncld.slv.2023011218.000.surface.grb2
>
> dset ^IDY25001.APS3.syncld.slv.2023011218.000.surface.grb2
>
> index ^IDY25001.APS3.syncld.slv.2023011218.000.surface.grb2.idx
>
> undef 9.999E+20
>
> title IDY25001.APS3.syncld.slv.2023011218.000.surface.grb2
>
> * produced by alt_g2ctl v1.0.8, use alt_gmp to make idx file
>
> * command line options: -short
> IDY25001.APS3.syncld.slv.2023011218.000.surface.grb2
>
> * alt_gmp options: update=0
>
> * alt_gmp options: nthreads=4
>
> * alt_gmp options: big=0
>
> * alt_gmp options: match=
>
> * wgrib2 inventory flags: -npts -set_ext_name 1 -end_FT -ext_name -lev
>
> * wgrib2 inv suffix: .invd01
>
> * griddef=1:0:(427 x 512):grid_template=0:winds(N/S): lat-lon grid:(427 x
> 512) units 1e-06 input WE:SN output WE:SN res 48 lat -54.902344 to 4.980469
> by 0.117187 lon 95.009765 to 169.892578 by 0.175781
> #points=218624:winds(N/S)
>
>
>
> dtype grib2
>
> ydef 512 linear -54.902344 0.117187
>
> xdef 427 linear 95.009765 0.175781
>
> tdef 1 linear 18Z12jan2023 1mo
>
> zdef 1 levels 1
>
> vars 2
>
> v1 0 0 "BRTEMP.GMS_5_VISSR_(GMS-5)_10.79_um_:no_level"
>
> v2 0 0 "BRTEMP.MTSAT-2_IMAGER/MTSAT-2_6.80_um_:no_level"
>
> endvars
>
>
>
>
>
> This was without the -short flag:
>
> BRTEMPdGMS_5_VISSR_(GMS5)_10d79_um_ 0 0
> "BRTEMP.GMS_5_VISSR_(GMS-5)_10.79_um_:no_level" *
> BRTEMP.GMS_5_VISSR_(GMS-5)_10.79_um_:no_level
>
> BRTEMPdMTSAT2_IMAGER/MTSAT2_6d80_um_ 0 0
> "BRTEMP.MTSAT-2_IMAGER/MTSAT-2_6.80_um_:no_level" *
> BRTEMP.MTSAT-2_IMAGER/MTSAT-2_6.80_um_:no_level
>
>
>
>
>
>
>
> I have another question with ICON data now from DWD.  I’ve been following
> this article here -
> https://www.ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/_README.ICON.DWD and can
> get wgrib2 to display an inventory all good:
>
>
>
> wgrib2 icon.grb -if "^(1|2):" -grid_def -else -s -lon 95.00 0 -lon 180.00
> -55 -endif
>
> 1:0
>
> 2:5898409
>
> 3:11796818:d=2023011218:TPRATE:surface:0-3 hour acc
> fcst::lon=94.945312,lat=-0.002945,val=0:lon=180.000000,lat=-55.049820,val=0.0371094
>
>
>
> However when I try and run alt_g2ctl I get an undefined grid so was
> wondering what needs to be changed in alt_g2ctl to allow me to make a
> control file for ICON data.  I’d like to be able to use templating if
> possible as well:
>
>
>
> alt_g2ctl icon.grb
>
>
>
> dset ^icon.grb
>
> index ^icon.grb.idx
>
> undef 9.999E+20
>
> title icon.grb
>
> * produced by alt_g2ctl v1.0.8, use alt_gmp to make idx file
>
> * command line options: icon.grb
>
> * alt_gmp options: update=0
>
> * alt_gmp options: nthreads=4
>
> * alt_gmp options: big=0
>
> * alt_gmp options: match=
>
> * wgrib2 inventory flags: -npts -set_ext_name 1 -end_FT -ext_name -lev
>
> * wgrib2 inv suffix: .invd01
>
> * griddef=1:0:(2949120 x 1):grid_template=101: General Unstructured Grid
> grid=26 ref_grid=1 uuid=a27b8de6-18c4-11e4-820a-b5b098c6a5c0 #points=2949120
>
> :
>
>
>
> dtype grib2
>
> *** script needs to be modified ***
>
> *** unknown grid ***
>
>
>
> Cheers, Mike
>
>
>
> *From:* gradsusr <gradsusr-bounces at gradsusr.org> *On Behalf Of *Wesley
> Ebisuzaki - NOAA Federal
> *Sent:* Saturday, 7 January 2023 12:44 AM
> *To:* GrADS Users Forum <gradsusr at gradsusr.org>
> *Subject:* Re: [gradsusr] Issue with duplicate Varible in GRIB Data
>
>
>
> Mike,
>
>
>
> I would like to see if alt_g2clt can handle the  offending file.
>
> Can you send me a small version of the data?
>
>
>
> wgrib2 IN.grib -rpn 0 -grib_out ZERO.grb
>
>
>
> Wesley
>
>
>
> On Sat, Dec 10, 2022 at 1:44 AM <mike at weatherwatch.net.au> wrote:
>
> Hi all,
>
>
>
> I am currently dealing with an issue with GRIB2 data from The Australian
> BOM – They offer GRIB2 data for forecast satellite imagery.  The issue
> though is that the 2 different satellite images both use the same variable
> name of BRTEMP (see below):
>
>
>
> wgrib2 IDY25001.APS3.syncld.slv.2022120612.120.surface.grb2
>
> 1:0:d=2022120612:BRTEMP:no_level:120 hour fcst:GMS 5 VISSR (GMS-5) 10.79 um
>
> 2:171088:d=2022120612:BRTEMP:no_level:120 hour fcst:MTSAT-2 IMAGER/MTSAT-2
> 6.80 um
>
>
>
> How do I deal with this when creating a control file with g2ctl? Is this
> data unusable with GRADS because of this situation? If I just run g2ctl
> with no config, it produces a control file as follows:
>
> dtype grib2
>
> options template
>
> ydef 512 linear -54.902344 0.117187
>
> xdef 427 linear 95.009765 0.175781
>
> tdef 5 linear 00Z10dec2022 3hr
>
> zdef 1 linear 1 1
>
> vars 1
>
> BRTEMP  0,-1 0,5,7 ** no_level Brightness Temperature [K]
>
> ENDVARS
>
>
>
> Then when I try and display that data I get all undefined data.
>
>
>
> Cheers, Mike
>
> _______________________________________________
> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20230118/a01e40ad/attachment.html>


More information about the gradsusr mailing list