[Met_help] MET v1.1 test case error

Tue Jul 15 20:48:28 MDT 2008

I downloaded the METv1.1 code to a 32-bit Red Hat Linux system.  I was able to successfully compile the code using the gnu g77, gcc, and g++ compilers.  I ran through the test suite calling the test_all.sh script. The tests ran sucessfully until it tried to read the sample prepbufr data.  The log from the test is attached.  The error is found in line 400 of the log and says that the BUFR file is not valid.  

I looked at the data file using the od command and the output is attached.  It appears the prepbufr file has an 8-byte header before the BUFR string.  On my 32-bit machine, the BUFR-reading programs will not read this file.  (I tried to read the file using readpb and got the same error.)

I did a quick check of a couple of prepbuf files on the NCEP web site.  The "blocked" files have four-byte headers, and the unblocked files have no header.

I speculate that either the sample data was blocked twice, or was blocked on a 64-bit machine (which makes files unreadable on at least 32-bit machines and maybe all of them).  I had similar issues when I first tried to get MET to run on a 64-bit machine.

I am at the WRF tutorial in Boulder this week and can show you my setup if desired.

Steve Masters  masters.steve at ensco.com
ENSCO, Inc.
Melbourne, FL

......................................................................
The information contained in this email message is intended only for the use of the individuals to whom it is
addressed and may contain information that is privileged and sensitive. If the reader of this message is not
the intended recipient, you are hereby notified that any dissemination, distribution or copying of this
communication is strictly prohibited. If you have received this communication in error, please notify the
sender immediately by email at the above referenced address. Thank you.
-------------- next part --------------

*** Testing the Model Evaluation Tools Project ***

*** Testing PCP_COMBINE application ***

*** Running PCP_COMBINE to combine 3-hourly APCP accumulation forecasts into a single 12 hour accumulation forecast ***
Performing sum command: init_time/in_accum/valid_time/out_accum Times = 1123372800/10800/1123416000/43200
Searching for 4 files with accumulation times of 3 hours to sum to a total accumulation time of 12 hours in ../data/sample_fcst/2005080700
[1] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
[2] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_09.tm00_G212
[3] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_06.tm00_G212
[4] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_03.tm00_G212
Writing output file: ../out/pcp_combine/sample_fcst_12L_2005080712V_12A.nc

*** Running PCP_COMBINE to combine 3-hourly APCP accumulation forecasts into a single 24 hour accumulation forecast ***
Performing sum command: init_time/in_accum/valid_time/out_accum Times = 1123372800/10800/1123459200/86400
Searching for 8 files with accumulation times of 3 hours to sum to a total accumulation time of 24 hours in ../data/sample_fcst/2005080700
[1] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_24.tm00_G212
[2] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_21.tm00_G212
[3] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_18.tm00_G212
[4] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_15.tm00_G212
[5] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
[6] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_09.tm00_G212
[7] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_06.tm00_G212
[8] Reading input file: ../data/sample_fcst/2005080700/wrfprs_ruc13_03.tm00_G212
Writing output file: ../out/pcp_combine/sample_fcst_24L_2005080800V_24A.nc

*** Running PCP_COMBINE to combine 1-hourly APCP accumulation observations into a single 12 hour accumulation observation ***
Performing sum command: init_time/in_accum/valid_time/out_accum Times = 0/3600/1123416000/43200
Searching for 12 files with accumulation times of 1 hours to sum to a total accumulation time of 12 hours in ../data/sample_obs/ST2ml
[1] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080712.Grb_G212
[2] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080711.Grb_G212
[3] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080710.Grb_G212
[4] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080709.Grb_G212
[5] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080708.Grb_G212
[6] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080707.Grb_G212
[7] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080706.Grb_G212
[8] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080705.Grb_G212
[9] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080704.Grb_G212
[10] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080703.Grb_G212
[11] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080702.Grb_G212
[12] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080701.Grb_G212
Writing output file: ../out/pcp_combine/sample_obs_2005080712V_12A.nc

