[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