[Met_help] Grid_Stat error using NARR to do verification of WRF output

[David Andrade]

Hi,

I am running into some trouble when I try to run gridstat using NARR and wrf
output.
Essentially, I get the following error:

Forecast File: wrfout_regrid.grb
Observation File: narr_regrid.grb
Configuration File: GridStatConfig_default
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072124576585

ERROR: compute_cntinfo() -> the forecast and observation arrays must have
the same non-zero length!

terminate called after throwing an instance of 'int'
Aborted

where wrfout_regrid.grb is the WPS output regridded to the NARR grid using
copygb and narr_regrid.grb is a NARR grib file
that has been formed by concatenating several single-time NARR files
together. I think the error may have something to do
with the time labeling in each file. In the wrfoutput, subsequent times are
referred to as a forecast for +X hours from the initial time,
whereas, in the concatenated NARR data, each time has its own label. Am I
doing something wrong? Any help would be
greatly appreciated.


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

The following config_file was used:

model = "WRF";
vx_grib_code[] = [ "11/P500" ];
thresholds[] = [ "gt0.0 ge5.0" ];
mask_grids[] = [ "FULL" ];
mask_polys[] = [ "MET_BASE/data/poly/EAST.poly" ];
ci_alpha[] = [ 0.05 ];
boot_interval = 1;
boot_rep_prop = 1.0;
n_boot_rep = 1000;
boot_rng = "mt19937";
boot_seed = "";
interp_method[] = [ "UW_MEAN" ];
interp_width[] = [ 1 ];
interp_threshold = 1.0;
nbr_width[] = [ 3, 5 ];
nbr_threshold = 1.0;
nbr_frac_threshold[] = [ "ge0.5" ];
output_flag[] = [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1 ];
rank_corr_flag = 1;
ncep_defaults = 1;
tmp_dir = "./tmp";
version = "V1.1";

-- 
David Andrade
andraded at eas.slu.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/met_help/attachments/20081201/60e2bdba/attachment.html