*** Running PCP_COMBINE to combine 1-hourly APCP accumulation observations into a single 24 hour accumulation observation ***
Performing sum command: init_time/in_accum/valid_time/out_accum Times = 0/3600/1123459200/86400
Searching for 24 files with accumulation times of 1 hours to sum to a total accumulation time of 24 hours in ../data/sample_obs/ST2ml
[1] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080800.Grb_G212
[2] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080723.Grb_G212
[3] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080722.Grb_G212
[4] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080721.Grb_G212
[5] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080720.Grb_G212
[6] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080719.Grb_G212
[7] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080718.Grb_G212
[8] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080717.Grb_G212
[9] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080716.Grb_G212
[10] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080715.Grb_G212
[11] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080714.Grb_G212
[12] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080713.Grb_G212
[13] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080712.Grb_G212
[14] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080711.Grb_G212
[15] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080710.Grb_G212
[16] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080709.Grb_G212
[17] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080708.Grb_G212
[18] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080707.Grb_G212
[19] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080706.Grb_G212
[20] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080705.Grb_G212
[21] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080704.Grb_G212
[22] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080703.Grb_G212
[23] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080702.Grb_G212
[24] Reading input file: ../data/sample_obs/ST2ml/ST2ml2005080701.Grb_G212
Writing output file: ../out/pcp_combine/sample_obs_2005080800V_24A.nc

*** Testing MODE application ***

*** Running MODE on APCP using netCDF input for both forecast and observation ***
Forecast File: ../out/pcp_combine/sample_fcst_12L_2005080712V_12A.nc
Observation File: ../out/pcp_combine/sample_obs_2005080712V_12A.nc
Match Config File: config/WrfModeConfig_APCP_12
Merge Config File: config/WrfModeConfig_APCP_12
NetCDF APCP Variable (min, max) = (0, 38.8841)
NetCDF APCP Variable (min, max) = (0, 104.1)
Forecast Field: APCP_12 at SFC
Observation Field: APCP_12 at SFC
Identifying objects in the forecast and observation fields...
Computing contingency table statistics...
Identified: 7 forecast objects and 4 observation objects.
Performing merging (threshold merging) in the forecast field.
Performing merging (threshold merging) in the observation field.
Remaining: 7 forecast objects and 4 observation objects.
Performing matching between the forecast and observation fields.
Creating Contingency Table Statistics file: ../out/mode/mode_APCP_12_SFC_vs_APCP_12_SFC_120000L_20050807_120000V_120000A_cts.txt
Creating Fcst-Obs Object Statistics file: ../out/mode/mode_APCP_12_SFC_vs_APCP_12_SFC_120000L_20050807_120000V_120000A_obj.txt
Creating Object NetCDF file: ../out/mode/mode_APCP_12_SFC_vs_APCP_12_SFC_120000L_20050807_120000V_120000A_obj.nc
Loading forecast raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Loading observation raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Creating postscript file: ../out/mode/mode_APCP_12_SFC_vs_APCP_12_SFC_120000L_20050807_120000V_120000A.ps

*** Running MODE on APCP using a Grib forecast and netCDF observation ***
Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_24.tm00_G212
Observation File: ../out/pcp_combine/sample_obs_2005080800V_24A.nc
Match Config File: config/WrfModeConfig_APCP_24
Merge Config File: config/WrfModeConfig_APCP_24
Searching grib forecast file for grib code 61 with accumulation/level indicator of A24.
NetCDF APCP Variable (min, max) = (0, 93.1002)
Forecast Field: APCP_24 at SFC
Observation Field: APCP_24 at SFC
Identifying objects in the forecast and observation fields...
Computing contingency table statistics...
Identified: 4 forecast objects and 6 observation objects.
Performing merging (no merging) in the forecast field.
Performing merging (no merging) in the observation field.
Remaining: 4 forecast objects and 6 observation objects.
Performing matching between the forecast and observation fields.
Creating Contingency Table Statistics file: ../out/mode/mode_APCP_24_SFC_vs_APCP_24_SFC_240000L_20050808_000000V_240000A_cts.txt
Creating Fcst-Obs Object Statistics file: ../out/mode/mode_APCP_24_SFC_vs_APCP_24_SFC_240000L_20050808_000000V_240000A_obj.txt
Creating Object NetCDF file: ../out/mode/mode_APCP_24_SFC_vs_APCP_24_SFC_240000L_20050808_000000V_240000A_obj.nc
Loading forecast raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Loading observation raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Creating postscript file: ../out/mode/mode_APCP_24_SFC_vs_APCP_24_SFC_240000L_20050808_000000V_240000A.ps

