[Met_help] [rt.rap.ucar.edu #82980] History for grib_code rejection with multiple variables evaluation of point_stat mode
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:04:07 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Zhang, Mandy has shared OneDrive for Business files with you. To view them, click the links below.
<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1a926be162db74a49821d4185c60d4e88&authkey=AZj15sdHUqPB0g2JURrWLQM&expiration=2017-12-06T22%3a07%3a48.000Z>
[https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1a926be162db74a49821d4185c60d4e88&authkey=AZj15sdHUqPB0g2JURrWLQM&expiration=2017-12-06T22%3a07%3a48.000Z>
201711012000_height_met.nc<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1a926be162db74a49821d4185c60d4e88&authkey=AZj15sdHUqPB0g2JURrWLQM&expiration=2017-12-06T22%3a07%3a48.000Z>
<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1dddaebb76fdd400aa8fa20bf3e53e873&authkey=AellTpArmhPUzSiExMe7i74&expiration=2017-12-06T22%3a08%3a13.000Z>
[https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1dddaebb76fdd400aa8fa20bf3e53e873&authkey=AellTpArmhPUzSiExMe7i74&expiration=2017-12-06T22%3a08%3a13.000Z>
mettmp.nc<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1dddaebb76fdd400aa8fa20bf3e53e873&authkey=AellTpArmhPUzSiExMe7i74&expiration=2017-12-06T22%3a08%3a13.000Z>
<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1c2dc3eacce1343a7af093da549d6682c&authkey=AerQ14m-xvL6ApcS-E1BPWs&expiration=2017-12-06T22%3a08%3a38.000Z>
[https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1c2dc3eacce1343a7af093da549d6682c&authkey=AerQ14m-xvL6ApcS-E1BPWs&expiration=2017-12-06T22%3a08%3a38.000Z>
PointStatConfig_levs<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=1c2dc3eacce1343a7af093da549d6682c&authkey=AerQ14m-xvL6ApcS-E1BPWs&expiration=2017-12-06T22%3a08%3a38.000Z>
<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=17ef4e61c73654a96bee2e15b31cd9ded&authkey=Ae9wRURZblgucVgbuDf7wvI&expiration=2017-12-06T22%3a11%3a29.000Z>
[https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=17ef4e61c73654a96bee2e15b31cd9ded&authkey=Ae9wRURZblgucVgbuDf7wvI&expiration=2017-12-06T22%3a11%3a29.000Z>
run_point_stat(3).log<https://sooners-my.sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=17ef4e61c73654a96bee2e15b31cd9ded&authkey=Ae9wRURZblgucVgbuDf7wvI&expiration=2017-12-06T22%3a11%3a29.000Z>
Hi All,
I try to use mode point_stat to compare the forecast to observation data.
I'm not sure how to define the attribution information of multiple level forecast data.
my forecast data information is the belowing:
dimensions:
lev = 5 ;
lat = 363 ;
lon = 423 ;
variables:
float lev(lev) ;
lev:_FillValue = 9.96921e+36f ;
float lat(lat, lon) ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
float lon(lat, lon) ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
float TMP(lev, lat, lon) ;
TMP:_FillValue = -9999.f ;
TMP:_name = "temp" ;
TMP:long_name = "TMPC" ;
TMP:level = "p925" ;
TMP:units = "C" ;
TMP:init_time = "20171019_120000" ;
TMP:init_time_ut = "1508432400" ;
TMP:valid_time = "20171019_120000" ;
TMP:valid_time_ut = "1508432400" ;
TMP:accum_time = "010000" ;
TMP:accum_time_sec = 3600 ;
******
then I run point_stat, the error is:
*** Running Point-Stat on sample NAM data ***
DEBUG 1: Default Config File: /home1/04865/tg841649/tools/met/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: PointStatConfig_levs
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1
DEBUG 1: Forecast File: /scratch/04865/tg841649/met_predata/toZhangman/data2d/mettmp.nc
DEBUG 1: Observation File: 201711012000_height_met.nc
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMPP200.
ERROR :
ERROR : MetNcFile::data(NcVar *, const LongArray &, DataPlane &) -> needed 3 arguments for variable TMP, got 2
ERROR :
##################
I don't whether the error comes from the definition of forecast data.
Does anyone know what cause this issue?
please find my attach data and configure files!
Thank you very much!
Mandy
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: grib_code rejection with multiple variables evaluation of point_stat mode
From: John Halley Gotway
Time: Tue Nov 14 13:34:29 2017
Hi Mandy,
I see that you're wondering why you're getting 0 matched for when
evaluating the U and V components for wind. Thanks for sending your
sample
data files.
Using the data you sent, along with a sample forecast file, I was able
to
run met-6.0 point_stat and produce 5 matched pairs for TMP, UGRD, and
VGRD.
Looking at the log file from point_stat, I notice an important
difference.
Here's a selection from your log file:
DEBUG 2: Processing UGRDP925 versus UGRD/P925, for observation type
ADPUPA,
over region FULL, for interpolation method MEDIAN(9), using 0 pairs.
DEBUG 3: Number of matched pairs = 0
DEBUG 3: Observations processed = 85
DEBUG 3: Rejected: SID exclusion = 9
DEBUG 3: Rejected: GRIB code = 76
And here's a selection from mine:
DEBUG 2: Processing UGRD/P925 versus UGRD/P925, for observation type
ADPUPA, over region FULL, for interpolation method NEAREST(1), using 5
pairs.
DEBUG 3: Number of matched pairs = 5
DEBUG 3: Observations processed = 340
DEBUG 3: Rejected: SID exclusion = 36
DEBUG 3: Rejected: GRIB code = 228
The big question is, why are you processing 85 observations while I'm
processing 340 of them?
So let's step backwards... back to ascii2nc. I saved the point
observation
data you sent me to an ascii file named met.txt. Here's what I see
when I
run it through ascii2nc:
met-6.0/bin/ascii2nc met.txt met.nc
DEBUG 1: Config File:
/usr/local/met-6.0/share/met/config/Ascii2NcConfig_default
DEBUG 2: Rejected 0 observations off the masking grid.
DEBUG 2: Rejected 0 observations outside the masking polyline.
DEBUG 2: Rejected 0 observations not matched with station ID's.
DEBUG 2: Processing observations for 85 headers.
DEBUG 1: Creating NetCDF Observation file: met.nc
DEBUG 2: Finished processing 340 observations for 85 headers.
Notice that I'm getting 340 observations from 85 locations. What do
you
see when you run it?
Thanks,
John
On Fri, Nov 10, 2017 at 12:14 PM, Zhang, Mandy via RT
<met_help at ucar.edu>
wrote:
>
> Fri Nov 10 12:14:07 2017: Request 82980 was acted upon.
> Transaction: Ticket created by Man.Zhang-1 at ou.edu
> Queue: met_help
> Subject: grib_code rejection with multiple variables evaluation
of
> point_stat mode
> Owner: Nobody
> Requestors: Man.Zhang-1 at ou.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=82980 >
>
>
> Zhang, Mandy has shared OneDrive for Business files with you. To
view
> them, click the links below.
>
>
> <https://sooners-my.sharepoint.com/personal/man_
> zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=
> 1182a80e2655e47c4bcae5a6dff393fcf&authkey=AXGsbQHtrRqWMb3HW5Lt8Ag&
> expiration=2017-12-10T19%3a13%3a45.000Z>
> [https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<
> https://sooners-my.sharepoint.com/personal/man_zhang-
1_ou_edu/_layouts/15/
> guestaccess.aspx?docid=1182a80e2655e47c4bcae5a6dff393fcf&authkey=
> AXGsbQHtrRqWMb3HW5Lt8Ag&expiration=2017-12-10T19%3a13%3a45.000Z>
>
> PointStatConfig_levs(2)<https://sooners-my.sharepoint.com/
> personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=
> 1182a80e2655e47c4bcae5a6dff393fcf&authkey=AXGsbQHtrRqWMb3HW5Lt8Ag&
> expiration=2017-12-10T19%3a13%3a45.000Z>
>
> <https://sooners-my.sharepoint.com/personal/man_
> zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=
> 1f31fab5a4d5c45b4af71b8d8bc8c1530&authkey=ATrhucSN8Mb7ldFUgtP32uk&
> expiration=2017-12-10T19%3a13%3a50.000Z>
> [https://r1.res.office365.com/owa/prem/images/dc-generic_20.png]<
> https://sooners-my.sharepoint.com/personal/man_zhang-
1_ou_edu/_layouts/15/
> guestaccess.aspx?docid=1f31fab5a4d5c45b4af71b8d8bc8c1530&authkey=
> ATrhucSN8Mb7ldFUgtP32uk&expiration=2017-12-10T19%3a13%3a50.000Z>
>
> run_point_stat(5).log<https://sooners-my.sharepoint.com/
> personal/man_zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=
> 1f31fab5a4d5c45b4af71b8d8bc8c1530&authkey=ATrhucSN8Mb7ldFUgtP32uk&
> expiration=2017-12-10T19%3a13%3a50.000Z>
>
> <https://sooners-my.sharepoint.com/personal/man_
> zhang-1_ou_edu/_layouts/15/guestaccess.aspx?docid=
> 1e5b29d066b83488cb97375c833947b87&authkey=ASR5CawchNhWEYAivEtTIIs&
> expiration=2017-12-10T19%3a13%3a54.000Z>
> [https://r1.res.office365.com/owa/prem/images/dc-txt_20.png]<
> https://sooners-my.sharepoint.com/personal/man_zhang-
1_ou_edu/_layouts/15/
> guestaccess.aspx?docid=1e5b29d066b83488cb97375c833947b87&authkey=
> ASR5CawchNhWEYAivEtTIIs&expiration=2017-12-10T19%3a13%3a54.000Z>
>
> 201711012000_height925_met.txt<https://sooners-my.
> sharepoint.com/personal/man_zhang-1_ou_edu/_layouts/15/
> guestaccess.aspx?docid=1e5b29d066b83488cb97375c833947b87&authkey=
> ASR5CawchNhWEYAivEtTIIs&expiration=2017-12-10T19%3a13%3a54.000Z>
>
>
>
>
> Hi John,
>
> When I evaluate the single variable of temperature, the mode
> point_stat works well. But when I add other variables as ugrd and
vgrd,
> then run point_stat mode, the error of u and v is DEBUG 3: Rejected:
GRIB
> code = 76
>
> . Can you help me find any clue of what cause this issue?
>
> please find my observation data in ascii format and running script,
output
> information in attached files.
>
> Thank you very much!
>
> Mandy
>
>
>
------------------------------------------------
More information about the Met_help
mailing list