[Met_news] Troubles with point_stat

Model Evaluation Tools (MET) Announcements met_news at mailman.ucar.edu
Wed Feb 3 03:34:53 MST 2016


Dear MET users,

I have a problem with, I think, the point_stat tool.

Here a brief history of what I did, just to check if I made some former
mistake.

I made a 24 h simulation with ARW initialised with GFS and starting
2016020100.
I setted the  io_form_history=5 to get a grib1 wrfout file (called
wrfout.grb)

I created an ascii file containing the accumulated precipitation in three
hours for
the same dates filling the eleven fields described in table 4.4 of
Met_user_guide_v5.1
just for two stations.
Here is the file (ascii_obs.txt)
NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA  2.6
NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA  2.6
NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA  2.9
NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA  3.6
NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA  7.1
NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA 12.0
NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA 12.5
NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA 14.2
NA AA000A002 20160201_030000 8.832 40.941   6 061 3   6 NA  4.2
NA AA000A002 20160201_060000 8.832 40.941   6 061 3   6 NA  5.6
NA AA000A002 20160201_090000 8.832 40.941   6 061 3   6 NA 12.9
NA AA000A002 20160201_120000 8.832 40.941   6 061 3   6 NA 13.1
NA AA000A002 20160201_150000 8.832 40.941   6 061 3   6 NA 13.1
NA AA000A002 20160201_180000 8.832 40.941   6 061 3   6 NA 14.0
NA AA000A002 20160201_210000 8.832 40.941   6 061 3   6 NA 15.5
NA AA000A002 20160202_000000 8.832 40.941   6 061 3   6 NA 16.1

Afterward I created an nc file using the tool
./ascii2nc ascii_obs.txt ascii_obs.nc -v 2

The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.

I modified PointStatConfig_default like this

fcst = {
   sid_exc      = [];

   field = [
      {
        name       = "APCP";
        level      = [ "A03" ];
        cat_thresh = [ >1.0 ];
        file_type  = [ 1 ] ;
        init_time  = [ "20160201_000000" ];
        valid_time = [ "20160202_000000" ];
      }
   ];

}
obs = fcst;

and using  "./point_stat wrfout.grb ascii_obs.nc PointStatConfig_modified
-outdir ./  -v 2"
the program just give me the headers of the statistic files and the
following message error

DEBUG 1: User Config File: PointStatConfig_modified
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1087746390
DEBUG 1: Forecast File: wrfout_gfs.grib1
DEBUG 1: Observation File: prova_ascii_grib1.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for APCP/A03.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching APCP/A03 found in
file: wrfout_gfs.grib1
WARNING:
.....


I also tried to modify the input file like this

fcst = {
   sid_exc      = [];

   field = [
      {
        name       = "APCP_03";
        level      = [ "(*,*)" ];
        cat_thresh = [ >1.0 ];
        file_type  = [ 1 ] ;
        init_time  = [ "20160201_000000" ];
        valid_time = [ "20160202_000000" ];
      }
   ];

}

obs = {
   message_type = [ "ADPSFC" ];
   sid_exc      = [];

   field = [
      {
        name       = "APCP";
        level      = [ "A03" ];
        cat_thresh = [ >1.0 ];
        file_type  = [ 1 ] ;
        init_time  = [ "20160201_000000" ];
        valid_time = [ "20160202_000000" ];
      }
   ];

}

but in this case the message error is even worse and I get

DEBUG 1: User Config File: PointStatConfig_modified
ERROR  :
ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field abbreviation
'APCP_03' for table version 2
ERROR  :

without any statistic output file.


Any advise is really welcome.

Thanks,
Samuele
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/met_news/attachments/20160203/540162f7/attachment.html 


More information about the Met_news mailing list