*** Running MODE on RH at 500mb using Grib input for both forecast and observation ***
Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
Observation File: ../data/sample_fcst/2005080712/wrfprs_ruc13_00.tm00_G212
Match Config File: config/WrfModeConfig_RH
Merge Config File: config/WrfModeConfig_RH
Searching grib forecast file for grib code 52 with accumulation/level indicator of P500-500.
Searching grib observation file for grib code 52 with accumulation/level indicator of P500-500.
Forecast Field: RH at ISBL_500
Observation Field: RH at ISBL_500
Identifying objects in the forecast and observation fields...
Computing contingency table statistics...
Identified: 8 forecast objects and 12 observation objects.
Performing merging (no merging) in the forecast field.
Performing merging (no merging) in the observation field.
Remaining: 8 forecast objects and 12 observation objects.
Performing matching between the forecast and observation fields.
Creating Contingency Table Statistics file: ../out/mode/mode_RH_ISBL_500_vs_RH_ISBL_500_120000L_20050807_120000V_000000A_cts.txt
Creating Fcst-Obs Object Statistics file: ../out/mode/mode_RH_ISBL_500_vs_RH_ISBL_500_120000L_20050807_120000V_000000A_obj.txt
Creating Object NetCDF file: ../out/mode/mode_RH_ISBL_500_vs_RH_ISBL_500_120000L_20050807_120000V_000000A_obj.nc
Loading forecast raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Loading observation raw color table: /home/stevem/sw/METv1.1/data/colortables/mode_raw.ctable
Creating postscript file: ../out/mode/mode_RH_ISBL_500_vs_RH_ISBL_500_120000L_20050807_120000V_000000A.ps

*** Testing GRID_STAT application ***

*** Running GRID_STAT using netCDF input for both forecast and observation ***
Forecast File: ../out/pcp_combine/sample_fcst_12L_2005080712V_12A.nc
Observation File: ../out/pcp_combine/sample_obs_2005080712V_12A.nc
Configuration File: config/GridStatConfig_APCP_12
NetCDF APCP Variable (min, max) = (0, 38.8841)
NetCDF APCP Variable (min, max) = (0, 104.1)
Creating Output VSDB file: ../out/grid_stat/grid_stat_120000L_20050807_120000V.vsdb
Creating Output FHO file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_fho.txt
Creating Output CTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctc.txt
Creating Output CTP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctp.txt
Creating Output CFP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cfp.txt
Creating Output COP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cop.txt
Creating Output CTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cts.txt
Creating Output CNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cnt.txt
Creating Output SL1L2 file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_sl1l2.txt
Creating Output NBRCTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrctc.txt
Creating Output NBRCTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcts.txt
Creating Output NBRCNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcnt.txt
Creating Output NetCDF file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_pairs.nc

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

Processing Grib Code 61 with accumulation/level indicator of A12
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A12, for interpolation method UW_MEAN(1), over region DTC165, using 2650 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A12, for interpolation method UW_MEAN(1), over region DTC166, using 3807 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A12, for interpolation method UW_MEAN(1), over region MDW, using 914 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A12, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >0.000, over region DTC165 using 2360 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >0.000, over region DTC166 using 3501 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >0.000, over region MDW using 878 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >0.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >=5.000, over region DTC165 using 2360 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >=5.000, over region DTC166 using 3501 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >=5.000, over region MDW using 878 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 3, raw threshold of >=5.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >0.000, over region DTC165 using 2070 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >0.000, over region DTC166 using 3208 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >0.000, over region MDW using 829 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >0.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >=5.000, over region DTC165 using 2070 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >=5.000, over region DTC166 using 3208 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >=5.000, over region MDW using 829 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A12, for neighborhood width of 5, raw threshold of >=5.000, over region LMV using 521 pairs

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

