[gradsusr] Issue with duplicate Varible in GRIB Data

Fri Jan 27 09:14:42 EST 2023

Mike,

I went to
https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/

and the only grib files that I saw were like

uv.t12z.grbf01.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf01.grib2>
    25-Jan-2023 17:35   92K  uv.t12z.grbf02.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf02.grib2>
    25-Jan-2023 17:35   93K  uv.t12z.grbf03.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf03.grib2>
    25-Jan-2023 17:35   93K  uv.t12z.grbf04.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf04.grib2>
    25-Jan-2023 17:34   92K  uv.t12z.grbf05.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf05.grib2>
    25-Jan-2023 17:34   92K  uv.t12z.grbf06.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf06.grib2>
    25-Jan-2023 17:34   93K  uv.t12z.grbf07.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf07.grib2>
    25-Jan-2023 17:35   93K  uv.t12z.grbf08.grib2
<https://ftpprd.ncep.noaa.gov/data/nccf/com/uvi/prod/uvi.20230125/uv.t12z.grbf08.grib2>
    25-Jan-2023 17:34   93K

I download one, and it was 720x360 lat-lon grid.

I had no problems displaying it after converting it from jpeg2000
to simple packing.  (The GrADS that comes with Ubuntu 20.04LTS
neglected include the jpeg2000 library.)

Wesley

On Tue, Jan 24, 2023 at 12:35 AM <mike at weatherwatch.net.au> wrote:

> 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> 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
>
> _______________________________________________
> 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/20230127/23c9f31e/attachment-0001.html>