[Met_help] [rt.rap.ucar.edu #52023] History for Problem in Point Stat
Paul Oldenburg via RT
met_help at ucar.edu
Thu Jan 5 09:37:11 MST 2012
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear MET Help
Actually I've gone through the Ticket
51020<https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51020> and
I found I am also facing same kind of problems.
In my case Point-Stat is working for a bigger domain but not for smaller
one.
for smaller domain it's giving as follows:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=360310862
Forecast File: ../DATA/FCST/wrfprs_d01_2011-06-25-000-f1-26-29_76-79
Climatology File: none
Configuration File: ../POINT_STAT/PointStatConfig_tmp
Observation File: ../DATA/OBS/GTS_OBS_TMP_2011-06-25-000-26-29_76-79.nc
--------------------------------------------------------------------------------
Reading records for TMP/Z2.
GRIB Record 244: Init = 20110624_000000, Valid = 20110625_000000, Accum =
000000, Min = 267.302, Max = 304.324
For TMP/Z2 found 1 forecast levels and 0 climatology levels.
--------------------------------------------------------------------------------
Searching 4 observations from 4 PrepBufr messages.
--------------------------------------------------------------------------------
Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over region
FULL, for interpolation method UW_MEAN(1), using 0 pairs.
Number of matched pairs = 0
Observations processed = 4
Rejected: GRIB code = 0
Rejected: valid time = 4
Rejected: bad obs value = 0
Rejected: off the grid = 0
Rejected: level mismatch = 0
Rejected: message type = 0
Rejected: masking region = 0
Rejected: bad fcst value = 0
and for bigger domain as follows:
--------------------------------------------------------------------------------
Searching 195 observations from 195 PrepBufr messages.
--------------------------------------------------------------------------------
Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over region
FULL, for interpolation method UW_MEAN(1), using 26 pairs.
Number of matched pairs = 26
Observations processed = 195
Rejected: GRIB code = 0
Rejected: valid time = 0
Rejected: bad obs value = 0
Rejected: off the grid = 144
Rejected: level mismatch = 0
Rejected: message type = 0
Rejected: masking region = 0
Rejected: bad fcst value = 25
Computing Categorical Statistics.
Computing Continuous Statistics.
Computing Scalar Partial Sums.
So, I've a few questions regarding this:
Does Point-Stat work for smaller regions, as it was working for bigger one?
Can we get the list of matched pairs lat/lon?
All the data I have transferred in incoming/irap/met_help/hvs_data for your
verification purpose.
--
* Regards*
*********************************
* HARVIR SINGH <http://www.harvir007.webs.com/>*
* **
***********************************
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Problem in Point Stat
From: Paul Oldenburg
Time: Wed Dec 14 15:38:09 2011
Harvir,
There are two problems with your small domain that I can see:
1. The obs for the small domain are valid at 20110601_00 and the
forecast data is valid at 20110625_00. For the obs
data, you can see the valid time using the ncdump utility:
$ ncdump -v hdr_vld GTS_OBS_TMP_2011-06-25-000-26-29_76-79.nc
The discrepancy in valid time can be overcome, if you wish to do so,
using the obs_valid_beg/end options, like you had
tried:
$MET_BASE/bin/point_stat \
wrfprs_d01_2011-06-25-000-f1-26-29_76-79 \
GTS_OBS_TMP_2011-06-25-000-26-29_76-79.nc \
PointStatConfig_tmp \
-obs_valid_beg 20110515 \
-obs_valid_end 20110701 \
-outdir . -v 3
2. The second problem is more important. It appears that the domain
of your model data is over Egypt, while your
observation points are in northern India. I attached two plots
showing the location of the data:
GTS_OBS_TMP_2011-06-25-000-26-29_76-79.png - created using the MET
utility plot_point_obs, showing the obs points (in
red) over northern India, I think
wrfprs_d01_2011-06-25-000-f1-26-29_76-79_domain.png - created using
IDV, showing the model data over southern Egypt, I think
If you have any other questions, please let me know.
Thanks,
Paul
On 12/14/2011 04:19 AM, Harvir Singh via RT wrote:
>
> Wed Dec 14 04:19:11 2011: Request 52023 was acted upon.
> Transaction: Ticket created by harvirduphy at gmail.com
> Queue: met_help
> Subject: Problem in Point Stat
> Owner: Nobody
> Requestors: harvirduphy at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=52023 >
>
>
> Dear MET Help
>
> Actually I've gone through the Ticket
> 51020<https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51020> and
> I found I am also facing same kind of problems.
> In my case Point-Stat is working for a bigger domain but not for
smaller
> one.
>
> for smaller domain it's giving as follows:
>
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=360310862
> Forecast File: ../DATA/FCST/wrfprs_d01_2011-06-25-000-f1-26-29_76-79
> Climatology File: none
> Configuration File: ../POINT_STAT/PointStatConfig_tmp
> Observation File: ../DATA/OBS/GTS_OBS_TMP_2011-06-25-000-26-29_76-
79.nc
>
>
--------------------------------------------------------------------------------
>
> Reading records for TMP/Z2.
> GRIB Record 244: Init = 20110624_000000, Valid = 20110625_000000,
Accum =
> 000000, Min = 267.302, Max = 304.324
> For TMP/Z2 found 1 forecast levels and 0 climatology levels.
>
>
--------------------------------------------------------------------------------
>
> Searching 4 observations from 4 PrepBufr messages.
>
>
--------------------------------------------------------------------------------
>
> Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over
region
> FULL, for interpolation method UW_MEAN(1), using 0 pairs.
> Number of matched pairs = 0
> Observations processed = 4
> Rejected: GRIB code = 0
> Rejected: valid time = 4
> Rejected: bad obs value = 0
> Rejected: off the grid = 0
> Rejected: level mismatch = 0
> Rejected: message type = 0
> Rejected: masking region = 0
> Rejected: bad fcst value = 0
>
>
> and for bigger domain as follows:
>
>
--------------------------------------------------------------------------------
>
> Searching 195 observations from 195 PrepBufr messages.
>
>
--------------------------------------------------------------------------------
>
> Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over
region
> FULL, for interpolation method UW_MEAN(1), using 26 pairs.
> Number of matched pairs = 26
> Observations processed = 195
> Rejected: GRIB code = 0
> Rejected: valid time = 0
> Rejected: bad obs value = 0
> Rejected: off the grid = 144
> Rejected: level mismatch = 0
> Rejected: message type = 0
> Rejected: masking region = 0
> Rejected: bad fcst value = 25
> Computing Categorical Statistics.
> Computing Continuous Statistics.
> Computing Scalar Partial Sums.
>
>
> So, I've a few questions regarding this:
>
> Does Point-Stat work for smaller regions, as it was working for
bigger one?
> Can we get the list of matched pairs lat/lon?
>
> All the data I have transferred in incoming/irap/met_help/hvs_data
for your
> verification purpose.
>
------------------------------------------------
Subject: Problem in Point Stat
From: Harvir Singh
Time: Fri Dec 16 00:00:54 2011
Paul,
I would like to thank you for your prompt reply.
Actually its working now.
Thank you again!!!
On 15 December 2011 04:08, Paul Oldenburg via RT <met_help at ucar.edu>
wrote:
> Harvir,
>
> There are two problems with your small domain that I can see:
>
> 1. The obs for the small domain are valid at 20110601_00 and the
forecast
> data is valid at 20110625_00. For the obs
> data, you can see the valid time using the ncdump utility:
>
> $ ncdump -v hdr_vld GTS_OBS_TMP_2011-06-25-000-26-29_76-79.nc
>
> The discrepancy in valid time can be overcome, if you wish to do so,
using
> the obs_valid_beg/end options, like you had
> tried:
>
> $MET_BASE/bin/point_stat \
> wrfprs_d01_2011-06-25-000-f1-26-29_76-79 \
> GTS_OBS_TMP_2011-06-25-000-26-29_76-79.nc \
> PointStatConfig_tmp \
> -obs_valid_beg 20110515 \
> -obs_valid_end 20110701 \
> -outdir . -v 3
>
>
> 2. The second problem is more important. It appears that the
domain of
> your model data is over Egypt, while your
> observation points are in northern India. I attached two plots
showing
> the location of the data:
>
> GTS_OBS_TMP_2011-06-25-000-26-29_76-79.png - created using the MET
utility
> plot_point_obs, showing the obs points (in
> red) over northern India, I think
> wrfprs_d01_2011-06-25-000-f1-26-29_76-79_domain.png - created using
IDV,
> showing the model data over southern Egypt, I think
>
> If you have any other questions, please let me know.
>
> Thanks,
>
> Paul
>
>
> On 12/14/2011 04:19 AM, Harvir Singh via RT wrote:
> >
> > Wed Dec 14 04:19:11 2011: Request 52023 was acted upon.
> > Transaction: Ticket created by harvirduphy at gmail.com
> > Queue: met_help
> > Subject: Problem in Point Stat
> > Owner: Nobody
> > Requestors: harvirduphy at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=52023 >
> >
> >
> > Dear MET Help
> >
> > Actually I've gone through the Ticket
> > 51020<https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51020> and
> > I found I am also facing same kind of problems.
> > In my case Point-Stat is working for a bigger domain but not for
smaller
> > one.
> >
> > for smaller domain it's giving as follows:
> >
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=360310862
> > Forecast File: ../DATA/FCST/wrfprs_d01_2011-06-25-000-f1-26-29_76-
79
> > Climatology File: none
> > Configuration File: ../POINT_STAT/PointStatConfig_tmp
> > Observation File: ../DATA/OBS/GTS_OBS_TMP_2011-06-25-000-26-29_76-
79.nc
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Reading records for TMP/Z2.
> > GRIB Record 244: Init = 20110624_000000, Valid = 20110625_000000,
Accum =
> > 000000, Min = 267.302, Max = 304.324
> > For TMP/Z2 found 1 forecast levels and 0 climatology levels.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Searching 4 observations from 4 PrepBufr messages.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over
region
> > FULL, for interpolation method UW_MEAN(1), using 0 pairs.
> > Number of matched pairs = 0
> > Observations processed = 4
> > Rejected: GRIB code = 0
> > Rejected: valid time = 4
> > Rejected: bad obs value = 0
> > Rejected: off the grid = 0
> > Rejected: level mismatch = 0
> > Rejected: message type = 0
> > Rejected: masking region = 0
> > Rejected: bad fcst value = 0
> >
> >
> > and for bigger domain as follows:
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Searching 195 observations from 195 PrepBufr messages.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Processing TMP/Z2 versus TMP/Z2, for observation type ADPSFC, over
region
> > FULL, for interpolation method UW_MEAN(1), using 26 pairs.
> > Number of matched pairs = 26
> > Observations processed = 195
> > Rejected: GRIB code = 0
> > Rejected: valid time = 0
> > Rejected: bad obs value = 0
> > Rejected: off the grid = 144
> > Rejected: level mismatch = 0
> > Rejected: message type = 0
> > Rejected: masking region = 0
> > Rejected: bad fcst value = 25
> > Computing Categorical Statistics.
> > Computing Continuous Statistics.
> > Computing Scalar Partial Sums.
> >
> >
> > So, I've a few questions regarding this:
> >
> > Does Point-Stat work for smaller regions, as it was working for
bigger
> one?
> > Can we get the list of matched pairs lat/lon?
> >
> > All the data I have transferred in incoming/irap/met_help/hvs_data
for
> your
> > verification purpose.
> >
>
>
>
--
* Regards*
*********************************
* HARVIR SINGH <http://www.harvir007.webs.com/>*
* Project Scientist**
HCL @ NCMRWF
***********************************
------------------------------------------------
More information about the Met_help
mailing list