[gradsusr] Issue with duplicate Varible in GRIB Data

Fri Jan 13 02:45:06 EST 2023

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 <mailto: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 <mailto:gradsusr at gradsusr.org> 
http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20230113/a0e6d62e/attachment.html>