[Met_help] [rt.rap.ucar.edu #55790] History for Geeta met_help problem

Tue Apr 24 12:13:00 MDT 2012

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Dear administrator, 
I wish to seek help from the met_help, and I am sending my message for the last 3-4 days. But everytime I am told that there is some error that has occurred. 
I donot know whtheer my meage is reaching across or not? I have already registered but my query is not answered. 

kindly help me
Regards
geeta

Dr G. Agnihotri, 
Met Centre, IMD, Palace Road, 
Bangalore, India. 
080-22244419

----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help problem
From: Paul Oldenburg
Time: Wed Apr 11 09:24:51 2012

Geeta,

I created a met_help ticket for you.  Please respond to this email and
describe the MET problem you are having in detail.

Thanks,

Paul

On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
> 	"Geeta met_help problem",
> a summary of which appears below.
>
> There is no need to reply to this message right now.  Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #55790].
>
> Please include the string:
>
>           [rt.rap.ucar.edu #55790]
>
> in the subject line of all future correspondence about this issue.
To do so,
> you may reply to this message.
>
>                          Thank you,
>                          met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> Dear administrator,
> I wish to seek help from the met_help, and I am sending my message
for the last 3-4 days. But everytime I am told that there is some
error that has occurred.
> I donot know whtheer my meage is reaching across or not? I have
already registered but my query is not answered.
>
> kindly help me
> Regards
> geeta
>
> Dr G. Agnihotri,
> Met Centre, IMD, Palace Road,
> Bangalore, India.
> 080-22244419
>

------------------------------------------------
Subject: Geeta met_help problem
From: Paul Oldenburg
Time: Thu Apr 12 09:03:20 2012

Geeta,

First of all, if you want to send data files to me that are larger
than 2MB, please upload them to our FTP site using
the instructions found here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp.

I found several problems, but they were all easy to find and fix.
There was a config file syntax error in the version
of PointStatConfig_rain that you sent me.  You must use double quotes
around values in the obs_field setting:

obs_field[]  = [ "61/A24" ];

Next, I recommend using a verbosity setting of 3 or higher when you
start using point_stat.  This will cause point_stat
to print out the reason why points are rejected.  One reason is that
there is a mis-match between the valid time of your
model data and the valid time of your observation data.  The following
ncdump commands show the discrepancy:

[pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
...
                 APCP_24:valid_time = "20110610_030000" ;

[pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
precip24_2011061000.nc
...
  hdr_vld =
   "20110610_083000",
   ...

You can tell point_stat to allow this discrepancy using the following
two settings in PointStatConfig_rain:

beg_ds = -21600;
end_ds =  21600;

Next, the message type of your observations is ADPSFC, which you can
see using the following ncdump command:

[pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
precip24_2011061000.nc
...
  hdr_typ =
   "ADPSFC",
   ...

In PointStatConfig_rain, the following setting must be used to specify
this message type:

message_type[] = [ "ADPSFC" ];

Once these changes were made, I ran point_stat with your model and
observation data and the result was 26 matched pairs.
  I attached the output below and the PointStatConfig_rain file that I
used.  Please let me know if you have any more
questions.

One last thing, I think you might have a date typo in your
precip24_2011061000.ascii file on line 2.  Should the date
string "20110613_083000" actually be "20110610_083000" like all the
others?

Thanks,

Paul

$ point_stat model24_day.nc precip24_2011061000.nc
PointStatConfig_rain -v 3
DEBUG 3: It's a Mercator grid...
DEBUG 3: Mercator Data:
DEBUG 3: lat_ll = 9.097
DEBUG 3: lon_ll = -73.943
DEBUG 3: lat_ur = 21.711
DEBUG 3: lon_ur = -87.057
DEBUG 3: ny = 163
DEBUG 3: nx = 163
DEBUG 1: Default Config File:
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
DEBUG 1: User Config File: PointStatConfig_rain
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072433272399
DEBUG 1: Forecast File: model24_day.nc
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: precip24_2011061000.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for APCP_24(*,*).
DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0 climatology
levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 27 observations from 27 messages.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation type
ADPSFC, over region FULL, for interpolation method
UW_MEAN(1), using 26 pairs.
DEBUG 3: Number of matched pairs  = 26
DEBUG 3: Observations processed   = 27
DEBUG 3: Rejected: GRIB code      = 0
DEBUG 3: Rejected: valid time     = 1
DEBUG 3: Rejected: bad obs value  = 0
DEBUG 3: Rejected: off the grid   = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: message type   = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_fho.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_ctc.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_cts.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_cnt.txt
DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_sl1l2.txt
DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_sal1l2.txt
DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_vl1l2.txt
DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_val1l2.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_pct.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_pstd.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_pjc.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_prc.txt
DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V_mpr.txt

On 04/11/2012 11:53 PM, Geeta Geeta wrote:
>
> Dear met_help,
> I am using point stat tool for Rainfall verification in my region
Western part of Indian peninsula (9-21N/74-87E).
>
>
>   1. for this, I have converted the ascii observations to netcdf
using the ascii2nc tool (O/p file precip24_2011061000.nc)
>   2. I have converted the model o/p in .nc format to grib1 format
using WPP.
>   3. I have then used pcpcombine to get the desired RF/output ie. by
subtracting the RF of 27th hour from the 4th hour.  (model24_day.nc is
attached)
>
>
>   4. now i am running
>   -bash-3.2$ cd ..
>   -bash-3.2$ cd verf-met/
>   /oprn/model/wrf3/utils/met/METv3.0/verf-met
>   -bash-3.2$ ../bin/point_stat model24_day.nc precip24_2011061000.nc
~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
>   GSL_RNG_TYPE=mt19937
>   GSL_RNG_SEED=801800783
>   Forecast File: model24_day.nc
>   Climatology File: none
>   Configuration File:
/oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
>   Observation File: precip24_2011061000.nc
>
>
--------------------------------------------------------------------------------
>
>   Reading records for APCP_24(*,*).
>   For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.
>
>
--------------------------------------------------------------------------------
>
>   Searching 27 observations from 27 PrepBufr messages.
>
>
--------------------------------------------------------------------------------
>
>   Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using 0
pairs.
>
--------------------------------------------------------------------------------
>
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
>   Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
>   -bash-3.2$
>
>
>
>   My problem is that there is no data in the output files, when i
open them, it is only headers.
>   Kindly help.
>   I have plotted the grib1 files (after WPP by using grib2ctl.pl) to
see the domain and if all the stations are falling within this domain,
So all looks okay to me
>   geeta
>
>
>
>> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
problem
>> From: met_help at ucar.edu
>> To: geeta124 at hotmail.com
>> Date: Wed, 11 Apr 2012 09:24:51 -0600
>>
>> Geeta,
>>
>> I created a met_help ticket for you.  Please respond to this email
and describe the MET problem you are having in detail.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
>>> Greetings,
>>>
>>> This message has been automatically generated in response to the
>>> creation of a trouble ticket regarding:
>>> 	"Geeta met_help problem",
>>> a summary of which appears below.
>>>
>>> There is no need to reply to this message right now.  Your ticket
has been
>>> assigned an ID of [rt.rap.ucar.edu #55790].
>>>
>>> Please include the string:
>>>
>>>            [rt.rap.ucar.edu #55790]
>>>
>>> in the subject line of all future correspondence about this issue.
To do so,
>>> you may reply to this message.
>>>
>>>                           Thank you,
>>>                           met_help at ucar.edu
>>>
>>>
-------------------------------------------------------------------------
>>> Dear administrator,
>>> I wish to seek help from the met_help, and I am sending my message
for the last 3-4 days. But everytime I am told that there is some
error that has occurred.
>>> I donot know whtheer my meage is reaching across or not? I have
already registered but my query is not answered.
>>>
>>> kindly help me
>>> Regards
>>> geeta
>>>
>>> Dr G. Agnihotri,
>>> Met Centre, IMD, Palace Road,
>>> Bangalore, India.
>>> 080-22244419
>>>
>>
>>
>

------------------------------------------------
Subject: Geeta met_help problem
From: Paul Oldenburg
Time: Thu Apr 12 09:03:20 2012

////////////////////////////////////////////////////////////////////////////////
//
// Default point_stat configuration file
//
////////////////////////////////////////////////////////////////////////////////

//
// Specify a name to designate the model being verified.  This name
will be
// written to the second column of the ASCII output generated.
//
model = "WRF";

//
// Beginning and ending time offset values in seconds for observations
// to be used.  These time offsets are defined in reference to the
// forecast valid time, v.  Observations with a valid time falling in
the
// window [v+beg_ds, v+end_ds] will be used.
// These selections are overridden by the command line arguments
// -obs_valid_beg and -obs_valid_end.
//
beg_ds = -21600;
end_ds =  21600;

//
// Specify a comma-separated list of fields to be verified.  The
forecast and
// observation fields may be specified separately.  If the obs_field
parameter
// is left blank, it will default to the contents of fcst_field.
//
// Each field is specified as a GRIB code or abbreviation followed by
an
// accumulation or vertical level indicator for GRIB files or as a
variable name
// followed by a list of dimensions for NetCDF files output from
p_interp or MET.
//
// Specifying verification fields for GRIB files:
//    GC/ANNN for accumulation interval NNN
//    GC/ZNNN for vertical level NNN
//    GC/ZNNN-NNN for a range of vertical levels (MSL or AGL)
//    GC/PNNN for pressure level NNN in hPa
//    GC/PNNN-NNN for a range of pressure levels in hPa
//    GC/LNNN for a generic level type
//    GC/RNNN for a specific GRIB record number
//    Where GC is the number of or abbreviation for the grib code
//    to be verified.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
// Specifying verification fields for NetCDF files:
//    var_name(i,...,j,*,*) for a single field
//    var_name(i-j,*,*) for a range of fields
//    Where var_name is the name of the NetCDF variable,
//    and i,...,j specifies fixed dimension values,
//    and i-j specifies a range of values for a single dimension,
//    and *,* specifies the two dimensions for the gridded field.
//
//    NOTE: To verify winds as vectors rather than scalars,
//          specify UGRD (or 33) followed by VGRD (or 34) with the
//          same level values.
//
//    NOTE: To process a probability field, add "/PROB", such as
"POP/Z0/PROB".
//
// e.g. fcst_field[] = [ "SPFH/P500", "TMP/P500" ]; for a GRIB input
// e.g. fcst_field[] = [ "QVAPOR(0,5,*,*)", "TT(0,5,*,*)" ]; for
NetCDF input
//
fcst_field[] = [ "APCP_24(*,*)"];
obs_field[]  = [ "61/A24" ];

//
// Specify a comma-separated list of groups of thresholds to be
applied to the
// fields listed above.  Thresholds for the forecast and observation
fields
// may be specified separately.  If the obs_thresh parameter is left
blank,
// it will default to the contents of fcst_thresh.
//
// At least one threshold must be provided for each field listed
above.  The
// lengths of the "fcst_field" and "fcst_thresh" arrays must match, as
must
// lengths of the "obs_field" and "obs_thresh" arrays.  To apply
multiple
// thresholds to a field, separate the threshold values with a space.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
//    'lt' for less than     'le' for less than or equal to
//    'eq' for equal to      'ne' for not equal to
//    'gt' for greater than  'ge' for greater than or equal to
//
// NOTE: Thresholds for probabilities must begin with 0.0, end with
1.0,
//       and be preceeded by "ge".
//
// e.g. fcst_thresh[] = [ "gt80", "gt273" ];
//
fcst_thresh[] = [ "gt20"];
obs_thresh[]  = [ "gt20"];

//
// Specify a comma-separated list of thresholds to be used when
computing
// VL1L2 and VAL1L2 partial sums for winds.  The thresholds are
applied to the
// wind speed values derived from each U/V pair.  Only those U/V pairs
which meet
// the wind speed threshold criteria are retained.  If the
obs_wind_thresh
// parameter is left blank, it will default to the contents of
fcst_wind_thresh.
//
// To apply multiple wind speed thresholds, separate the threshold
values with a
// space.  Use "NA" to indicate that no wind speed threshold should be
applied.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
//    'lt' for less than     'le' for less than or equal to
//    'eq' for equal to      'ne' for not equal to
//    'gt' for greater than  'ge' for greater than or equal to
//    'NA' for no threshold
//
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
//
fcst_wind_thresh[] = [ "NA" ];
obs_wind_thresh[]  = [];

//
// Specify a comma-separated list of PrepBufr message types with which
// to perform the verification.  Statistics will be computed
separately
// for each message type specified.  At least one PrepBufr message
type
// must be provided.
// List of valid message types:
//    ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
//    MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
//    SFCSHP SPSSMI SYNDAT VADWND
//    ANYAIR (= AIRCAR, AIRCFT)
//    ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
//    ONLYSF (= ADPSFC, SFCSHP)
//
http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
//
// e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
//
message_type[] = [ "ADPSFC" ];

//
// Specify a comma-separated list of grids to be used in masking the
data over
// which to perform scoring.  An empty list indicates that no masking
grid
// should be performed.  The standard NCEP grids are named "GNNN"
where NNN
// indicates the three digit grid number.  Enter "FULL" to score over
the
// entire domain.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
//
// e.g. mask_grid[] = [ "FULL" ];
//
mask_grid[] = [ "FULL" ];

//
// Specify a comma-separated list of masking regions to be applied.
// An empty list indicates that no additional masks should be used.
// The masking regions may be defined in one of 4 ways:
//
// (1) An ASCII file containing a lat/lon polygon.
//     Latitude in degrees north and longitude in degrees east.
//     By default, the first and last polygon points are connected.
//     e.g. "MET_BASE/data/poly/EAST.poly" which consists of n points:
//          "poly_name lat1 lon1 lat2 lon2... latn lonn"
//
// (2) The NetCDF output of the gen_poly_mask tool.
//
// (3) A NetCDF data file, followed by the name of the NetCDF variable
//     to be used, and optionally, a threshold to be applied to the
field.
//     e.g. "sample.nc var_name gt0.00"
//
// (4) A GRIB data file, followed by a description of the field
//     to be used, and optionally, a threshold to be applied to the
field.
//     e.g. "sample.grb APCP/A3 gt0.00"
//
// Any NetCDF or GRIB file used must have the same grid dimensions as
the
// data being verified.
//
// MET_BASE may be used in the path for the files above.
//
// e.g. mask_poly[] = [ "MET_BASE/data/poly/EAST.poly",
//                      "poly_mask.ncf",
//                      "sample.nc APCP",
//                      "sample.grb HGT/Z0 gt100.0" ];
//
mask_poly[] = [];

//
// Specify the name of an ASCII file containing a space-separated list
of
// station ID's at which to perform verification.  Each station ID
specified
// is treated as an individual masking region.
//
// An empty list file name indicates that no station ID masks should
be used.
//
// MET_BASE may be used in the path for the station ID mask file name.
//
// e.g. mask_sid = "CONUS.stations";
//
mask_sid = "";

//
// Specify a comma-separated list of values for alpha to be used when
computing
// confidence intervals.  Values of alpha must be between 0 and 1.
//
// e.g. ci_alpha[] = [ 0.05, 0.10 ];
//
ci_alpha[] = [ 0.05 ];

//
// Specify the method to be used for computing bootstrap confidence
intervals.
// The value for this is interpreted as follows:
//    (0) Use the BCa interval method (computationally intensive)
//    (1) Use the percentile interval method
//
boot_interval = 1;

//
// Specify a proportion between 0 and 1 to define the replicate sample
size
// to be used when computing percentile intervals.  The replicate
sample
// size is set to boot_rep_prop * n, where n is the number of raw data
points.
//
// e.g boot_rep_prop = 0.80;
//
boot_rep_prop = 1.0;

//
// Specify the number of times each set of matched pair data should be
// resampled when computing bootstrap confidence intervals.  A value
of
// zero disables the computation of bootstrap condifence intervals.
//
// e.g. n_boot_rep = 1000;
//
n_boot_rep = 1000;

//
// Specify the name of the random number generator to be used.  See
the MET
// Users Guide for a list of possible random number generators.
//
boot_rng = "mt19937";

//
// Specify the seed value to be used when computing bootstrap
confidence
// intervals.  If left unspecified, the seed will change for each run
and
// the computed bootstrap confidence intervals will not be
reproducable.
//
boot_seed = "";

//
// Specify a comma-separated list of interpolation method(s) to be
used
// for comparing the forecast grid to the observation points.  String
values
// are interpreted as follows:
//    MIN     = Minimum in the neighborhood
//    MAX     = Maximum in the neighborhood
//    MEDIAN  = Median in the neighborhood
//    UW_MEAN = Unweighted mean in the neighborhood
//    DW_MEAN = Distance-weighted mean in the neighborhood
//    LS_FIT  = Least-squares fit in the neighborhood
//
// In all cases, vertical interpolation is performed in the natural
log
// of pressure of the levels above and below the observation.
//
// e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
//
interp_method[] = [ "UW_MEAN" ];

//
// Specify a comma-separated list of box widths to be used by the
// interpolation techniques listed above.  A value of 1 indicates that
// the nearest neighbor approach should be used.  For a value of n
// greater than 1, the n*n grid points closest to the observation
define
// the neighborhood.
//
// e.g. interp_width = [ 1, 3, 5 ];
//
interp_width[] = [ 1 ];

//
// When interpolating, compute a ratio of the number of valid data
points
// to the total number of points in the neighborhood.  If that ratio
is
// less than this threshold, do not include the observation.  This
// threshold must be between 0 and 1.  Setting this threshold to 1
will
// require that each observation be surrounded by n*n valid forecast
// points.
//
// e.g. interp_thresh = 1.0;
//
interp_thresh = 1.0;

//
// Specify flags to indicate the type of data to be output:
//    (1) STAT and FHO Text Files, Forecast, Hit, Observation Rates:
//           Total (TOTAL),
//           Forecast Rate (F_RATE),
//           Hit Rate (H_RATE),
//           Observation Rate (O_RATE)
//
//    (2) STAT and CTC Text Files, Contingency Table Counts:
//           Total (TOTAL),
//           Forecast Yes and Observation Yes Count (FY_OY),
//           Forecast Yes and Observation No Count (FY_ON),
//           Forecast No and Observation Yes Count (FN_OY),
//           Forecast No and Observation No Count (FN_ON)
//
//    (3) STAT and CTS Text Files, Contingency Table Scores:
//           Total (TOTAL),
//           Base Rate (BASER),
//           Forecast Mean (FMEAN),
//           Accuracy (ACC),
//           Frequency Bias (FBIAS),
//           Probability of Detecting Yes (PODY),
//           Probability of Detecting No (PODN),
//           Probability of False Detection (POFD),
//           False Alarm Ratio (FAR),
//           Critical Success Index (CSI),
//           Gilbert Skill Score (GSS),
//           Hanssen and Kuipers Discriminant (HK),
//           Heidke Skill Score (HSS),
//           Odds Ratio (ODDS),
//           NOTE: All statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (4) STAT and MCTC Text Files, NxN Multi-Category Contingency
Table Counts:
//           Total (TOTAL),
//           Number of Categories (N_CAT),
//           Contingency Table Count columns repeated N_CAT*N_CAT
times
//
//    (5) STAT and MCTS Text Files, NxN Multi-Category Contingency
Table Scores:
//           Total (TOTAL),
//           Number of Categories (N_CAT),
//           Accuracy (ACC),
//           Hanssen and Kuipers Discriminant (HK),
//           Heidke Skill Score (HSS),
//           Gerrity Score (GER),
//           NOTE: All statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (6) STAT and CNT Text Files, Statistics of Continuous Variables:
//           Total (TOTAL),
//           Forecast Mean (FBAR),
//           Forecast Standard Deviation (FSTDEV),
//           Observation Mean (OBAR),
//           Observation Standard Deviation (OSTDEV),
//           Pearson's Correlation Coefficient (PR_CORR),
//           Spearman's Rank Correlation Coefficient (SP_CORR),
//           Kendall Tau Rank Correlation Coefficient (KT_CORR),
//           Number of ranks compared (RANKS),
//           Number of tied ranks in the forecast field (FRANK_TIES),
//           Number of tied ranks in the observation field
(ORANK_TIES),
//           Mean Error (ME),
//           Standard Deviation of the Error (ESTDEV),
//           Multiplicative Bias (MBIAS = FBAR - OBAR),
//           Mean Absolute Error (MAE),
//           Mean Squared Error (MSE),
//           Bias-Corrected Mean Squared Error (BCMSE),
//           Root Mean Squared Error (RMSE),
//           Percentiles of the Error (E10, E25, E50, E75, E90)
//           NOTE: Most statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (7) STAT and SL1L2 Text Files, Scalar Partial Sums:
//           Total (TOTAL),
//           Forecast Mean (FBAR),
//              = mean(f)
//           Observation Mean (OBAR),
//              = mean(o)
//           Forecast*Observation Product Mean (FOBAR),
//              = mean(f*o)
//           Forecast Squared Mean (FFBAR),
//              = mean(f^2)
//           Observation Squared Mean (OOBAR)
//              = mean(o^2)
//
//    (8) STAT and SAL1L2 Text Files, Scalar Anomaly Partial Sums:
//           Total (TOTAL),
//           Forecast Anomaly Mean (FABAR),
//              = mean(f-c)
//           Observation Anomaly Mean (OABAR),
//              = mean(o-c)
//           Product of Forecast and Observation Anomalies Mean
(FOABAR),
//              = mean((f-c)*(o-c))
//           Forecast Anomaly Squared Mean (FFABAR),
//              = mean((f-c)^2)
//           Observation Anomaly Squared Mean (OOABAR)
//              = mean((o-c)^2)
//
//    (9) STAT and VL1L2 Text Files, Vector Partial Sums:
//           Total (TOTAL),
//           U-Forecast Mean (UFBAR),
//              = mean(uf)
//           V-Forecast Mean (VFBAR),
//              = mean(vf)
//           U-Observation Mean (UOBAR),
//              = mean(uo)
//           V-Observation Mean (VOBAR),
//              = mean(vo)
//           U-Product Plus V-Product (UVFOBAR),
//              = mean(uf*uo+vf*vo)
//           U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
//              = mean(uf^2+vf^2)
//           U-Observation Squared Plus V-Observation Squared
(UVOOBAR)
//              = mean(uo^2+vo^2)
//
//   (10) STAT and VAL1L2 Text Files, Vector Anomaly Partial Sums:
//           U-Forecast Anomaly Mean (UFABAR),
//              = mean(uf-uc)
//           V-Forecast Anomaly Mean (VFABAR),
//              = mean(vf-vc)
//           U-Observation Anomaly Mean (UOABAR),
//              = mean(uo-uc)
//           V-Observation Anomaly Mean (VOABAR),
//              = mean(vo-vc)
//           U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
//              = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
//           U-Forecast Anomaly Squared Plus V-Forecast Anomaly
Squared (UVFFABAR),
//              = mean((uf-uc)^2+(vf-vc)^2)
//           U-Observation Anomaly Squared Plus V-Observation Anomaly
Squared (UVOOABAR)
//              = mean((uo-uc)^2+(vo-vc)^2)
//
//   (11) STAT and PCT Text Files, Nx2 Probability Contingency Table
Counts:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Row Observation Yes Count (OY_i),
//           Row Observation No Count (ON_i),
//           NOTE: Previous 3 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (12) STAT and PSTD Text Files, Nx2 Probability Contingency Table
Scores:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Base Rate (BASER) with confidence interval limits,
//           Reliability (RELIABILITY),
//           Resolution (RESOLUTION),
//           Uncertainty (UNCERTAINTY),
//           Area Under the ROC Curve (ROC_AUC),
//           Brier Score (BRIER) with confidence interval limits,
//           Probability Threshold Value (THRESH_i)
//           NOTE: Previous column repeated for each probability
threshold.
//
//   (13) STAT and PJC Text Files, Joint/Continuous Statistics of
//                                 Probabilistic Variables:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Observation Yes Count Divided by Total (OY_TP_i),
//           Observation No Count Divided by Total (ON_TP_i),
//           Calibration (CALIBRATION_i),
//           Refinement (REFINEMENT_i),
//           Likelikhood (LIKELIHOOD_i),
//           Base Rate (BASER_i),
//           NOTE: Previous 7 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (14) STAT and PRC Text Files, ROC Curve Points for
//                                 Probabilistic Variables:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Probability of Detecting Yes (PODY_i),
//           Probability of False Detection (POFD_i),
//           NOTE: Previous 3 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (15) STAT and MPR Text Files, Matched Pair Data:
//           Total (TOTAL),
//           Index (INDEX),
//           Observation Station ID (OBS_SID),
//           Observation Latitude (OBS_LAT),
//           Observation Longitude (OBS_LON),
//           Observation Level (OBS_LVL),
//           Observation Elevation (OBS_ELV),
//           Forecast Value (FCST),
//           Observation Value (OBS),
//           Climatological Value (CLIMO)
//
//   In the expressions above, f are forecast values, o are observed
values,
//   and c are climatological values.
//
// Values for these flags are interpreted as follows:
//    (0) Do not generate output of this type
//    (1) Write output to a STAT file
//    (2) Write output to a STAT file and a text file
//
output_flag[] = [ 2, 2, 2, 0, 0, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];

//
// Flag to indicate whether Kendall's Tau and Spearman's Rank
Correlation
// Coefficients should be computed.  Computing them over large
datasets is
// computationally intensive and slows down the runtime execution
significantly.
//    (0) Do not compute these correlation coefficients
//    (1) Compute these correlation coefficients
//
rank_corr_flag = 1;

//
// Specify the GRIB Table 2 parameter table version number to be used
// for interpreting GRIB codes.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
grib_ptv = 2;

//
// Directory where temporary files should be written.
//
tmp_dir = "/tmp";

//
// Prefix to be used for the output file names.
//
output_prefix = "";

//
// Indicate a version number for the contents of this configuration
file.
// The value should generally not be modified.
//
version = "V3.0";

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Mon Apr 16 03:08:18 2012

Dear Met_help,
I have done the corrections in my configuration file and it works fine
now.
Thanks for Ur support.

geeta

> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
problem, error with point-stat
> From: met_help at ucar.edu
> To: geeta124 at hotmail.com
> Date: Thu, 12 Apr 2012 09:03:20 -0600
>
> Geeta,
>
> First of all, if you want to send data files to me that are larger
than 2MB, please upload them to our FTP site using
> the instructions found here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp.
>
> I found several problems, but they were all easy to find and fix.
There was a config file syntax error in the version
> of PointStatConfig_rain that you sent me. You must use double quotes
around values in the obs_field setting:
>
> obs_field[] = [ "61/A24" ];
>
> Next, I recommend using a verbosity setting of 3 or higher when you
start using point_stat. This will cause point_stat
> to print out the reason why points are rejected. One reason is that
there is a mis-match between the valid time of your
> model data and the valid time of your observation data. The
following ncdump commands show the discrepancy:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
> ...
> APCP_24:valid_time = "20110610_030000" ;
>
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
precip24_2011061000.nc
> ...
> hdr_vld =
> "20110610_083000",
> ...
>
> You can tell point_stat to allow this discrepancy using the
following two settings in PointStatConfig_rain:
>
> beg_ds = -21600;
> end_ds = 21600;
>
> Next, the message type of your observations is ADPSFC, which you can
see using the following ncdump command:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
precip24_2011061000.nc
> ...
> hdr_typ =
> "ADPSFC",
> ...
>
> In PointStatConfig_rain, the following setting must be used to
specify this message type:
>
> message_type[] = [ "ADPSFC" ];
>
> Once these changes were made, I ran point_stat with your model and
observation data and the result was 26 matched pairs.
> I attached the output below and the PointStatConfig_rain file that I
used. Please let me know if you have any more
> questions.
>
> One last thing, I think you might have a date typo in your
precip24_2011061000.ascii file on line 2. Should the date
> string "20110613_083000" actually be "20110610_083000" like all the
others?
>
> Thanks,
>
> Paul
>
>
> $ point_stat model24_day.nc precip24_2011061000.nc
PointStatConfig_rain -v 3
> DEBUG 3: It's a Mercator grid...
> DEBUG 3: Mercator Data:
> DEBUG 3: lat_ll = 9.097
> DEBUG 3: lon_ll = -73.943
> DEBUG 3: lat_ur = 21.711
> DEBUG 3: lon_ur = -87.057
> DEBUG 3: ny = 163
> DEBUG 3: nx = 163
> DEBUG 1: Default Config File:
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
> DEBUG 1: User Config File: PointStatConfig_rain
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072433272399
> DEBUG 1: Forecast File: model24_day.nc
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File: precip24_2011061000.nc
> DEBUG 2:
> DEBUG 2:
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for APCP_24(*,*).
> DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0 climatology
levels.
> DEBUG 2:
> DEBUG 2:
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 27 observations from 27 messages.
> DEBUG 2:
> DEBUG 2:
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation type
ADPSFC, over region FULL, for interpolation method
> UW_MEAN(1), using 26 pairs.
> DEBUG 3: Number of matched pairs = 26
> DEBUG 3: Observations processed = 27
> DEBUG 3: Rejected: GRIB code = 0
> DEBUG 3: Rejected: valid time = 1
> DEBUG 3: Rejected: bad obs value = 0
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 0
> DEBUG 3: Rejected: message type = 0
> DEBUG 3: Rejected: masking region = 0
> DEBUG 3: Rejected: bad fcst value = 0
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Computing Scalar Partial Sums.
> DEBUG 2:
> DEBUG 2:
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_fho.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_ctc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cts.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cnt.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_sl1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_sal1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_vl1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_val1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pct.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pstd.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pjc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_prc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_mpr.txt
>
>
>
> On 04/11/2012 11:53 PM, Geeta Geeta wrote:
> >
> > Dear met_help,
> > I am using point stat tool for Rainfall verification in my region
Western part of Indian peninsula (9-21N/74-87E).
> >
> >
> > 1. for this, I have converted the ascii observations to netcdf
using the ascii2nc tool (O/p file precip24_2011061000.nc)
> > 2. I have converted the model o/p in .nc format to grib1 format
using WPP.
> > 3. I have then used pcpcombine to get the desired RF/output ie. by
subtracting the RF of 27th hour from the 4th hour. (model24_day.nc is
attached)
> >
> >
> > 4. now i am running
> > -bash-3.2$ cd ..
> > -bash-3.2$ cd verf-met/
> > /oprn/model/wrf3/utils/met/METv3.0/verf-met
> > -bash-3.2$ ../bin/point_stat model24_day.nc precip24_2011061000.nc
~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=801800783
> > Forecast File: model24_day.nc
> > Climatology File: none
> > Configuration File:
/oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > Observation File: precip24_2011061000.nc
> >
> >
--------------------------------------------------------------------------------
> >
> > Reading records for APCP_24(*,*).
> > For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.
> >
> >
--------------------------------------------------------------------------------
> >
> > Searching 27 observations from 27 PrepBufr messages.
> >
> >
--------------------------------------------------------------------------------
> >
> > Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using 0
pairs.
> >
--------------------------------------------------------------------------------
> >
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
> > Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
> > -bash-3.2$
> >
> >
> >
> > My problem is that there is no data in the output files, when i
open them, it is only headers.
> > Kindly help.
> > I have plotted the grib1 files (after WPP by using grib2ctl.pl) to
see the domain and if all the stations are falling within this domain,
So all looks okay to me
> > geeta
> >
> >
> >
> >> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
problem
> >> From: met_help at ucar.edu
> >> To: geeta124 at hotmail.com
> >> Date: Wed, 11 Apr 2012 09:24:51 -0600
> >>
> >> Geeta,
> >>
> >> I created a met_help ticket for you. Please respond to this email
and describe the MET problem you are having in detail.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> >>> Greetings,
> >>>
> >>> This message has been automatically generated in response to the
> >>> creation of a trouble ticket regarding:
> >>> "Geeta met_help problem",
> >>> a summary of which appears below.
> >>>
> >>> There is no need to reply to this message right now. Your ticket
has been
> >>> assigned an ID of [rt.rap.ucar.edu #55790].
> >>>
> >>> Please include the string:
> >>>
> >>> [rt.rap.ucar.edu #55790]
> >>>
> >>> in the subject line of all future correspondence about this
issue. To do so,
> >>> you may reply to this message.
> >>>
> >>> Thank you,
> >>> met_help at ucar.edu
> >>>
> >>>
-------------------------------------------------------------------------
> >>> Dear administrator,
> >>> I wish to seek help from the met_help, and I am sending my
message for the last 3-4 days. But everytime I am told that there is
some error that has occurred.
> >>> I donot know whtheer my meage is reaching across or not? I have
already registered but my query is not answered.
> >>>
> >>> kindly help me
> >>> Regards
> >>> geeta
> >>>
> >>> Dr G. Agnihotri,
> >>> Met Centre, IMD, Palace Road,
> >>> Bangalore, India.
> >>> 080-22244419
> >>>
> >>
> >>
> >
>
>

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

Dear met_help,
I have a small problem which goes like this. I am unablle to make
point_stat write the output in the folder that I want . It writes in
the default folder, METv3.0/out/point-stat.
There is no option in the configuration file specifying so,.
I am attaching the files.
Thanks.
geeta

> To: met_help_archive at mailman.ucar.edu
> From: met_help at ucar.edu
> Date: Mon, 16 Apr 2012 08:52:43 -0600
> CC: met_help at ucar.edu
> Subject: [Met_help] [rt.rap.ucar.edu #55790] History for Geeta
met_help problem
>
> ----------------------------------------------------------------
> Initial Request
> ----------------------------------------------------------------
>
> Dear administrator,
> I wish to seek help from the met_help, and I am sending my message
for the last 3-4 days. But everytime I am told that there is some
error that has occurred.
> I donot know whtheer my meage is reaching across or not? I have
already registered but my query is not answered.
>
> kindly help me
> Regards
> geeta
>
> Dr G. Agnihotri,
> Met Centre, IMD, Palace Road,
> Bangalore, India.
> 080-22244419
>
>
>
> ----------------------------------------------------------------
> Complete Ticket History
> ----------------------------------------------------------------
>
> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
problem
> From: Paul Oldenburg
> Time: Wed Apr 11 09:24:51 2012
>
> Geeta,
>
> I created a met_help ticket for you. Please respond to this email
and
> describe the MET problem you are having in detail.
>
> Thanks,
>
> Paul
>
>
> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> > Greetings,
> >
> > This message has been automatically generated in response to the
> > creation of a trouble ticket regarding:
> > "Geeta met_help problem",
> > a summary of which appears below.
> >
> > There is no need to reply to this message right now. Your ticket
> has been
> > assigned an ID of [rt.rap.ucar.edu #55790].
> >
> > Please include the string:
> >
> > [rt.rap.ucar.edu #55790]
> >
> > in the subject line of all future correspondence about this issue.
> To do so,
> > you may reply to this message.
> >
> > Thank you,
> > met_help at ucar.edu
> >
> >
>
-------------------------------------------------------------------------
> > Dear administrator,
> > I wish to seek help from the met_help, and I am sending my message
> for the last 3-4 days. But everytime I am told that there is some
> error that has occurred.
> > I donot know whtheer my meage is reaching across or not? I have
> already registered but my query is not answered.
> >
> > kindly help me
> > Regards
> > geeta
> >
> > Dr G. Agnihotri,
> > Met Centre, IMD, Palace Road,
> > Bangalore, India.
> > 080-22244419
> >
>
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Paul Oldenburg
> Time: Thu Apr 12 09:03:20 2012
>
> Geeta,
>
> First of all, if you want to send data files to me that are larger
> than 2MB, please upload them to our FTP site using
> the instructions found here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
>
> I found several problems, but they were all easy to find and fix.
> There was a config file syntax error in the version
> of PointStatConfig_rain that you sent me. You must use double quotes
> around values in the obs_field setting:
>
> obs_field[] = [ "61/A24" ];
>
> Next, I recommend using a verbosity setting of 3 or higher when you
> start using point_stat. This will cause point_stat
> to print out the reason why points are rejected. One reason is that
> there is a mis-match between the valid time of your
> model data and the valid time of your observation data. The
following
> ncdump commands show the discrepancy:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
> ...
> APCP_24:valid_time = "20110610_030000" ;
>
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
> precip24_2011061000.nc
> ...
> hdr_vld =
> "20110610_083000",
> ...
>
> You can tell point_stat to allow this discrepancy using the
following
> two settings in PointStatConfig_rain:
>
> beg_ds = -21600;
> end_ds = 21600;
>
> Next, the message type of your observations is ADPSFC, which you can
> see using the following ncdump command:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
> precip24_2011061000.nc
> ...
> hdr_typ =
> "ADPSFC",
> ...
>
> In PointStatConfig_rain, the following setting must be used to
specify
> this message type:
>
> message_type[] = [ "ADPSFC" ];
>
> Once these changes were made, I ran point_stat with your model and
> observation data and the result was 26 matched pairs.
> I attached the output below and the PointStatConfig_rain file that I
> used. Please let me know if you have any more
> questions.
>
> One last thing, I think you might have a date typo in your
> precip24_2011061000.ascii file on line 2. Should the date
> string "20110613_083000" actually be "20110610_083000" like all the
> others?
>
> Thanks,
>
> Paul
>
>
> $ point_stat model24_day.nc precip24_2011061000.nc
> PointStatConfig_rain -v 3
> DEBUG 3: It's a Mercator grid...
> DEBUG 3: Mercator Data:
> DEBUG 3: lat_ll = 9.097
> DEBUG 3: lon_ll = -73.943
> DEBUG 3: lat_ur = 21.711
> DEBUG 3: lon_ur = -87.057
> DEBUG 3: ny = 163
> DEBUG 3: nx = 163
> DEBUG 1: Default Config File:
>
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
> DEBUG 1: User Config File: PointStatConfig_rain
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072433272399
> DEBUG 1: Forecast File: model24_day.nc
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File: precip24_2011061000.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for APCP_24(*,*).
> DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0 climatology
> levels.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 27 observations from 27 messages.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation type
> ADPSFC, over region FULL, for interpolation method
> UW_MEAN(1), using 26 pairs.
> DEBUG 3: Number of matched pairs = 26
> DEBUG 3: Observations processed = 27
> DEBUG 3: Rejected: GRIB code = 0
> DEBUG 3: Rejected: valid time = 1
> DEBUG 3: Rejected: bad obs value = 0
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 0
> DEBUG 3: Rejected: message type = 0
> DEBUG 3: Rejected: masking region = 0
> DEBUG 3: Rejected: bad fcst value = 0
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Computing Scalar Partial Sums.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_fho.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_ctc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cts.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cnt.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sl1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sal1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_vl1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_val1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pct.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pstd.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pjc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_prc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_mpr.txt
>
>
>
> On 04/11/2012 11:53 PM, Geeta Geeta wrote:
> >
> > Dear met_help,
> > I am using point stat tool for Rainfall verification in my region
> Western part of Indian peninsula (9-21N/74-87E).
> >
> >
> > 1. for this, I have converted the ascii observations to netcdf
> using the ascii2nc tool (O/p file precip24_2011061000.nc)
> > 2. I have converted the model o/p in .nc format to grib1 format
> using WPP.
> > 3. I have then used pcpcombine to get the desired RF/output ie. by
> subtracting the RF of 27th hour from the 4th hour. (model24_day.nc
is
> attached)
> >
> >
> > 4. now i am running
> > -bash-3.2$ cd ..
> > -bash-3.2$ cd verf-met/
> > /oprn/model/wrf3/utils/met/METv3.0/verf-met
> > -bash-3.2$ ../bin/point_stat model24_day.nc precip24_2011061000.nc
> ~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=801800783
> > Forecast File: model24_day.nc
> > Climatology File: none
> > Configuration File:
> /oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > Observation File: precip24_2011061000.nc
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Reading records for APCP_24(*,*).
> > For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Searching 27 observations from 27 PrepBufr messages.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Processing APCP_24(*,*) versus APCP_24/A24, for observation type
> ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
0
> pairs.
> >
>
--------------------------------------------------------------------------------
> >
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
> > -bash-3.2$
> >
> >
> >
> > My problem is that there is no data in the output files, when i
> open them, it is only headers.
> > Kindly help.
> > I have plotted the grib1 files (after WPP by using grib2ctl.pl) to
> see the domain and if all the stations are falling within this
domain,
> So all looks okay to me
> > geeta
> >
> >
> >
> >> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem
> >> From: met_help at ucar.edu
> >> To: geeta124 at hotmail.com
> >> Date: Wed, 11 Apr 2012 09:24:51 -0600
> >>
> >> Geeta,
> >>
> >> I created a met_help ticket for you. Please respond to this email
> and describe the MET problem you are having in detail.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> >>> Greetings,
> >>>
> >>> This message has been automatically generated in response to the
> >>> creation of a trouble ticket regarding:
> >>> "Geeta met_help problem",
> >>> a summary of which appears below.
> >>>
> >>> There is no need to reply to this message right now. Your ticket
> has been
> >>> assigned an ID of [rt.rap.ucar.edu #55790].
> >>>
> >>> Please include the string:
> >>>
> >>> [rt.rap.ucar.edu #55790]
> >>>
> >>> in the subject line of all future correspondence about this
issue.
> To do so,
> >>> you may reply to this message.
> >>>
> >>> Thank you,
> >>> met_help at ucar.edu
> >>>
> >>>
>
-------------------------------------------------------------------------
> >>> Dear administrator,
> >>> I wish to seek help from the met_help, and I am sending my
message
> for the last 3-4 days. But everytime I am told that there is some
> error that has occurred.
> >>> I donot know whtheer my meage is reaching across or not? I have
> already registered but my query is not answered.
> >>>
> >>> kindly help me
> >>> Regards
> >>> geeta
> >>>
> >>> Dr G. Agnihotri,
> >>> Met Centre, IMD, Palace Road,
> >>> Bangalore, India.
> >>> 080-22244419
> >>>
> >>
> >>
> >
>
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Paul Oldenburg
> Time: Thu Apr 12 09:03:20 2012
>
>
////////////////////////////////////////////////////////////////////////////////
> //
> // Default point_stat configuration file
> //
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Specify a name to designate the model being verified. This name
> will be
> // written to the second column of the ASCII output generated.
> //
> model = "WRF";
>
> //
> // Beginning and ending time offset values in seconds for
observations
> // to be used. These time offsets are defined in reference to the
> // forecast valid time, v. Observations with a valid time falling in
> the
> // window [v+beg_ds, v+end_ds] will be used.
> // These selections are overridden by the command line arguments
> // -obs_valid_beg and -obs_valid_end.
> //
> beg_ds = -21600;
> end_ds = 21600;
>
> //
> // Specify a comma-separated list of fields to be verified. The
> forecast and
> // observation fields may be specified separately. If the obs_field
> parameter
> // is left blank, it will default to the contents of fcst_field.
> //
> // Each field is specified as a GRIB code or abbreviation followed
by
> an
> // accumulation or vertical level indicator for GRIB files or as a
> variable name
> // followed by a list of dimensions for NetCDF files output from
> p_interp or MET.
> //
> // Specifying verification fields for GRIB files:
> // GC/ANNN for accumulation interval NNN
> // GC/ZNNN for vertical level NNN
> // GC/ZNNN-NNN for a range of vertical levels (MSL or AGL)
> // GC/PNNN for pressure level NNN in hPa
> // GC/PNNN-NNN for a range of pressure levels in hPa
> // GC/LNNN for a generic level type
> // GC/RNNN for a specific GRIB record number
> // Where GC is the number of or abbreviation for the grib code
> // to be verified.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> //
> // Specifying verification fields for NetCDF files:
> // var_name(i,...,j,*,*) for a single field
> // var_name(i-j,*,*) for a range of fields
> // Where var_name is the name of the NetCDF variable,
> // and i,...,j specifies fixed dimension values,
> // and i-j specifies a range of values for a single dimension,
> // and *,* specifies the two dimensions for the gridded field.
> //
> // NOTE: To verify winds as vectors rather than scalars,
> // specify UGRD (or 33) followed by VGRD (or 34) with the
> // same level values.
> //
> // NOTE: To process a probability field, add "/PROB", such as
> "POP/Z0/PROB".
> //
> // e.g. fcst_field[] = [ "SPFH/P500", "TMP/P500" ]; for a GRIB input
> // e.g. fcst_field[] = [ "QVAPOR(0,5,*,*)", "TT(0,5,*,*)" ]; for
> NetCDF input
> //
> fcst_field[] = [ "APCP_24(*,*)"];
> obs_field[] = [ "61/A24" ];
>
> //
> // Specify a comma-separated list of groups of thresholds to be
> applied to the
> // fields listed above. Thresholds for the forecast and observation
> fields
> // may be specified separately. If the obs_thresh parameter is left
> blank,
> // it will default to the contents of fcst_thresh.
> //
> // At least one threshold must be provided for each field listed
> above. The
> // lengths of the "fcst_field" and "fcst_thresh" arrays must match,
as
> must
> // lengths of the "obs_field" and "obs_thresh" arrays. To apply
> multiple
> // thresholds to a field, separate the threshold values with a
space.
> //
> // Each threshold must be preceded by a two letter indicator for the
> type of
> // thresholding to be performed:
> // 'lt' for less than 'le' for less than or equal to
> // 'eq' for equal to 'ne' for not equal to
> // 'gt' for greater than 'ge' for greater than or equal to
> //
> // NOTE: Thresholds for probabilities must begin with 0.0, end with
> 1.0,
> // and be preceeded by "ge".
> //
> // e.g. fcst_thresh[] = [ "gt80", "gt273" ];
> //
> fcst_thresh[] = [ "gt20"];
> obs_thresh[] = [ "gt20"];
>
> //
> // Specify a comma-separated list of thresholds to be used when
> computing
> // VL1L2 and VAL1L2 partial sums for winds. The thresholds are
> applied to the
> // wind speed values derived from each U/V pair. Only those U/V
pairs
> which meet
> // the wind speed threshold criteria are retained. If the
> obs_wind_thresh
> // parameter is left blank, it will default to the contents of
> fcst_wind_thresh.
> //
> // To apply multiple wind speed thresholds, separate the threshold
> values with a
> // space. Use "NA" to indicate that no wind speed threshold should
be
> applied.
> //
> // Each threshold must be preceded by a two letter indicator for the
> type of
> // thresholding to be performed:
> // 'lt' for less than 'le' for less than or equal to
> // 'eq' for equal to 'ne' for not equal to
> // 'gt' for greater than 'ge' for greater than or equal to
> // 'NA' for no threshold
> //
> // e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
> //
> fcst_wind_thresh[] = [ "NA" ];
> obs_wind_thresh[] = [];
>
> //
> // Specify a comma-separated list of PrepBufr message types with
which
> // to perform the verification. Statistics will be computed
> separately
> // for each message type specified. At least one PrepBufr message
> type
> // must be provided.
> // List of valid message types:
> // ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
> // MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
> // SFCSHP SPSSMI SYNDAT VADWND
> // ANYAIR (= AIRCAR, AIRCFT)
> // ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
> // ONLYSF (= ADPSFC, SFCSHP)
> //
>
http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
> //
> // e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
> //
> message_type[] = [ "ADPSFC" ];
>
> //
> // Specify a comma-separated list of grids to be used in masking the
> data over
> // which to perform scoring. An empty list indicates that no masking
> grid
> // should be performed. The standard NCEP grids are named "GNNN"
> where NNN
> // indicates the three digit grid number. Enter "FULL" to score over
> the
> // entire domain.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
> //
> // e.g. mask_grid[] = [ "FULL" ];
> //
> mask_grid[] = [ "FULL" ];
>
> //
> // Specify a comma-separated list of masking regions to be applied.
> // An empty list indicates that no additional masks should be used.
> // The masking regions may be defined in one of 4 ways:
> //
> // (1) An ASCII file containing a lat/lon polygon.
> // Latitude in degrees north and longitude in degrees east.
> // By default, the first and last polygon points are connected.
> // e.g. "MET_BASE/data/poly/EAST.poly" which consists of n points:
> // "poly_name lat1 lon1 lat2 lon2... latn lonn"
> //
> // (2) The NetCDF output of the gen_poly_mask tool.
> //
> // (3) A NetCDF data file, followed by the name of the NetCDF
variable
> // to be used, and optionally, a threshold to be applied to the
> field.
> // e.g. "sample.nc var_name gt0.00"
> //
> // (4) A GRIB data file, followed by a description of the field
> // to be used, and optionally, a threshold to be applied to the
> field.
> // e.g. "sample.grb APCP/A3 gt0.00"
> //
> // Any NetCDF or GRIB file used must have the same grid dimensions
as
> the
> // data being verified.
> //
> // MET_BASE may be used in the path for the files above.
> //
> // e.g. mask_poly[] = [ "MET_BASE/data/poly/EAST.poly",
> // "poly_mask.ncf",
> // "sample.nc APCP",
> // "sample.grb HGT/Z0 gt100.0" ];
> //
> mask_poly[] = [];
>
> //
> // Specify the name of an ASCII file containing a space-separated
list
> of
> // station ID's at which to perform verification. Each station ID
> specified
> // is treated as an individual masking region.
> //
> // An empty list file name indicates that no station ID masks should
> be used.
> //
> // MET_BASE may be used in the path for the station ID mask file
name.
> //
> // e.g. mask_sid = "CONUS.stations";
> //
> mask_sid = "";
>
> //
> // Specify a comma-separated list of values for alpha to be used
when
> computing
> // confidence intervals. Values of alpha must be between 0 and 1.
> //
> // e.g. ci_alpha[] = [ 0.05, 0.10 ];
> //
> ci_alpha[] = [ 0.05 ];
>
> //
> // Specify the method to be used for computing bootstrap confidence
> intervals.
> // The value for this is interpreted as follows:
> // (0) Use the BCa interval method (computationally intensive)
> // (1) Use the percentile interval method
> //
> boot_interval = 1;
>
> //
> // Specify a proportion between 0 and 1 to define the replicate
sample
> size
> // to be used when computing percentile intervals. The replicate
> sample
> // size is set to boot_rep_prop * n, where n is the number of raw
data
> points.
> //
> // e.g boot_rep_prop = 0.80;
> //
> boot_rep_prop = 1.0;
>
> //
> // Specify the number of times each set of matched pair data should
be
> // resampled when computing bootstrap confidence intervals. A value
> of
> // zero disables the computation of bootstrap condifence intervals.
> //
> // e.g. n_boot_rep = 1000;
> //
> n_boot_rep = 1000;
>
> //
> // Specify the name of the random number generator to be used. See
> the MET
> // Users Guide for a list of possible random number generators.
> //
> boot_rng = "mt19937";
>
> //
> // Specify the seed value to be used when computing bootstrap
> confidence
> // intervals. If left unspecified, the seed will change for each run
> and
> // the computed bootstrap confidence intervals will not be
> reproducable.
> //
> boot_seed = "";
>
> //
> // Specify a comma-separated list of interpolation method(s) to be
> used
> // for comparing the forecast grid to the observation points. String
> values
> // are interpreted as follows:
> // MIN = Minimum in the neighborhood
> // MAX = Maximum in the neighborhood
> // MEDIAN = Median in the neighborhood
> // UW_MEAN = Unweighted mean in the neighborhood
> // DW_MEAN = Distance-weighted mean in the neighborhood
> // LS_FIT = Least-squares fit in the neighborhood
> //
> // In all cases, vertical interpolation is performed in the natural
> log
> // of pressure of the levels above and below the observation.
> //
> // e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
> //
> interp_method[] = [ "UW_MEAN" ];
>
> //
> // Specify a comma-separated list of box widths to be used by the
> // interpolation techniques listed above. A value of 1 indicates
that
> // the nearest neighbor approach should be used. For a value of n
> // greater than 1, the n*n grid points closest to the observation
> define
> // the neighborhood.
> //
> // e.g. interp_width = [ 1, 3, 5 ];
> //
> interp_width[] = [ 1 ];
>
> //
> // When interpolating, compute a ratio of the number of valid data
> points
> // to the total number of points in the neighborhood. If that ratio
> is
> // less than this threshold, do not include the observation. This
> // threshold must be between 0 and 1. Setting this threshold to 1
> will
> // require that each observation be surrounded by n*n valid forecast
> // points.
> //
> // e.g. interp_thresh = 1.0;
> //
> interp_thresh = 1.0;
>
> //
> // Specify flags to indicate the type of data to be output:
> // (1) STAT and FHO Text Files, Forecast, Hit, Observation Rates:
> // Total (TOTAL),
> // Forecast Rate (F_RATE),
> // Hit Rate (H_RATE),
> // Observation Rate (O_RATE)
> //
> // (2) STAT and CTC Text Files, Contingency Table Counts:
> // Total (TOTAL),
> // Forecast Yes and Observation Yes Count (FY_OY),
> // Forecast Yes and Observation No Count (FY_ON),
> // Forecast No and Observation Yes Count (FN_OY),
> // Forecast No and Observation No Count (FN_ON)
> //
> // (3) STAT and CTS Text Files, Contingency Table Scores:
> // Total (TOTAL),
> // Base Rate (BASER),
> // Forecast Mean (FMEAN),
> // Accuracy (ACC),
> // Frequency Bias (FBIAS),
> // Probability of Detecting Yes (PODY),
> // Probability of Detecting No (PODN),
> // Probability of False Detection (POFD),
> // False Alarm Ratio (FAR),
> // Critical Success Index (CSI),
> // Gilbert Skill Score (GSS),
> // Hanssen and Kuipers Discriminant (HK),
> // Heidke Skill Score (HSS),
> // Odds Ratio (ODDS),
> // NOTE: All statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (4) STAT and MCTC Text Files, NxN Multi-Category Contingency
> Table Counts:
> // Total (TOTAL),
> // Number of Categories (N_CAT),
> // Contingency Table Count columns repeated N_CAT*N_CAT
> times
> //
> // (5) STAT and MCTS Text Files, NxN Multi-Category Contingency
> Table Scores:
> // Total (TOTAL),
> // Number of Categories (N_CAT),
> // Accuracy (ACC),
> // Hanssen and Kuipers Discriminant (HK),
> // Heidke Skill Score (HSS),
> // Gerrity Score (GER),
> // NOTE: All statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (6) STAT and CNT Text Files, Statistics of Continuous Variables:
> // Total (TOTAL),
> // Forecast Mean (FBAR),
> // Forecast Standard Deviation (FSTDEV),
> // Observation Mean (OBAR),
> // Observation Standard Deviation (OSTDEV),
> // Pearson's Correlation Coefficient (PR_CORR),
> // Spearman's Rank Correlation Coefficient (SP_CORR),
> // Kendall Tau Rank Correlation Coefficient (KT_CORR),
> // Number of ranks compared (RANKS),
> // Number of tied ranks in the forecast field (FRANK_TIES),
> // Number of tied ranks in the observation field
> (ORANK_TIES),
> // Mean Error (ME),
> // Standard Deviation of the Error (ESTDEV),
> // Multiplicative Bias (MBIAS = FBAR - OBAR),
> // Mean Absolute Error (MAE),
> // Mean Squared Error (MSE),
> // Bias-Corrected Mean Squared Error (BCMSE),
> // Root Mean Squared Error (RMSE),
> // Percentiles of the Error (E10, E25, E50, E75, E90)
> // NOTE: Most statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (7) STAT and SL1L2 Text Files, Scalar Partial Sums:
> // Total (TOTAL),
> // Forecast Mean (FBAR),
> // = mean(f)
> // Observation Mean (OBAR),
> // = mean(o)
> // Forecast*Observation Product Mean (FOBAR),
> // = mean(f*o)
> // Forecast Squared Mean (FFBAR),
> // = mean(f^2)
> // Observation Squared Mean (OOBAR)
> // = mean(o^2)
> //
> // (8) STAT and SAL1L2 Text Files, Scalar Anomaly Partial Sums:
> // Total (TOTAL),
> // Forecast Anomaly Mean (FABAR),
> // = mean(f-c)
> // Observation Anomaly Mean (OABAR),
> // = mean(o-c)
> // Product of Forecast and Observation Anomalies Mean
> (FOABAR),
> // = mean((f-c)*(o-c))
> // Forecast Anomaly Squared Mean (FFABAR),
> // = mean((f-c)^2)
> // Observation Anomaly Squared Mean (OOABAR)
> // = mean((o-c)^2)
> //
> // (9) STAT and VL1L2 Text Files, Vector Partial Sums:
> // Total (TOTAL),
> // U-Forecast Mean (UFBAR),
> // = mean(uf)
> // V-Forecast Mean (VFBAR),
> // = mean(vf)
> // U-Observation Mean (UOBAR),
> // = mean(uo)
> // V-Observation Mean (VOBAR),
> // = mean(vo)
> // U-Product Plus V-Product (UVFOBAR),
> // = mean(uf*uo+vf*vo)
> // U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
> // = mean(uf^2+vf^2)
> // U-Observation Squared Plus V-Observation Squared
> (UVOOBAR)
> // = mean(uo^2+vo^2)
> //
> // (10) STAT and VAL1L2 Text Files, Vector Anomaly Partial Sums:
> // U-Forecast Anomaly Mean (UFABAR),
> // = mean(uf-uc)
> // V-Forecast Anomaly Mean (VFABAR),
> // = mean(vf-vc)
> // U-Observation Anomaly Mean (UOABAR),
> // = mean(uo-uc)
> // V-Observation Anomaly Mean (VOABAR),
> // = mean(vo-vc)
> // U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
> // = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
> // U-Forecast Anomaly Squared Plus V-Forecast Anomaly
> Squared (UVFFABAR),
> // = mean((uf-uc)^2+(vf-vc)^2)
> // U-Observation Anomaly Squared Plus V-Observation Anomaly
> Squared (UVOOABAR)
> // = mean((uo-uc)^2+(vo-vc)^2)
> //
> // (11) STAT and PCT Text Files, Nx2 Probability Contingency Table
> Counts:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Row Observation Yes Count (OY_i),
> // Row Observation No Count (ON_i),
> // NOTE: Previous 3 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (12) STAT and PSTD Text Files, Nx2 Probability Contingency Table
> Scores:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Base Rate (BASER) with confidence interval limits,
> // Reliability (RELIABILITY),
> // Resolution (RESOLUTION),
> // Uncertainty (UNCERTAINTY),
> // Area Under the ROC Curve (ROC_AUC),
> // Brier Score (BRIER) with confidence interval limits,
> // Probability Threshold Value (THRESH_i)
> // NOTE: Previous column repeated for each probability
> threshold.
> //
> // (13) STAT and PJC Text Files, Joint/Continuous Statistics of
> // Probabilistic Variables:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Observation Yes Count Divided by Total (OY_TP_i),
> // Observation No Count Divided by Total (ON_TP_i),
> // Calibration (CALIBRATION_i),
> // Refinement (REFINEMENT_i),
> // Likelikhood (LIKELIHOOD_i),
> // Base Rate (BASER_i),
> // NOTE: Previous 7 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (14) STAT and PRC Text Files, ROC Curve Points for
> // Probabilistic Variables:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Probability of Detecting Yes (PODY_i),
> // Probability of False Detection (POFD_i),
> // NOTE: Previous 3 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (15) STAT and MPR Text Files, Matched Pair Data:
> // Total (TOTAL),
> // Index (INDEX),
> // Observation Station ID (OBS_SID),
> // Observation Latitude (OBS_LAT),
> // Observation Longitude (OBS_LON),
> // Observation Level (OBS_LVL),
> // Observation Elevation (OBS_ELV),
> // Forecast Value (FCST),
> // Observation Value (OBS),
> // Climatological Value (CLIMO)
> //
> // In the expressions above, f are forecast values, o are observed
> values,
> // and c are climatological values.
> //
> // Values for these flags are interpreted as follows:
> // (0) Do not generate output of this type
> // (1) Write output to a STAT file
> // (2) Write output to a STAT file and a text file
> //
> output_flag[] = [ 2, 2, 2, 0, 0, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
>
> //
> // Flag to indicate whether Kendall's Tau and Spearman's Rank
> Correlation
> // Coefficients should be computed. Computing them over large
> datasets is
> // computationally intensive and slows down the runtime execution
> significantly.
> // (0) Do not compute these correlation coefficients
> // (1) Compute these correlation coefficients
> //
> rank_corr_flag = 1;
>
> //
> // Specify the GRIB Table 2 parameter table version number to be
used
> // for interpreting GRIB codes.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> //
> grib_ptv = 2;
>
> //
> // Directory where temporary files should be written.
> //
> tmp_dir = "/tmp";
>
> //
> // Prefix to be used for the output file names.
> //
> output_prefix = "";
>
> //
> // Indicate a version number for the contents of this configuration
> file.
> // The value should generally not be modified.
> //
> version = "V3.0";
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Geeta Geeta
> Time: Mon Apr 16 03:08:18 2012
>
>
> Dear Met_help,
> I have done the corrections in my configuration file and it works
fine
> now.
> Thanks for Ur support.
>
> geeta
>
>
> > Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem, error with point-stat
> > From: met_help at ucar.edu
> > To: geeta124 at hotmail.com
> > Date: Thu, 12 Apr 2012 09:03:20 -0600
> >
> > Geeta,
> >
> > First of all, if you want to send data files to me that are larger
> than 2MB, please upload them to our FTP site using
> > the instructions found here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
> >
> > I found several problems, but they were all easy to find and fix.
> There was a config file syntax error in the version
> > of PointStatConfig_rain that you sent me. You must use double
quotes
> around values in the obs_field setting:
> >
> > obs_field[] = [ "61/A24" ];
> >
> > Next, I recommend using a verbosity setting of 3 or higher when
you
> start using point_stat. This will cause point_stat
> > to print out the reason why points are rejected. One reason is
that
> there is a mis-match between the valid time of your
> > model data and the valid time of your observation data. The
> following ncdump commands show the discrepancy:
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
> > ...
> > APCP_24:valid_time = "20110610_030000" ;
> >
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
> precip24_2011061000.nc
> > ...
> > hdr_vld =
> > "20110610_083000",
> > ...
> >
> > You can tell point_stat to allow this discrepancy using the
> following two settings in PointStatConfig_rain:
> >
> > beg_ds = -21600;
> > end_ds = 21600;
> >
> > Next, the message type of your observations is ADPSFC, which you
can
> see using the following ncdump command:
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
> precip24_2011061000.nc
> > ...
> > hdr_typ =
> > "ADPSFC",
> > ...
> >
> > In PointStatConfig_rain, the following setting must be used to
> specify this message type:
> >
> > message_type[] = [ "ADPSFC" ];
> >
> > Once these changes were made, I ran point_stat with your model and
> observation data and the result was 26 matched pairs.
> > I attached the output below and the PointStatConfig_rain file that
I
> used. Please let me know if you have any more
> > questions.
> >
> > One last thing, I think you might have a date typo in your
> precip24_2011061000.ascii file on line 2. Should the date
> > string "20110613_083000" actually be "20110610_083000" like all
the
> others?
> >
> > Thanks,
> >
> > Paul
> >
> >
> > $ point_stat model24_day.nc precip24_2011061000.nc
> PointStatConfig_rain -v 3
> > DEBUG 3: It's a Mercator grid...
> > DEBUG 3: Mercator Data:
> > DEBUG 3: lat_ll = 9.097
> > DEBUG 3: lon_ll = -73.943
> > DEBUG 3: lat_ur = 21.711
> > DEBUG 3: lon_ur = -87.057
> > DEBUG 3: ny = 163
> > DEBUG 3: nx = 163
> > DEBUG 1: Default Config File:
>
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
> > DEBUG 1: User Config File: PointStatConfig_rain
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744072433272399
> > DEBUG 1: Forecast File: model24_day.nc
> > DEBUG 1: Climatology File: none
> > DEBUG 1: Observation File: precip24_2011061000.nc
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for APCP_24(*,*).
> > DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0
climatology
> levels.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Searching 27 observations from 27 messages.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation
type
> ADPSFC, over region FULL, for interpolation method
> > UW_MEAN(1), using 26 pairs.
> > DEBUG 3: Number of matched pairs = 26
> > DEBUG 3: Observations processed = 27
> > DEBUG 3: Rejected: GRIB code = 0
> > DEBUG 3: Rejected: valid time = 1
> > DEBUG 3: Rejected: bad obs value = 0
> > DEBUG 3: Rejected: off the grid = 0
> > DEBUG 3: Rejected: level mismatch = 0
> > DEBUG 3: Rejected: message type = 0
> > DEBUG 3: Rejected: masking region = 0
> > DEBUG 3: Rejected: bad fcst value = 0
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Computing Scalar Partial Sums.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_fho.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_ctc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_cts.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_cnt.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sl1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sal1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_vl1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_val1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pct.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pstd.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pjc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_prc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_mpr.txt
> >
> >
> >
> > On 04/11/2012 11:53 PM, Geeta Geeta wrote:
> > >
> > > Dear met_help,
> > > I am using point stat tool for Rainfall verification in my
region
> Western part of Indian peninsula (9-21N/74-87E).
> > >
> > >
> > > 1. for this, I have converted the ascii observations to netcdf
> using the ascii2nc tool (O/p file precip24_2011061000.nc)
> > > 2. I have converted the model o/p in .nc format to grib1 format
> using WPP.
> > > 3. I have then used pcpcombine to get the desired RF/output ie.
by
> subtracting the RF of 27th hour from the 4th hour. (model24_day.nc
is
> attached)
> > >
> > >
> > > 4. now i am running
> > > -bash-3.2$ cd ..
> > > -bash-3.2$ cd verf-met/
> > > /oprn/model/wrf3/utils/met/METv3.0/verf-met
> > > -bash-3.2$ ../bin/point_stat model24_day.nc
precip24_2011061000.nc
> ~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=801800783
> > > Forecast File: model24_day.nc
> > > Climatology File: none
> > > Configuration File:
> /oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > > Observation File: precip24_2011061000.nc
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Reading records for APCP_24(*,*).
> > > For APCP_24(*,*) found 1 forecast levels and 0 climatology
levels.
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Searching 27 observations from 27 PrepBufr messages.
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Processing APCP_24(*,*) versus APCP_24/A24, for observation type
> ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
0
> pairs.
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
> > > -bash-3.2$
> > >
> > >
> > >
> > > My problem is that there is no data in the output files, when i
> open them, it is only headers.
> > > Kindly help.
> > > I have plotted the grib1 files (after WPP by using grib2ctl.pl)
to
> see the domain and if all the stations are falling within this
domain,
> So all looks okay to me
> > > geeta
> > >
> > >
> > >
> > >> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem
> > >> From: met_help at ucar.edu
> > >> To: geeta124 at hotmail.com
> > >> Date: Wed, 11 Apr 2012 09:24:51 -0600
> > >>
> > >> Geeta,
> > >>
> > >> I created a met_help ticket for you. Please respond to this
email
> and describe the MET problem you are having in detail.
> > >>
> > >> Thanks,
> > >>
> > >> Paul
> > >>
> > >>
> > >> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> > >>> Greetings,
> > >>>
> > >>> This message has been automatically generated in response to
the
> > >>> creation of a trouble ticket regarding:
> > >>> "Geeta met_help problem",
> > >>> a summary of which appears below.
> > >>>
> > >>> There is no need to reply to this message right now. Your
ticket
> has been
> > >>> assigned an ID of [rt.rap.ucar.edu #55790].
> > >>>
> > >>> Please include the string:
> > >>>
> > >>> [rt.rap.ucar.edu #55790]
> > >>>
> > >>> in the subject line of all future correspondence about this
> issue. To do so,
> > >>> you may reply to this message.
> > >>>
> > >>> Thank you,
> > >>> met_help at ucar.edu
> > >>>
> > >>>
>
-------------------------------------------------------------------------
> > >>> Dear administrator,
> > >>> I wish to seek help from the met_help, and I am sending my
> message for the last 3-4 days. But everytime I am told that there is
> some error that has occurred.
> > >>> I donot know whtheer my meage is reaching across or not? I
have
> already registered but my query is not answered.
> > >>>
> > >>> kindly help me
> > >>> Regards
> > >>> geeta
> > >>>
> > >>> Dr G. Agnihotri,
> > >>> Met Centre, IMD, Palace Road,
> > >>> Bangalore, India.
> > >>> 080-22244419
> > >>>
> > >>
> > >>
> > >
> >
> >
>
> ------------------------------------------------
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

Creating NetCDF Observation file: ../myascii2nc/precip24_2011061000.nc
Reading ASCII Observation file: precip24_2011061000.ascii
20110600
/oprn/model/wrf3/WRFV3/postprd
-pcpdir                  WRFPRS_d01.000           WRFPRS_d01.001
WRFPRS_d01.001.idx       WRFPRS_d01.002           WRFPRS_d01.003
WRFPRS_d01.004           WRFPRS_d01.005           WRFPRS_d01.006
WRFPRS_d01.007           WRFPRS_d01.008           WRFPRS_d01.009
WRFPRS_d01.010           WRFPRS_d01.011           WRFPRS_d01.012
WRFPRS_d01.013           WRFPRS_d01.014           WRFPRS_d01.015
WRFPRS_d01.016           WRFPRS_d01.017           WRFPRS_d01.018
WRFPRS_d01.019           WRFPRS_d01.020           WRFPRS_d01.021
WRFPRS_d01.022           WRFPRS_d01.023           WRFPRS_d01.024
WRFPRS_d01.025           WRFPRS_d01.026           WRFPRS_d01.027
WRFPRS_d01.028           WRFPRS_d01.029           WRFPRS_d01.030
WRFPRS_d01.031           WRFPRS_d01.032           WRFPRS_d01.033
WRFPRS_d01.034           WRFPRS_d01.035           WRFPRS_d01.036
WRFPRS_d01.037           WRFPRS_d01.038           WRFPRS_d01.039
WRFPRS_d01.040           WRFPRS_d01.041           WRFPRS_d01.042
WRFPRS_d01.043           WRFPRS_d01.044           WRFPRS_d01.045
WRFPRS_d01.046           WRFPRS_d01.047           WRFPRS_d01.048
WRFPRS_d01.049           WRFPRS_d01.050           WRFPRS_d01.051
copygb_hwrf.txt          eta_micro_lookup.dat     fort.110
fort.14                  griddef.out              itag
out1                     run_wrfpost              run_wrfpost-
donotchange
temp.ctl                 tempfile                 testfile
wrf_cntrl.parm           wrfpost_d01.000.out      wrfpost_d01.001.out
wrfpost_d01.002.out      wrfpost_d01.003.out      wrfpost_d01.004.out
wrfpost_d01.005.out      wrfpost_d01.006.out      wrfpost_d01.007.out
wrfpost_d01.008.out      wrfpost_d01.009.out      wrfpost_d01.010.out
wrfpost_d01.011.out      wrfpost_d01.012.out      wrfpost_d01.013.out
wrfpost_d01.014.out      wrfpost_d01.015.out      wrfpost_d01.016.out
wrfpost_d01.017.out      wrfpost_d01.018.out      wrfpost_d01.019.out
wrfpost_d01.020.out      wrfpost_d01.021.out      wrfpost_d01.022.out
wrfpost_d01.023.out      wrfpost_d01.024.out      wrfpost_d01.025.out
wrfpost_d01.026.out      wrfpost_d01.027.out      wrfpost_d01.028.out
wrfpost_d01.029.out      wrfpost_d01.030.out      wrfpost_d01.031.out
wrfpost_d01.032.out      wrfpost_d01.033.out      wrfpost_d01.034.out
wrfpost_d01.035.out      wrfpost_d01.036.out      wrfpost_d01.037.out
wrfpost_d01.038.out      wrfpost_d01.039.out      wrfpost_d01.040.out
wrfpost_d01.041.out      wrfpost_d01.042.out      wrfpost_d01.043.out
wrfpost_d01.044.out      wrfpost_d01.045.out      wrfpost_d01.046.out
wrfpost_d01.047.out      wrfpost_d01.048.out      wrfpost_d01.049.out
wrfpost_d01.050.out      wrfpost_d01.051.out      wrfpost_d01.052.out
wrfprs_d01.000           wrfprs_d01.001           wrfprs_d01.002
wrfprs_d01.003           wrfprs_d01.004           wrfprs_d01.005
wrfprs_d01.006           wrfprs_d01.007           wrfprs_d01.008
wrfprs_d01.009           wrfprs_d01.010           wrfprs_d01.011
wrfprs_d01.012           wrfprs_d01.013           wrfprs_d01.014
wrfprs_d01.015           wrfprs_d01.016           wrfprs_d01.017
wrfprs_d01.018           wrfprs_d01.019           wrfprs_d01.020
wrfprs_d01.021           wrfprs_d01.022           wrfprs_d01.023
wrfprs_d01.024           wrfprs_d01.025           wrfprs_d01.026
wrfprs_d01.027           wrfprs_d01.028           wrfprs_d01.029
wrfprs_d01.030           wrfprs_d01.031           wrfprs_d01.032
wrfprs_d01.033           wrfprs_d01.034           wrfprs_d01.035
wrfprs_d01.036           wrfprs_d01.037           wrfprs_d01.038
wrfprs_d01.039           wrfprs_d01.040           wrfprs_d01.041
wrfprs_d01.042           wrfprs_d01.043           wrfprs_d01.044
wrfprs_d01.045           wrfprs_d01.046           wrfprs_d01.047
wrfprs_d01.048           wrfprs_d01.049           wrfprs_d01.050
wrfprs_d01.051
NEWDATE 2011060900
YY 2011
-rwxr-xr-x    1 model    staff      11408030 Apr 19 14:37
WRFPRS_d01.000
NEWDATE 2011060901
YY 2011
-rwxr-xr-x    1 model    staff      12643498 Apr 19 14:37
WRFPRS_d01.001
NEWDATE 2011060902
YY 2011
-rwxr-xr-x    1 model    staff      12643500 Apr 19 14:37
WRFPRS_d01.002
NEWDATE 2011060903
YY 2011
-rwxr-xr-x    1 model    staff      12600320 Apr 19 14:37
WRFPRS_d01.003
NEWDATE 2011060904
YY 2011
-rwxr-xr-x    1 model    staff      12590344 Apr 19 14:38
WRFPRS_d01.004
NEWDATE 2011060905
YY 2011
-rwxr-xr-x    1 model    staff      12623562 Apr 19 14:38
WRFPRS_d01.005
NEWDATE 2011060906
YY 2011
-rwxr-xr-x    1 model    staff      12583710 Apr 19 14:38
WRFPRS_d01.006
NEWDATE 2011060907
YY 2011
-rwxr-xr-x    1 model    staff      12537212 Apr 19 14:38
WRFPRS_d01.007
NEWDATE 2011060908
YY 2011
-rwxr-xr-x    1 model    staff      12587026 Apr 19 14:38
WRFPRS_d01.008
NEWDATE 2011060909
YY 2011
-rwxr-xr-x    1 model    staff      12630192 Apr 19 14:39
WRFPRS_d01.009
NEWDATE 2011060910
YY 2011
-rwxr-xr-x    1 model    staff      12643474 Apr 19 14:39
WRFPRS_d01.010
NEWDATE 2011060911
YY 2011
-rwxr-xr-x    1 model    staff      12580376 Apr 19 14:39
WRFPRS_d01.011
NEWDATE 2011060912
YY 2011
-rwxr-xr-x    1 model    staff      12593662 Apr 19 14:39
WRFPRS_d01.012
NEWDATE 2011060913
YY 2011
-rwxr-xr-x    1 model    staff      12613590 Apr 19 14:39
WRFPRS_d01.013
NEWDATE 2011060914
YY 2011
-rwxr-xr-x    1 model    staff      12520594 Apr 19 14:40
WRFPRS_d01.014
NEWDATE 2011060915
YY 2011
-rwxr-xr-x    1 model    staff      12530568 Apr 19 14:40
WRFPRS_d01.015
NEWDATE 2011060916
YY 2011
-rwxr-xr-x    1 model    staff      12587014 Apr 19 14:40
WRFPRS_d01.016
NEWDATE 2011060917
YY 2011
-rwxr-xr-x    1 model    staff      12606944 Apr 19 14:40
WRFPRS_d01.017
NEWDATE 2011060918
YY 2011
-rwxr-xr-x    1 model    staff      12567086 Apr 19 14:40
WRFPRS_d01.018
NEWDATE 2011060919
YY 2011
-rwxr-xr-x    1 model    staff      12543850 Apr 19 14:41
WRFPRS_d01.019
NEWDATE 2011060920
YY 2011
-rwxr-xr-x    1 model    staff      12504000 Apr 19 14:41
WRFPRS_d01.020
NEWDATE 2011060921
YY 2011
-rwxr-xr-x    1 model    staff      12487396 Apr 19 14:41
WRFPRS_d01.021
NEWDATE 2011060922
YY 2011
-rwxr-xr-x    1 model    staff      12563782 Apr 19 14:41
WRFPRS_d01.022
NEWDATE 2011060923
YY 2011
-rwxr-xr-x    1 model    staff      12490712 Apr 19 14:41
WRFPRS_d01.023
NEWDATE 2011061000
YY 2011
-rwxr-xr-x    1 model    staff      12560452 Apr 19 14:42
WRFPRS_d01.024
NEWDATE 2011061001
YY 2011
-rwxr-xr-x    1 model    staff      12550486 Apr 19 14:42
WRFPRS_d01.025
NEWDATE 2011061002
YY 2011
-rwxr-xr-x    1 model    staff      12596992 Apr 19 14:42
WRFPRS_d01.026
NEWDATE 2011061003
YY 2011
-rwxr-xr-x    1 model    staff      12623568 Apr 19 14:42
WRFPRS_d01.027
NEWDATE 2011061004
YY 2011
-rwxr-xr-x    1 model    staff      12620248 Apr 19 14:42
WRFPRS_d01.028
NEWDATE 2011061005
YY 2011
-rwxr-xr-x    1 model    staff      12626896 Apr 19 14:43
WRFPRS_d01.029
NEWDATE 2011061006
YY 2011
-rwxr-xr-x    1 model    staff      12613612 Apr 19 14:43
WRFPRS_d01.030
NEWDATE 2011061007
YY 2011
-rwxr-xr-x    1 model    staff      12583714 Apr 19 14:43
WRFPRS_d01.031
NEWDATE 2011061008
YY 2011
-rwxr-xr-x    1 model    staff      12573740 Apr 19 14:43
WRFPRS_d01.032
NEWDATE 2011061009
YY 2011
-rwxr-xr-x    1 model    staff      12603626 Apr 19 14:43
WRFPRS_d01.033
NEWDATE 2011061010
YY 2011
-rwxr-xr-x    1 model    staff      12613596 Apr 19 14:44
WRFPRS_d01.034
NEWDATE 2011061011
YY 2011
-rwxr-xr-x    1 model    staff      12577060 Apr 19 14:44
WRFPRS_d01.035
NEWDATE 2011061012
YY 2011
-rwxr-xr-x    1 model    staff      12616916 Apr 19 14:44
WRFPRS_d01.036
NEWDATE 2011061013
YY 2011
-rwxr-xr-x    1 model    staff      12610266 Apr 19 14:44
WRFPRS_d01.037
NEWDATE 2011061014
YY 2011
-rwxr-xr-x    1 model    staff      12550486 Apr 19 14:44
WRFPRS_d01.038
NEWDATE 2011061015
YY 2011
-rwxr-xr-x    1 model    staff      12527246 Apr 19 14:45
WRFPRS_d01.039
NEWDATE 2011061016
YY 2011
-rwxr-xr-x    1 model    staff      12530566 Apr 19 14:45
WRFPRS_d01.040
NEWDATE 2011061017
YY 2011
-rwxr-xr-x    1 model    staff      12567092 Apr 19 14:45
WRFPRS_d01.041
NEWDATE 2011061018
YY 2011
-rwxr-xr-x    1 model    staff      12517278 Apr 19 14:45
WRFPRS_d01.042
NEWDATE 2011061019
YY 2011
-rwxr-xr-x    1 model    staff      12500672 Apr 19 14:45
WRFPRS_d01.043
NEWDATE 2011061020
YY 2011
-rwxr-xr-x    1 model    staff      12527242 Apr 19 14:46
WRFPRS_d01.044
NEWDATE 2011061021
YY 2011
-rwxr-xr-x    1 model    staff      12550488 Apr 19 14:46
WRFPRS_d01.045
NEWDATE 2011061022
YY 2011
-rwxr-xr-x    1 model    staff      12583708 Apr 19 14:46
WRFPRS_d01.046
NEWDATE 2011061023
YY 2011
-rwxr-xr-x    1 model    staff      12547174 Apr 19 14:46
WRFPRS_d01.047
NEWDATE 2011061100
YY 2011
-rwxr-xr-x    1 model    staff      12670052 Apr 19 14:46
WRFPRS_d01.048
NEWDATE 2011061101
YY 2011
-rwxr-xr-x    1 model    staff      12656766 Apr 19 14:47
WRFPRS_d01.049
NEWDATE 2011061102
YY 2011
-rwxr-xr-x    1 model    staff      12650132 Apr 19 14:47
WRFPRS_d01.050
NEWDATE 2011061103
YY 2011
-rwxr-xr-x    1 model    staff      12610284 Apr 19 14:47
WRFPRS_d01.051
Thu Apr 19 14:47:36 IST 2012
End of Output Job
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.027
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.003
Performing subtraction command.
Writing output file: ./myoutput/pcp_combine/24hr/model24_2011060900.nc
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.051
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.027
Performing subtraction command.
Writing output file: ./myoutput/pcp_combine/48hr/model48_2011060900.nc
Forecast File:
/oprn/model/wrf3/utils/met/METv3.0/myoutput/pcp_combine/24hr/model24_2011060900.nc
Climatology File: none
Configuration File:
/oprn/model/wrf3/utils/met/METv3.0/myscripts/PointStatConfig_rain
Observation File:
/oprn/model/wrf3/utils/met/METv3.0/myascii2nc/precip24_2011061000.nc

--------------------------------------------------------------------------------

Reading records for APCP_24(*,*).
For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.

--------------------------------------------------------------------------------

Searching 29 observations from 29 PrepBufr messages.

--------------------------------------------------------------------------------

Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
27 pairs.
Computing Categorical Statistics.
Computing Continuous Statistics.
Computing Scalar Partial Sums.

--------------------------------------------------------------------------------

Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V.stat
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_fho.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_ctc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_cts.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_cnt.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_sl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_sal1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_vl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_val1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pct.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pstd.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pjc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_prc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_mpr.txt
Forecast File: ./myoutput/pcp_combine/48hr/model48_2011060900.nc
Climatology File: none
Configuration File: ./myscripts/PointStatConfig_rain
Observation File: ./myascii2nc/precip24_2011061100.nc

--------------------------------------------------------------------------------

Reading records for APCP_24(*,*).
For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.

--------------------------------------------------------------------------------

Searching 27 observations from 27 PrepBufr messages.

--------------------------------------------------------------------------------

Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
27 pairs.
Computing Categorical Statistics.
Computing Continuous Statistics.
Computing Scalar Partial Sums.

--------------------------------------------------------------------------------

Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V.stat
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_fho.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_ctc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_cts.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_cnt.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_sl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_sal1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_vl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_val1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pct.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pstd.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pjc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_prc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_mpr.txt

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

////////////////////////////////////////////////////////////////////////////////
//
// Default point_stat configuration file
//
////////////////////////////////////////////////////////////////////////////////

//
// Specify a name to designate the model being verified.  This name
will be
// written to the second column of the ASCII output generated.
//
model = "WRF";

//
// Beginning and ending time offset values in seconds for observations
// to be used.  These time offsets are defined in reference to the
// forecast valid time, v.  Observations with a valid time falling in
the
// window [v+beg_ds, v+end_ds] will be used.
// These selections are overridden by the command line arguments
// -obs_valid_beg and -obs_valid_end.
//
beg_ds = -21600;
end_ds =  21600;

//
// Specify a comma-separated list of fields to be verified.  The
forecast and
// observation fields may be specified separately.  If the obs_field
parameter
// is left blank, it will default to the contents of fcst_field.
//
// Each field is specified as a GRIB code or abbreviation followed by
an
// accumulation or vertical level indicator for GRIB files or as a
variable name
// followed by a list of dimensions for NetCDF files output from
p_interp or MET.
//
// Specifying verification fields for GRIB files:
//    GC/ANNN for accumulation interval NNN
//    GC/ZNNN for vertical level NNN
//    GC/ZNNN-NNN for a range of vertical levels (MSL or AGL)
//    GC/PNNN for pressure level NNN in hPa
//    GC/PNNN-NNN for a range of pressure levels in hPa
//    GC/LNNN for a generic level type
//    GC/RNNN for a specific GRIB record number
//    Where GC is the number of or abbreviation for the grib code
//    to be verified.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
// Specifying verification fields for NetCDF files:
//    var_name(i,...,j,*,*) for a single field
//    var_name(i-j,*,*) for a range of fields
//    Where var_name is the name of the NetCDF variable,
//    and i,...,j specifies fixed dimension values,
//    and i-j specifies a range of values for a single dimension,
//    and *,* specifies the two dimensions for the gridded field.
//
//    NOTE: To verify winds as vectors rather than scalars,
//          specify UGRD (or 33) followed by VGRD (or 34) with the
//          same level values.
//
//    NOTE: To process a probability field, add "/PROB", such as
"POP/Z0/PROB".
//
// e.g. fcst_field[] = [ "SPFH/P500", "TMP/P500" ]; for a GRIB input
// e.g. fcst_field[] = [ "QVAPOR(0,5,*,*)", "TT(0,5,*,*)" ]; for
NetCDF input
//
fcst_field[] = [ "APCP_24(*,*)"];
obs_field[]  = ["61/A24"];

//
// Specify a comma-separated list of groups of thresholds to be
applied to the
// fields listed above.  Thresholds for the forecast and observation
fields
// may be specified separately.  If the obs_thresh parameter is left
blank,
// it will default to the contents of fcst_thresh.
//
// At least one threshold must be provided for each field listed
above.  The
// lengths of the "fcst_field" and "fcst_thresh" arrays must match, as
must
// lengths of the "obs_field" and "obs_thresh" arrays.  To apply
multiple
// thresholds to a field, separate the threshold values with a space.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
//    'lt' for less than     'le' for less than or equal to
//    'eq' for equal to      'ne' for not equal to
//    'gt' for greater than  'ge' for greater than or equal to
//
// NOTE: Thresholds for probabilities must begin with 0.0, end with
1.0,
//       and be preceeded by "ge".
//
// e.g. fcst_thresh[] = [ "gt80", "gt273" ];
//
fcst_thresh[] = [ "gt0.0"];
obs_thresh[]  = [ "gt0.0"];

//
// Specify a comma-separated list of thresholds to be used when
computing
// VL1L2 and VAL1L2 partial sums for winds.  The thresholds are
applied to the
// wind speed values derived from each U/V pair.  Only those U/V pairs
which meet
// the wind speed threshold criteria are retained.  If the
obs_wind_thresh
// parameter is left blank, it will default to the contents of
fcst_wind_thresh.
//
// To apply multiple wind speed thresholds, separate the threshold
values with a
// space.  Use "NA" to indicate that no wind speed threshold should be
applied.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
//    'lt' for less than     'le' for less than or equal to
//    'eq' for equal to      'ne' for not equal to
//    'gt' for greater than  'ge' for greater than or equal to
//    'NA' for no threshold
//
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
//
fcst_wind_thresh[] = [ "NA" ];
obs_wind_thresh[]  = [];

//
// Specify a comma-separated list of PrepBufr message types with which
// to perform the verification.  Statistics will be computed
separately
// for each message type specified.  At least one PrepBufr message
type
// must be provided.
// List of valid message types:
//    ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
//    MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
//    SFCSHP SPSSMI SYNDAT VADWND
//    ANYAIR (= AIRCAR, AIRCFT)
//    ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
//    ONLYSF (= ADPSFC, SFCSHP)
//
http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
//
// e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
//
message_type[] = [ "ADPSFC" ];

//
// Specify a comma-separated list of grids to be used in masking the
data over
// which to perform scoring.  An empty list indicates that no masking
grid
// should be performed.  The standard NCEP grids are named "GNNN"
where NNN
// indicates the three digit grid number.  Enter "FULL" to score over
the
// entire domain.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
//
// e.g. mask_grid[] = [ "FULL" ];
//
mask_grid[] = [ "FULL" ];

//
// Specify a comma-separated list of masking regions to be applied.
// An empty list indicates that no additional masks should be used.
// The masking regions may be defined in one of 4 ways:
//
// (1) An ASCII file containing a lat/lon polygon.
//     Latitude in degrees north and longitude in degrees east.
//     By default, the first and last polygon points are connected.
//     e.g. "MET_BASE/data/poly/EAST.poly" which consists of n points:
//          "poly_name lat1 lon1 lat2 lon2... latn lonn"
//
// (2) The NetCDF output of the gen_poly_mask tool.
//
// (3) A NetCDF data file, followed by the name of the NetCDF variable
//     to be used, and optionally, a threshold to be applied to the
field.
//     e.g. "sample.nc var_name gt0.00"
//
// (4) A GRIB data file, followed by a description of the field
//     to be used, and optionally, a threshold to be applied to the
field.
//     e.g. "sample.grb APCP/A3 gt0.00"
//
// Any NetCDF or GRIB file used must have the same grid dimensions as
the
// data being verified.
//
// MET_BASE may be used in the path for the files above.
//
// e.g. mask_poly[] = [ "MET_BASE/data/poly/EAST.poly",
//                      "poly_mask.ncf",
//                      "sample.nc APCP",
//                      "sample.grb HGT/Z0 gt100.0" ];
//
mask_poly[] = [];

//
// Specify the name of an ASCII file containing a space-separated list
of
// station ID's at which to perform verification.  Each station ID
specified
// is treated as an individual masking region.
//
// An empty list file name indicates that no station ID masks should
be used.
//
// MET_BASE may be used in the path for the station ID mask file name.
//
// e.g. mask_sid = "CONUS.stations";
//
mask_sid = "";

//
// Specify a comma-separated list of values for alpha to be used when
computing
// confidence intervals.  Values of alpha must be between 0 and 1.
//
// e.g. ci_alpha[] = [ 0.05, 0.10 ];
//
ci_alpha[] = [ 0.05 ];

//
// Specify the method to be used for computing bootstrap confidence
intervals.
// The value for this is interpreted as follows:
//    (0) Use the BCa interval method (computationally intensive)
//    (1) Use the percentile interval method
//
boot_interval = 1;

//
// Specify a proportion between 0 and 1 to define the replicate sample
size
// to be used when computing percentile intervals.  The replicate
sample
// size is set to boot_rep_prop * n, where n is the number of raw data
points.
//
// e.g boot_rep_prop = 0.80;
//
boot_rep_prop = 1.0;

//
// Specify the number of times each set of matched pair data should be
// resampled when computing bootstrap confidence intervals.  A value
of
// zero disables the computation of bootstrap condifence intervals.
//
// e.g. n_boot_rep = 1000;
//
n_boot_rep = 1000;

//
// Specify the name of the random number generator to be used.  See
the MET
// Users Guide for a list of possible random number generators.
//
boot_rng = "mt19937";

//
// Specify the seed value to be used when computing bootstrap
confidence
// intervals.  If left unspecified, the seed will change for each run
and
// the computed bootstrap confidence intervals will not be
reproducable.
//
boot_seed = "";

//
// Specify a comma-separated list of interpolation method(s) to be
used
// for comparing the forecast grid to the observation points.  String
values
// are interpreted as follows:
//    MIN     = Minimum in the neighborhood
//    MAX     = Maximum in the neighborhood
//    MEDIAN  = Median in the neighborhood
//    UW_MEAN = Unweighted mean in the neighborhood
//    DW_MEAN = Distance-weighted mean in the neighborhood
//    LS_FIT  = Least-squares fit in the neighborhood
//
// In all cases, vertical interpolation is performed in the natural
log
// of pressure of the levels above and below the observation.
//
// e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
//
interp_method[] = [ "UW_MEAN" ];

//
// Specify a comma-separated list of box widths to be used by the
// interpolation techniques listed above.  A value of 1 indicates that
// the nearest neighbor approach should be used.  For a value of n
// greater than 1, the n*n grid points closest to the observation
define
// the neighborhood.
//
// e.g. interp_width = [ 1, 3, 5 ];
//
interp_width[] = [ 1 ];

//
// When interpolating, compute a ratio of the number of valid data
points
// to the total number of points in the neighborhood.  If that ratio
is
// less than this threshold, do not include the observation.  This
// threshold must be between 0 and 1.  Setting this threshold to 1
will
// require that each observation be surrounded by n*n valid forecast
// points.
//
// e.g. interp_thresh = 1.0;
//
interp_thresh = 1.0;

//
// Specify flags to indicate the type of data to be output:
//    (1) STAT and FHO Text Files, Forecast, Hit, Observation Rates:
//           Total (TOTAL),
//           Forecast Rate (F_RATE),
//           Hit Rate (H_RATE),
//           Observation Rate (O_RATE)
//
//    (2) STAT and CTC Text Files, Contingency Table Counts:
//           Total (TOTAL),
//           Forecast Yes and Observation Yes Count (FY_OY),
//           Forecast Yes and Observation No Count (FY_ON),
//           Forecast No and Observation Yes Count (FN_OY),
//           Forecast No and Observation No Count (FN_ON)
//
//    (3) STAT and CTS Text Files, Contingency Table Scores:
//           Total (TOTAL),
//           Base Rate (BASER),
//           Forecast Mean (FMEAN),
//           Accuracy (ACC),
//           Frequency Bias (FBIAS),
//           Probability of Detecting Yes (PODY),
//           Probability of Detecting No (PODN),
//           Probability of False Detection (POFD),
//           False Alarm Ratio (FAR),
//           Critical Success Index (CSI),
//           Gilbert Skill Score (GSS),
//           Hanssen and Kuipers Discriminant (HK),
//           Heidke Skill Score (HSS),
//           Odds Ratio (ODDS),
//           NOTE: All statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (4) STAT and MCTC Text Files, NxN Multi-Category Contingency
Table Counts:
//           Total (TOTAL),
//           Number of Categories (N_CAT),
//           Contingency Table Count columns repeated N_CAT*N_CAT
times
//
//    (5) STAT and MCTS Text Files, NxN Multi-Category Contingency
Table Scores:
//           Total (TOTAL),
//           Number of Categories (N_CAT),
//           Accuracy (ACC),
//           Hanssen and Kuipers Discriminant (HK),
//           Heidke Skill Score (HSS),
//           Gerrity Score (GER),
//           NOTE: All statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (6) STAT and CNT Text Files, Statistics of Continuous Variables:
//           Total (TOTAL),
//           Forecast Mean (FBAR),
//           Forecast Standard Deviation (FSTDEV),
//           Observation Mean (OBAR),
//           Observation Standard Deviation (OSTDEV),
//           Pearson's Correlation Coefficient (PR_CORR),
//           Spearman's Rank Correlation Coefficient (SP_CORR),
//           Kendall Tau Rank Correlation Coefficient (KT_CORR),
//           Number of ranks compared (RANKS),
//           Number of tied ranks in the forecast field (FRANK_TIES),
//           Number of tied ranks in the observation field
(ORANK_TIES),
//           Mean Error (ME),
//           Standard Deviation of the Error (ESTDEV),
//           Multiplicative Bias (MBIAS = FBAR - OBAR),
//           Mean Absolute Error (MAE),
//           Mean Squared Error (MSE),
//           Bias-Corrected Mean Squared Error (BCMSE),
//           Root Mean Squared Error (RMSE),
//           Percentiles of the Error (E10, E25, E50, E75, E90)
//           NOTE: Most statistics listed above contain parametric
and/or
//                 non-parametric confidence interval limits.
//
//    (7) STAT and SL1L2 Text Files, Scalar Partial Sums:
//           Total (TOTAL),
//           Forecast Mean (FBAR),
//              = mean(f)
//           Observation Mean (OBAR),
//              = mean(o)
//           Forecast*Observation Product Mean (FOBAR),
//              = mean(f*o)
//           Forecast Squared Mean (FFBAR),
//              = mean(f^2)
//           Observation Squared Mean (OOBAR)
//              = mean(o^2)
//
//    (8) STAT and SAL1L2 Text Files, Scalar Anomaly Partial Sums:
//           Total (TOTAL),
//           Forecast Anomaly Mean (FABAR),
//              = mean(f-c)
//           Observation Anomaly Mean (OABAR),
//              = mean(o-c)
//           Product of Forecast and Observation Anomalies Mean
(FOABAR),
//              = mean((f-c)*(o-c))
//           Forecast Anomaly Squared Mean (FFABAR),
//              = mean((f-c)^2)
//           Observation Anomaly Squared Mean (OOABAR)
//              = mean((o-c)^2)
//
//    (9) STAT and VL1L2 Text Files, Vector Partial Sums:
//           Total (TOTAL),
//           U-Forecast Mean (UFBAR),
//              = mean(uf)
//           V-Forecast Mean (VFBAR),
//              = mean(vf)
//           U-Observation Mean (UOBAR),
//              = mean(uo)
//           V-Observation Mean (VOBAR),
//              = mean(vo)
//           U-Product Plus V-Product (UVFOBAR),
//              = mean(uf*uo+vf*vo)
//           U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
//              = mean(uf^2+vf^2)
//           U-Observation Squared Plus V-Observation Squared
(UVOOBAR)
//              = mean(uo^2+vo^2)
//
//   (10) STAT and VAL1L2 Text Files, Vector Anomaly Partial Sums:
//           U-Forecast Anomaly Mean (UFABAR),
//              = mean(uf-uc)
//           V-Forecast Anomaly Mean (VFABAR),
//              = mean(vf-vc)
//           U-Observation Anomaly Mean (UOABAR),
//              = mean(uo-uc)
//           V-Observation Anomaly Mean (VOABAR),
//              = mean(vo-vc)
//           U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
//              = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
//           U-Forecast Anomaly Squared Plus V-Forecast Anomaly
Squared (UVFFABAR),
//              = mean((uf-uc)^2+(vf-vc)^2)
//           U-Observation Anomaly Squared Plus V-Observation Anomaly
Squared (UVOOABAR)
//              = mean((uo-uc)^2+(vo-vc)^2)
//
//   (11) STAT and PCT Text Files, Nx2 Probability Contingency Table
Counts:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Row Observation Yes Count (OY_i),
//           Row Observation No Count (ON_i),
//           NOTE: Previous 3 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (12) STAT and PSTD Text Files, Nx2 Probability Contingency Table
Scores:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Base Rate (BASER) with confidence interval limits,
//           Reliability (RELIABILITY),
//           Resolution (RESOLUTION),
//           Uncertainty (UNCERTAINTY),
//           Area Under the ROC Curve (ROC_AUC),
//           Brier Score (BRIER) with confidence interval limits,
//           Probability Threshold Value (THRESH_i)
//           NOTE: Previous column repeated for each probability
threshold.
//
//   (13) STAT and PJC Text Files, Joint/Continuous Statistics of
//                                 Probabilistic Variables:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Observation Yes Count Divided by Total (OY_TP_i),
//           Observation No Count Divided by Total (ON_TP_i),
//           Calibration (CALIBRATION_i),
//           Refinement (REFINEMENT_i),
//           Likelikhood (LIKELIHOOD_i),
//           Base Rate (BASER_i),
//           NOTE: Previous 7 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (14) STAT and PRC Text Files, ROC Curve Points for
//                                 Probabilistic Variables:
//           Total (TOTAL),
//           Number of Forecast Probability Thresholds (N_THRESH),
//           Probability Threshold Value (THRESH_i),
//           Probability of Detecting Yes (PODY_i),
//           Probability of False Detection (POFD_i),
//           NOTE: Previous 3 columns repeated for each row in the
table.
//           Last Probability Threshold Value (THRESH_n)
//
//   (15) STAT and MPR Text Files, Matched Pair Data:
//           Total (TOTAL),
//           Index (INDEX),
//           Observation Station ID (OBS_SID),
//           Observation Latitude (OBS_LAT),
//           Observation Longitude (OBS_LON),
//           Observation Level (OBS_LVL),
//           Observation Elevation (OBS_ELV),
//           Forecast Value (FCST),
//           Observation Value (OBS),
//           Climatological Value (CLIMO)
//
//   In the expressions above, f are forecast values, o are observed
values,
//   and c are climatological values.
//
// Values for these flags are interpreted as follows:
//    (0) Do not generate output of this type
//    (1) Write output to a STAT file
//    (2) Write output to a STAT file and a text file
//
output_flag[] = [ 2, 2, 2, 0, 0, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];

//
// Flag to indicate whether Kendall's Tau and Spearman's Rank
Correlation
// Coefficients should be computed.  Computing them over large
datasets is
// computationally intensive and slows down the runtime execution
significantly.
//    (0) Do not compute these correlation coefficients
//    (1) Compute these correlation coefficients
//
rank_corr_flag = 1;

//
// Specify the GRIB Table 2 parameter table version number to be used
// for interpreting GRIB codes.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
grib_ptv = 2;

//
// Directory where temporary files should be written.
//
tmp_dir = "/tmp";

//
// Prefix to be used for the output file names.
//
output_prefix = "24hr";

//
// Indicate a version number for the contents of this configuration
file.
// The value should generally not be modified.
//
version = "V3.0";

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

#!/bin/sh
cd ~/wrf3/utils/met/METv3.0/myascii
../bin/ascii2nc precip24_2011061000.ascii
../myascii2nc/precip24_2011061000.nc
cd ~/wrf3/WRFV3/wrfprd
ls -lt wrfout_d01_*00*.nc | tail -1 >out1
cp out1 ~/wrf3/WRFV3/postprd
export yr=`cat out1 | cut -c70-73`
export mn=`cat out1 | cut -c75-76`
export dy=`cat out1 | cut -c78-79`
export hr=`cat out1 | cut -c81-82`
cd ~/wrf3/WRFV3/postprd
./run_wrfpost
## checking the no of WRFPRS files
export chk1=`ls -l WRFPRS* | wc -l`
cd ~/wrf3/utils/met/METv3.0/
./bin/pcp_combine -subtract ~/wrf3/WRFV3/postprd/WRFPRS_d01.027 27
~/wrf3/WRFV3/postprd/WRFPRS_d01.003 3
./myoutput/pcp_combine/24hr/model24_$yr$mn$dy$hr.nc
./bin/pcp_combine -subtract ~/wrf3/WRFV3/postprd/WRFPRS_d01.051 51
~/wrf3/WRFV3/postprd/WRFPRS_d01.027 27
./myoutput/pcp_combine/48hr/model48_$yr$mn$dy$hr.nc
# ############################
# step 3 using point-stat tool
# ##################################################################
./bin/point_stat
~/wrf3/utils/met/METv3.0/myoutput/pcp_combine/24hr/model24_$yr$mn$dy$hr.nc
~/wrf3/utils/met/METv3.0/myascii2nc/precip24_2011061000.nc
~/wrf3/utils/met/METv3.0/myscripts/PointStatConfig_rain
~/wrf3/utils/met/METv3.0/myoutput/point_stat/24hr/
./bin/point_stat ./myoutput/pcp_combine/48hr/model48_$yr$mn$dy$hr.nc
./myascii2nc/precip24_2011061100.nc ./myscripts/PointStatConfig_rain
./myoutput/point_stat/48hr/
rm ~/wrf3/WRFV3/wrfprd/wrfout_d*.nc
rm ~/wrf3/WRFV3/postprd/WRFPRS_d*
rm ~/wrf3/WRFV3/postprd/wrfpost*
exit

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

Dear met_help,
I have a small problem which goes like this. I am unablle to make
point_stat write the output in the folder that I want . It writes in
the default folder, METv3.0/out/point-stat.
There is no option in the configuration file specifying so,.
I am attaching the files.
Thanks.
geeta

> To: met_help_archive at mailman.ucar.edu
> From: met_help at ucar.edu
> Date: Mon, 16 Apr 2012 08:52:43 -0600
> CC: met_help at ucar.edu
> Subject: [Met_help] [rt.rap.ucar.edu #55790] History for Geeta
met_help problem
>
> ----------------------------------------------------------------
> Initial Request
> ----------------------------------------------------------------
>
> Dear administrator,
> I wish to seek help from the met_help, and I am sending my message
for the last 3-4 days. But everytime I am told that there is some
error that has occurred.
> I donot know whtheer my meage is reaching across or not? I have
already registered but my query is not answered.
>
> kindly help me
> Regards
> geeta
>
> Dr G. Agnihotri,
> Met Centre, IMD, Palace Road,
> Bangalore, India.
> 080-22244419
>
>
>
> ----------------------------------------------------------------
> Complete Ticket History
> ----------------------------------------------------------------
>
> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
problem
> From: Paul Oldenburg
> Time: Wed Apr 11 09:24:51 2012
>
> Geeta,
>
> I created a met_help ticket for you. Please respond to this email
and
> describe the MET problem you are having in detail.
>
> Thanks,
>
> Paul
>
>
> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> > Greetings,
> >
> > This message has been automatically generated in response to the
> > creation of a trouble ticket regarding:
> > "Geeta met_help problem",
> > a summary of which appears below.
> >
> > There is no need to reply to this message right now. Your ticket
> has been
> > assigned an ID of [rt.rap.ucar.edu #55790].
> >
> > Please include the string:
> >
> > [rt.rap.ucar.edu #55790]
> >
> > in the subject line of all future correspondence about this issue.
> To do so,
> > you may reply to this message.
> >
> > Thank you,
> > met_help at ucar.edu
> >
> >
>
-------------------------------------------------------------------------
> > Dear administrator,
> > I wish to seek help from the met_help, and I am sending my message
> for the last 3-4 days. But everytime I am told that there is some
> error that has occurred.
> > I donot know whtheer my meage is reaching across or not? I have
> already registered but my query is not answered.
> >
> > kindly help me
> > Regards
> > geeta
> >
> > Dr G. Agnihotri,
> > Met Centre, IMD, Palace Road,
> > Bangalore, India.
> > 080-22244419
> >
>
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Paul Oldenburg
> Time: Thu Apr 12 09:03:20 2012
>
> Geeta,
>
> First of all, if you want to send data files to me that are larger
> than 2MB, please upload them to our FTP site using
> the instructions found here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
>
> I found several problems, but they were all easy to find and fix.
> There was a config file syntax error in the version
> of PointStatConfig_rain that you sent me. You must use double quotes
> around values in the obs_field setting:
>
> obs_field[] = [ "61/A24" ];
>
> Next, I recommend using a verbosity setting of 3 or higher when you
> start using point_stat. This will cause point_stat
> to print out the reason why points are rejected. One reason is that
> there is a mis-match between the valid time of your
> model data and the valid time of your observation data. The
following
> ncdump commands show the discrepancy:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
> ...
> APCP_24:valid_time = "20110610_030000" ;
>
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
> precip24_2011061000.nc
> ...
> hdr_vld =
> "20110610_083000",
> ...
>
> You can tell point_stat to allow this discrepancy using the
following
> two settings in PointStatConfig_rain:
>
> beg_ds = -21600;
> end_ds = 21600;
>
> Next, the message type of your observations is ADPSFC, which you can
> see using the following ncdump command:
>
> [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
> precip24_2011061000.nc
> ...
> hdr_typ =
> "ADPSFC",
> ...
>
> In PointStatConfig_rain, the following setting must be used to
specify
> this message type:
>
> message_type[] = [ "ADPSFC" ];
>
> Once these changes were made, I ran point_stat with your model and
> observation data and the result was 26 matched pairs.
> I attached the output below and the PointStatConfig_rain file that I
> used. Please let me know if you have any more
> questions.
>
> One last thing, I think you might have a date typo in your
> precip24_2011061000.ascii file on line 2. Should the date
> string "20110613_083000" actually be "20110610_083000" like all the
> others?
>
> Thanks,
>
> Paul
>
>
> $ point_stat model24_day.nc precip24_2011061000.nc
> PointStatConfig_rain -v 3
> DEBUG 3: It's a Mercator grid...
> DEBUG 3: Mercator Data:
> DEBUG 3: lat_ll = 9.097
> DEBUG 3: lon_ll = -73.943
> DEBUG 3: lat_ur = 21.711
> DEBUG 3: lon_ur = -87.057
> DEBUG 3: ny = 163
> DEBUG 3: nx = 163
> DEBUG 1: Default Config File:
>
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
> DEBUG 1: User Config File: PointStatConfig_rain
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072433272399
> DEBUG 1: Forecast File: model24_day.nc
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File: precip24_2011061000.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for APCP_24(*,*).
> DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0 climatology
> levels.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 27 observations from 27 messages.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation type
> ADPSFC, over region FULL, for interpolation method
> UW_MEAN(1), using 26 pairs.
> DEBUG 3: Number of matched pairs = 26
> DEBUG 3: Observations processed = 27
> DEBUG 3: Rejected: GRIB code = 0
> DEBUG 3: Rejected: valid time = 1
> DEBUG 3: Rejected: bad obs value = 0
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 0
> DEBUG 3: Rejected: message type = 0
> DEBUG 3: Rejected: masking region = 0
> DEBUG 3: Rejected: bad fcst value = 0
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Computing Scalar Partial Sums.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_fho.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_ctc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cts.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_cnt.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sl1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sal1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_vl1l2.txt
> DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_val1l2.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pct.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pstd.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_pjc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_prc.txt
> DEBUG 1: Output file:
out/point_stat_270000L_20110610_030000V_mpr.txt
>
>
>
> On 04/11/2012 11:53 PM, Geeta Geeta wrote:
> >
> > Dear met_help,
> > I am using point stat tool for Rainfall verification in my region
> Western part of Indian peninsula (9-21N/74-87E).
> >
> >
> > 1. for this, I have converted the ascii observations to netcdf
> using the ascii2nc tool (O/p file precip24_2011061000.nc)
> > 2. I have converted the model o/p in .nc format to grib1 format
> using WPP.
> > 3. I have then used pcpcombine to get the desired RF/output ie. by
> subtracting the RF of 27th hour from the 4th hour. (model24_day.nc
is
> attached)
> >
> >
> > 4. now i am running
> > -bash-3.2$ cd ..
> > -bash-3.2$ cd verf-met/
> > /oprn/model/wrf3/utils/met/METv3.0/verf-met
> > -bash-3.2$ ../bin/point_stat model24_day.nc precip24_2011061000.nc
> ~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=801800783
> > Forecast File: model24_day.nc
> > Climatology File: none
> > Configuration File:
> /oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > Observation File: precip24_2011061000.nc
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Reading records for APCP_24(*,*).
> > For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Searching 27 observations from 27 PrepBufr messages.
> >
> >
>
--------------------------------------------------------------------------------
> >
> > Processing APCP_24(*,*) versus APCP_24/A24, for observation type
> ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
0
> pairs.
> >
>
--------------------------------------------------------------------------------
> >
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
> > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
> > -bash-3.2$
> >
> >
> >
> > My problem is that there is no data in the output files, when i
> open them, it is only headers.
> > Kindly help.
> > I have plotted the grib1 files (after WPP by using grib2ctl.pl) to
> see the domain and if all the stations are falling within this
domain,
> So all looks okay to me
> > geeta
> >
> >
> >
> >> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem
> >> From: met_help at ucar.edu
> >> To: geeta124 at hotmail.com
> >> Date: Wed, 11 Apr 2012 09:24:51 -0600
> >>
> >> Geeta,
> >>
> >> I created a met_help ticket for you. Please respond to this email
> and describe the MET problem you are having in detail.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> >>> Greetings,
> >>>
> >>> This message has been automatically generated in response to the
> >>> creation of a trouble ticket regarding:
> >>> "Geeta met_help problem",
> >>> a summary of which appears below.
> >>>
> >>> There is no need to reply to this message right now. Your ticket
> has been
> >>> assigned an ID of [rt.rap.ucar.edu #55790].
> >>>
> >>> Please include the string:
> >>>
> >>> [rt.rap.ucar.edu #55790]
> >>>
> >>> in the subject line of all future correspondence about this
issue.
> To do so,
> >>> you may reply to this message.
> >>>
> >>> Thank you,
> >>> met_help at ucar.edu
> >>>
> >>>
>
-------------------------------------------------------------------------
> >>> Dear administrator,
> >>> I wish to seek help from the met_help, and I am sending my
message
> for the last 3-4 days. But everytime I am told that there is some
> error that has occurred.
> >>> I donot know whtheer my meage is reaching across or not? I have
> already registered but my query is not answered.
> >>>
> >>> kindly help me
> >>> Regards
> >>> geeta
> >>>
> >>> Dr G. Agnihotri,
> >>> Met Centre, IMD, Palace Road,
> >>> Bangalore, India.
> >>> 080-22244419
> >>>
> >>
> >>
> >
>
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Paul Oldenburg
> Time: Thu Apr 12 09:03:20 2012
>
>
////////////////////////////////////////////////////////////////////////////////
> //
> // Default point_stat configuration file
> //
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Specify a name to designate the model being verified. This name
> will be
> // written to the second column of the ASCII output generated.
> //
> model = "WRF";
>
> //
> // Beginning and ending time offset values in seconds for
observations
> // to be used. These time offsets are defined in reference to the
> // forecast valid time, v. Observations with a valid time falling in
> the
> // window [v+beg_ds, v+end_ds] will be used.
> // These selections are overridden by the command line arguments
> // -obs_valid_beg and -obs_valid_end.
> //
> beg_ds = -21600;
> end_ds = 21600;
>
> //
> // Specify a comma-separated list of fields to be verified. The
> forecast and
> // observation fields may be specified separately. If the obs_field
> parameter
> // is left blank, it will default to the contents of fcst_field.
> //
> // Each field is specified as a GRIB code or abbreviation followed
by
> an
> // accumulation or vertical level indicator for GRIB files or as a
> variable name
> // followed by a list of dimensions for NetCDF files output from
> p_interp or MET.
> //
> // Specifying verification fields for GRIB files:
> // GC/ANNN for accumulation interval NNN
> // GC/ZNNN for vertical level NNN
> // GC/ZNNN-NNN for a range of vertical levels (MSL or AGL)
> // GC/PNNN for pressure level NNN in hPa
> // GC/PNNN-NNN for a range of pressure levels in hPa
> // GC/LNNN for a generic level type
> // GC/RNNN for a specific GRIB record number
> // Where GC is the number of or abbreviation for the grib code
> // to be verified.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> //
> // Specifying verification fields for NetCDF files:
> // var_name(i,...,j,*,*) for a single field
> // var_name(i-j,*,*) for a range of fields
> // Where var_name is the name of the NetCDF variable,
> // and i,...,j specifies fixed dimension values,
> // and i-j specifies a range of values for a single dimension,
> // and *,* specifies the two dimensions for the gridded field.
> //
> // NOTE: To verify winds as vectors rather than scalars,
> // specify UGRD (or 33) followed by VGRD (or 34) with the
> // same level values.
> //
> // NOTE: To process a probability field, add "/PROB", such as
> "POP/Z0/PROB".
> //
> // e.g. fcst_field[] = [ "SPFH/P500", "TMP/P500" ]; for a GRIB input
> // e.g. fcst_field[] = [ "QVAPOR(0,5,*,*)", "TT(0,5,*,*)" ]; for
> NetCDF input
> //
> fcst_field[] = [ "APCP_24(*,*)"];
> obs_field[] = [ "61/A24" ];
>
> //
> // Specify a comma-separated list of groups of thresholds to be
> applied to the
> // fields listed above. Thresholds for the forecast and observation
> fields
> // may be specified separately. If the obs_thresh parameter is left
> blank,
> // it will default to the contents of fcst_thresh.
> //
> // At least one threshold must be provided for each field listed
> above. The
> // lengths of the "fcst_field" and "fcst_thresh" arrays must match,
as
> must
> // lengths of the "obs_field" and "obs_thresh" arrays. To apply
> multiple
> // thresholds to a field, separate the threshold values with a
space.
> //
> // Each threshold must be preceded by a two letter indicator for the
> type of
> // thresholding to be performed:
> // 'lt' for less than 'le' for less than or equal to
> // 'eq' for equal to 'ne' for not equal to
> // 'gt' for greater than 'ge' for greater than or equal to
> //
> // NOTE: Thresholds for probabilities must begin with 0.0, end with
> 1.0,
> // and be preceeded by "ge".
> //
> // e.g. fcst_thresh[] = [ "gt80", "gt273" ];
> //
> fcst_thresh[] = [ "gt20"];
> obs_thresh[] = [ "gt20"];
>
> //
> // Specify a comma-separated list of thresholds to be used when
> computing
> // VL1L2 and VAL1L2 partial sums for winds. The thresholds are
> applied to the
> // wind speed values derived from each U/V pair. Only those U/V
pairs
> which meet
> // the wind speed threshold criteria are retained. If the
> obs_wind_thresh
> // parameter is left blank, it will default to the contents of
> fcst_wind_thresh.
> //
> // To apply multiple wind speed thresholds, separate the threshold
> values with a
> // space. Use "NA" to indicate that no wind speed threshold should
be
> applied.
> //
> // Each threshold must be preceded by a two letter indicator for the
> type of
> // thresholding to be performed:
> // 'lt' for less than 'le' for less than or equal to
> // 'eq' for equal to 'ne' for not equal to
> // 'gt' for greater than 'ge' for greater than or equal to
> // 'NA' for no threshold
> //
> // e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
> //
> fcst_wind_thresh[] = [ "NA" ];
> obs_wind_thresh[] = [];
>
> //
> // Specify a comma-separated list of PrepBufr message types with
which
> // to perform the verification. Statistics will be computed
> separately
> // for each message type specified. At least one PrepBufr message
> type
> // must be provided.
> // List of valid message types:
> // ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
> // MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
> // SFCSHP SPSSMI SYNDAT VADWND
> // ANYAIR (= AIRCAR, AIRCFT)
> // ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
> // ONLYSF (= ADPSFC, SFCSHP)
> //
>
http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
> //
> // e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
> //
> message_type[] = [ "ADPSFC" ];
>
> //
> // Specify a comma-separated list of grids to be used in masking the
> data over
> // which to perform scoring. An empty list indicates that no masking
> grid
> // should be performed. The standard NCEP grids are named "GNNN"
> where NNN
> // indicates the three digit grid number. Enter "FULL" to score over
> the
> // entire domain.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
> //
> // e.g. mask_grid[] = [ "FULL" ];
> //
> mask_grid[] = [ "FULL" ];
>
> //
> // Specify a comma-separated list of masking regions to be applied.
> // An empty list indicates that no additional masks should be used.
> // The masking regions may be defined in one of 4 ways:
> //
> // (1) An ASCII file containing a lat/lon polygon.
> // Latitude in degrees north and longitude in degrees east.
> // By default, the first and last polygon points are connected.
> // e.g. "MET_BASE/data/poly/EAST.poly" which consists of n points:
> // "poly_name lat1 lon1 lat2 lon2... latn lonn"
> //
> // (2) The NetCDF output of the gen_poly_mask tool.
> //
> // (3) A NetCDF data file, followed by the name of the NetCDF
variable
> // to be used, and optionally, a threshold to be applied to the
> field.
> // e.g. "sample.nc var_name gt0.00"
> //
> // (4) A GRIB data file, followed by a description of the field
> // to be used, and optionally, a threshold to be applied to the
> field.
> // e.g. "sample.grb APCP/A3 gt0.00"
> //
> // Any NetCDF or GRIB file used must have the same grid dimensions
as
> the
> // data being verified.
> //
> // MET_BASE may be used in the path for the files above.
> //
> // e.g. mask_poly[] = [ "MET_BASE/data/poly/EAST.poly",
> // "poly_mask.ncf",
> // "sample.nc APCP",
> // "sample.grb HGT/Z0 gt100.0" ];
> //
> mask_poly[] = [];
>
> //
> // Specify the name of an ASCII file containing a space-separated
list
> of
> // station ID's at which to perform verification. Each station ID
> specified
> // is treated as an individual masking region.
> //
> // An empty list file name indicates that no station ID masks should
> be used.
> //
> // MET_BASE may be used in the path for the station ID mask file
name.
> //
> // e.g. mask_sid = "CONUS.stations";
> //
> mask_sid = "";
>
> //
> // Specify a comma-separated list of values for alpha to be used
when
> computing
> // confidence intervals. Values of alpha must be between 0 and 1.
> //
> // e.g. ci_alpha[] = [ 0.05, 0.10 ];
> //
> ci_alpha[] = [ 0.05 ];
>
> //
> // Specify the method to be used for computing bootstrap confidence
> intervals.
> // The value for this is interpreted as follows:
> // (0) Use the BCa interval method (computationally intensive)
> // (1) Use the percentile interval method
> //
> boot_interval = 1;
>
> //
> // Specify a proportion between 0 and 1 to define the replicate
sample
> size
> // to be used when computing percentile intervals. The replicate
> sample
> // size is set to boot_rep_prop * n, where n is the number of raw
data
> points.
> //
> // e.g boot_rep_prop = 0.80;
> //
> boot_rep_prop = 1.0;
>
> //
> // Specify the number of times each set of matched pair data should
be
> // resampled when computing bootstrap confidence intervals. A value
> of
> // zero disables the computation of bootstrap condifence intervals.
> //
> // e.g. n_boot_rep = 1000;
> //
> n_boot_rep = 1000;
>
> //
> // Specify the name of the random number generator to be used. See
> the MET
> // Users Guide for a list of possible random number generators.
> //
> boot_rng = "mt19937";
>
> //
> // Specify the seed value to be used when computing bootstrap
> confidence
> // intervals. If left unspecified, the seed will change for each run
> and
> // the computed bootstrap confidence intervals will not be
> reproducable.
> //
> boot_seed = "";
>
> //
> // Specify a comma-separated list of interpolation method(s) to be
> used
> // for comparing the forecast grid to the observation points. String
> values
> // are interpreted as follows:
> // MIN = Minimum in the neighborhood
> // MAX = Maximum in the neighborhood
> // MEDIAN = Median in the neighborhood
> // UW_MEAN = Unweighted mean in the neighborhood
> // DW_MEAN = Distance-weighted mean in the neighborhood
> // LS_FIT = Least-squares fit in the neighborhood
> //
> // In all cases, vertical interpolation is performed in the natural
> log
> // of pressure of the levels above and below the observation.
> //
> // e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
> //
> interp_method[] = [ "UW_MEAN" ];
>
> //
> // Specify a comma-separated list of box widths to be used by the
> // interpolation techniques listed above. A value of 1 indicates
that
> // the nearest neighbor approach should be used. For a value of n
> // greater than 1, the n*n grid points closest to the observation
> define
> // the neighborhood.
> //
> // e.g. interp_width = [ 1, 3, 5 ];
> //
> interp_width[] = [ 1 ];
>
> //
> // When interpolating, compute a ratio of the number of valid data
> points
> // to the total number of points in the neighborhood. If that ratio
> is
> // less than this threshold, do not include the observation. This
> // threshold must be between 0 and 1. Setting this threshold to 1
> will
> // require that each observation be surrounded by n*n valid forecast
> // points.
> //
> // e.g. interp_thresh = 1.0;
> //
> interp_thresh = 1.0;
>
> //
> // Specify flags to indicate the type of data to be output:
> // (1) STAT and FHO Text Files, Forecast, Hit, Observation Rates:
> // Total (TOTAL),
> // Forecast Rate (F_RATE),
> // Hit Rate (H_RATE),
> // Observation Rate (O_RATE)
> //
> // (2) STAT and CTC Text Files, Contingency Table Counts:
> // Total (TOTAL),
> // Forecast Yes and Observation Yes Count (FY_OY),
> // Forecast Yes and Observation No Count (FY_ON),
> // Forecast No and Observation Yes Count (FN_OY),
> // Forecast No and Observation No Count (FN_ON)
> //
> // (3) STAT and CTS Text Files, Contingency Table Scores:
> // Total (TOTAL),
> // Base Rate (BASER),
> // Forecast Mean (FMEAN),
> // Accuracy (ACC),
> // Frequency Bias (FBIAS),
> // Probability of Detecting Yes (PODY),
> // Probability of Detecting No (PODN),
> // Probability of False Detection (POFD),
> // False Alarm Ratio (FAR),
> // Critical Success Index (CSI),
> // Gilbert Skill Score (GSS),
> // Hanssen and Kuipers Discriminant (HK),
> // Heidke Skill Score (HSS),
> // Odds Ratio (ODDS),
> // NOTE: All statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (4) STAT and MCTC Text Files, NxN Multi-Category Contingency
> Table Counts:
> // Total (TOTAL),
> // Number of Categories (N_CAT),
> // Contingency Table Count columns repeated N_CAT*N_CAT
> times
> //
> // (5) STAT and MCTS Text Files, NxN Multi-Category Contingency
> Table Scores:
> // Total (TOTAL),
> // Number of Categories (N_CAT),
> // Accuracy (ACC),
> // Hanssen and Kuipers Discriminant (HK),
> // Heidke Skill Score (HSS),
> // Gerrity Score (GER),
> // NOTE: All statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (6) STAT and CNT Text Files, Statistics of Continuous Variables:
> // Total (TOTAL),
> // Forecast Mean (FBAR),
> // Forecast Standard Deviation (FSTDEV),
> // Observation Mean (OBAR),
> // Observation Standard Deviation (OSTDEV),
> // Pearson's Correlation Coefficient (PR_CORR),
> // Spearman's Rank Correlation Coefficient (SP_CORR),
> // Kendall Tau Rank Correlation Coefficient (KT_CORR),
> // Number of ranks compared (RANKS),
> // Number of tied ranks in the forecast field (FRANK_TIES),
> // Number of tied ranks in the observation field
> (ORANK_TIES),
> // Mean Error (ME),
> // Standard Deviation of the Error (ESTDEV),
> // Multiplicative Bias (MBIAS = FBAR - OBAR),
> // Mean Absolute Error (MAE),
> // Mean Squared Error (MSE),
> // Bias-Corrected Mean Squared Error (BCMSE),
> // Root Mean Squared Error (RMSE),
> // Percentiles of the Error (E10, E25, E50, E75, E90)
> // NOTE: Most statistics listed above contain parametric
> and/or
> // non-parametric confidence interval limits.
> //
> // (7) STAT and SL1L2 Text Files, Scalar Partial Sums:
> // Total (TOTAL),
> // Forecast Mean (FBAR),
> // = mean(f)
> // Observation Mean (OBAR),
> // = mean(o)
> // Forecast*Observation Product Mean (FOBAR),
> // = mean(f*o)
> // Forecast Squared Mean (FFBAR),
> // = mean(f^2)
> // Observation Squared Mean (OOBAR)
> // = mean(o^2)
> //
> // (8) STAT and SAL1L2 Text Files, Scalar Anomaly Partial Sums:
> // Total (TOTAL),
> // Forecast Anomaly Mean (FABAR),
> // = mean(f-c)
> // Observation Anomaly Mean (OABAR),
> // = mean(o-c)
> // Product of Forecast and Observation Anomalies Mean
> (FOABAR),
> // = mean((f-c)*(o-c))
> // Forecast Anomaly Squared Mean (FFABAR),
> // = mean((f-c)^2)
> // Observation Anomaly Squared Mean (OOABAR)
> // = mean((o-c)^2)
> //
> // (9) STAT and VL1L2 Text Files, Vector Partial Sums:
> // Total (TOTAL),
> // U-Forecast Mean (UFBAR),
> // = mean(uf)
> // V-Forecast Mean (VFBAR),
> // = mean(vf)
> // U-Observation Mean (UOBAR),
> // = mean(uo)
> // V-Observation Mean (VOBAR),
> // = mean(vo)
> // U-Product Plus V-Product (UVFOBAR),
> // = mean(uf*uo+vf*vo)
> // U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
> // = mean(uf^2+vf^2)
> // U-Observation Squared Plus V-Observation Squared
> (UVOOBAR)
> // = mean(uo^2+vo^2)
> //
> // (10) STAT and VAL1L2 Text Files, Vector Anomaly Partial Sums:
> // U-Forecast Anomaly Mean (UFABAR),
> // = mean(uf-uc)
> // V-Forecast Anomaly Mean (VFABAR),
> // = mean(vf-vc)
> // U-Observation Anomaly Mean (UOABAR),
> // = mean(uo-uc)
> // V-Observation Anomaly Mean (VOABAR),
> // = mean(vo-vc)
> // U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
> // = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
> // U-Forecast Anomaly Squared Plus V-Forecast Anomaly
> Squared (UVFFABAR),
> // = mean((uf-uc)^2+(vf-vc)^2)
> // U-Observation Anomaly Squared Plus V-Observation Anomaly
> Squared (UVOOABAR)
> // = mean((uo-uc)^2+(vo-vc)^2)
> //
> // (11) STAT and PCT Text Files, Nx2 Probability Contingency Table
> Counts:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Row Observation Yes Count (OY_i),
> // Row Observation No Count (ON_i),
> // NOTE: Previous 3 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (12) STAT and PSTD Text Files, Nx2 Probability Contingency Table
> Scores:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Base Rate (BASER) with confidence interval limits,
> // Reliability (RELIABILITY),
> // Resolution (RESOLUTION),
> // Uncertainty (UNCERTAINTY),
> // Area Under the ROC Curve (ROC_AUC),
> // Brier Score (BRIER) with confidence interval limits,
> // Probability Threshold Value (THRESH_i)
> // NOTE: Previous column repeated for each probability
> threshold.
> //
> // (13) STAT and PJC Text Files, Joint/Continuous Statistics of
> // Probabilistic Variables:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Observation Yes Count Divided by Total (OY_TP_i),
> // Observation No Count Divided by Total (ON_TP_i),
> // Calibration (CALIBRATION_i),
> // Refinement (REFINEMENT_i),
> // Likelikhood (LIKELIHOOD_i),
> // Base Rate (BASER_i),
> // NOTE: Previous 7 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (14) STAT and PRC Text Files, ROC Curve Points for
> // Probabilistic Variables:
> // Total (TOTAL),
> // Number of Forecast Probability Thresholds (N_THRESH),
> // Probability Threshold Value (THRESH_i),
> // Probability of Detecting Yes (PODY_i),
> // Probability of False Detection (POFD_i),
> // NOTE: Previous 3 columns repeated for each row in the
> table.
> // Last Probability Threshold Value (THRESH_n)
> //
> // (15) STAT and MPR Text Files, Matched Pair Data:
> // Total (TOTAL),
> // Index (INDEX),
> // Observation Station ID (OBS_SID),
> // Observation Latitude (OBS_LAT),
> // Observation Longitude (OBS_LON),
> // Observation Level (OBS_LVL),
> // Observation Elevation (OBS_ELV),
> // Forecast Value (FCST),
> // Observation Value (OBS),
> // Climatological Value (CLIMO)
> //
> // In the expressions above, f are forecast values, o are observed
> values,
> // and c are climatological values.
> //
> // Values for these flags are interpreted as follows:
> // (0) Do not generate output of this type
> // (1) Write output to a STAT file
> // (2) Write output to a STAT file and a text file
> //
> output_flag[] = [ 2, 2, 2, 0, 0, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 ];
>
> //
> // Flag to indicate whether Kendall's Tau and Spearman's Rank
> Correlation
> // Coefficients should be computed. Computing them over large
> datasets is
> // computationally intensive and slows down the runtime execution
> significantly.
> // (0) Do not compute these correlation coefficients
> // (1) Compute these correlation coefficients
> //
> rank_corr_flag = 1;
>
> //
> // Specify the GRIB Table 2 parameter table version number to be
used
> // for interpreting GRIB codes.
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> //
> grib_ptv = 2;
>
> //
> // Directory where temporary files should be written.
> //
> tmp_dir = "/tmp";
>
> //
> // Prefix to be used for the output file names.
> //
> output_prefix = "";
>
> //
> // Indicate a version number for the contents of this configuration
> file.
> // The value should generally not be modified.
> //
> version = "V3.0";
>
> ------------------------------------------------
> Subject: Geeta met_help problem
> From: Geeta Geeta
> Time: Mon Apr 16 03:08:18 2012
>
>
> Dear Met_help,
> I have done the corrections in my configuration file and it works
fine
> now.
> Thanks for Ur support.
>
> geeta
>
>
> > Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem, error with point-stat
> > From: met_help at ucar.edu
> > To: geeta124 at hotmail.com
> > Date: Thu, 12 Apr 2012 09:03:20 -0600
> >
> > Geeta,
> >
> > First of all, if you want to send data files to me that are larger
> than 2MB, please upload them to our FTP site using
> > the instructions found here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp.
> >
> > I found several problems, but they were all easy to find and fix.
> There was a config file syntax error in the version
> > of PointStatConfig_rain that you sent me. You must use double
quotes
> around values in the obs_field setting:
> >
> > obs_field[] = [ "61/A24" ];
> >
> > Next, I recommend using a verbosity setting of 3 or higher when
you
> start using point_stat. This will cause point_stat
> > to print out the reason why points are rejected. One reason is
that
> there is a mis-match between the valid time of your
> > model data and the valid time of your observation data. The
> following ncdump commands show the discrepancy:
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -h model24_day.nc
> > ...
> > APCP_24:valid_time = "20110610_030000" ;
> >
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_vld
> precip24_2011061000.nc
> > ...
> > hdr_vld =
> > "20110610_083000",
> > ...
> >
> > You can tell point_stat to allow this discrepancy using the
> following two settings in PointStatConfig_rain:
> >
> > beg_ds = -21600;
> > end_ds = 21600;
> >
> > Next, the message type of your observations is ADPSFC, which you
can
> see using the following ncdump command:
> >
> > [pgoldenb at orval 20120412.geeta]$ ncdump -v hdr_typ
> precip24_2011061000.nc
> > ...
> > hdr_typ =
> > "ADPSFC",
> > ...
> >
> > In PointStatConfig_rain, the following setting must be used to
> specify this message type:
> >
> > message_type[] = [ "ADPSFC" ];
> >
> > Once these changes were made, I ran point_stat with your model and
> observation data and the result was 26 matched pairs.
> > I attached the output below and the PointStatConfig_rain file that
I
> used. Please let me know if you have any more
> > questions.
> >
> > One last thing, I think you might have a date typo in your
> precip24_2011061000.ascii file on line 2. Should the date
> > string "20110613_083000" actually be "20110610_083000" like all
the
> others?
> >
> > Thanks,
> >
> > Paul
> >
> >
> > $ point_stat model24_day.nc precip24_2011061000.nc
> PointStatConfig_rain -v 3
> > DEBUG 3: It's a Mercator grid...
> > DEBUG 3: Mercator Data:
> > DEBUG 3: lat_ll = 9.097
> > DEBUG 3: lon_ll = -73.943
> > DEBUG 3: lat_ur = 21.711
> > DEBUG 3: lon_ur = -87.057
> > DEBUG 3: ny = 163
> > DEBUG 3: nx = 163
> > DEBUG 1: Default Config File:
>
/d1/pgoldenb/opt/MET_builds/METv3.1_gnu4/data/config/PointStatConfig_default
> > DEBUG 1: User Config File: PointStatConfig_rain
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744072433272399
> > DEBUG 1: Forecast File: model24_day.nc
> > DEBUG 1: Climatology File: none
> > DEBUG 1: Observation File: precip24_2011061000.nc
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for APCP_24(*,*).
> > DEBUG 2: For APCP_24(*,*) found 1 forecast levels and 0
climatology
> levels.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Searching 27 observations from 27 messages.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Processing APCP_24(*,*) versus 61/A24, for observation
type
> ADPSFC, over region FULL, for interpolation method
> > UW_MEAN(1), using 26 pairs.
> > DEBUG 3: Number of matched pairs = 26
> > DEBUG 3: Observations processed = 27
> > DEBUG 3: Rejected: GRIB code = 0
> > DEBUG 3: Rejected: valid time = 1
> > DEBUG 3: Rejected: bad obs value = 0
> > DEBUG 3: Rejected: off the grid = 0
> > DEBUG 3: Rejected: level mismatch = 0
> > DEBUG 3: Rejected: message type = 0
> > DEBUG 3: Rejected: masking region = 0
> > DEBUG 3: Rejected: bad fcst value = 0
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Computing Scalar Partial Sums.
> > DEBUG 2:
> > DEBUG 2:
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 1: Output file: out/point_stat_270000L_20110610_030000V.stat
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_fho.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_ctc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_cts.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_cnt.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sl1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_sal1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_vl1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_val1l2.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pct.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pstd.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_pjc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_prc.txt
> > DEBUG 1: Output file:
> out/point_stat_270000L_20110610_030000V_mpr.txt
> >
> >
> >
> > On 04/11/2012 11:53 PM, Geeta Geeta wrote:
> > >
> > > Dear met_help,
> > > I am using point stat tool for Rainfall verification in my
region
> Western part of Indian peninsula (9-21N/74-87E).
> > >
> > >
> > > 1. for this, I have converted the ascii observations to netcdf
> using the ascii2nc tool (O/p file precip24_2011061000.nc)
> > > 2. I have converted the model o/p in .nc format to grib1 format
> using WPP.
> > > 3. I have then used pcpcombine to get the desired RF/output ie.
by
> subtracting the RF of 27th hour from the 4th hour. (model24_day.nc
is
> attached)
> > >
> > >
> > > 4. now i am running
> > > -bash-3.2$ cd ..
> > > -bash-3.2$ cd verf-met/
> > > /oprn/model/wrf3/utils/met/METv3.0/verf-met
> > > -bash-3.2$ ../bin/point_stat model24_day.nc
precip24_2011061000.nc
> ~/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=801800783
> > > Forecast File: model24_day.nc
> > > Climatology File: none
> > > Configuration File:
> /oprn/model/wrf3/utils/met/METv3.0/data/config/PointStatConfig_rain
> > > Observation File: precip24_2011061000.nc
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Reading records for APCP_24(*,*).
> > > For APCP_24(*,*) found 1 forecast levels and 0 climatology
levels.
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Searching 27 observations from 27 PrepBufr messages.
> > >
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Processing APCP_24(*,*) versus APCP_24/A24, for observation type
> ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
0
> pairs.
> > >
>
--------------------------------------------------------------------------------
> > >
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V.stat
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_fho.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_ctc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cts.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_cnt.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sl1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_sal1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_vl1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_val1l2.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pct.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pstd.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_pjc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_prc.txt
> > > Output file:
>
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_270000L_20110610_030000V_mpr.txt
> > > -bash-3.2$
> > >
> > >
> > >
> > > My problem is that there is no data in the output files, when i
> open them, it is only headers.
> > > Kindly help.
> > > I have plotted the grib1 files (after WPP by using grib2ctl.pl)
to
> see the domain and if all the stations are falling within this
domain,
> So all looks okay to me
> > > geeta
> > >
> > >
> > >
> > >> Subject: Re: [rt.rap.ucar.edu #55790] AutoReply: Geeta met_help
> problem
> > >> From: met_help at ucar.edu
> > >> To: geeta124 at hotmail.com
> > >> Date: Wed, 11 Apr 2012 09:24:51 -0600
> > >>
> > >> Geeta,
> > >>
> > >> I created a met_help ticket for you. Please respond to this
email
> and describe the MET problem you are having in detail.
> > >>
> > >> Thanks,
> > >>
> > >> Paul
> > >>
> > >>
> > >> On 04/11/2012 09:21 AM, met_help at ucar.edu via RT wrote:
> > >>> Greetings,
> > >>>
> > >>> This message has been automatically generated in response to
the
> > >>> creation of a trouble ticket regarding:
> > >>> "Geeta met_help problem",
> > >>> a summary of which appears below.
> > >>>
> > >>> There is no need to reply to this message right now. Your
ticket
> has been
> > >>> assigned an ID of [rt.rap.ucar.edu #55790].
> > >>>
> > >>> Please include the string:
> > >>>
> > >>> [rt.rap.ucar.edu #55790]
> > >>>
> > >>> in the subject line of all future correspondence about this
> issue. To do so,
> > >>> you may reply to this message.
> > >>>
> > >>> Thank you,
> > >>> met_help at ucar.edu
> > >>>
> > >>>
>
-------------------------------------------------------------------------
> > >>> Dear administrator,
> > >>> I wish to seek help from the met_help, and I am sending my
> message for the last 3-4 days. But everytime I am told that there is
> some error that has occurred.
> > >>> I donot know whtheer my meage is reaching across or not? I
have
> already registered but my query is not answered.
> > >>>
> > >>> kindly help me
> > >>> Regards
> > >>> geeta
> > >>>
> > >>> Dr G. Agnihotri,
> > >>> Met Centre, IMD, Palace Road,
> > >>> Bangalore, India.
> > >>> 080-22244419
> > >>>
> > >>
> > >>
> > >
> >
> >
>
> ------------------------------------------------
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

Creating NetCDF Observation file: ../myascii2nc/precip24_2011061000.nc
Reading ASCII Observation file: precip24_2011061000.ascii
20110600
/oprn/model/wrf3/WRFV3/postprd
-pcpdir                  WRFPRS_d01.000           WRFPRS_d01.001
WRFPRS_d01.001.idx       WRFPRS_d01.002           WRFPRS_d01.003
WRFPRS_d01.004           WRFPRS_d01.005           WRFPRS_d01.006
WRFPRS_d01.007           WRFPRS_d01.008           WRFPRS_d01.009
WRFPRS_d01.010           WRFPRS_d01.011           WRFPRS_d01.012
WRFPRS_d01.013           WRFPRS_d01.014           WRFPRS_d01.015
WRFPRS_d01.016           WRFPRS_d01.017           WRFPRS_d01.018
WRFPRS_d01.019           WRFPRS_d01.020           WRFPRS_d01.021
WRFPRS_d01.022           WRFPRS_d01.023           WRFPRS_d01.024
WRFPRS_d01.025           WRFPRS_d01.026           WRFPRS_d01.027
WRFPRS_d01.028           WRFPRS_d01.029           WRFPRS_d01.030
WRFPRS_d01.031           WRFPRS_d01.032           WRFPRS_d01.033
WRFPRS_d01.034           WRFPRS_d01.035           WRFPRS_d01.036
WRFPRS_d01.037           WRFPRS_d01.038           WRFPRS_d01.039
WRFPRS_d01.040           WRFPRS_d01.041           WRFPRS_d01.042
WRFPRS_d01.043           WRFPRS_d01.044           WRFPRS_d01.045
WRFPRS_d01.046           WRFPRS_d01.047           WRFPRS_d01.048
WRFPRS_d01.049           WRFPRS_d01.050           WRFPRS_d01.051
copygb_hwrf.txt          eta_micro_lookup.dat     fort.110
fort.14                  griddef.out              itag
out1                     run_wrfpost              run_wrfpost-
donotchange
temp.ctl                 tempfile                 testfile
wrf_cntrl.parm           wrfpost_d01.000.out      wrfpost_d01.001.out
wrfpost_d01.002.out      wrfpost_d01.003.out      wrfpost_d01.004.out
wrfpost_d01.005.out      wrfpost_d01.006.out      wrfpost_d01.007.out
wrfpost_d01.008.out      wrfpost_d01.009.out      wrfpost_d01.010.out
wrfpost_d01.011.out      wrfpost_d01.012.out      wrfpost_d01.013.out
wrfpost_d01.014.out      wrfpost_d01.015.out      wrfpost_d01.016.out
wrfpost_d01.017.out      wrfpost_d01.018.out      wrfpost_d01.019.out
wrfpost_d01.020.out      wrfpost_d01.021.out      wrfpost_d01.022.out
wrfpost_d01.023.out      wrfpost_d01.024.out      wrfpost_d01.025.out
wrfpost_d01.026.out      wrfpost_d01.027.out      wrfpost_d01.028.out
wrfpost_d01.029.out      wrfpost_d01.030.out      wrfpost_d01.031.out
wrfpost_d01.032.out      wrfpost_d01.033.out      wrfpost_d01.034.out
wrfpost_d01.035.out      wrfpost_d01.036.out      wrfpost_d01.037.out
wrfpost_d01.038.out      wrfpost_d01.039.out      wrfpost_d01.040.out
wrfpost_d01.041.out      wrfpost_d01.042.out      wrfpost_d01.043.out
wrfpost_d01.044.out      wrfpost_d01.045.out      wrfpost_d01.046.out
wrfpost_d01.047.out      wrfpost_d01.048.out      wrfpost_d01.049.out
wrfpost_d01.050.out      wrfpost_d01.051.out      wrfpost_d01.052.out
wrfprs_d01.000           wrfprs_d01.001           wrfprs_d01.002
wrfprs_d01.003           wrfprs_d01.004           wrfprs_d01.005
wrfprs_d01.006           wrfprs_d01.007           wrfprs_d01.008
wrfprs_d01.009           wrfprs_d01.010           wrfprs_d01.011
wrfprs_d01.012           wrfprs_d01.013           wrfprs_d01.014
wrfprs_d01.015           wrfprs_d01.016           wrfprs_d01.017
wrfprs_d01.018           wrfprs_d01.019           wrfprs_d01.020
wrfprs_d01.021           wrfprs_d01.022           wrfprs_d01.023
wrfprs_d01.024           wrfprs_d01.025           wrfprs_d01.026
wrfprs_d01.027           wrfprs_d01.028           wrfprs_d01.029
wrfprs_d01.030           wrfprs_d01.031           wrfprs_d01.032
wrfprs_d01.033           wrfprs_d01.034           wrfprs_d01.035
wrfprs_d01.036           wrfprs_d01.037           wrfprs_d01.038
wrfprs_d01.039           wrfprs_d01.040           wrfprs_d01.041
wrfprs_d01.042           wrfprs_d01.043           wrfprs_d01.044
wrfprs_d01.045           wrfprs_d01.046           wrfprs_d01.047
wrfprs_d01.048           wrfprs_d01.049           wrfprs_d01.050
wrfprs_d01.051
NEWDATE 2011060900
YY 2011
-rwxr-xr-x    1 model    staff      11408030 Apr 19 14:37
WRFPRS_d01.000
NEWDATE 2011060901
YY 2011
-rwxr-xr-x    1 model    staff      12643498 Apr 19 14:37
WRFPRS_d01.001
NEWDATE 2011060902
YY 2011
-rwxr-xr-x    1 model    staff      12643500 Apr 19 14:37
WRFPRS_d01.002
NEWDATE 2011060903
YY 2011
-rwxr-xr-x    1 model    staff      12600320 Apr 19 14:37
WRFPRS_d01.003
NEWDATE 2011060904
YY 2011
-rwxr-xr-x    1 model    staff      12590344 Apr 19 14:38
WRFPRS_d01.004
NEWDATE 2011060905
YY 2011
-rwxr-xr-x    1 model    staff      12623562 Apr 19 14:38
WRFPRS_d01.005
NEWDATE 2011060906
YY 2011
-rwxr-xr-x    1 model    staff      12583710 Apr 19 14:38
WRFPRS_d01.006
NEWDATE 2011060907
YY 2011
-rwxr-xr-x    1 model    staff      12537212 Apr 19 14:38
WRFPRS_d01.007
NEWDATE 2011060908
YY 2011
-rwxr-xr-x    1 model    staff      12587026 Apr 19 14:38
WRFPRS_d01.008
NEWDATE 2011060909
YY 2011
-rwxr-xr-x    1 model    staff      12630192 Apr 19 14:39
WRFPRS_d01.009
NEWDATE 2011060910
YY 2011
-rwxr-xr-x    1 model    staff      12643474 Apr 19 14:39
WRFPRS_d01.010
NEWDATE 2011060911
YY 2011
-rwxr-xr-x    1 model    staff      12580376 Apr 19 14:39
WRFPRS_d01.011
NEWDATE 2011060912
YY 2011
-rwxr-xr-x    1 model    staff      12593662 Apr 19 14:39
WRFPRS_d01.012
NEWDATE 2011060913
YY 2011
-rwxr-xr-x    1 model    staff      12613590 Apr 19 14:39
WRFPRS_d01.013
NEWDATE 2011060914
YY 2011
-rwxr-xr-x    1 model    staff      12520594 Apr 19 14:40
WRFPRS_d01.014
NEWDATE 2011060915
YY 2011
-rwxr-xr-x    1 model    staff      12530568 Apr 19 14:40
WRFPRS_d01.015
NEWDATE 2011060916
YY 2011
-rwxr-xr-x    1 model    staff      12587014 Apr 19 14:40
WRFPRS_d01.016
NEWDATE 2011060917
YY 2011
-rwxr-xr-x    1 model    staff      12606944 Apr 19 14:40
WRFPRS_d01.017
NEWDATE 2011060918
YY 2011
-rwxr-xr-x    1 model    staff      12567086 Apr 19 14:40
WRFPRS_d01.018
NEWDATE 2011060919
YY 2011
-rwxr-xr-x    1 model    staff      12543850 Apr 19 14:41
WRFPRS_d01.019
NEWDATE 2011060920
YY 2011
-rwxr-xr-x    1 model    staff      12504000 Apr 19 14:41
WRFPRS_d01.020
NEWDATE 2011060921
YY 2011
-rwxr-xr-x    1 model    staff      12487396 Apr 19 14:41
WRFPRS_d01.021
NEWDATE 2011060922
YY 2011
-rwxr-xr-x    1 model    staff      12563782 Apr 19 14:41
WRFPRS_d01.022
NEWDATE 2011060923
YY 2011
-rwxr-xr-x    1 model    staff      12490712 Apr 19 14:41
WRFPRS_d01.023
NEWDATE 2011061000
YY 2011
-rwxr-xr-x    1 model    staff      12560452 Apr 19 14:42
WRFPRS_d01.024
NEWDATE 2011061001
YY 2011
-rwxr-xr-x    1 model    staff      12550486 Apr 19 14:42
WRFPRS_d01.025
NEWDATE 2011061002
YY 2011
-rwxr-xr-x    1 model    staff      12596992 Apr 19 14:42
WRFPRS_d01.026
NEWDATE 2011061003
YY 2011
-rwxr-xr-x    1 model    staff      12623568 Apr 19 14:42
WRFPRS_d01.027
NEWDATE 2011061004
YY 2011
-rwxr-xr-x    1 model    staff      12620248 Apr 19 14:42
WRFPRS_d01.028
NEWDATE 2011061005
YY 2011
-rwxr-xr-x    1 model    staff      12626896 Apr 19 14:43
WRFPRS_d01.029
NEWDATE 2011061006
YY 2011
-rwxr-xr-x    1 model    staff      12613612 Apr 19 14:43
WRFPRS_d01.030
NEWDATE 2011061007
YY 2011
-rwxr-xr-x    1 model    staff      12583714 Apr 19 14:43
WRFPRS_d01.031
NEWDATE 2011061008
YY 2011
-rwxr-xr-x    1 model    staff      12573740 Apr 19 14:43
WRFPRS_d01.032
NEWDATE 2011061009
YY 2011
-rwxr-xr-x    1 model    staff      12603626 Apr 19 14:43
WRFPRS_d01.033
NEWDATE 2011061010
YY 2011
-rwxr-xr-x    1 model    staff      12613596 Apr 19 14:44
WRFPRS_d01.034
NEWDATE 2011061011
YY 2011
-rwxr-xr-x    1 model    staff      12577060 Apr 19 14:44
WRFPRS_d01.035
NEWDATE 2011061012
YY 2011
-rwxr-xr-x    1 model    staff      12616916 Apr 19 14:44
WRFPRS_d01.036
NEWDATE 2011061013
YY 2011
-rwxr-xr-x    1 model    staff      12610266 Apr 19 14:44
WRFPRS_d01.037
NEWDATE 2011061014
YY 2011
-rwxr-xr-x    1 model    staff      12550486 Apr 19 14:44
WRFPRS_d01.038
NEWDATE 2011061015
YY 2011
-rwxr-xr-x    1 model    staff      12527246 Apr 19 14:45
WRFPRS_d01.039
NEWDATE 2011061016
YY 2011
-rwxr-xr-x    1 model    staff      12530566 Apr 19 14:45
WRFPRS_d01.040
NEWDATE 2011061017
YY 2011
-rwxr-xr-x    1 model    staff      12567092 Apr 19 14:45
WRFPRS_d01.041
NEWDATE 2011061018
YY 2011
-rwxr-xr-x    1 model    staff      12517278 Apr 19 14:45
WRFPRS_d01.042
NEWDATE 2011061019
YY 2011
-rwxr-xr-x    1 model    staff      12500672 Apr 19 14:45
WRFPRS_d01.043
NEWDATE 2011061020
YY 2011
-rwxr-xr-x    1 model    staff      12527242 Apr 19 14:46
WRFPRS_d01.044
NEWDATE 2011061021
YY 2011
-rwxr-xr-x    1 model    staff      12550488 Apr 19 14:46
WRFPRS_d01.045
NEWDATE 2011061022
YY 2011
-rwxr-xr-x    1 model    staff      12583708 Apr 19 14:46
WRFPRS_d01.046
NEWDATE 2011061023
YY 2011
-rwxr-xr-x    1 model    staff      12547174 Apr 19 14:46
WRFPRS_d01.047
NEWDATE 2011061100
YY 2011
-rwxr-xr-x    1 model    staff      12670052 Apr 19 14:46
WRFPRS_d01.048
NEWDATE 2011061101
YY 2011
-rwxr-xr-x    1 model    staff      12656766 Apr 19 14:47
WRFPRS_d01.049
NEWDATE 2011061102
YY 2011
-rwxr-xr-x    1 model    staff      12650132 Apr 19 14:47
WRFPRS_d01.050
NEWDATE 2011061103
YY 2011
-rwxr-xr-x    1 model    staff      12610284 Apr 19 14:47
WRFPRS_d01.051
Thu Apr 19 14:47:36 IST 2012
End of Output Job
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.027
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.003
Performing subtraction command.
Writing output file: ./myoutput/pcp_combine/24hr/model24_2011060900.nc
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.051
Reading input file: /oprn/model/wrf3/WRFV3/postprd/WRFPRS_d01.027
Performing subtraction command.
Writing output file: ./myoutput/pcp_combine/48hr/model48_2011060900.nc
Forecast File:
/oprn/model/wrf3/utils/met/METv3.0/myoutput/pcp_combine/24hr/model24_2011060900.nc
Climatology File: none
Configuration File:
/oprn/model/wrf3/utils/met/METv3.0/myscripts/PointStatConfig_rain
Observation File:
/oprn/model/wrf3/utils/met/METv3.0/myascii2nc/precip24_2011061000.nc

--------------------------------------------------------------------------------

Reading records for APCP_24(*,*).
For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.

--------------------------------------------------------------------------------

Searching 29 observations from 29 PrepBufr messages.

--------------------------------------------------------------------------------

Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
27 pairs.
Computing Categorical Statistics.
Computing Continuous Statistics.
Computing Scalar Partial Sums.

--------------------------------------------------------------------------------

Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V.stat
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_fho.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_ctc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_cts.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_cnt.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_sl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_sal1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_vl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_val1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pct.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pstd.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_pjc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_prc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_270000L_20110610_030000V_mpr.txt
Forecast File: ./myoutput/pcp_combine/48hr/model48_2011060900.nc
Climatology File: none
Configuration File: ./myscripts/PointStatConfig_rain
Observation File: ./myascii2nc/precip24_2011061100.nc

--------------------------------------------------------------------------------

Reading records for APCP_24(*,*).
For APCP_24(*,*) found 1 forecast levels and 0 climatology levels.

--------------------------------------------------------------------------------

Searching 27 observations from 27 PrepBufr messages.

--------------------------------------------------------------------------------

Processing APCP_24(*,*) versus APCP_24/A24, for observation type
ADPSFC, over region FULL, for interpolation method UW_MEAN(1), using
27 pairs.
Computing Categorical Statistics.
Computing Continuous Statistics.
Computing Scalar Partial Sums.

--------------------------------------------------------------------------------

Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V.stat
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_fho.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_ctc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_cts.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_cnt.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_sl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_sal1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_vl1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_val1l2.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pct.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pstd.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_pjc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_prc.txt
Output file:
/oprn/model/wrf3/utils/met/METv3.0/out/point_stat/point_stat_24hr_510000L_20110611_030000V_mpr.txt

------------------------------------------------
Subject: Geeta met_help problem
From: Geeta Geeta
Time: Thu Apr 19 23:09:07 2012

////////////////////////////////////////////////////////////////////////////////
//
// Default point_stat configuration file
//
////////////////////////////////////////////////////////////////////////////////