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

[John Halley Gotway]

David,

It'd probably be easiest for you to just send me a sample forecast file, a sample observation file, and the config file you're using.

I'll run it through Grid-Stat in debug mode and figure out what's going on.  I'm guessing that for whatever reason, no matched pairs were found.  But I'll take a look.

You can either email the files to me or post them on our anonymous ftp site:
ftp ftp.rap.ucar.edu
username = anonymous
password = "your email address"
cd incoming/irap/johnhg
put "names of files"

Thanks,
John Halley-Gotway
johnhg at ucar.edu

David Andrade wrote:
> 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";
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help