[gradsusr] Issue with duplicate Varible in GRIB Data

mike at weatherwatch.net.au mike at weatherwatch.net.au
Mon Jan 23 23:54:14 EST 2023


Thanks Wesley, I got that working with the DWD info you provided.

 

One more question I promise – I’m downloading the UVIndex GRIB data files from the NOAA site - https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230123/

 

When I try and view any data at all using grads, I get entire grid undefined.  Now I did a text dump of the data (wgrib2 gfs.uv.t12z.pgrb2.0p25.f014 -text tmp.pout) and when viewing data, I see the following for example:

720 361

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

0.0450439

 

I assume I have to *100 however grads doesn’t want a bar of it still which is quite odd:

 

/grads2/opengrads-2.2.1.oga.1/Contents/grads -blc

ga-> open gfs.uv.12z.ctl

Scanning description file: gfs.uv.12z.ctl

Data file gfs.uv.t12z.pgrb2.0p25.f%f3 is open as file 1

LON set to 0 360

LAT set to -90 90

LEV set to 1 1

Time values set: 2023:1:23:13 2023:1:23:13

E set to 1 1

ga-> q vars

uvisfc

description: ** surface Ultra Violet Index [W/m^2]

levels: 0

units: 0    7    196  -999 -999 -999 -999 -999 1    0    -999 -999 -999 -999 -999 -999

vecpair=-999

isu=0

offset=0

recoff=0

dfrm=0

var_t=0

 

ga-> d uvisfc*100

Cannot contour grid - all undefined values

ga-> quit

GX Package Terminated

 

Is there something that needs to be done differently here for the UV Index?  I have tried all different time steps and they all give all undefined errors. 

 

 

Cheers, Mike 

 

From: gradsusr <gradsusr-bounces at gradsusr.org> On Behalf Of Wesley Ebisuzaki - NOAA Federal
Sent: Thursday, 19 January 2023 7:02 AM
To: GrADS Users Forum <gradsusr at gradsusr.org>
Subject: Re: [gradsusr] Issue with duplicate Varible in GRIB Data

 

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 <mailto: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 <mailto: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 <mailto: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

_______________________________________________
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/20230124/4a2fd383/attachment-0001.html>


More information about the gradsusr mailing list