[Met_help] [rt.rap.ucar.edu #56795] History for Question about PointStat problems
Paul Oldenburg via RT
met_help at ucar.edu
Tue Jul 10 08:28:27 MDT 2012
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello
a) I run the ncks utility to obtain selected forecast hours from my
WRF-ARW wrfout file (09Z to 15Z).
b) I run p_interp to obtain the pressure levels from the a) file and I got
wrfout_d01_2011-04-07_09:00:00_PLEV
c) I create an observation file with ADPUPA levels and convert it to netcdf
using the ascii2nc program.
d) I run PointStat
/u/charo/METv3.1/bin/point_stat
/u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV SalidaDeUnDia.nc
PointStatConfig_tutorial_run1 -v 2
and I got the following error message:
DEBUG 1: Default Config File:
/u/charo/METv3.1/data/config/PointStatConfig_default
DEBUG 1: User Config File: PointStatConfig_tutorial_run1
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=2605698639
DEBUG 1: Forecast File:
/u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: SalidaDeUnDia.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TT(0,1,*,*).
WARNING:
WARNING: MetNcPinterpDataFile::data_plane_array() -> no range flag found in
magic string "TT(0,1,*,*)"
WARNING:
ERROR :
ERROR : process_fcst_climo_files() -> no fields matching TT(0,1,*,*) found
in file: /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
e) I check the TT fields on the _PLEV file and they look OK. I send you the
Config file
Thank you very much for your help
Gisela
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #56795] Question about PointStat problems
From: Paul Oldenburg
Time: Tue Jun 05 14:21:48 2012
Gisela,
The syntax for specifying the level information in a p_interp field is
a little confusing. The level index in the
second position only accepts a range, specified like 1-1. The syntax
in general is: {lower}-{upper}, but in this case
both the lower and upper level is 1. Try the following setting in
your point_stat config file
PointStatConfig_tutorial_run1 to see if that will fix the problem:
fcst_field [] = [ "TT(0,1-1,*,*)" ];
Sorry for the confusion. Please let me know if you have any
questions.
Thanks,
Paul
On 06/05/2012 02:08 PM, Gisela Daniela Charó via RT wrote:
>
> Tue Jun 05 14:08:19 2012: Request 56795 was acted upon.
> Transaction: Ticket created by charogiselad at gmail.com
> Queue: met_help
> Subject: Question about PointStat problems
> Owner: Nobody
> Requestors: charogiselad at gmail.com
> Status: new
> Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795>
>
>
> Hello
>
> a) I run the ncks utility to obtain selected forecast hours from my
> WRF-ARW wrfout file (09Z to 15Z).
>
> b) I run p_interp to obtain the pressure levels from the a) file
and I got
> wrfout_d01_2011-04-07_09:00:00_PLEV
>
> c) I create an observation file with ADPUPA levels and convert it to
netcdf
> using the ascii2nc program.
>
> d) I run PointStat
> /u/charo/METv3.1/bin/point_stat
> /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
SalidaDeUnDia.nc
> PointStatConfig_tutorial_run1 -v 2
>
> and I got the following error message:
> DEBUG 1: Default Config File:
> /u/charo/METv3.1/data/config/PointStatConfig_default
> DEBUG 1: User Config File: PointStatConfig_tutorial_run1
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=2605698639
> DEBUG 1: Forecast File:
> /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File: SalidaDeUnDia.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for TT(0,1,*,*).
> WARNING:
> WARNING: MetNcPinterpDataFile::data_plane_array() -> no range flag
found in
> magic string "TT(0,1,*,*)"
> WARNING:
> ERROR :
> ERROR : process_fcst_climo_files() -> no fields matching
TT(0,1,*,*) found
> in file: /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
>
> e) I check the TT fields on the _PLEV file and they look OK. I send
you the
> Config file
>
>
> Thank you very much for your help
>
> Gisela
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Tue Jun 05 16:55:56 2012
Paul,
I will do that.If I have any doubts I will write to you again.
Thank you very much for your help.
Gisela
2012/6/5 Paul Oldenburg via RT <met_help at ucar.edu>
> Gisela,
>
> The syntax for specifying the level information in a p_interp field
is a
> little confusing. The level index in the
> second position only accepts a range, specified like 1-1. The
syntax in
> general is: {lower}-{upper}, but in this case
> both the lower and upper level is 1. Try the following setting in
your
> point_stat config file
> PointStatConfig_tutorial_run1 to see if that will fix the problem:
>
> fcst_field [] = [ "TT(0,1-1,*,*)" ];
>
> Sorry for the confusion. Please let me know if you have any
questions.
>
> Thanks,
>
> Paul
>
>
> On 06/05/2012 02:08 PM, Gisela Daniela Charó via RT wrote:
> >
> > Tue Jun 05 14:08:19 2012: Request 56795 was acted upon.
> > Transaction: Ticket created by charogiselad at gmail.com
> > Queue: met_help
> > Subject: Question about PointStat problems
> > Owner: Nobody
> > Requestors: charogiselad at gmail.com
> > Status: new
> > Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795>
> >
> >
> > Hello
> >
> > a) I run the ncks utility to obtain selected forecast hours from
my
> > WRF-ARW wrfout file (09Z to 15Z).
> >
> > b) I run p_interp to obtain the pressure levels from the a) file
and I
> got
> > wrfout_d01_2011-04-07_09:00:00_PLEV
> >
> > c) I create an observation file with ADPUPA levels and convert it
to
> netcdf
> > using the ascii2nc program.
> >
> > d) I run PointStat
> > /u/charo/METv3.1/bin/point_stat
> > /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
SalidaDeUnDia.nc
> > PointStatConfig_tutorial_run1 -v 2
> >
> > and I got the following error message:
> > DEBUG 1: Default Config File:
> > /u/charo/METv3.1/data/config/PointStatConfig_default
> > DEBUG 1: User Config File: PointStatConfig_tutorial_run1
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=2605698639
> > DEBUG 1: Forecast File:
> > /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
> > DEBUG 1: Climatology File: none
> > DEBUG 1: Observation File: SalidaDeUnDia.nc
> > DEBUG 2:
> > DEBUG 2:
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for TT(0,1,*,*).
> > WARNING:
> > WARNING: MetNcPinterpDataFile::data_plane_array() -> no range
flag
> found in
> > magic string "TT(0,1,*,*)"
> > WARNING:
> > ERROR :
> > ERROR : process_fcst_climo_files() -> no fields matching
TT(0,1,*,*)
> found
> > in file: /u/charo/P_INTERP/wrfout_d01_2011-04-07_09:00:00_PLEV
> >
> > e) I check the TT fields on the _PLEV file and they look OK. I
send you
> the
> > Config file
> >
> >
> > Thank you very much for your help
> >
> > Gisela
>
>
>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Mon Jun 11 14:00:51 2012
Paul,
Following your advice, PointStat works, but there was no results, only
the
titles.
Searching for the errors, we look at the observation file in nc format
resulting from the ASCII2NC. We are using ADPUPA observations from
different Argentine Meteorological Stations. When we made ncdump -v
-obs_arr, we see differences compared with the Sesion2.pdf from the
online
TUTORIAl.
Tutorial:
obs_arr =
0, 7, 837, 1618, 1618,
1, 11, 837, 1618, 273.05,
2, 17, 837, 1618, 271.85,
3, 52, 837, 1618, 92,
4, 53, 837, 1618, 0.00417,
5, 7, 826, 1724, 1724,
6, 11, 826, 1724, 274.55,
7, 17, 826, 1724, 272.15,
8, 52, 826, 1724, 84,
9, 53, 826, 1724, 0.00432,
10, 7, 815.3, 1829, 1829,
11, 11, 815.3, 1829, 276.45,
12, 17, 815.3, 1829, 265.75,
13, 52, 815.3, 1829, 45,
14, 53, 815.3, 1829, 0.0027,
15, 7, 815, 1832, 1832,
16, 11, 815, 1832, 276.55,
17, 17, 815, 1832, 265.55,
18, 52, 815, 1832, 44,
19, 53, 815, 1832, 0.00266,
20, 7, 784.7, 2134, 2134,
21, 11, 784.7, 2134, 274.05,
22, 17, 784.7, 2134, 264.15,
23, 52, 784.7, 2134, 47
Our output:
obs_arr =
0, 7, 20, -9999, 26460,
0, 11, 20, -9999, 222.1,
0, 17, 20, -9999, -9999,
0, 31, 20, -9999, 80,
0, 32, 20, -9999, 13.3744,
0, 33, 20, -9999, 2.32244,
0, 34, 20, -9999, 13.1712,
0, 7, 30, -9999, 23840,
0, 11, 30, -9999, 218.1,
0, 17, 30, -9999, -9999,
0, 31, 30, -9999, 40,
0, 32, 30, -9999, 7.2016,
0, 33, 30, -9999, 5.51675,
0, 34, 30, -9999, 4.6291,
0, 7, 50, -9999, 20640,
0, 11, 50, -9999, 208.5,
0, 17, 50, -9999, -9999,
0, 31, 50, -9999, 40,
0, 32, 50, -9999, 6.1728,
0, 33, 50, -9999, 4.72864,
0, 34, 50, -9999, 3.9678,
0, 7, 70, -9999, 18610,
0, 11, 70, -9999, 199.1,
0, 17, 70, -9999, -9999,
0, 31, 70, -9999, 245,
0, 32, 70, -9999, 4.6296,
0, 33, 70, -9999, -1.95655,
0, 34, 70, -9999, -4.19584,
0, 7, 100, -9999, 16540,
0, 11, 100, -9999, 201.3,
0, 17, 100, -9999, -9999,
0, 31, 100, -9999, 255,
0, 32, 100, -9999, 14.9176,
0, 33, 100, -9999, -3.86096,
0, 34, 100, -9999, -14.4093,
0, 7, 150, -9999, 14140,
0, 11, 150, -9999, 204.1,
0, 17, 150, -9999, -9999,
0, 31, 150, -9999, 275,
0, 32, 150, -9999, 41.152,
0, 33, 150, -9999, 3.58663,
0, 34, 150, -9999, -40.9954,
0, 7, 200, -9999, 12350,
0, 11, 200, -9999, 217.9,
0, 17, 200, -9999, -9999,
0, 31, 200, -9999, 275,
0, 32, 200, -9999, 36.008,
0, 33, 200, -9999, 3.1383,
0, 34, 200, -9999, -35.871,
0, 7, 250, -9999, 10890,
0, 11, 250, -9999, 228.9,
0, 17, 250, -9999, -9999,
0, 31, 250, -9999, 280,
0, 32, 250, -9999, 29.3208,
0, 33, 250, -9999, 5.0915,
0, 34, 250, -9999, -28.8754,
0, 7, 300, -9999, 9640,
0, 11, 300, -9999, 238.3,
0, 17, 300, -9999, 214.3,
0, 31, 300, -9999, 290,
0, 32, 300, -9999, 22.6336,
0, 33, 300, -9999, 7.74115,
0, 34, 300, -9999, -21.2686,
0, 7, 400, -9999, 7570,
0, 11, 400, -9999, 253.1,
0, 17, 400, -9999, 229.1,
0, 31, 400, -9999, 280,
0, 32, 400, -9999, 17.4896,
0, 33, 400, -9999, 3.03704,
0, 34, 400, -9999, -17.2239,
0, 7, 500, -9999, 5880,
0, 11, 500, -9999, 265.7,
0, 17, 500, -9999, 257.7,
0, 31, 500, -9999, 305,
0, 32, 500, -9999, 5.144,
0, 33, 500, -9999, 2.95048,
0, 34, 500, -9999, -4.21372,
0, 7, 700, -9999, 3145,
0, 11, 700, -9999, 279.6,
0, 17, 700, -9999, 279.6,
0, 31, 700, -9999, 325,
0, 32, 700, -9999, 2.0576,
0, 33, 700, -9999, 1.68549,
0, 34, 700, -9999, -1.18019,
0, 7, 850, -9999, 1518,
0, 11, 850, -9999, 288.6,
0, 17, 850, -9999, 287.1,
0, 31, 850, -9999, 40,
0, 32, 850, -9999, 1.5432,
0, 33, 850, -9999, 1.18216,
0, 34, 850, -9999, 0.99195,
0, 7, 925, -9999, -9999,
0, 11, 925, -9999, -9999,
0, 17, 925, -9999, -9999,
0, 31, 925, -9999, -9999,
0, 32, 925, -9999, -9999,
0, 33, 925, -9999, -9999,
0, 34, 925, -9999, -9999,
0, 7, 1000, -9999, -9999,
0, 11, 1000, -9999, -9999,
0, 17, 1000, -9999, -9999,
0, 31, 1000, -9999, -9999,
0, 32, 1000, -9999, -9999,
0, 33, 1000, -9999, -9999,
0, 34, 1000, -9999, -9999,
Thanks for your help
Gisela
2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #56795] Resolved: Question about PointStat problems
From: Paul Oldenburg
Time: Mon Jun 11 14:06:04 2012
Gisela,
I think it would be easiest for you to send me the data that you are
trying to feed point_stat so that I can reproduce
your problem. Can you please send the following files?
* pinterp model_data file
* observation ASCII data file
* point_stat config file
If the data files are larger than about 2MB, please upload these files
to our FTP site using the instructions here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp. Usually,
it is easiest to create a gzipped tar file that
contains all the files and email it or put that on the FTP site.
Please let me know if you have any questions.
Thanks,
Paul
On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>
> Paul,
>
> Following your advice, PointStat works, but there was no results,
only the
> titles.
>
> Searching for the errors, we look at the observation file in nc
format
> resulting from the ASCII2NC. We are using ADPUPA observations from
> different Argentine Meteorological Stations. When we made ncdump -v
> -obs_arr, we see differences compared with the Sesion2.pdf from the
online
> TUTORIAl.
>
> Tutorial:
> obs_arr =
> 0, 7, 837, 1618, 1618,
> 1, 11, 837, 1618, 273.05,
> 2, 17, 837, 1618, 271.85,
> 3, 52, 837, 1618, 92,
> 4, 53, 837, 1618, 0.00417,
> 5, 7, 826, 1724, 1724,
> 6, 11, 826, 1724, 274.55,
> 7, 17, 826, 1724, 272.15,
> 8, 52, 826, 1724, 84,
> 9, 53, 826, 1724, 0.00432,
> 10, 7, 815.3, 1829, 1829,
> 11, 11, 815.3, 1829, 276.45,
> 12, 17, 815.3, 1829, 265.75,
> 13, 52, 815.3, 1829, 45,
> 14, 53, 815.3, 1829, 0.0027,
> 15, 7, 815, 1832, 1832,
> 16, 11, 815, 1832, 276.55,
> 17, 17, 815, 1832, 265.55,
> 18, 52, 815, 1832, 44,
> 19, 53, 815, 1832, 0.00266,
> 20, 7, 784.7, 2134, 2134,
> 21, 11, 784.7, 2134, 274.05,
> 22, 17, 784.7, 2134, 264.15,
> 23, 52, 784.7, 2134, 47
>
> Our output:
>
> obs_arr =
> 0, 7, 20, -9999, 26460,
> 0, 11, 20, -9999, 222.1,
> 0, 17, 20, -9999, -9999,
> 0, 31, 20, -9999, 80,
> 0, 32, 20, -9999, 13.3744,
> 0, 33, 20, -9999, 2.32244,
> 0, 34, 20, -9999, 13.1712,
> 0, 7, 30, -9999, 23840,
> 0, 11, 30, -9999, 218.1,
> 0, 17, 30, -9999, -9999,
> 0, 31, 30, -9999, 40,
> 0, 32, 30, -9999, 7.2016,
> 0, 33, 30, -9999, 5.51675,
> 0, 34, 30, -9999, 4.6291,
> 0, 7, 50, -9999, 20640,
> 0, 11, 50, -9999, 208.5,
> 0, 17, 50, -9999, -9999,
> 0, 31, 50, -9999, 40,
> 0, 32, 50, -9999, 6.1728,
> 0, 33, 50, -9999, 4.72864,
> 0, 34, 50, -9999, 3.9678,
> 0, 7, 70, -9999, 18610,
> 0, 11, 70, -9999, 199.1,
> 0, 17, 70, -9999, -9999,
> 0, 31, 70, -9999, 245,
> 0, 32, 70, -9999, 4.6296,
> 0, 33, 70, -9999, -1.95655,
> 0, 34, 70, -9999, -4.19584,
> 0, 7, 100, -9999, 16540,
> 0, 11, 100, -9999, 201.3,
> 0, 17, 100, -9999, -9999,
> 0, 31, 100, -9999, 255,
> 0, 32, 100, -9999, 14.9176,
> 0, 33, 100, -9999, -3.86096,
> 0, 34, 100, -9999, -14.4093,
> 0, 7, 150, -9999, 14140,
> 0, 11, 150, -9999, 204.1,
> 0, 17, 150, -9999, -9999,
> 0, 31, 150, -9999, 275,
> 0, 32, 150, -9999, 41.152,
> 0, 33, 150, -9999, 3.58663,
> 0, 34, 150, -9999, -40.9954,
> 0, 7, 200, -9999, 12350,
> 0, 11, 200, -9999, 217.9,
> 0, 17, 200, -9999, -9999,
> 0, 31, 200, -9999, 275,
> 0, 32, 200, -9999, 36.008,
> 0, 33, 200, -9999, 3.1383,
> 0, 34, 200, -9999, -35.871,
> 0, 7, 250, -9999, 10890,
> 0, 11, 250, -9999, 228.9,
> 0, 17, 250, -9999, -9999,
> 0, 31, 250, -9999, 280,
> 0, 32, 250, -9999, 29.3208,
> 0, 33, 250, -9999, 5.0915,
> 0, 34, 250, -9999, -28.8754,
> 0, 7, 300, -9999, 9640,
> 0, 11, 300, -9999, 238.3,
> 0, 17, 300, -9999, 214.3,
> 0, 31, 300, -9999, 290,
> 0, 32, 300, -9999, 22.6336,
> 0, 33, 300, -9999, 7.74115,
> 0, 34, 300, -9999, -21.2686,
> 0, 7, 400, -9999, 7570,
> 0, 11, 400, -9999, 253.1,
> 0, 17, 400, -9999, 229.1,
> 0, 31, 400, -9999, 280,
> 0, 32, 400, -9999, 17.4896,
> 0, 33, 400, -9999, 3.03704,
> 0, 34, 400, -9999, -17.2239,
> 0, 7, 500, -9999, 5880,
> 0, 11, 500, -9999, 265.7,
> 0, 17, 500, -9999, 257.7,
> 0, 31, 500, -9999, 305,
> 0, 32, 500, -9999, 5.144,
> 0, 33, 500, -9999, 2.95048,
> 0, 34, 500, -9999, -4.21372,
> 0, 7, 700, -9999, 3145,
> 0, 11, 700, -9999, 279.6,
> 0, 17, 700, -9999, 279.6,
> 0, 31, 700, -9999, 325,
> 0, 32, 700, -9999, 2.0576,
> 0, 33, 700, -9999, 1.68549,
> 0, 34, 700, -9999, -1.18019,
> 0, 7, 850, -9999, 1518,
> 0, 11, 850, -9999, 288.6,
> 0, 17, 850, -9999, 287.1,
> 0, 31, 850, -9999, 40,
> 0, 32, 850, -9999, 1.5432,
> 0, 33, 850, -9999, 1.18216,
> 0, 34, 850, -9999, 0.99195,
> 0, 7, 925, -9999, -9999,
> 0, 11, 925, -9999, -9999,
> 0, 17, 925, -9999, -9999,
> 0, 31, 925, -9999, -9999,
> 0, 32, 925, -9999, -9999,
> 0, 33, 925, -9999, -9999,
> 0, 34, 925, -9999, -9999,
> 0, 7, 1000, -9999, -9999,
> 0, 11, 1000, -9999, -9999,
> 0, 17, 1000, -9999, -9999,
> 0, 31, 1000, -9999, -9999,
> 0, 32, 1000, -9999, -9999,
> 0, 33, 1000, -9999, -9999,
> 0, 34, 1000, -9999, -9999,
>
> Thanks for your help
>
> Gisela
>
>
>
>
>
> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> According to our records, your request has been resolved. If you
have any
>> further questions or concerns, please respond to this message.
>>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Mon Jun 11 21:14:44 2012
Paul,
I put the files on the FTP site.If you need another file, please tell
me.
Thanks for your help!!
Gisela
2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> Gisela,
>
> I think it would be easiest for you to send me the data that you are
> trying to feed point_stat so that I can reproduce
> your problem. Can you please send the following files?
>
> * pinterp model_data file
> * observation ASCII data file
> * point_stat config file
>
> If the data files are larger than about 2MB, please upload these
files to
> our FTP site using the instructions here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually, it
> is easiest to create a gzipped tar file that
> contains all the files and email it or put that on the FTP site.
Please
> let me know if you have any questions.
>
> Thanks,
>
> Paul
>
>
> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >
> > Paul,
> >
> > Following your advice, PointStat works, but there was no results,
only
> the
> > titles.
> >
> > Searching for the errors, we look at the observation file in nc
format
> > resulting from the ASCII2NC. We are using ADPUPA observations from
> > different Argentine Meteorological Stations. When we made ncdump
-v
> > -obs_arr, we see differences compared with the Sesion2.pdf from
the
> online
> > TUTORIAl.
> >
> > Tutorial:
> > obs_arr =
> > 0, 7, 837, 1618, 1618,
> > 1, 11, 837, 1618, 273.05,
> > 2, 17, 837, 1618, 271.85,
> > 3, 52, 837, 1618, 92,
> > 4, 53, 837, 1618, 0.00417,
> > 5, 7, 826, 1724, 1724,
> > 6, 11, 826, 1724, 274.55,
> > 7, 17, 826, 1724, 272.15,
> > 8, 52, 826, 1724, 84,
> > 9, 53, 826, 1724, 0.00432,
> > 10, 7, 815.3, 1829, 1829,
> > 11, 11, 815.3, 1829, 276.45,
> > 12, 17, 815.3, 1829, 265.75,
> > 13, 52, 815.3, 1829, 45,
> > 14, 53, 815.3, 1829, 0.0027,
> > 15, 7, 815, 1832, 1832,
> > 16, 11, 815, 1832, 276.55,
> > 17, 17, 815, 1832, 265.55,
> > 18, 52, 815, 1832, 44,
> > 19, 53, 815, 1832, 0.00266,
> > 20, 7, 784.7, 2134, 2134,
> > 21, 11, 784.7, 2134, 274.05,
> > 22, 17, 784.7, 2134, 264.15,
> > 23, 52, 784.7, 2134, 47
> >
> > Our output:
> >
> > obs_arr =
> > 0, 7, 20, -9999, 26460,
> > 0, 11, 20, -9999, 222.1,
> > 0, 17, 20, -9999, -9999,
> > 0, 31, 20, -9999, 80,
> > 0, 32, 20, -9999, 13.3744,
> > 0, 33, 20, -9999, 2.32244,
> > 0, 34, 20, -9999, 13.1712,
> > 0, 7, 30, -9999, 23840,
> > 0, 11, 30, -9999, 218.1,
> > 0, 17, 30, -9999, -9999,
> > 0, 31, 30, -9999, 40,
> > 0, 32, 30, -9999, 7.2016,
> > 0, 33, 30, -9999, 5.51675,
> > 0, 34, 30, -9999, 4.6291,
> > 0, 7, 50, -9999, 20640,
> > 0, 11, 50, -9999, 208.5,
> > 0, 17, 50, -9999, -9999,
> > 0, 31, 50, -9999, 40,
> > 0, 32, 50, -9999, 6.1728,
> > 0, 33, 50, -9999, 4.72864,
> > 0, 34, 50, -9999, 3.9678,
> > 0, 7, 70, -9999, 18610,
> > 0, 11, 70, -9999, 199.1,
> > 0, 17, 70, -9999, -9999,
> > 0, 31, 70, -9999, 245,
> > 0, 32, 70, -9999, 4.6296,
> > 0, 33, 70, -9999, -1.95655,
> > 0, 34, 70, -9999, -4.19584,
> > 0, 7, 100, -9999, 16540,
> > 0, 11, 100, -9999, 201.3,
> > 0, 17, 100, -9999, -9999,
> > 0, 31, 100, -9999, 255,
> > 0, 32, 100, -9999, 14.9176,
> > 0, 33, 100, -9999, -3.86096,
> > 0, 34, 100, -9999, -14.4093,
> > 0, 7, 150, -9999, 14140,
> > 0, 11, 150, -9999, 204.1,
> > 0, 17, 150, -9999, -9999,
> > 0, 31, 150, -9999, 275,
> > 0, 32, 150, -9999, 41.152,
> > 0, 33, 150, -9999, 3.58663,
> > 0, 34, 150, -9999, -40.9954,
> > 0, 7, 200, -9999, 12350,
> > 0, 11, 200, -9999, 217.9,
> > 0, 17, 200, -9999, -9999,
> > 0, 31, 200, -9999, 275,
> > 0, 32, 200, -9999, 36.008,
> > 0, 33, 200, -9999, 3.1383,
> > 0, 34, 200, -9999, -35.871,
> > 0, 7, 250, -9999, 10890,
> > 0, 11, 250, -9999, 228.9,
> > 0, 17, 250, -9999, -9999,
> > 0, 31, 250, -9999, 280,
> > 0, 32, 250, -9999, 29.3208,
> > 0, 33, 250, -9999, 5.0915,
> > 0, 34, 250, -9999, -28.8754,
> > 0, 7, 300, -9999, 9640,
> > 0, 11, 300, -9999, 238.3,
> > 0, 17, 300, -9999, 214.3,
> > 0, 31, 300, -9999, 290,
> > 0, 32, 300, -9999, 22.6336,
> > 0, 33, 300, -9999, 7.74115,
> > 0, 34, 300, -9999, -21.2686,
> > 0, 7, 400, -9999, 7570,
> > 0, 11, 400, -9999, 253.1,
> > 0, 17, 400, -9999, 229.1,
> > 0, 31, 400, -9999, 280,
> > 0, 32, 400, -9999, 17.4896,
> > 0, 33, 400, -9999, 3.03704,
> > 0, 34, 400, -9999, -17.2239,
> > 0, 7, 500, -9999, 5880,
> > 0, 11, 500, -9999, 265.7,
> > 0, 17, 500, -9999, 257.7,
> > 0, 31, 500, -9999, 305,
> > 0, 32, 500, -9999, 5.144,
> > 0, 33, 500, -9999, 2.95048,
> > 0, 34, 500, -9999, -4.21372,
> > 0, 7, 700, -9999, 3145,
> > 0, 11, 700, -9999, 279.6,
> > 0, 17, 700, -9999, 279.6,
> > 0, 31, 700, -9999, 325,
> > 0, 32, 700, -9999, 2.0576,
> > 0, 33, 700, -9999, 1.68549,
> > 0, 34, 700, -9999, -1.18019,
> > 0, 7, 850, -9999, 1518,
> > 0, 11, 850, -9999, 288.6,
> > 0, 17, 850, -9999, 287.1,
> > 0, 31, 850, -9999, 40,
> > 0, 32, 850, -9999, 1.5432,
> > 0, 33, 850, -9999, 1.18216,
> > 0, 34, 850, -9999, 0.99195,
> > 0, 7, 925, -9999, -9999,
> > 0, 11, 925, -9999, -9999,
> > 0, 17, 925, -9999, -9999,
> > 0, 31, 925, -9999, -9999,
> > 0, 32, 925, -9999, -9999,
> > 0, 33, 925, -9999, -9999,
> > 0, 34, 925, -9999, -9999,
> > 0, 7, 1000, -9999, -9999,
> > 0, 11, 1000, -9999, -9999,
> > 0, 17, 1000, -9999, -9999,
> > 0, 31, 1000, -9999, -9999,
> > 0, 32, 1000, -9999, -9999,
> > 0, 33, 1000, -9999, -9999,
> > 0, 34, 1000, -9999, -9999,
> >
> > Thanks for your help
> >
> > Gisela
> >
> >
> >
> >
> >
> > 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >
> >> According to our records, your request has been resolved. If you
have
> any
> >> further questions or concerns, please respond to this message.
> >>
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #56795] Resolved: Question about PointStat problems
From: Paul Oldenburg
Time: Tue Jun 12 11:06:25 2012
Gisela,
If you use the verbosity option in point_stat, you will see the
following output when you run it:
DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for observation
type ADPUPA, over region FULL, for interpolation
method UW_MEAN(1), using 0 pairs.
DEBUG 3: Number of matched pairs = 0
DEBUG 3: Observations processed = 770
DEBUG 3: Rejected: GRIB code = 660
DEBUG 3: Rejected: valid time = 110
...
Using the utility ncdump with "-v hdr_vld", I see that your
observations are valid at 20110407_120000. I used the NCAR
Integrated Data Viewer (IDV -
http://www.unidata.ucar.edu/software/idv/) to view the model data and
saw that there are
three forecast times present in the temperature variable:
20110407_090000
20110407_120000
20110407_150000
So, I changed the following setting in your configuration file (using
a 1 in the first index, to denote the second time
in the list above):
fcst_field [] = [ "TT(1,1-1,*,*)" ];
With this change, I get 3 matched pairs:
DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for observation
type ADPUPA, over region FULL, for interpolation
method UW_MEAN(1), using 3 pairs.
DEBUG 3: Number of matched pairs = 3
DEBUG 3: Observations processed = 770
DEBUG 3: Rejected: GRIB code = 660
DEBUG 3: Rejected: valid time = 0
DEBUG 3: Rejected: bad obs value = 7
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 100
Another option, if you want to compare the 20110407_090000 forecast
against your observations is to increase the valid
time window using the beg_ds and end_ds settings. Please let me know
if you have any questions.
Paul
On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>
> Paul,
> I put the files on the FTP site.If you need another file, please
tell me.
> Thanks for your help!!
> Gisela
>
> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> Gisela,
>>
>> I think it would be easiest for you to send me the data that you
are
>> trying to feed point_stat so that I can reproduce
>> your problem. Can you please send the following files?
>>
>> * pinterp model_data file
>> * observation ASCII data file
>> * point_stat config file
>>
>> If the data files are larger than about 2MB, please upload these
files to
>> our FTP site using the instructions here:
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually, it
>> is easiest to create a gzipped tar file that
>> contains all the files and email it or put that on the FTP site.
Please
>> let me know if you have any questions.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>>>
>>> Paul,
>>>
>>> Following your advice, PointStat works, but there was no results,
only
>> the
>>> titles.
>>>
>>> Searching for the errors, we look at the observation file in nc
format
>>> resulting from the ASCII2NC. We are using ADPUPA observations from
>>> different Argentine Meteorological Stations. When we made ncdump
-v
>>> -obs_arr, we see differences compared with the Sesion2.pdf from
the
>> online
>>> TUTORIAl.
>>>
>>> Tutorial:
>>> obs_arr =
>>> 0, 7, 837, 1618, 1618,
>>> 1, 11, 837, 1618, 273.05,
>>> 2, 17, 837, 1618, 271.85,
>>> 3, 52, 837, 1618, 92,
>>> 4, 53, 837, 1618, 0.00417,
>>> 5, 7, 826, 1724, 1724,
>>> 6, 11, 826, 1724, 274.55,
>>> 7, 17, 826, 1724, 272.15,
>>> 8, 52, 826, 1724, 84,
>>> 9, 53, 826, 1724, 0.00432,
>>> 10, 7, 815.3, 1829, 1829,
>>> 11, 11, 815.3, 1829, 276.45,
>>> 12, 17, 815.3, 1829, 265.75,
>>> 13, 52, 815.3, 1829, 45,
>>> 14, 53, 815.3, 1829, 0.0027,
>>> 15, 7, 815, 1832, 1832,
>>> 16, 11, 815, 1832, 276.55,
>>> 17, 17, 815, 1832, 265.55,
>>> 18, 52, 815, 1832, 44,
>>> 19, 53, 815, 1832, 0.00266,
>>> 20, 7, 784.7, 2134, 2134,
>>> 21, 11, 784.7, 2134, 274.05,
>>> 22, 17, 784.7, 2134, 264.15,
>>> 23, 52, 784.7, 2134, 47
>>>
>>> Our output:
>>>
>>> obs_arr =
>>> 0, 7, 20, -9999, 26460,
>>> 0, 11, 20, -9999, 222.1,
>>> 0, 17, 20, -9999, -9999,
>>> 0, 31, 20, -9999, 80,
>>> 0, 32, 20, -9999, 13.3744,
>>> 0, 33, 20, -9999, 2.32244,
>>> 0, 34, 20, -9999, 13.1712,
>>> 0, 7, 30, -9999, 23840,
>>> 0, 11, 30, -9999, 218.1,
>>> 0, 17, 30, -9999, -9999,
>>> 0, 31, 30, -9999, 40,
>>> 0, 32, 30, -9999, 7.2016,
>>> 0, 33, 30, -9999, 5.51675,
>>> 0, 34, 30, -9999, 4.6291,
>>> 0, 7, 50, -9999, 20640,
>>> 0, 11, 50, -9999, 208.5,
>>> 0, 17, 50, -9999, -9999,
>>> 0, 31, 50, -9999, 40,
>>> 0, 32, 50, -9999, 6.1728,
>>> 0, 33, 50, -9999, 4.72864,
>>> 0, 34, 50, -9999, 3.9678,
>>> 0, 7, 70, -9999, 18610,
>>> 0, 11, 70, -9999, 199.1,
>>> 0, 17, 70, -9999, -9999,
>>> 0, 31, 70, -9999, 245,
>>> 0, 32, 70, -9999, 4.6296,
>>> 0, 33, 70, -9999, -1.95655,
>>> 0, 34, 70, -9999, -4.19584,
>>> 0, 7, 100, -9999, 16540,
>>> 0, 11, 100, -9999, 201.3,
>>> 0, 17, 100, -9999, -9999,
>>> 0, 31, 100, -9999, 255,
>>> 0, 32, 100, -9999, 14.9176,
>>> 0, 33, 100, -9999, -3.86096,
>>> 0, 34, 100, -9999, -14.4093,
>>> 0, 7, 150, -9999, 14140,
>>> 0, 11, 150, -9999, 204.1,
>>> 0, 17, 150, -9999, -9999,
>>> 0, 31, 150, -9999, 275,
>>> 0, 32, 150, -9999, 41.152,
>>> 0, 33, 150, -9999, 3.58663,
>>> 0, 34, 150, -9999, -40.9954,
>>> 0, 7, 200, -9999, 12350,
>>> 0, 11, 200, -9999, 217.9,
>>> 0, 17, 200, -9999, -9999,
>>> 0, 31, 200, -9999, 275,
>>> 0, 32, 200, -9999, 36.008,
>>> 0, 33, 200, -9999, 3.1383,
>>> 0, 34, 200, -9999, -35.871,
>>> 0, 7, 250, -9999, 10890,
>>> 0, 11, 250, -9999, 228.9,
>>> 0, 17, 250, -9999, -9999,
>>> 0, 31, 250, -9999, 280,
>>> 0, 32, 250, -9999, 29.3208,
>>> 0, 33, 250, -9999, 5.0915,
>>> 0, 34, 250, -9999, -28.8754,
>>> 0, 7, 300, -9999, 9640,
>>> 0, 11, 300, -9999, 238.3,
>>> 0, 17, 300, -9999, 214.3,
>>> 0, 31, 300, -9999, 290,
>>> 0, 32, 300, -9999, 22.6336,
>>> 0, 33, 300, -9999, 7.74115,
>>> 0, 34, 300, -9999, -21.2686,
>>> 0, 7, 400, -9999, 7570,
>>> 0, 11, 400, -9999, 253.1,
>>> 0, 17, 400, -9999, 229.1,
>>> 0, 31, 400, -9999, 280,
>>> 0, 32, 400, -9999, 17.4896,
>>> 0, 33, 400, -9999, 3.03704,
>>> 0, 34, 400, -9999, -17.2239,
>>> 0, 7, 500, -9999, 5880,
>>> 0, 11, 500, -9999, 265.7,
>>> 0, 17, 500, -9999, 257.7,
>>> 0, 31, 500, -9999, 305,
>>> 0, 32, 500, -9999, 5.144,
>>> 0, 33, 500, -9999, 2.95048,
>>> 0, 34, 500, -9999, -4.21372,
>>> 0, 7, 700, -9999, 3145,
>>> 0, 11, 700, -9999, 279.6,
>>> 0, 17, 700, -9999, 279.6,
>>> 0, 31, 700, -9999, 325,
>>> 0, 32, 700, -9999, 2.0576,
>>> 0, 33, 700, -9999, 1.68549,
>>> 0, 34, 700, -9999, -1.18019,
>>> 0, 7, 850, -9999, 1518,
>>> 0, 11, 850, -9999, 288.6,
>>> 0, 17, 850, -9999, 287.1,
>>> 0, 31, 850, -9999, 40,
>>> 0, 32, 850, -9999, 1.5432,
>>> 0, 33, 850, -9999, 1.18216,
>>> 0, 34, 850, -9999, 0.99195,
>>> 0, 7, 925, -9999, -9999,
>>> 0, 11, 925, -9999, -9999,
>>> 0, 17, 925, -9999, -9999,
>>> 0, 31, 925, -9999, -9999,
>>> 0, 32, 925, -9999, -9999,
>>> 0, 33, 925, -9999, -9999,
>>> 0, 34, 925, -9999, -9999,
>>> 0, 7, 1000, -9999, -9999,
>>> 0, 11, 1000, -9999, -9999,
>>> 0, 17, 1000, -9999, -9999,
>>> 0, 31, 1000, -9999, -9999,
>>> 0, 32, 1000, -9999, -9999,
>>> 0, 33, 1000, -9999, -9999,
>>> 0, 34, 1000, -9999, -9999,
>>>
>>> Thanks for your help
>>>
>>> Gisela
>>>
>>>
>>>
>>>
>>>
>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>>>
>>>> According to our records, your request has been resolved. If you
have
>> any
>>>> further questions or concerns, please respond to this message.
>>>>
>>
>>
>>
>>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Tue Jun 12 12:10:19 2012
Paul,
I tested and it works!
The idea was just to check 12 Z forecast.
I would like to know where I can find an updated version of the
documentation about the different MET modules.Because what I found
online
does not agree with the latest METv3.1 version.
Thank you very much for your help!!
Best regards,
Gisela
2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>
> Gisela,
>
> If you use the verbosity option in point_stat, you will see the
following
> output when you run it:
>
> DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for observation
type
> ADPUPA, over region FULL, for interpolation
> method UW_MEAN(1), using 0 pairs.
> DEBUG 3: Number of matched pairs = 0
> DEBUG 3: Observations processed = 770
> DEBUG 3: Rejected: GRIB code = 660
> DEBUG 3: Rejected: valid time = 110
> ...
>
> Using the utility ncdump with "-v hdr_vld", I see that your
observations
> are valid at 20110407_120000. I used the NCAR
> Integrated Data Viewer (IDV -
http://www.unidata.ucar.edu/software/idv/)
> to view the model data and saw that there are
> three forecast times present in the temperature variable:
>
> 20110407_090000
> 20110407_120000
> 20110407_150000
>
> So, I changed the following setting in your configuration file
(using a 1
> in the first index, to denote the second time
> in the list above):
>
> fcst_field [] = [ "TT(1,1-1,*,*)" ];
>
> With this change, I get 3 matched pairs:
>
> DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for observation
type
> ADPUPA, over region FULL, for interpolation
> method UW_MEAN(1), using 3 pairs.
> DEBUG 3: Number of matched pairs = 3
> DEBUG 3: Observations processed = 770
> DEBUG 3: Rejected: GRIB code = 660
> DEBUG 3: Rejected: valid time = 0
> DEBUG 3: Rejected: bad obs value = 7
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 100
>
> Another option, if you want to compare the 20110407_090000 forecast
> against your observations is to increase the valid
> time window using the beg_ds and end_ds settings. Please let me
know if
> you have any questions.
>
> Paul
>
>
>
> On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >
> > Paul,
> > I put the files on the FTP site.If you need another file, please
tell me.
> > Thanks for your help!!
> > Gisela
> >
> > 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >
> >> Gisela,
> >>
> >> I think it would be easiest for you to send me the data that you
are
> >> trying to feed point_stat so that I can reproduce
> >> your problem. Can you please send the following files?
> >>
> >> * pinterp model_data file
> >> * observation ASCII data file
> >> * point_stat config file
> >>
> >> If the data files are larger than about 2MB, please upload these
files
> to
> >> our FTP site using the instructions here:
> >> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually,
> it
> >> is easiest to create a gzipped tar file that
> >> contains all the files and email it or put that on the FTP site.
Please
> >> let me know if you have any questions.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >>>
> >>> Paul,
> >>>
> >>> Following your advice, PointStat works, but there was no
results, only
> >> the
> >>> titles.
> >>>
> >>> Searching for the errors, we look at the observation file in nc
format
> >>> resulting from the ASCII2NC. We are using ADPUPA observations
from
> >>> different Argentine Meteorological Stations. When we made
ncdump -v
> >>> -obs_arr, we see differences compared with the Sesion2.pdf from
the
> >> online
> >>> TUTORIAl.
> >>>
> >>> Tutorial:
> >>> obs_arr =
> >>> 0, 7, 837, 1618, 1618,
> >>> 1, 11, 837, 1618, 273.05,
> >>> 2, 17, 837, 1618, 271.85,
> >>> 3, 52, 837, 1618, 92,
> >>> 4, 53, 837, 1618, 0.00417,
> >>> 5, 7, 826, 1724, 1724,
> >>> 6, 11, 826, 1724, 274.55,
> >>> 7, 17, 826, 1724, 272.15,
> >>> 8, 52, 826, 1724, 84,
> >>> 9, 53, 826, 1724, 0.00432,
> >>> 10, 7, 815.3, 1829, 1829,
> >>> 11, 11, 815.3, 1829, 276.45,
> >>> 12, 17, 815.3, 1829, 265.75,
> >>> 13, 52, 815.3, 1829, 45,
> >>> 14, 53, 815.3, 1829, 0.0027,
> >>> 15, 7, 815, 1832, 1832,
> >>> 16, 11, 815, 1832, 276.55,
> >>> 17, 17, 815, 1832, 265.55,
> >>> 18, 52, 815, 1832, 44,
> >>> 19, 53, 815, 1832, 0.00266,
> >>> 20, 7, 784.7, 2134, 2134,
> >>> 21, 11, 784.7, 2134, 274.05,
> >>> 22, 17, 784.7, 2134, 264.15,
> >>> 23, 52, 784.7, 2134, 47
> >>>
> >>> Our output:
> >>>
> >>> obs_arr =
> >>> 0, 7, 20, -9999, 26460,
> >>> 0, 11, 20, -9999, 222.1,
> >>> 0, 17, 20, -9999, -9999,
> >>> 0, 31, 20, -9999, 80,
> >>> 0, 32, 20, -9999, 13.3744,
> >>> 0, 33, 20, -9999, 2.32244,
> >>> 0, 34, 20, -9999, 13.1712,
> >>> 0, 7, 30, -9999, 23840,
> >>> 0, 11, 30, -9999, 218.1,
> >>> 0, 17, 30, -9999, -9999,
> >>> 0, 31, 30, -9999, 40,
> >>> 0, 32, 30, -9999, 7.2016,
> >>> 0, 33, 30, -9999, 5.51675,
> >>> 0, 34, 30, -9999, 4.6291,
> >>> 0, 7, 50, -9999, 20640,
> >>> 0, 11, 50, -9999, 208.5,
> >>> 0, 17, 50, -9999, -9999,
> >>> 0, 31, 50, -9999, 40,
> >>> 0, 32, 50, -9999, 6.1728,
> >>> 0, 33, 50, -9999, 4.72864,
> >>> 0, 34, 50, -9999, 3.9678,
> >>> 0, 7, 70, -9999, 18610,
> >>> 0, 11, 70, -9999, 199.1,
> >>> 0, 17, 70, -9999, -9999,
> >>> 0, 31, 70, -9999, 245,
> >>> 0, 32, 70, -9999, 4.6296,
> >>> 0, 33, 70, -9999, -1.95655,
> >>> 0, 34, 70, -9999, -4.19584,
> >>> 0, 7, 100, -9999, 16540,
> >>> 0, 11, 100, -9999, 201.3,
> >>> 0, 17, 100, -9999, -9999,
> >>> 0, 31, 100, -9999, 255,
> >>> 0, 32, 100, -9999, 14.9176,
> >>> 0, 33, 100, -9999, -3.86096,
> >>> 0, 34, 100, -9999, -14.4093,
> >>> 0, 7, 150, -9999, 14140,
> >>> 0, 11, 150, -9999, 204.1,
> >>> 0, 17, 150, -9999, -9999,
> >>> 0, 31, 150, -9999, 275,
> >>> 0, 32, 150, -9999, 41.152,
> >>> 0, 33, 150, -9999, 3.58663,
> >>> 0, 34, 150, -9999, -40.9954,
> >>> 0, 7, 200, -9999, 12350,
> >>> 0, 11, 200, -9999, 217.9,
> >>> 0, 17, 200, -9999, -9999,
> >>> 0, 31, 200, -9999, 275,
> >>> 0, 32, 200, -9999, 36.008,
> >>> 0, 33, 200, -9999, 3.1383,
> >>> 0, 34, 200, -9999, -35.871,
> >>> 0, 7, 250, -9999, 10890,
> >>> 0, 11, 250, -9999, 228.9,
> >>> 0, 17, 250, -9999, -9999,
> >>> 0, 31, 250, -9999, 280,
> >>> 0, 32, 250, -9999, 29.3208,
> >>> 0, 33, 250, -9999, 5.0915,
> >>> 0, 34, 250, -9999, -28.8754,
> >>> 0, 7, 300, -9999, 9640,
> >>> 0, 11, 300, -9999, 238.3,
> >>> 0, 17, 300, -9999, 214.3,
> >>> 0, 31, 300, -9999, 290,
> >>> 0, 32, 300, -9999, 22.6336,
> >>> 0, 33, 300, -9999, 7.74115,
> >>> 0, 34, 300, -9999, -21.2686,
> >>> 0, 7, 400, -9999, 7570,
> >>> 0, 11, 400, -9999, 253.1,
> >>> 0, 17, 400, -9999, 229.1,
> >>> 0, 31, 400, -9999, 280,
> >>> 0, 32, 400, -9999, 17.4896,
> >>> 0, 33, 400, -9999, 3.03704,
> >>> 0, 34, 400, -9999, -17.2239,
> >>> 0, 7, 500, -9999, 5880,
> >>> 0, 11, 500, -9999, 265.7,
> >>> 0, 17, 500, -9999, 257.7,
> >>> 0, 31, 500, -9999, 305,
> >>> 0, 32, 500, -9999, 5.144,
> >>> 0, 33, 500, -9999, 2.95048,
> >>> 0, 34, 500, -9999, -4.21372,
> >>> 0, 7, 700, -9999, 3145,
> >>> 0, 11, 700, -9999, 279.6,
> >>> 0, 17, 700, -9999, 279.6,
> >>> 0, 31, 700, -9999, 325,
> >>> 0, 32, 700, -9999, 2.0576,
> >>> 0, 33, 700, -9999, 1.68549,
> >>> 0, 34, 700, -9999, -1.18019,
> >>> 0, 7, 850, -9999, 1518,
> >>> 0, 11, 850, -9999, 288.6,
> >>> 0, 17, 850, -9999, 287.1,
> >>> 0, 31, 850, -9999, 40,
> >>> 0, 32, 850, -9999, 1.5432,
> >>> 0, 33, 850, -9999, 1.18216,
> >>> 0, 34, 850, -9999, 0.99195,
> >>> 0, 7, 925, -9999, -9999,
> >>> 0, 11, 925, -9999, -9999,
> >>> 0, 17, 925, -9999, -9999,
> >>> 0, 31, 925, -9999, -9999,
> >>> 0, 32, 925, -9999, -9999,
> >>> 0, 33, 925, -9999, -9999,
> >>> 0, 34, 925, -9999, -9999,
> >>> 0, 7, 1000, -9999, -9999,
> >>> 0, 11, 1000, -9999, -9999,
> >>> 0, 17, 1000, -9999, -9999,
> >>> 0, 31, 1000, -9999, -9999,
> >>> 0, 32, 1000, -9999, -9999,
> >>> 0, 33, 1000, -9999, -9999,
> >>> 0, 34, 1000, -9999, -9999,
> >>>
> >>> Thanks for your help
> >>>
> >>> Gisela
> >>>
> >>>
> >>>
> >>>
> >>>
> >>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>
> >>>> According to our records, your request has been resolved. If
you have
> >> any
> >>>> further questions or concerns, please respond to this message.
> >>>>
> >>
> >>
> >>
> >>
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #56795] Resolved: Question about PointStat problems
From: Paul Oldenburg
Time: Tue Jun 12 12:14:14 2012
Gisela,
The most up-to-date documentation that we have is posted on the MET
website here
http://www.dtcenter.org/met/users/docs/overview.php. If you find
discrepancies between the documentation posted there
and the corresponding instance of MET that you are using, please
report them to us and we will fix them accordingly.
Sorry for the trouble.
Paul
On 06/12/2012 12:10 PM, Gisela Daniela Charó via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>
> Paul,
> I tested and it works!
> The idea was just to check 12 Z forecast.
>
> I would like to know where I can find an updated version of the
> documentation about the different MET modules.Because what I found
online
> does not agree with the latest METv3.1 version.
> Thank you very much for your help!!
> Best regards,
> Gisela
>
> 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> Gisela,
>>
>> If you use the verbosity option in point_stat, you will see the
following
>> output when you run it:
>>
>> DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for observation
type
>> ADPUPA, over region FULL, for interpolation
>> method UW_MEAN(1), using 0 pairs.
>> DEBUG 3: Number of matched pairs = 0
>> DEBUG 3: Observations processed = 770
>> DEBUG 3: Rejected: GRIB code = 660
>> DEBUG 3: Rejected: valid time = 110
>> ...
>>
>> Using the utility ncdump with "-v hdr_vld", I see that your
observations
>> are valid at 20110407_120000. I used the NCAR
>> Integrated Data Viewer (IDV -
http://www.unidata.ucar.edu/software/idv/)
>> to view the model data and saw that there are
>> three forecast times present in the temperature variable:
>>
>> 20110407_090000
>> 20110407_120000
>> 20110407_150000
>>
>> So, I changed the following setting in your configuration file
(using a 1
>> in the first index, to denote the second time
>> in the list above):
>>
>> fcst_field [] = [ "TT(1,1-1,*,*)" ];
>>
>> With this change, I get 3 matched pairs:
>>
>> DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for observation
type
>> ADPUPA, over region FULL, for interpolation
>> method UW_MEAN(1), using 3 pairs.
>> DEBUG 3: Number of matched pairs = 3
>> DEBUG 3: Observations processed = 770
>> DEBUG 3: Rejected: GRIB code = 660
>> DEBUG 3: Rejected: valid time = 0
>> DEBUG 3: Rejected: bad obs value = 7
>> DEBUG 3: Rejected: off the grid = 0
>> DEBUG 3: Rejected: level mismatch = 100
>>
>> Another option, if you want to compare the 20110407_090000 forecast
>> against your observations is to increase the valid
>> time window using the beg_ds and end_ds settings. Please let me
know if
>> you have any questions.
>>
>> Paul
>>
>>
>>
>> On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>>>
>>> Paul,
>>> I put the files on the FTP site.If you need another file, please
tell me.
>>> Thanks for your help!!
>>> Gisela
>>>
>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>>>
>>>> Gisela,
>>>>
>>>> I think it would be easiest for you to send me the data that you
are
>>>> trying to feed point_stat so that I can reproduce
>>>> your problem. Can you please send the following files?
>>>>
>>>> * pinterp model_data file
>>>> * observation ASCII data file
>>>> * point_stat config file
>>>>
>>>> If the data files are larger than about 2MB, please upload these
files
>> to
>>>> our FTP site using the instructions here:
>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually,
>> it
>>>> is easiest to create a gzipped tar file that
>>>> contains all the files and email it or put that on the FTP site.
Please
>>>> let me know if you have any questions.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>>>>>
>>>>> Paul,
>>>>>
>>>>> Following your advice, PointStat works, but there was no
results, only
>>>> the
>>>>> titles.
>>>>>
>>>>> Searching for the errors, we look at the observation file in nc
format
>>>>> resulting from the ASCII2NC. We are using ADPUPA observations
from
>>>>> different Argentine Meteorological Stations. When we made
ncdump -v
>>>>> -obs_arr, we see differences compared with the Sesion2.pdf from
the
>>>> online
>>>>> TUTORIAl.
>>>>>
>>>>> Tutorial:
>>>>> obs_arr =
>>>>> 0, 7, 837, 1618, 1618,
>>>>> 1, 11, 837, 1618, 273.05,
>>>>> 2, 17, 837, 1618, 271.85,
>>>>> 3, 52, 837, 1618, 92,
>>>>> 4, 53, 837, 1618, 0.00417,
>>>>> 5, 7, 826, 1724, 1724,
>>>>> 6, 11, 826, 1724, 274.55,
>>>>> 7, 17, 826, 1724, 272.15,
>>>>> 8, 52, 826, 1724, 84,
>>>>> 9, 53, 826, 1724, 0.00432,
>>>>> 10, 7, 815.3, 1829, 1829,
>>>>> 11, 11, 815.3, 1829, 276.45,
>>>>> 12, 17, 815.3, 1829, 265.75,
>>>>> 13, 52, 815.3, 1829, 45,
>>>>> 14, 53, 815.3, 1829, 0.0027,
>>>>> 15, 7, 815, 1832, 1832,
>>>>> 16, 11, 815, 1832, 276.55,
>>>>> 17, 17, 815, 1832, 265.55,
>>>>> 18, 52, 815, 1832, 44,
>>>>> 19, 53, 815, 1832, 0.00266,
>>>>> 20, 7, 784.7, 2134, 2134,
>>>>> 21, 11, 784.7, 2134, 274.05,
>>>>> 22, 17, 784.7, 2134, 264.15,
>>>>> 23, 52, 784.7, 2134, 47
>>>>>
>>>>> Our output:
>>>>>
>>>>> obs_arr =
>>>>> 0, 7, 20, -9999, 26460,
>>>>> 0, 11, 20, -9999, 222.1,
>>>>> 0, 17, 20, -9999, -9999,
>>>>> 0, 31, 20, -9999, 80,
>>>>> 0, 32, 20, -9999, 13.3744,
>>>>> 0, 33, 20, -9999, 2.32244,
>>>>> 0, 34, 20, -9999, 13.1712,
>>>>> 0, 7, 30, -9999, 23840,
>>>>> 0, 11, 30, -9999, 218.1,
>>>>> 0, 17, 30, -9999, -9999,
>>>>> 0, 31, 30, -9999, 40,
>>>>> 0, 32, 30, -9999, 7.2016,
>>>>> 0, 33, 30, -9999, 5.51675,
>>>>> 0, 34, 30, -9999, 4.6291,
>>>>> 0, 7, 50, -9999, 20640,
>>>>> 0, 11, 50, -9999, 208.5,
>>>>> 0, 17, 50, -9999, -9999,
>>>>> 0, 31, 50, -9999, 40,
>>>>> 0, 32, 50, -9999, 6.1728,
>>>>> 0, 33, 50, -9999, 4.72864,
>>>>> 0, 34, 50, -9999, 3.9678,
>>>>> 0, 7, 70, -9999, 18610,
>>>>> 0, 11, 70, -9999, 199.1,
>>>>> 0, 17, 70, -9999, -9999,
>>>>> 0, 31, 70, -9999, 245,
>>>>> 0, 32, 70, -9999, 4.6296,
>>>>> 0, 33, 70, -9999, -1.95655,
>>>>> 0, 34, 70, -9999, -4.19584,
>>>>> 0, 7, 100, -9999, 16540,
>>>>> 0, 11, 100, -9999, 201.3,
>>>>> 0, 17, 100, -9999, -9999,
>>>>> 0, 31, 100, -9999, 255,
>>>>> 0, 32, 100, -9999, 14.9176,
>>>>> 0, 33, 100, -9999, -3.86096,
>>>>> 0, 34, 100, -9999, -14.4093,
>>>>> 0, 7, 150, -9999, 14140,
>>>>> 0, 11, 150, -9999, 204.1,
>>>>> 0, 17, 150, -9999, -9999,
>>>>> 0, 31, 150, -9999, 275,
>>>>> 0, 32, 150, -9999, 41.152,
>>>>> 0, 33, 150, -9999, 3.58663,
>>>>> 0, 34, 150, -9999, -40.9954,
>>>>> 0, 7, 200, -9999, 12350,
>>>>> 0, 11, 200, -9999, 217.9,
>>>>> 0, 17, 200, -9999, -9999,
>>>>> 0, 31, 200, -9999, 275,
>>>>> 0, 32, 200, -9999, 36.008,
>>>>> 0, 33, 200, -9999, 3.1383,
>>>>> 0, 34, 200, -9999, -35.871,
>>>>> 0, 7, 250, -9999, 10890,
>>>>> 0, 11, 250, -9999, 228.9,
>>>>> 0, 17, 250, -9999, -9999,
>>>>> 0, 31, 250, -9999, 280,
>>>>> 0, 32, 250, -9999, 29.3208,
>>>>> 0, 33, 250, -9999, 5.0915,
>>>>> 0, 34, 250, -9999, -28.8754,
>>>>> 0, 7, 300, -9999, 9640,
>>>>> 0, 11, 300, -9999, 238.3,
>>>>> 0, 17, 300, -9999, 214.3,
>>>>> 0, 31, 300, -9999, 290,
>>>>> 0, 32, 300, -9999, 22.6336,
>>>>> 0, 33, 300, -9999, 7.74115,
>>>>> 0, 34, 300, -9999, -21.2686,
>>>>> 0, 7, 400, -9999, 7570,
>>>>> 0, 11, 400, -9999, 253.1,
>>>>> 0, 17, 400, -9999, 229.1,
>>>>> 0, 31, 400, -9999, 280,
>>>>> 0, 32, 400, -9999, 17.4896,
>>>>> 0, 33, 400, -9999, 3.03704,
>>>>> 0, 34, 400, -9999, -17.2239,
>>>>> 0, 7, 500, -9999, 5880,
>>>>> 0, 11, 500, -9999, 265.7,
>>>>> 0, 17, 500, -9999, 257.7,
>>>>> 0, 31, 500, -9999, 305,
>>>>> 0, 32, 500, -9999, 5.144,
>>>>> 0, 33, 500, -9999, 2.95048,
>>>>> 0, 34, 500, -9999, -4.21372,
>>>>> 0, 7, 700, -9999, 3145,
>>>>> 0, 11, 700, -9999, 279.6,
>>>>> 0, 17, 700, -9999, 279.6,
>>>>> 0, 31, 700, -9999, 325,
>>>>> 0, 32, 700, -9999, 2.0576,
>>>>> 0, 33, 700, -9999, 1.68549,
>>>>> 0, 34, 700, -9999, -1.18019,
>>>>> 0, 7, 850, -9999, 1518,
>>>>> 0, 11, 850, -9999, 288.6,
>>>>> 0, 17, 850, -9999, 287.1,
>>>>> 0, 31, 850, -9999, 40,
>>>>> 0, 32, 850, -9999, 1.5432,
>>>>> 0, 33, 850, -9999, 1.18216,
>>>>> 0, 34, 850, -9999, 0.99195,
>>>>> 0, 7, 925, -9999, -9999,
>>>>> 0, 11, 925, -9999, -9999,
>>>>> 0, 17, 925, -9999, -9999,
>>>>> 0, 31, 925, -9999, -9999,
>>>>> 0, 32, 925, -9999, -9999,
>>>>> 0, 33, 925, -9999, -9999,
>>>>> 0, 34, 925, -9999, -9999,
>>>>> 0, 7, 1000, -9999, -9999,
>>>>> 0, 11, 1000, -9999, -9999,
>>>>> 0, 17, 1000, -9999, -9999,
>>>>> 0, 31, 1000, -9999, -9999,
>>>>> 0, 32, 1000, -9999, -9999,
>>>>> 0, 33, 1000, -9999, -9999,
>>>>> 0, 34, 1000, -9999, -9999,
>>>>>
>>>>> Thanks for your help
>>>>>
>>>>> Gisela
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>>>>>
>>>>>> According to our records, your request has been resolved. If
you have
>>>> any
>>>>>> further questions or concerns, please respond to this message.
>>>>>>
>>>>
>>>>
>>>>
>>>>
>>
>>
>>
>>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Tue Jun 26 11:43:18 2012
Hi Paul,
Finally I got the Poinstat working with soundings of pressure levels
1000
hPa throughout 70 hPa, but It didn't work for levels higher than 70
hPa.
Have you any suggestion about that?
Thanks in advanced,
Gisela
2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>t
> Gisela,
>
> The most up-to-date documentation that we have is posted on the MET
> website here
> http://www.dtcenter.org/met/users/docs/overview.php. If you find
> discrepancies between the documentation posted there
> and the corresponding instance of MET that you are using, please
report
> them to us and we will fix them accordingly.
> Sorry for the trouble.
>
> Paul
>
>
> On 06/12/2012 12:10 PM, Gisela Daniela Charó via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >
> > Paul,
> > I tested and it works!
> > The idea was just to check 12 Z forecast.
> >
> > I would like to know where I can find an updated version of the
> > documentation about the different MET modules.Because what I found
online
> > does not agree with the latest METv3.1 version.
> > Thank you very much for your help!!
> > Best regards,
> > Gisela
> >
> > 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>
> >
> >> Gisela,
> >>
> >> If you use the verbosity option in point_stat, you will see the
> following
> >> output when you run it:
> >>
> >> DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for
observation type
> >> ADPUPA, over region FULL, for interpolation
> >> method UW_MEAN(1), using 0 pairs.
> >> DEBUG 3: Number of matched pairs = 0
> >> DEBUG 3: Observations processed = 770
> >> DEBUG 3: Rejected: GRIB code = 660
> >> DEBUG 3: Rejected: valid time = 110
> >> ...
> >>
> >> Using the utility ncdump with "-v hdr_vld", I see that your
observations
> >> are valid at 20110407_120000. I used the NCAR
> >> Integrated Data Viewer (IDV -
http://www.unidata.ucar.edu/software/idv/
> )
> >> to view the model data and saw that there are
> >> three forecast times present in the temperature variable:
> >>
> >> 20110407_090000
> >> 20110407_120000
> >> 20110407_150000
> >>
> >> So, I changed the following setting in your configuration file
(using a
> 1
> >> in the first index, to denote the second time
> >> in the list above):
> >>
> >> fcst_field [] = [ "TT(1,1-1,*,*)" ];
> >>
> >> With this change, I get 3 matched pairs:
> >>
> >> DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for
observation type
> >> ADPUPA, over region FULL, for interpolation
> >> method UW_MEAN(1), using 3 pairs.
> >> DEBUG 3: Number of matched pairs = 3
> >> DEBUG 3: Observations processed = 770
> >> DEBUG 3: Rejected: GRIB code = 660
> >> DEBUG 3: Rejected: valid time = 0
> >> DEBUG 3: Rejected: bad obs value = 7
> >> DEBUG 3: Rejected: off the grid = 0
> >> DEBUG 3: Rejected: level mismatch = 100
> >>
> >> Another option, if you want to compare the 20110407_090000
forecast
> >> against your observations is to increase the valid
> >> time window using the beg_ds and end_ds settings. Please let me
know if
> >> you have any questions.
> >>
> >> Paul
> >>
> >>
> >>
> >> On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >>>
> >>> Paul,
> >>> I put the files on the FTP site.If you need another file, please
tell
> me.
> >>> Thanks for your help!!
> >>> Gisela
> >>>
> >>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>
> >>>> Gisela,
> >>>>
> >>>> I think it would be easiest for you to send me the data that
you are
> >>>> trying to feed point_stat so that I can reproduce
> >>>> your problem. Can you please send the following files?
> >>>>
> >>>> * pinterp model_data file
> >>>> * observation ASCII data file
> >>>> * point_stat config file
> >>>>
> >>>> If the data files are larger than about 2MB, please upload
these files
> >> to
> >>>> our FTP site using the instructions here:
> >>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually,
> >> it
> >>>> is easiest to create a gzipped tar file that
> >>>> contains all the files and email it or put that on the FTP
site.
> Please
> >>>> let me know if you have any questions.
> >>>>
> >>>> Thanks,
> >>>>
> >>>> Paul
> >>>>
> >>>>
> >>>> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795
>
> >>>>>
> >>>>> Paul,
> >>>>>
> >>>>> Following your advice, PointStat works, but there was no
results,
> only
> >>>> the
> >>>>> titles.
> >>>>>
> >>>>> Searching for the errors, we look at the observation file in
nc
> format
> >>>>> resulting from the ASCII2NC. We are using ADPUPA observations
from
> >>>>> different Argentine Meteorological Stations. When we made
ncdump -v
> >>>>> -obs_arr, we see differences compared with the Sesion2.pdf
from the
> >>>> online
> >>>>> TUTORIAl.
> >>>>>
> >>>>> Tutorial:
> >>>>> obs_arr =
> >>>>> 0, 7, 837, 1618, 1618,
> >>>>> 1, 11, 837, 1618, 273.05,
> >>>>> 2, 17, 837, 1618, 271.85,
> >>>>> 3, 52, 837, 1618, 92,
> >>>>> 4, 53, 837, 1618, 0.00417,
> >>>>> 5, 7, 826, 1724, 1724,
> >>>>> 6, 11, 826, 1724, 274.55,
> >>>>> 7, 17, 826, 1724, 272.15,
> >>>>> 8, 52, 826, 1724, 84,
> >>>>> 9, 53, 826, 1724, 0.00432,
> >>>>> 10, 7, 815.3, 1829, 1829,
> >>>>> 11, 11, 815.3, 1829, 276.45,
> >>>>> 12, 17, 815.3, 1829, 265.75,
> >>>>> 13, 52, 815.3, 1829, 45,
> >>>>> 14, 53, 815.3, 1829, 0.0027,
> >>>>> 15, 7, 815, 1832, 1832,
> >>>>> 16, 11, 815, 1832, 276.55,
> >>>>> 17, 17, 815, 1832, 265.55,
> >>>>> 18, 52, 815, 1832, 44,
> >>>>> 19, 53, 815, 1832, 0.00266,
> >>>>> 20, 7, 784.7, 2134, 2134,
> >>>>> 21, 11, 784.7, 2134, 274.05,
> >>>>> 22, 17, 784.7, 2134, 264.15,
> >>>>> 23, 52, 784.7, 2134, 47
> >>>>>
> >>>>> Our output:
> >>>>>
> >>>>> obs_arr =
> >>>>> 0, 7, 20, -9999, 26460,
> >>>>> 0, 11, 20, -9999, 222.1,
> >>>>> 0, 17, 20, -9999, -9999,
> >>>>> 0, 31, 20, -9999, 80,
> >>>>> 0, 32, 20, -9999, 13.3744,
> >>>>> 0, 33, 20, -9999, 2.32244,
> >>>>> 0, 34, 20, -9999, 13.1712,
> >>>>> 0, 7, 30, -9999, 23840,
> >>>>> 0, 11, 30, -9999, 218.1,
> >>>>> 0, 17, 30, -9999, -9999,
> >>>>> 0, 31, 30, -9999, 40,
> >>>>> 0, 32, 30, -9999, 7.2016,
> >>>>> 0, 33, 30, -9999, 5.51675,
> >>>>> 0, 34, 30, -9999, 4.6291,
> >>>>> 0, 7, 50, -9999, 20640,
> >>>>> 0, 11, 50, -9999, 208.5,
> >>>>> 0, 17, 50, -9999, -9999,
> >>>>> 0, 31, 50, -9999, 40,
> >>>>> 0, 32, 50, -9999, 6.1728,
> >>>>> 0, 33, 50, -9999, 4.72864,
> >>>>> 0, 34, 50, -9999, 3.9678,
> >>>>> 0, 7, 70, -9999, 18610,
> >>>>> 0, 11, 70, -9999, 199.1,
> >>>>> 0, 17, 70, -9999, -9999,
> >>>>> 0, 31, 70, -9999, 245,
> >>>>> 0, 32, 70, -9999, 4.6296,
> >>>>> 0, 33, 70, -9999, -1.95655,
> >>>>> 0, 34, 70, -9999, -4.19584,
> >>>>> 0, 7, 100, -9999, 16540,
> >>>>> 0, 11, 100, -9999, 201.3,
> >>>>> 0, 17, 100, -9999, -9999,
> >>>>> 0, 31, 100, -9999, 255,
> >>>>> 0, 32, 100, -9999, 14.9176,
> >>>>> 0, 33, 100, -9999, -3.86096,
> >>>>> 0, 34, 100, -9999, -14.4093,
> >>>>> 0, 7, 150, -9999, 14140,
> >>>>> 0, 11, 150, -9999, 204.1,
> >>>>> 0, 17, 150, -9999, -9999,
> >>>>> 0, 31, 150, -9999, 275,
> >>>>> 0, 32, 150, -9999, 41.152,
> >>>>> 0, 33, 150, -9999, 3.58663,
> >>>>> 0, 34, 150, -9999, -40.9954,
> >>>>> 0, 7, 200, -9999, 12350,
> >>>>> 0, 11, 200, -9999, 217.9,
> >>>>> 0, 17, 200, -9999, -9999,
> >>>>> 0, 31, 200, -9999, 275,
> >>>>> 0, 32, 200, -9999, 36.008,
> >>>>> 0, 33, 200, -9999, 3.1383,
> >>>>> 0, 34, 200, -9999, -35.871,
> >>>>> 0, 7, 250, -9999, 10890,
> >>>>> 0, 11, 250, -9999, 228.9,
> >>>>> 0, 17, 250, -9999, -9999,
> >>>>> 0, 31, 250, -9999, 280,
> >>>>> 0, 32, 250, -9999, 29.3208,
> >>>>> 0, 33, 250, -9999, 5.0915,
> >>>>> 0, 34, 250, -9999, -28.8754,
> >>>>> 0, 7, 300, -9999, 9640,
> >>>>> 0, 11, 300, -9999, 238.3,
> >>>>> 0, 17, 300, -9999, 214.3,
> >>>>> 0, 31, 300, -9999, 290,
> >>>>> 0, 32, 300, -9999, 22.6336,
> >>>>> 0, 33, 300, -9999, 7.74115,
> >>>>> 0, 34, 300, -9999, -21.2686,
> >>>>> 0, 7, 400, -9999, 7570,
> >>>>> 0, 11, 400, -9999, 253.1,
> >>>>> 0, 17, 400, -9999, 229.1,
> >>>>> 0, 31, 400, -9999, 280,
> >>>>> 0, 32, 400, -9999, 17.4896,
> >>>>> 0, 33, 400, -9999, 3.03704,
> >>>>> 0, 34, 400, -9999, -17.2239,
> >>>>> 0, 7, 500, -9999, 5880,
> >>>>> 0, 11, 500, -9999, 265.7,
> >>>>> 0, 17, 500, -9999, 257.7,
> >>>>> 0, 31, 500, -9999, 305,
> >>>>> 0, 32, 500, -9999, 5.144,
> >>>>> 0, 33, 500, -9999, 2.95048,
> >>>>> 0, 34, 500, -9999, -4.21372,
> >>>>> 0, 7, 700, -9999, 3145,
> >>>>> 0, 11, 700, -9999, 279.6,
> >>>>> 0, 17, 700, -9999, 279.6,
> >>>>> 0, 31, 700, -9999, 325,
> >>>>> 0, 32, 700, -9999, 2.0576,
> >>>>> 0, 33, 700, -9999, 1.68549,
> >>>>> 0, 34, 700, -9999, -1.18019,
> >>>>> 0, 7, 850, -9999, 1518,
> >>>>> 0, 11, 850, -9999, 288.6,
> >>>>> 0, 17, 850, -9999, 287.1,
> >>>>> 0, 31, 850, -9999, 40,
> >>>>> 0, 32, 850, -9999, 1.5432,
> >>>>> 0, 33, 850, -9999, 1.18216,
> >>>>> 0, 34, 850, -9999, 0.99195,
> >>>>> 0, 7, 925, -9999, -9999,
> >>>>> 0, 11, 925, -9999, -9999,
> >>>>> 0, 17, 925, -9999, -9999,
> >>>>> 0, 31, 925, -9999, -9999,
> >>>>> 0, 32, 925, -9999, -9999,
> >>>>> 0, 33, 925, -9999, -9999,
> >>>>> 0, 34, 925, -9999, -9999,
> >>>>> 0, 7, 1000, -9999, -9999,
> >>>>> 0, 11, 1000, -9999, -9999,
> >>>>> 0, 17, 1000, -9999, -9999,
> >>>>> 0, 31, 1000, -9999, -9999,
> >>>>> 0, 32, 1000, -9999, -9999,
> >>>>> 0, 33, 1000, -9999, -9999,
> >>>>> 0, 34, 1000, -9999, -9999,
> >>>>>
> >>>>> Thanks for your help
> >>>>>
> >>>>> Gisela
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>>>
> >>>>>> According to our records, your request has been resolved. If
you
> have
> >>>> any
> >>>>>> further questions or concerns, please respond to this
message.
> >>>>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>
> >>
> >>
> >>
>
>
>
>
------------------------------------------------
Subject: Question about PointStat problems
From: Paul Oldenburg
Time: Tue Jun 26 13:07:15 2012
Gisela,
I tried to run point_stat on the higher pressure levels for the data
that you sent me a couple weeks ago,
wrfout_d01_2011-04-07_09:00:00_PLEV. With the verbosity level set to
3 (-v 3), I saw the following output:
DEBUG 2: Processing TT(0,12-12,*,*) versus 11/P50, for observation
type ADPUPA, over region FULL, for interpolation
method UW_MEAN(1), using 0 pairs.
DEBUG 3: Number of matched pairs = 0
DEBUG 3: Observations processed = 770
DEBUG 3: Rejected: GRIB code = 660
DEBUG 3: Rejected: valid time = 0
DEBUG 3: Rejected: bad obs value = 7
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 96
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 7
I noticed, in particular, that there were 7 "bad fcst value" points
rejected. I plotted the data at both P70 and P50
(see attached files). It looks like the data at P50 is not valid. I
did not plot the data above P50, but I assume it
is bad as well. If MET notices that forecast or observation data
values are not "reasonable", meaning that they have
physically impossible values, it will not perform verification. That
is the case for these upper air temp fields.
Please let me know if you have any questions.
Thanks,
Paul
On 06/26/2012 11:43 AM, Gisela Daniela Charó via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>
> Hi Paul,
> Finally I got the Poinstat working with soundings of pressure levels
1000
> hPa throughout 70 hPa, but It didn't work for levels higher than
70 hPa.
> Have you any suggestion about that?
> Thanks in advanced,
> Gisela
>
> 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>t
>
>> Gisela,
>>
>> The most up-to-date documentation that we have is posted on the MET
>> website here
>> http://www.dtcenter.org/met/users/docs/overview.php. If you find
>> discrepancies between the documentation posted there
>> and the corresponding instance of MET that you are using, please
report
>> them to us and we will fix them accordingly.
>> Sorry for the trouble.
>>
>> Paul
>>
>>
>> On 06/12/2012 12:10 PM, Gisela Daniela Charó via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>>>
>>> Paul,
>>> I tested and it works!
>>> The idea was just to check 12 Z forecast.
>>>
>>> I would like to know where I can find an updated version of the
>>> documentation about the different MET modules.Because what I found
online
>>> does not agree with the latest METv3.1 version.
>>> Thank you very much for your help!!
>>> Best regards,
>>> Gisela
>>>
>>> 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>
>>>
>>>> Gisela,
>>>>
>>>> If you use the verbosity option in point_stat, you will see the
>> following
>>>> output when you run it:
>>>>
>>>> DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for
observation type
>>>> ADPUPA, over region FULL, for interpolation
>>>> method UW_MEAN(1), using 0 pairs.
>>>> DEBUG 3: Number of matched pairs = 0
>>>> DEBUG 3: Observations processed = 770
>>>> DEBUG 3: Rejected: GRIB code = 660
>>>> DEBUG 3: Rejected: valid time = 110
>>>> ...
>>>>
>>>> Using the utility ncdump with "-v hdr_vld", I see that your
observations
>>>> are valid at 20110407_120000. I used the NCAR
>>>> Integrated Data Viewer (IDV -
http://www.unidata.ucar.edu/software/idv/
>> )
>>>> to view the model data and saw that there are
>>>> three forecast times present in the temperature variable:
>>>>
>>>> 20110407_090000
>>>> 20110407_120000
>>>> 20110407_150000
>>>>
>>>> So, I changed the following setting in your configuration file
(using a
>> 1
>>>> in the first index, to denote the second time
>>>> in the list above):
>>>>
>>>> fcst_field [] = [ "TT(1,1-1,*,*)" ];
>>>>
>>>> With this change, I get 3 matched pairs:
>>>>
>>>> DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for
observation type
>>>> ADPUPA, over region FULL, for interpolation
>>>> method UW_MEAN(1), using 3 pairs.
>>>> DEBUG 3: Number of matched pairs = 3
>>>> DEBUG 3: Observations processed = 770
>>>> DEBUG 3: Rejected: GRIB code = 660
>>>> DEBUG 3: Rejected: valid time = 0
>>>> DEBUG 3: Rejected: bad obs value = 7
>>>> DEBUG 3: Rejected: off the grid = 0
>>>> DEBUG 3: Rejected: level mismatch = 100
>>>>
>>>> Another option, if you want to compare the 20110407_090000
forecast
>>>> against your observations is to increase the valid
>>>> time window using the beg_ds and end_ds settings. Please let me
know if
>>>> you have any questions.
>>>>
>>>> Paul
>>>>
>>>>
>>>>
>>>> On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>>>>>
>>>>> Paul,
>>>>> I put the files on the FTP site.If you need another file, please
tell
>> me.
>>>>> Thanks for your help!!
>>>>> Gisela
>>>>>
>>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>>>>>
>>>>>> Gisela,
>>>>>>
>>>>>> I think it would be easiest for you to send me the data that
you are
>>>>>> trying to feed point_stat so that I can reproduce
>>>>>> your problem. Can you please send the following files?
>>>>>>
>>>>>> * pinterp model_data file
>>>>>> * observation ASCII data file
>>>>>> * point_stat config file
>>>>>>
>>>>>> If the data files are larger than about 2MB, please upload
these files
>>>> to
>>>>>> our FTP site using the instructions here:
>>>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
Usually,
>>>> it
>>>>>> is easiest to create a gzipped tar file that
>>>>>> contains all the files and email it or put that on the FTP
site.
>> Please
>>>>>> let me know if you have any questions.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795
>
>>>>>>>
>>>>>>> Paul,
>>>>>>>
>>>>>>> Following your advice, PointStat works, but there was no
results,
>> only
>>>>>> the
>>>>>>> titles.
>>>>>>>
>>>>>>> Searching for the errors, we look at the observation file in
nc
>> format
>>>>>>> resulting from the ASCII2NC. We are using ADPUPA observations
from
>>>>>>> different Argentine Meteorological Stations. When we made
ncdump -v
>>>>>>> -obs_arr, we see differences compared with the Sesion2.pdf
from the
>>>>>> online
>>>>>>> TUTORIAl.
>>>>>>>
>>>>>>> Tutorial:
>>>>>>> obs_arr =
>>>>>>> 0, 7, 837, 1618, 1618,
>>>>>>> 1, 11, 837, 1618, 273.05,
>>>>>>> 2, 17, 837, 1618, 271.85,
>>>>>>> 3, 52, 837, 1618, 92,
>>>>>>> 4, 53, 837, 1618, 0.00417,
>>>>>>> 5, 7, 826, 1724, 1724,
>>>>>>> 6, 11, 826, 1724, 274.55,
>>>>>>> 7, 17, 826, 1724, 272.15,
>>>>>>> 8, 52, 826, 1724, 84,
>>>>>>> 9, 53, 826, 1724, 0.00432,
>>>>>>> 10, 7, 815.3, 1829, 1829,
>>>>>>> 11, 11, 815.3, 1829, 276.45,
>>>>>>> 12, 17, 815.3, 1829, 265.75,
>>>>>>> 13, 52, 815.3, 1829, 45,
>>>>>>> 14, 53, 815.3, 1829, 0.0027,
>>>>>>> 15, 7, 815, 1832, 1832,
>>>>>>> 16, 11, 815, 1832, 276.55,
>>>>>>> 17, 17, 815, 1832, 265.55,
>>>>>>> 18, 52, 815, 1832, 44,
>>>>>>> 19, 53, 815, 1832, 0.00266,
>>>>>>> 20, 7, 784.7, 2134, 2134,
>>>>>>> 21, 11, 784.7, 2134, 274.05,
>>>>>>> 22, 17, 784.7, 2134, 264.15,
>>>>>>> 23, 52, 784.7, 2134, 47
>>>>>>>
>>>>>>> Our output:
>>>>>>>
>>>>>>> obs_arr =
>>>>>>> 0, 7, 20, -9999, 26460,
>>>>>>> 0, 11, 20, -9999, 222.1,
>>>>>>> 0, 17, 20, -9999, -9999,
>>>>>>> 0, 31, 20, -9999, 80,
>>>>>>> 0, 32, 20, -9999, 13.3744,
>>>>>>> 0, 33, 20, -9999, 2.32244,
>>>>>>> 0, 34, 20, -9999, 13.1712,
>>>>>>> 0, 7, 30, -9999, 23840,
>>>>>>> 0, 11, 30, -9999, 218.1,
>>>>>>> 0, 17, 30, -9999, -9999,
>>>>>>> 0, 31, 30, -9999, 40,
>>>>>>> 0, 32, 30, -9999, 7.2016,
>>>>>>> 0, 33, 30, -9999, 5.51675,
>>>>>>> 0, 34, 30, -9999, 4.6291,
>>>>>>> 0, 7, 50, -9999, 20640,
>>>>>>> 0, 11, 50, -9999, 208.5,
>>>>>>> 0, 17, 50, -9999, -9999,
>>>>>>> 0, 31, 50, -9999, 40,
>>>>>>> 0, 32, 50, -9999, 6.1728,
>>>>>>> 0, 33, 50, -9999, 4.72864,
>>>>>>> 0, 34, 50, -9999, 3.9678,
>>>>>>> 0, 7, 70, -9999, 18610,
>>>>>>> 0, 11, 70, -9999, 199.1,
>>>>>>> 0, 17, 70, -9999, -9999,
>>>>>>> 0, 31, 70, -9999, 245,
>>>>>>> 0, 32, 70, -9999, 4.6296,
>>>>>>> 0, 33, 70, -9999, -1.95655,
>>>>>>> 0, 34, 70, -9999, -4.19584,
>>>>>>> 0, 7, 100, -9999, 16540,
>>>>>>> 0, 11, 100, -9999, 201.3,
>>>>>>> 0, 17, 100, -9999, -9999,
>>>>>>> 0, 31, 100, -9999, 255,
>>>>>>> 0, 32, 100, -9999, 14.9176,
>>>>>>> 0, 33, 100, -9999, -3.86096,
>>>>>>> 0, 34, 100, -9999, -14.4093,
>>>>>>> 0, 7, 150, -9999, 14140,
>>>>>>> 0, 11, 150, -9999, 204.1,
>>>>>>> 0, 17, 150, -9999, -9999,
>>>>>>> 0, 31, 150, -9999, 275,
>>>>>>> 0, 32, 150, -9999, 41.152,
>>>>>>> 0, 33, 150, -9999, 3.58663,
>>>>>>> 0, 34, 150, -9999, -40.9954,
>>>>>>> 0, 7, 200, -9999, 12350,
>>>>>>> 0, 11, 200, -9999, 217.9,
>>>>>>> 0, 17, 200, -9999, -9999,
>>>>>>> 0, 31, 200, -9999, 275,
>>>>>>> 0, 32, 200, -9999, 36.008,
>>>>>>> 0, 33, 200, -9999, 3.1383,
>>>>>>> 0, 34, 200, -9999, -35.871,
>>>>>>> 0, 7, 250, -9999, 10890,
>>>>>>> 0, 11, 250, -9999, 228.9,
>>>>>>> 0, 17, 250, -9999, -9999,
>>>>>>> 0, 31, 250, -9999, 280,
>>>>>>> 0, 32, 250, -9999, 29.3208,
>>>>>>> 0, 33, 250, -9999, 5.0915,
>>>>>>> 0, 34, 250, -9999, -28.8754,
>>>>>>> 0, 7, 300, -9999, 9640,
>>>>>>> 0, 11, 300, -9999, 238.3,
>>>>>>> 0, 17, 300, -9999, 214.3,
>>>>>>> 0, 31, 300, -9999, 290,
>>>>>>> 0, 32, 300, -9999, 22.6336,
>>>>>>> 0, 33, 300, -9999, 7.74115,
>>>>>>> 0, 34, 300, -9999, -21.2686,
>>>>>>> 0, 7, 400, -9999, 7570,
>>>>>>> 0, 11, 400, -9999, 253.1,
>>>>>>> 0, 17, 400, -9999, 229.1,
>>>>>>> 0, 31, 400, -9999, 280,
>>>>>>> 0, 32, 400, -9999, 17.4896,
>>>>>>> 0, 33, 400, -9999, 3.03704,
>>>>>>> 0, 34, 400, -9999, -17.2239,
>>>>>>> 0, 7, 500, -9999, 5880,
>>>>>>> 0, 11, 500, -9999, 265.7,
>>>>>>> 0, 17, 500, -9999, 257.7,
>>>>>>> 0, 31, 500, -9999, 305,
>>>>>>> 0, 32, 500, -9999, 5.144,
>>>>>>> 0, 33, 500, -9999, 2.95048,
>>>>>>> 0, 34, 500, -9999, -4.21372,
>>>>>>> 0, 7, 700, -9999, 3145,
>>>>>>> 0, 11, 700, -9999, 279.6,
>>>>>>> 0, 17, 700, -9999, 279.6,
>>>>>>> 0, 31, 700, -9999, 325,
>>>>>>> 0, 32, 700, -9999, 2.0576,
>>>>>>> 0, 33, 700, -9999, 1.68549,
>>>>>>> 0, 34, 700, -9999, -1.18019,
>>>>>>> 0, 7, 850, -9999, 1518,
>>>>>>> 0, 11, 850, -9999, 288.6,
>>>>>>> 0, 17, 850, -9999, 287.1,
>>>>>>> 0, 31, 850, -9999, 40,
>>>>>>> 0, 32, 850, -9999, 1.5432,
>>>>>>> 0, 33, 850, -9999, 1.18216,
>>>>>>> 0, 34, 850, -9999, 0.99195,
>>>>>>> 0, 7, 925, -9999, -9999,
>>>>>>> 0, 11, 925, -9999, -9999,
>>>>>>> 0, 17, 925, -9999, -9999,
>>>>>>> 0, 31, 925, -9999, -9999,
>>>>>>> 0, 32, 925, -9999, -9999,
>>>>>>> 0, 33, 925, -9999, -9999,
>>>>>>> 0, 34, 925, -9999, -9999,
>>>>>>> 0, 7, 1000, -9999, -9999,
>>>>>>> 0, 11, 1000, -9999, -9999,
>>>>>>> 0, 17, 1000, -9999, -9999,
>>>>>>> 0, 31, 1000, -9999, -9999,
>>>>>>> 0, 32, 1000, -9999, -9999,
>>>>>>> 0, 33, 1000, -9999, -9999,
>>>>>>> 0, 34, 1000, -9999, -9999,
>>>>>>>
>>>>>>> Thanks for your help
>>>>>>>
>>>>>>> Gisela
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
>>>>>>>
>>>>>>>> According to our records, your request has been resolved. If
you
>> have
>>>>>> any
>>>>>>>> further questions or concerns, please respond to this
message.
>>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>
>>>>
>>>>
>>>>
>>
>>
>>
>>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Tue Jun 26 14:10:47 2012
Paul,
I will check with my WRF configuration.
Thank you very much for your help,
Gisela
2012/6/26 Paul Oldenburg via RT <met_help at ucar.edu>
> Gisela,
>
> I tried to run point_stat on the higher pressure levels for the data
that
> you sent me a couple weeks ago,
> wrfout_d01_2011-04-07_09:00:00_PLEV. With the verbosity level set
to 3
> (-v 3), I saw the following output:
>
> DEBUG 2: Processing TT(0,12-12,*,*) versus 11/P50, for observation
type
> ADPUPA, over region FULL, for interpolation
> method UW_MEAN(1), using 0 pairs.
> DEBUG 3: Number of matched pairs = 0
> DEBUG 3: Observations processed = 770
> DEBUG 3: Rejected: GRIB code = 660
> DEBUG 3: Rejected: valid time = 0
> DEBUG 3: Rejected: bad obs value = 7
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 96
> DEBUG 3: Rejected: message type = 0
> DEBUG 3: Rejected: masking region = 0
> DEBUG 3: Rejected: bad fcst value = 7
>
> I noticed, in particular, that there were 7 "bad fcst value" points
> rejected. I plotted the data at both P70 and P50
> (see attached files). It looks like the data at P50 is not valid.
I did
> not plot the data above P50, but I assume it
> is bad as well. If MET notices that forecast or observation data
values
> are not "reasonable", meaning that they have
> physically impossible values, it will not perform verification.
That is
> the case for these upper air temp fields.
>
> Please let me know if you have any questions.
>
> Thanks,
>
> Paul
>
>
> On 06/26/2012 11:43 AM, Gisela Daniela Charó via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >
> > Hi Paul,
> > Finally I got the Poinstat working with soundings of pressure
levels 1000
> > hPa throughout 70 hPa, but It didn't work for levels higher than
70
> hPa.
> > Have you any suggestion about that?
> > Thanks in advanced,
> > Gisela
> >
> > 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>t
> >
> >> Gisela,
> >>
> >> The most up-to-date documentation that we have is posted on the
MET
> >> website here
> >> http://www.dtcenter.org/met/users/docs/overview.php. If you find
> >> discrepancies between the documentation posted there
> >> and the corresponding instance of MET that you are using, please
report
> >> them to us and we will fix them accordingly.
> >> Sorry for the trouble.
> >>
> >> Paul
> >>
> >>
> >> On 06/12/2012 12:10 PM, Gisela Daniela Charó via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >>>
> >>> Paul,
> >>> I tested and it works!
> >>> The idea was just to check 12 Z forecast.
> >>>
> >>> I would like to know where I can find an updated version of the
> >>> documentation about the different MET modules.Because what I
found
> online
> >>> does not agree with the latest METv3.1 version.
> >>> Thank you very much for your help!!
> >>> Best regards,
> >>> Gisela
> >>>
> >>> 2012/6/12 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>
> >>>> Gisela,
> >>>>
> >>>> If you use the verbosity option in point_stat, you will see the
> >> following
> >>>> output when you run it:
> >>>>
> >>>> DEBUG 2: Processing TT(0,1-1,*,*) versus 11/P1000, for
observation
> type
> >>>> ADPUPA, over region FULL, for interpolation
> >>>> method UW_MEAN(1), using 0 pairs.
> >>>> DEBUG 3: Number of matched pairs = 0
> >>>> DEBUG 3: Observations processed = 770
> >>>> DEBUG 3: Rejected: GRIB code = 660
> >>>> DEBUG 3: Rejected: valid time = 110
> >>>> ...
> >>>>
> >>>> Using the utility ncdump with "-v hdr_vld", I see that your
> observations
> >>>> are valid at 20110407_120000. I used the NCAR
> >>>> Integrated Data Viewer (IDV -
> http://www.unidata.ucar.edu/software/idv/
> >> )
> >>>> to view the model data and saw that there are
> >>>> three forecast times present in the temperature variable:
> >>>>
> >>>> 20110407_090000
> >>>> 20110407_120000
> >>>> 20110407_150000
> >>>>
> >>>> So, I changed the following setting in your configuration file
(using
> a
> >> 1
> >>>> in the first index, to denote the second time
> >>>> in the list above):
> >>>>
> >>>> fcst_field [] = [ "TT(1,1-1,*,*)" ];
> >>>>
> >>>> With this change, I get 3 matched pairs:
> >>>>
> >>>> DEBUG 2: Processing TT(1,1-1,*,*) versus 11/P1000, for
observation
> type
> >>>> ADPUPA, over region FULL, for interpolation
> >>>> method UW_MEAN(1), using 3 pairs.
> >>>> DEBUG 3: Number of matched pairs = 3
> >>>> DEBUG 3: Observations processed = 770
> >>>> DEBUG 3: Rejected: GRIB code = 660
> >>>> DEBUG 3: Rejected: valid time = 0
> >>>> DEBUG 3: Rejected: bad obs value = 7
> >>>> DEBUG 3: Rejected: off the grid = 0
> >>>> DEBUG 3: Rejected: level mismatch = 100
> >>>>
> >>>> Another option, if you want to compare the 20110407_090000
forecast
> >>>> against your observations is to increase the valid
> >>>> time window using the beg_ds and end_ds settings. Please let
me know
> if
> >>>> you have any questions.
> >>>>
> >>>> Paul
> >>>>
> >>>>
> >>>>
> >>>> On 06/11/2012 09:14 PM, Gisela Daniela Charó via RT wrote:
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795
>
> >>>>>
> >>>>> Paul,
> >>>>> I put the files on the FTP site.If you need another file,
please tell
> >> me.
> >>>>> Thanks for your help!!
> >>>>> Gisela
> >>>>>
> >>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>>>
> >>>>>> Gisela,
> >>>>>>
> >>>>>> I think it would be easiest for you to send me the data that
you are
> >>>>>> trying to feed point_stat so that I can reproduce
> >>>>>> your problem. Can you please send the following files?
> >>>>>>
> >>>>>> * pinterp model_data file
> >>>>>> * observation ASCII data file
> >>>>>> * point_stat config file
> >>>>>>
> >>>>>> If the data files are larger than about 2MB, please upload
these
> files
> >>>> to
> >>>>>> our FTP site using the instructions here:
> >>>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
> Usually,
> >>>> it
> >>>>>> is easiest to create a gzipped tar file that
> >>>>>> contains all the files and email it or put that on the FTP
site.
> >> Please
> >>>>>> let me know if you have any questions.
> >>>>>>
> >>>>>> Thanks,
> >>>>>>
> >>>>>> Paul
> >>>>>>
> >>>>>>
> >>>>>> On 06/11/2012 02:00 PM, Gisela Daniela Charó via RT wrote:
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >>>>>>>
> >>>>>>> Paul,
> >>>>>>>
> >>>>>>> Following your advice, PointStat works, but there was no
results,
> >> only
> >>>>>> the
> >>>>>>> titles.
> >>>>>>>
> >>>>>>> Searching for the errors, we look at the observation file in
nc
> >> format
> >>>>>>> resulting from the ASCII2NC. We are using ADPUPA
observations from
> >>>>>>> different Argentine Meteorological Stations. When we made
ncdump
> -v
> >>>>>>> -obs_arr, we see differences compared with the Sesion2.pdf
from
> the
> >>>>>> online
> >>>>>>> TUTORIAl.
> >>>>>>>
> >>>>>>> Tutorial:
> >>>>>>> obs_arr =
> >>>>>>> 0, 7, 837, 1618, 1618,
> >>>>>>> 1, 11, 837, 1618, 273.05,
> >>>>>>> 2, 17, 837, 1618, 271.85,
> >>>>>>> 3, 52, 837, 1618, 92,
> >>>>>>> 4, 53, 837, 1618, 0.00417,
> >>>>>>> 5, 7, 826, 1724, 1724,
> >>>>>>> 6, 11, 826, 1724, 274.55,
> >>>>>>> 7, 17, 826, 1724, 272.15,
> >>>>>>> 8, 52, 826, 1724, 84,
> >>>>>>> 9, 53, 826, 1724, 0.00432,
> >>>>>>> 10, 7, 815.3, 1829, 1829,
> >>>>>>> 11, 11, 815.3, 1829, 276.45,
> >>>>>>> 12, 17, 815.3, 1829, 265.75,
> >>>>>>> 13, 52, 815.3, 1829, 45,
> >>>>>>> 14, 53, 815.3, 1829, 0.0027,
> >>>>>>> 15, 7, 815, 1832, 1832,
> >>>>>>> 16, 11, 815, 1832, 276.55,
> >>>>>>> 17, 17, 815, 1832, 265.55,
> >>>>>>> 18, 52, 815, 1832, 44,
> >>>>>>> 19, 53, 815, 1832, 0.00266,
> >>>>>>> 20, 7, 784.7, 2134, 2134,
> >>>>>>> 21, 11, 784.7, 2134, 274.05,
> >>>>>>> 22, 17, 784.7, 2134, 264.15,
> >>>>>>> 23, 52, 784.7, 2134, 47
> >>>>>>>
> >>>>>>> Our output:
> >>>>>>>
> >>>>>>> obs_arr =
> >>>>>>> 0, 7, 20, -9999, 26460,
> >>>>>>> 0, 11, 20, -9999, 222.1,
> >>>>>>> 0, 17, 20, -9999, -9999,
> >>>>>>> 0, 31, 20, -9999, 80,
> >>>>>>> 0, 32, 20, -9999, 13.3744,
> >>>>>>> 0, 33, 20, -9999, 2.32244,
> >>>>>>> 0, 34, 20, -9999, 13.1712,
> >>>>>>> 0, 7, 30, -9999, 23840,
> >>>>>>> 0, 11, 30, -9999, 218.1,
> >>>>>>> 0, 17, 30, -9999, -9999,
> >>>>>>> 0, 31, 30, -9999, 40,
> >>>>>>> 0, 32, 30, -9999, 7.2016,
> >>>>>>> 0, 33, 30, -9999, 5.51675,
> >>>>>>> 0, 34, 30, -9999, 4.6291,
> >>>>>>> 0, 7, 50, -9999, 20640,
> >>>>>>> 0, 11, 50, -9999, 208.5,
> >>>>>>> 0, 17, 50, -9999, -9999,
> >>>>>>> 0, 31, 50, -9999, 40,
> >>>>>>> 0, 32, 50, -9999, 6.1728,
> >>>>>>> 0, 33, 50, -9999, 4.72864,
> >>>>>>> 0, 34, 50, -9999, 3.9678,
> >>>>>>> 0, 7, 70, -9999, 18610,
> >>>>>>> 0, 11, 70, -9999, 199.1,
> >>>>>>> 0, 17, 70, -9999, -9999,
> >>>>>>> 0, 31, 70, -9999, 245,
> >>>>>>> 0, 32, 70, -9999, 4.6296,
> >>>>>>> 0, 33, 70, -9999, -1.95655,
> >>>>>>> 0, 34, 70, -9999, -4.19584,
> >>>>>>> 0, 7, 100, -9999, 16540,
> >>>>>>> 0, 11, 100, -9999, 201.3,
> >>>>>>> 0, 17, 100, -9999, -9999,
> >>>>>>> 0, 31, 100, -9999, 255,
> >>>>>>> 0, 32, 100, -9999, 14.9176,
> >>>>>>> 0, 33, 100, -9999, -3.86096,
> >>>>>>> 0, 34, 100, -9999, -14.4093,
> >>>>>>> 0, 7, 150, -9999, 14140,
> >>>>>>> 0, 11, 150, -9999, 204.1,
> >>>>>>> 0, 17, 150, -9999, -9999,
> >>>>>>> 0, 31, 150, -9999, 275,
> >>>>>>> 0, 32, 150, -9999, 41.152,
> >>>>>>> 0, 33, 150, -9999, 3.58663,
> >>>>>>> 0, 34, 150, -9999, -40.9954,
> >>>>>>> 0, 7, 200, -9999, 12350,
> >>>>>>> 0, 11, 200, -9999, 217.9,
> >>>>>>> 0, 17, 200, -9999, -9999,
> >>>>>>> 0, 31, 200, -9999, 275,
> >>>>>>> 0, 32, 200, -9999, 36.008,
> >>>>>>> 0, 33, 200, -9999, 3.1383,
> >>>>>>> 0, 34, 200, -9999, -35.871,
> >>>>>>> 0, 7, 250, -9999, 10890,
> >>>>>>> 0, 11, 250, -9999, 228.9,
> >>>>>>> 0, 17, 250, -9999, -9999,
> >>>>>>> 0, 31, 250, -9999, 280,
> >>>>>>> 0, 32, 250, -9999, 29.3208,
> >>>>>>> 0, 33, 250, -9999, 5.0915,
> >>>>>>> 0, 34, 250, -9999, -28.8754,
> >>>>>>> 0, 7, 300, -9999, 9640,
> >>>>>>> 0, 11, 300, -9999, 238.3,
> >>>>>>> 0, 17, 300, -9999, 214.3,
> >>>>>>> 0, 31, 300, -9999, 290,
> >>>>>>> 0, 32, 300, -9999, 22.6336,
> >>>>>>> 0, 33, 300, -9999, 7.74115,
> >>>>>>> 0, 34, 300, -9999, -21.2686,
> >>>>>>> 0, 7, 400, -9999, 7570,
> >>>>>>> 0, 11, 400, -9999, 253.1,
> >>>>>>> 0, 17, 400, -9999, 229.1,
> >>>>>>> 0, 31, 400, -9999, 280,
> >>>>>>> 0, 32, 400, -9999, 17.4896,
> >>>>>>> 0, 33, 400, -9999, 3.03704,
> >>>>>>> 0, 34, 400, -9999, -17.2239,
> >>>>>>> 0, 7, 500, -9999, 5880,
> >>>>>>> 0, 11, 500, -9999, 265.7,
> >>>>>>> 0, 17, 500, -9999, 257.7,
> >>>>>>> 0, 31, 500, -9999, 305,
> >>>>>>> 0, 32, 500, -9999, 5.144,
> >>>>>>> 0, 33, 500, -9999, 2.95048,
> >>>>>>> 0, 34, 500, -9999, -4.21372,
> >>>>>>> 0, 7, 700, -9999, 3145,
> >>>>>>> 0, 11, 700, -9999, 279.6,
> >>>>>>> 0, 17, 700, -9999, 279.6,
> >>>>>>> 0, 31, 700, -9999, 325,
> >>>>>>> 0, 32, 700, -9999, 2.0576,
> >>>>>>> 0, 33, 700, -9999, 1.68549,
> >>>>>>> 0, 34, 700, -9999, -1.18019,
> >>>>>>> 0, 7, 850, -9999, 1518,
> >>>>>>> 0, 11, 850, -9999, 288.6,
> >>>>>>> 0, 17, 850, -9999, 287.1,
> >>>>>>> 0, 31, 850, -9999, 40,
> >>>>>>> 0, 32, 850, -9999, 1.5432,
> >>>>>>> 0, 33, 850, -9999, 1.18216,
> >>>>>>> 0, 34, 850, -9999, 0.99195,
> >>>>>>> 0, 7, 925, -9999, -9999,
> >>>>>>> 0, 11, 925, -9999, -9999,
> >>>>>>> 0, 17, 925, -9999, -9999,
> >>>>>>> 0, 31, 925, -9999, -9999,
> >>>>>>> 0, 32, 925, -9999, -9999,
> >>>>>>> 0, 33, 925, -9999, -9999,
> >>>>>>> 0, 34, 925, -9999, -9999,
> >>>>>>> 0, 7, 1000, -9999, -9999,
> >>>>>>> 0, 11, 1000, -9999, -9999,
> >>>>>>> 0, 17, 1000, -9999, -9999,
> >>>>>>> 0, 31, 1000, -9999, -9999,
> >>>>>>> 0, 32, 1000, -9999, -9999,
> >>>>>>> 0, 33, 1000, -9999, -9999,
> >>>>>>> 0, 34, 1000, -9999, -9999,
> >>>>>>>
> >>>>>>> Thanks for your help
> >>>>>>>
> >>>>>>> Gisela
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>> 2012/6/11 Paul Oldenburg via RT <met_help at ucar.edu>
> >>>>>>>
> >>>>>>>> According to our records, your request has been resolved.
If you
> >> have
> >>>>>> any
> >>>>>>>> further questions or concerns, please respond to this
message.
> >>>>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>
> >>
> >>
> >>
>
>
>
>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Fri Jul 06 13:38:38 2012
Paul,
We are processing succesfully poinstat comparing temperature from the
WRF-ARW vs soundings at standard levels,for October 2011.
Now, we want to process poinstat for U,V wind components , but it
fails.
We have the U fields from the p_interp applied to the wrfouts.. and
the U,V
components in the observation data file.Also, we follow your
instructions
to edit the configuration file e.g UU(0,0-0,*,*) and we get the
following
message:
DEBUG 2: Reading data for UU(0,0-0,*,*).
WARNING:
WARNING: PinterpFile::data(NcVar *, const LongArray &, DataPlane &,
double
&) const -> star found in bad slot
WARNING:
ERROR :
ERROR : process_fcst_climo_files() -> no fields matching UU(0,0-
0,*,*)
found in file: /u/charo/P_INTERP/wrfout_d01_2011-10-12_12:00:00_PLEV
Thanks in advanced for your help,
Gisela
2012/7/3 Paul Oldenburg via RT <met_help at ucar.edu>
> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #56795] Resolved: Question about PointStat problems
From: Paul Oldenburg
Time: Mon Jul 09 13:16:33 2012
Gisela,
MET support for p_interp fields is limited only to fields that are not
on a staggered grid, such as temperature:
float TT(Time, num_metgrid_levels, south_north, west_east) ;
Unfortunately, the p_interp software calculates and writes wind speed
components on the staggered grid:
float UU(Time, num_metgrid_levels, south_north, west_east_stag) ;
If you are interested in verifying wind speed components, you may want
to consider post-processing your model data using
the WRF Unified Post Processor (UPP - http://www.dtcenter.org/wrf-
nmm/users/overview/upp_overview.php). I'm sorry I
can't be more help.
Paul
On 07/06/2012 01:38 PM, Gisela Daniela Charó via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
>
> Paul,
> We are processing succesfully poinstat comparing temperature from
the
> WRF-ARW vs soundings at standard levels,for October 2011.
> Now, we want to process poinstat for U,V wind components , but it
fails.
>
> We have the U fields from the p_interp applied to the wrfouts.. and
the U,V
> components in the observation data file.Also, we follow your
instructions
> to edit the configuration file e.g UU(0,0-0,*,*) and we get the
following
> message:
>
> DEBUG 2: Reading data for UU(0,0-0,*,*).
> WARNING:
> WARNING: PinterpFile::data(NcVar *, const LongArray &, DataPlane &,
double
> &) const -> star found in bad slot
> WARNING:
> ERROR :
> ERROR : process_fcst_climo_files() -> no fields matching UU(0,0-
0,*,*)
> found in file: /u/charo/P_INTERP/wrfout_d01_2011-10-12_12:00:00_PLEV
>
> Thanks in advanced for your help,
> Gisela
> 2012/7/3 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> According to our records, your request has been resolved. If you
have any
>> further questions or concerns, please respond to this message.
>>
------------------------------------------------
Subject: Question about PointStat problems
From: Gisela Daniela Charó
Time: Tue Jul 10 08:08:24 2012
Paul,
Thanks for your help!
Gisela
2012/7/9 Paul Oldenburg via RT <met_help at ucar.edu>
> Gisela,
>
> MET support for p_interp fields is limited only to fields that are
not on
> a staggered grid, such as temperature:
>
> float TT(Time, num_metgrid_levels, south_north, west_east) ;
>
> Unfortunately, the p_interp software calculates and writes wind
speed
> components on the staggered grid:
>
> float UU(Time, num_metgrid_levels, south_north, west_east_stag) ;
>
> If you are interested in verifying wind speed components, you may
want to
> consider post-processing your model data using
> the WRF Unified Post Processor (UPP -
> http://www.dtcenter.org/wrf-nmm/users/overview/upp_overview.php).
I'm
> sorry I
> can't be more help.
>
> Paul
>
>
> On 07/06/2012 01:38 PM, Gisela Daniela Charó via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=56795 >
> >
> > Paul,
> > We are processing succesfully poinstat comparing temperature from
the
> > WRF-ARW vs soundings at standard levels,for October 2011.
> > Now, we want to process poinstat for U,V wind components , but it
fails.
> >
> > We have the U fields from the p_interp applied to the wrfouts..
and the
> U,V
> > components in the observation data file.Also, we follow your
> instructions
> > to edit the configuration file e.g UU(0,0-0,*,*) and we get the
following
> > message:
> >
> > DEBUG 2: Reading data for UU(0,0-0,*,*).
> > WARNING:
> > WARNING: PinterpFile::data(NcVar *, const LongArray &, DataPlane
&,
> double
> > &) const -> star found in bad slot
> > WARNING:
> > ERROR :
> > ERROR : process_fcst_climo_files() -> no fields matching UU(0,0-
0,*,*)
> > found in file: /u/charo/P_INTERP/wrfout_d01_2011-10-
12_12:00:00_PLEV
> >
> > Thanks in advanced for your help,
> > Gisela
> > 2012/7/3 Paul Oldenburg via RT <met_help at ucar.edu>
> >
> >> According to our records, your request has been resolved. If you
have
> any
> >> further questions or concerns, please respond to this message.
> >>
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list