Output VSDB file: ../out/grid_stat/grid_stat_120000L_20050807_120000V.vsdb
Output FHO file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_fho.txt
Output CTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctc.txt
Output CTP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctp.txt
Output CFP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cfp.txt
Output COP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cop.txt
Output CTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cts.txt
Output CNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cnt.txt
Output SL1L2 file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_sl1l2.txt
Output NBRCTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrctc.txt
Output NBRCTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcts.txt
Output NBRCNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcnt.txt
Output NetCDF file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_pairs.nc

*** Running GRID_STAT using a Grib forecast and netCDF observation ***
Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_24.tm00_G212
Observation File: ../out/pcp_combine/sample_obs_2005080800V_24A.nc
Configuration File: config/GridStatConfig_APCP_24
NetCDF APCP Variable (min, max) = (0, 93.1002)
Creating Output VSDB file: ../out/grid_stat/grid_stat_240000L_20050808_000000V.vsdb
Creating Output FHO file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_fho.txt
Creating Output CTC file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_ctc.txt
Creating Output CTP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_ctp.txt
Creating Output CFP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cfp.txt
Creating Output COP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cop.txt
Creating Output CTS file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cts.txt
Creating Output CNT file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cnt.txt
Creating Output SL1L2 file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_sl1l2.txt
Creating Output NBRCTC file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrctc.txt
Creating Output NBRCTS file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrcts.txt
Creating Output NBRCNT file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrcnt.txt
Creating Output NetCDF file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_pairs.nc

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

Processing Grib Code 61 with accumulation/level indicator of A24
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A24, for interpolation method UW_MEAN(1), over region DTC165, using 2645 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A24, for interpolation method UW_MEAN(1), over region DTC166, using 3764 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A24, for interpolation method UW_MEAN(1), over region MDW, using 914 pairs
Computing CNT, CTS, and SL1L2 for Grib Field APCP at A24, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >0.000, over region DTC165 using 2360 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >0.000, over region DTC166 using 3466 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >0.000, over region MDW using 878 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >0.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >=10.000, over region DTC165 using 2360 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >=10.000, over region DTC166 using 3466 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >=10.000, over region MDW using 878 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 3, raw threshold of >=10.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >0.000, over region DTC165 using 2070 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >0.000, over region DTC166 using 3181 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >0.000, over region MDW using 829 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >0.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >=10.000, over region DTC165 using 2070 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >=10.000, over region DTC166 using 3181 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >=10.000, over region MDW using 829 pairs
Computing NBRCNT and NBRCTS for Grib Field APCP at A24, for neighborhood width of 5, raw threshold of >=10.000, over region LMV using 521 pairs

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

Output VSDB file: ../out/grid_stat/grid_stat_240000L_20050808_000000V.vsdb
Output FHO file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_fho.txt
Output CTC file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_ctc.txt
Output CTP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_ctp.txt
Output CFP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cfp.txt
Output COP file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cop.txt
Output CTS file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cts.txt
Output CNT file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_cnt.txt
Output SL1L2 file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_sl1l2.txt
Output NBRCTC file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrctc.txt
Output NBRCTS file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrcts.txt
Output NBRCNT file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_nbrcnt.txt
Output NetCDF file: ../out/grid_stat/grid_stat_240000L_20050808_000000V_pairs.nc

*** Running GRID_STAT using Grib input for both forecast and observation ***
Forecast File: ../data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212
Observation File: ../data/sample_fcst/2005080712/wrfprs_ruc13_00.tm00_G212
Configuration File: config/GridStatConfig_all
Creating Output VSDB file: ../out/grid_stat/grid_stat_120000L_20050807_120000V.vsdb
Creating Output FHO file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_fho.txt
Creating Output CTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctc.txt
Creating Output CTP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctp.txt
Creating Output CFP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cfp.txt
Creating Output COP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cop.txt
Creating Output CTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cts.txt
Creating Output CNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cnt.txt
Creating Output SL1L2 file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_sl1l2.txt
Creating Output NBRCTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrctc.txt
Creating Output NBRCTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcts.txt
Creating Output NBRCNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcnt.txt
Creating Output NetCDF file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_pairs.nc

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

Processing Grib Code 11 with accumulation/level indicator of Z2
Computing CNT, CTS, and SL1L2 for Grib Field TMP at Z2, for interpolation method UW_MEAN(1), over region DTC165, using 5555 pairs
Computing CNT, CTS, and SL1L2 for Grib Field TMP at Z2, for interpolation method UW_MEAN(1), over region DTC166, using 6863 pairs
Computing CNT, CTS, and SL1L2 for Grib Field TMP at Z2, for interpolation method UW_MEAN(1), over region MDW, using 932 pairs
Computing CNT, CTS, and SL1L2 for Grib Field TMP at Z2, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of <300.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of <300.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of <300.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of <300.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of >=300.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of >=300.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of >=300.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 3, raw threshold of >=300.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of <300.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of <300.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of <300.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of <300.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of >=300.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of >=300.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of >=300.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field TMP at Z2, for neighborhood width of 5, raw threshold of >=300.000, over region LMV using 521 pairs

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

Processing Grib Code 52 with accumulation/level indicator of Z2
Computing CNT, CTS, and SL1L2 for Grib Field RH at Z2, for interpolation method UW_MEAN(1), over region DTC165, using 5555 pairs
Computing CNT, CTS, and SL1L2 for Grib Field RH at Z2, for interpolation method UW_MEAN(1), over region DTC166, using 6863 pairs
Computing CNT, CTS, and SL1L2 for Grib Field RH at Z2, for interpolation method UW_MEAN(1), over region MDW, using 932 pairs
Computing CNT, CTS, and SL1L2 for Grib Field RH at Z2, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 3, raw threshold of >=80.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 3, raw threshold of >=80.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 3, raw threshold of >=80.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 3, raw threshold of >=80.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 5, raw threshold of >=80.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 5, raw threshold of >=80.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 5, raw threshold of >=80.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field RH at Z2, for neighborhood width of 5, raw threshold of >=80.000, over region LMV using 521 pairs

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

Processing Grib Code 33 with accumulation/level indicator of Z10
Computing CNT, CTS, and SL1L2 for Grib Field UGRD at Z10, for interpolation method UW_MEAN(1), over region DTC165, using 5555 pairs
Computing CNT, CTS, and SL1L2 for Grib Field UGRD at Z10, for interpolation method UW_MEAN(1), over region DTC166, using 6863 pairs
Computing CNT, CTS, and SL1L2 for Grib Field UGRD at Z10, for interpolation method UW_MEAN(1), over region MDW, using 932 pairs
Computing CNT, CTS, and SL1L2 for Grib Field UGRD at Z10, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field UGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region LMV using 521 pairs

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

Processing Grib Code 34 with accumulation/level indicator of Z10
Computing CNT, CTS, and SL1L2 for Grib Field VGRD at Z10, for interpolation method UW_MEAN(1), over region DTC165, using 5555 pairs
Computing CNT, CTS, and SL1L2 for Grib Field VGRD at Z10, for interpolation method UW_MEAN(1), over region DTC166, using 6863 pairs
Computing CNT, CTS, and SL1L2 for Grib Field VGRD at Z10, for interpolation method UW_MEAN(1), over region MDW, using 932 pairs
Computing CNT, CTS, and SL1L2 for Grib Field VGRD at Z10, for interpolation method UW_MEAN(1), over region LMV, using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=2.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 3, raw threshold of >=5.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=2.000, over region LMV using 521 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region DTC165 using 5555 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region DTC166 using 6863 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region MDW using 932 pairs
Computing NBRCNT and NBRCTS for Grib Field VGRD at Z10, for neighborhood width of 5, raw threshold of >=5.000, over region LMV using 521 pairs

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

Output VSDB file: ../out/grid_stat/grid_stat_120000L_20050807_120000V.vsdb
Output FHO file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_fho.txt
Output CTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctc.txt
Output CTP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_ctp.txt
Output CFP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cfp.txt
Output COP file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cop.txt
Output CTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cts.txt
Output CNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_cnt.txt
Output SL1L2 file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_sl1l2.txt
Output NBRCTC file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrctc.txt
Output NBRCTS file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcts.txt
Output NBRCNT file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_nbrcnt.txt
Output NetCDF file: ../out/grid_stat/grid_stat_120000L_20050807_120000V_pairs.nc

*** Testing PB2NC application ***

*** Running PB2NC on a fortran-blocked PrepBufr file ***
Reading Config File:	config/PB2NCConfig_G212
Creating NetCDF File:	../out/pb2nc/sample_pb.nc
Reading PrepBufr File:	../data/sample_obs/prepbufr/ndas.t00z.prepbufr.tm12.20070401.nr.blk

 **************BUFR ARCHIVE LIBRARY ABORT*****************
 BUFRLIB: OPENBF - FIRST 4 BYTES READ FROM RECORD IN INPUT FILE CONNECTED TO UNIT  11 NOT 'BUFR', DOES NOT CONTAIN BUFR DATA     
 **************BUFR ARCHIVE LIBRARY ABORT*****************

*** Testing ASCII2NC application ***

*** Running ASCII2NC to reformat ASCII point observations into NetCDF ***
Creating NetCDF Observation file: ../out/ascii2nc/sample_ascii.nc
Reading ASCII Observation file: ../data/sample_obs/ascii/sample_ascii_obs.txt
Processing 2140 observations.
Finished processing 2140 observations.

*** Testing POINT_STAT application ***

*** Running POINT_STAT on sample NAM data ***
Forecast File:	../data/sample_fcst/2007033000/nam.t00z.awip1236.tm00.20070330.grb
Configuration File:	config/PointStatConfig
Climatology File:	none
Observation File:	../out/pb2nc/sample_pb.nc
Observation File:	../out/ascii2nc/sample_ascii.nc
Creating Output VSDB file:	../out/point_stat/point_stat_360000L_20070331_120000V.vsdb
Creating Output FHO file:	../out/point_stat/point_stat_360000L_20070331_120000V_fho.txt
Creating Output CTC file:	../out/point_stat/point_stat_360000L_20070331_120000V_ctc.txt
Creating Output CTP file:	../out/point_stat/point_stat_360000L_20070331_120000V_ctp.txt
Creating Output CFP file:	../out/point_stat/point_stat_360000L_20070331_120000V_cfp.txt
Creating Output COP file:	../out/point_stat/point_stat_360000L_20070331_120000V_cop.txt
Creating Output CTS file:	../out/point_stat/point_stat_360000L_20070331_120000V_cts.txt
Creating Output CNT file:	../out/point_stat/point_stat_360000L_20070331_120000V_cnt.txt
Creating Output SL1L2 file:	../out/point_stat/point_stat_360000L_20070331_120000V_sl1l2.txt
Creating Output SAL1L2 file:	../out/point_stat/point_stat_360000L_20070331_120000V_sal1l2.txt
Creating Output VL1L2 file:	../out/point_stat/point_stat_360000L_20070331_120000V_vl1l2.txt
Creating Output VAL1L2 file:	../out/point_stat/point_stat_360000L_20070331_120000V_val1l2.txt
Creating Output MPR file:	../out/point_stat/point_stat_360000L_20070331_120000V_mpr.txt

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

Reading records for Grib Field TMP at P900-750
For Grib Field TMP at P900-750 found 4 forecast levels and 0 climatology levels

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

Reading records for Grib Field UGRD at Z10
For Grib Field UGRD at Z10 found 1 forecast levels and 0 climatology levels

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

Reading records for Grib Field VGRD at Z10
For Grib Field VGRD at Z10 found 1 forecast levels and 0 climatology levels

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

*** Testing VSDB_ANALYSIS application ***

*** Running VSDB_ANALYSIS ***
Creating vsdb_analysis output file "../out/vsdb_analysis/vsdb_analysis.out"
VSDB Lines read     = 672
VSDB Lines retained = 672
Job 1: -job vsdb_job_filter -vx_mask G212 -line_type SL1L2 -dump_row ../out/vsdb_analysis/job_filter_SL1L2.vsdb
Creating VSDB output file "../out/vsdb_analysis/job_filter_SL1L2.vsdb"
Job 1 used 0 out of 672 VSDB lines.
Job 2: -job vsdb_job_summary -line_type CNT/0.050 -var TMP -column ME -dump_row ../out/vsdb_analysis/job_summary_ME.vsdb -alpha 0.05 -boot_interval 1 -boot_rep_prop 1.00 -n_boot_rep 1000 -boot_rng mt19937 -boot_seed 
Creating VSDB output file "../out/vsdb_analysis/job_summary_ME.vsdb"
Job 2 used 4 out of 672 VSDB lines.
Job 3: -job vsdb_job_aggr -vx_mask DTC165 -vx_mask DTC166 -line_type SL1L2 -var TMP -dump_row ../out/vsdb_analysis/job_aggr_SL1L2.vsdb
Creating VSDB output file "../out/vsdb_analysis/job_aggr_SL1L2.vsdb"
Job 3 used 2 out of 672 VSDB lines.
Job 4: -job vsdb_job_aggr_cnt -vx_mask DTC165 -vx_mask DTC166 -line_type SL1L2 -var TMP -dump_row ../out/vsdb_analysis/job_aggr_cnt_SL1L2.vsdb
Creating VSDB output file "../out/vsdb_analysis/job_aggr_cnt_SL1L2.vsdb"
Job 4 used 2 out of 672 VSDB lines.
Job 5: -job vsdb_job_aggr -vx_mask DTC165 -vx_mask DTC166 -line_type FHO<300.000 -var TMP -dump_row ../out/vsdb_analysis/job_aggr_FHO.vsdb
Creating VSDB output file "../out/vsdb_analysis/job_aggr_FHO.vsdb"
Job 5 used 2 out of 672 VSDB lines.
Job 6: -job vsdb_job_aggr_cts -vx_mask DTC165 -vx_mask DTC166 -line_type FHO<300.000 -var TMP -dump_row ../out/vsdb_analysis/job_aggr_cts_FHO.vsdb
Creating VSDB output file "../out/vsdb_analysis/job_aggr_cts_FHO.vsdb"
Job 6 used 2 out of 672 VSDB lines.

*** Testing MODE_ANALYSIS application ***

*** Running MODE_ANALYSIS to compute column summaries for simple forecast objects ***

Total mode lines read =  24
Total mode lines kept =   7

*** Running MODE_ANALYSIS to compute column summaries for simple observation objects ***

Total mode lines read =  24
Total mode lines kept =   4

*** Running MODE_ANALYSIS to summarize by case ***

Total mode lines read =  44
Total mode lines kept =  21

*** Finished Testing the Model Evaluation Tools Project ***
-------------- next part --------------
0000000 270   &  \0  \0  \0  \0  \0  \0   B   U   F   R  \0   & 262 003
0000020  \0  \0 022  \0 003  \a  \0  \0  \v 001  \f 001  \0  \0  \0  \0
0000040  \0  \0  \0  \0   &  \0  \0 001 200   C  \0 037 001  \0 001  \0
0000060 002  \0 003   A  \0 037 001 300 004   E  \0 037 001 300 003 205
0000100   @   A  \0 037 001  \0 036  \0  \0   &   n  \0 022   2   3   0