[Met_help] [rt.rap.ucar.edu #71731] History for Problem with point stat tool

John Halley Gotway via RT met_help at ucar.edu
Wed Aug 26 09:36:34 MDT 2015


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

Dear all,

I received the following error message while running the point_stat tool of met 5.0.I am trying to use Point-Stat to compare the NetCDF output of PCP-
Combine to point observations generated using the ASCII2NC tool. My configuration file set were:fcst_field  = [ "APCP_03(*,*)" ];
obs_field   = [ "APCP/A3" ];
 
when configuration fiel set 
fcst_field[]  = [ "APCP_03(*,*)" ];
obs_field[]   = [ "APCP/A3" ];
the error was 
ERROR  : yyerror() -> syntax error in file "./PointStatConfig_3hrpre.d02"
ERROR  : 
ERROR  :    line   = 20
ERROR  : 
ERROR  :    column = 11
ERROR  : 
ERROR  :    text   = "["
ERROR  : 
ERROR  : 
ERROR  : fcst_field[] = [ "APCP_03(*,*)" ];
ERROR  : __________^_______________________
ERROR  : 
 
Kindly advise.

Thanks in advance.
 
Yuanpu
 
 
DEBUG 1: Default Config File: /data1/home/liuyp/software/met-5.0/met-install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_3hrpre.d02
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1
DEBUG 1: Forecast File: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: /data1/home/liuyp/software/met-5.0/ascii2nc/output/SURF_obs_ascii_201504.nc
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFHP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching SPFHP500 found in file: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMPP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching TMPP500 found in file: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for HGTP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching HGTP500 found in file: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for UGRDP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching UGRDP500 found in file: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for VGRDP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching VGRDP500 found in file: /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 523850 observations from 39657 messages.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V.stat
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V_fho.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V_ctc.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V_cts.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V_cnt.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-39L_19700101_000000V_mpr.txt







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

Subject: Problem with point stat tool
From: 柳媛普
Time: Sun May 03 21:49:31 2015





Dear all,

I received the following error message while running the point_stat
tool of met 5.0.
I am trying to use Point-Stat to compare the NetCDF output of PCP-
Combine to point observations generated using the ASCII2NC tool.

My configuration file set were:
fcst_field  = [ "APCP_03(*,*)" ];
obs_field   = [ "APCP/A3" ];

when configuration fiel set
fcst_field[]  = [ "APCP_03(*,*)" ];
obs_field[]   = [ "APCP/A3" ];
the error was
ERROR  : yyerror() -> syntax error in file
"./PointStatConfig_3hrpre.d02"
ERROR  :
ERROR  :    line   = 20
ERROR  :
ERROR  :    column = 11
ERROR  :
ERROR  :    text   = "["
ERROR  :
ERROR  :
ERROR  : fcst_field[] = [ "APCP_03(*,*)" ];
ERROR  : __________^_______________________
ERROR  :

Kindly advise.

Thanks in advance.

Yuanpu


DEBUG 1: Default Config File: /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_3hrpre.d02
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1
DEBUG 1: Forecast File: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: /data1/home/liuyp/software/met-
5.0/ascii2nc/output/SURF_obs_ascii_201504.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFHP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching SPFHP500
found in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMPP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching TMPP500
found in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for HGTP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching HGTP500
found in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for UGRDP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching UGRDP500
found in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for VGRDP500.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching VGRDP500
found in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/WRFPRS_APCP03_d02.03.nc
WARNING:
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 523850 observations from 39657 messages.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V.stat
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_fho.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_ctc.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cts.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cnt.txt
DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_mpr.txt




------------------------------------------------
Subject: Problem with point stat tool
From: Julie Prestopnik
Time: Mon May 04 09:56:13 2015

I'm sorry to hear you're having problems running point_stat.  There
appears
to be a problem with the way your configuration file is set up.

I can try to help some with the information you provided.  You can
also
look at the default config file:

/data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default

to see how it is set up.

Without seeing your configuration file, I can only refer to the fields
you
include in your email.  You mentioned that you have:

fcst_field[]  = [ "APCP_03(*,*)" ];
obs_field[]   = [ "APCP/A3" ];

This isn't the correct syntax for the point stat configuration.
Instead,
you want to have something like this with the wind_thresh and
cat_thresh
set however you'd like, but following this sort of format:

fcst = {
     wind_thresh  = [ NA ];

     field = [
        {
          name       = "APCP_03";
          level      = [ "(*,*)" ];
          cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0, >=30.0
];
        }
     ];
};
obs = {
     wind_thresh  = [ NA ];
     message_type = ["ADPSFC"];

     field = [
        {
          name       = "APCP";
          level      = [ "A03" ];
          cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0, >=30.0
];
        }
     ];
};


If setting the fields this way doesn't solve your problem, please let
me
know and please include your configuration files as well.

Thanks,
Julie

On Sun, May 3, 2015 at 9:33 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> Sun May 03 21:33:42 2015: Request 71731 was acted upon.
> Transaction: Ticket created by lyp at lzb.ac.cn
>        Queue: met_help
>      Subject: Problem with point stat tool
>        Owner: Nobody
>   Requestors: lyp at lzb.ac.cn
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
>
> Dear all,
>
> I received the following error message while running the point_stat
tool
> of met 5.0.I am trying to use Point-Stat to compare the NetCDF
output of
> PCP-
> Combine to point observations generated using the ASCII2NC tool. My
> configuration file set were:fcst_field  = [ "APCP_03(*,*)" ];
> obs_field   = [ "APCP/A3" ];
>
> when configuration fiel set
> fcst_field[]  = [ "APCP_03(*,*)" ];
> obs_field[]   = [ "APCP/A3" ];
> the error was
> ERROR  : yyerror() -> syntax error in file
"./PointStatConfig_3hrpre.d02"
> ERROR  :
> ERROR  :    line   = 20
> ERROR  :
> ERROR  :    column = 11
> ERROR  :
> ERROR  :    text   = "["
> ERROR  :
> ERROR  :
> ERROR  : fcst_field[] = [ "APCP_03(*,*)" ];
> ERROR  : __________^_______________________
> ERROR  :
>
> Kindly advise.
>
> Thanks in advance.
>
> Yuanpu
>
>
> DEBUG 1: Default Config File:
> /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_3hrpre.d02
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1
> DEBUG 1: Forecast File:
> /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File:
> /data1/home/liuyp/software/met-
5.0/ascii2nc/output/SURF_obs_ascii_201504.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SPFHP500.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching SPFHP500
found
> in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for TMPP500.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching TMPP500
found in
> file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for HGTP500.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching HGTP500
found in
> file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for UGRDP500.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching UGRDP500
found
> in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for VGRDP500.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching VGRDP500
found
> in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> WRFPRS_APCP03_d02.03.nc
> WARNING:
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 523850 observations from 39657 messages.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V.stat
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_fho.txt
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_ctc.txt
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cts.txt
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cnt.txt
> DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_mpr.txt
>
>
>
>
>
>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Mon May 04 20:40:32 2015


Hi Julie

When I set the fcst and obs fields like the sort of format you
suggested . It comes the correct results.

Thank you for your reply and answers.

I appreciate .

Have a good day!

Yuanpu

> -----原始邮件-----
> 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> 发送时间: 2015年5月4日 星期一
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #71731] Problem with point stat tool
>
> I'm sorry to hear you're having problems running point_stat.  There
appears
> to be a problem with the way your configuration file is set up.
>
> I can try to help some with the information you provided.  You can
also
> look at the default config file:
>
> /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
>
> to see how it is set up.
>
> Without seeing your configuration file, I can only refer to the
fields you
> include in your email.  You mentioned that you have:
>
> fcst_field[]  = [ "APCP_03(*,*)" ];
> obs_field[]   = [ "APCP/A3" ];
>
> This isn't the correct syntax for the point stat configuration.
Instead,
> you want to have something like this with the wind_thresh and
cat_thresh
> set however you'd like, but following this sort of format:
>
> fcst = {
>      wind_thresh  = [ NA ];
>
>      field = [
>         {
>           name       = "APCP_03";
>           level      = [ "(*,*)" ];
>           cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0,
>=30.0 ];
>         }
>      ];
> };
> obs = {
>      wind_thresh  = [ NA ];
>      message_type = ["ADPSFC"];
>
>      field = [
>         {
>           name       = "APCP";
>           level      = [ "A03" ];
>           cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0,
>=30.0 ];
>         }
>      ];
> };
>
>
> If setting the fields this way doesn't solve your problem, please
let me
> know and please include your configuration files as well.
>
> Thanks,
> Julie
>
> On Sun, May 3, 2015 at 9:33 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > Sun May 03 21:33:42 2015: Request 71731 was acted upon.
> > Transaction: Ticket created by lyp at lzb.ac.cn
> >        Queue: met_help
> >      Subject: Problem with point stat tool
> >        Owner: Nobody
> >   Requestors: lyp at lzb.ac.cn
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> >
> > Dear all,
> >
> > I received the following error message while running the
point_stat tool
> > of met 5.0.I am trying to use Point-Stat to compare the NetCDF
output of
> > PCP-
> > Combine to point observations generated using the ASCII2NC tool.
My
> > configuration file set were:fcst_field  = [ "APCP_03(*,*)" ];
> > obs_field   = [ "APCP/A3" ];
> >
> > when configuration fiel set
> > fcst_field[]  = [ "APCP_03(*,*)" ];
> > obs_field[]   = [ "APCP/A3" ];
> > the error was
> > ERROR  : yyerror() -> syntax error in file
"./PointStatConfig_3hrpre.d02"
> > ERROR  :
> > ERROR  :    line   = 20
> > ERROR  :
> > ERROR  :    column = 11
> > ERROR  :
> > ERROR  :    text   = "["
> > ERROR  :
> > ERROR  :
> > ERROR  : fcst_field[] = [ "APCP_03(*,*)" ];
> > ERROR  : __________^_______________________
> > ERROR  :
> >
> > Kindly advise.
> >
> > Thanks in advance.
> >
> > Yuanpu
> >
> >
> > DEBUG 1: Default Config File:
> > /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: ./PointStatConfig_3hrpre.d02
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=1
> > DEBUG 1: Forecast File:
> > /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > DEBUG 1: Climatology File: none
> > DEBUG 1: Observation File:
> > /data1/home/liuyp/software/met-
5.0/ascii2nc/output/SURF_obs_ascii_201504.nc
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for SPFHP500.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching SPFHP500
found
> > in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for TMPP500.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching TMPP500
found in
> > file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for HGTP500.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching HGTP500
found in
> > file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for UGRDP500.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching UGRDP500
found
> > in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for VGRDP500.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching VGRDP500
found
> > in file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > WRFPRS_APCP03_d02.03.nc
> > WARNING:
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Searching 523850 observations from 39657 messages.
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V.stat
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_fho.txt
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_ctc.txt
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cts.txt
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_cnt.txt
> > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V_mpr.txt
> >
> >
> >
> >
> >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>




------------------------------------------------
Subject: Problem with point stat tool
From: Julie Prestopnik
Time: Tue May 05 09:37:08 2015

Great!  I'm so glad to hear it is working for you now.

I hope you'll have a good day too!

Julie

On Mon, May 4, 2015 at 8:40 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
>
> Hi Julie
>
> When I set the fcst and obs fields like the sort of format you
suggested .
> It comes the correct results.
>
> Thank you for your reply and answers.
>
> I appreciate .
>
> Have a good day!
>
> Yuanpu
>
> > -----原始邮件-----
> > 发件人: "Julie Prestopnik via RT" <met_help at ucar.edu>
> > 发送时间: 2015年5月4日 星期一
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: [rt.rap.ucar.edu #71731] Problem with point stat tool
> >
> > I'm sorry to hear you're having problems running point_stat.
There
> appears
> > to be a problem with the way your configuration file is set up.
> >
> > I can try to help some with the information you provided.  You can
also
> > look at the default config file:
> >
> >
> /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
> >
> > to see how it is set up.
> >
> > Without seeing your configuration file, I can only refer to the
fields
> you
> > include in your email.  You mentioned that you have:
> >
> > fcst_field[]  = [ "APCP_03(*,*)" ];
> > obs_field[]   = [ "APCP/A3" ];
> >
> > This isn't the correct syntax for the point stat configuration.
Instead,
> > you want to have something like this with the wind_thresh and
cat_thresh
> > set however you'd like, but following this sort of format:
> >
> > fcst = {
> >      wind_thresh  = [ NA ];
> >
> >      field = [
> >         {
> >           name       = "APCP_03";
> >           level      = [ "(*,*)" ];
> >           cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0,
>=30.0 ];
> >         }
> >      ];
> > };
> > obs = {
> >      wind_thresh  = [ NA ];
> >      message_type = ["ADPSFC"];
> >
> >      field = [
> >         {
> >           name       = "APCP";
> >           level      = [ "A03" ];
> >           cat_thresh = [ >=5.0, >=10.0, >=15.0, >=20.0, >=25.0,
>=30.0 ];
> >         }
> >      ];
> > };
> >
> >
> > If setting the fields this way doesn't solve your problem, please
let me
> > know and please include your configuration files as well.
> >
> > Thanks,
> > Julie
> >
> > On Sun, May 3, 2015 at 9:33 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > Sun May 03 21:33:42 2015: Request 71731 was acted upon.
> > > Transaction: Ticket created by lyp at lzb.ac.cn
> > >        Queue: met_help
> > >      Subject: Problem with point stat tool
> > >        Owner: Nobody
> > >   Requestors: lyp at lzb.ac.cn
> > >       Status: new
> > >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> >
> > >
> > >
> > > Dear all,
> > >
> > > I received the following error message while running the
point_stat
> tool
> > > of met 5.0.I am trying to use Point-Stat to compare the NetCDF
output
> of
> > > PCP-
> > > Combine to point observations generated using the ASCII2NC tool.
My
> > > configuration file set were:fcst_field  = [ "APCP_03(*,*)" ];
> > > obs_field   = [ "APCP/A3" ];
> > >
> > > when configuration fiel set
> > > fcst_field[]  = [ "APCP_03(*,*)" ];
> > > obs_field[]   = [ "APCP/A3" ];
> > > the error was
> > > ERROR  : yyerror() -> syntax error in file
> "./PointStatConfig_3hrpre.d02"
> > > ERROR  :
> > > ERROR  :    line   = 20
> > > ERROR  :
> > > ERROR  :    column = 11
> > > ERROR  :
> > > ERROR  :    text   = "["
> > > ERROR  :
> > > ERROR  :
> > > ERROR  : fcst_field[] = [ "APCP_03(*,*)" ];
> > > ERROR  : __________^_______________________
> > > ERROR  :
> > >
> > > Kindly advise.
> > >
> > > Thanks in advance.
> > >
> > > Yuanpu
> > >
> > >
> > > DEBUG 1: Default Config File:
> > >
> /data1/home/liuyp/software/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: ./PointStatConfig_3hrpre.d02
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=1
> > > DEBUG 1: Forecast File:
> > > /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > DEBUG 1: Climatology File: none
> > > DEBUG 1: Observation File:
> > >
> /data1/home/liuyp/software/met-
5.0/ascii2nc/output/SURF_obs_ascii_201504.nc
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for SPFHP500.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
SPFHP500
> found
> > > in file:
> /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for TMPP500.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
TMPP500
> found in
> > > file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for HGTP500.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
HGTP500
> found in
> > > file: /data1/home/liuyp/software/met-
5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for UGRDP500.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
UGRDP500
> found
> > > in file:
> /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for VGRDP500.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
VGRDP500
> found
> > > in file:
> /data1/home/liuyp/software/met-5.0/pcp_combine/output/2015040500/
> > > WRFPRS_APCP03_d02.03.nc
> > > WARNING:
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Searching 523850 observations from 39657 messages.
> > > DEBUG 2:
> > > DEBUG 2:
> > >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 1: Output file: ./point_stat_wrf_-02-46-
39L_19700101_000000V.stat
> > > DEBUG 1: Output file:
> ./point_stat_wrf_-02-46-39L_19700101_000000V_fho.txt
> > > DEBUG 1: Output file:
> ./point_stat_wrf_-02-46-39L_19700101_000000V_ctc.txt
> > > DEBUG 1: Output file:
> ./point_stat_wrf_-02-46-39L_19700101_000000V_cts.txt
> > > DEBUG 1: Output file:
> ./point_stat_wrf_-02-46-39L_19700101_000000V_cnt.txt
> > > DEBUG 1: Output file:
> ./point_stat_wrf_-02-46-39L_19700101_000000V_mpr.txt
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
>
>
>
>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Tue Jun 02 19:08:57 2015



Hello met_help members,

I have some questions about mode output plots. First, what graph
software was used to plot the PostScript plots ? Second, which
directory was the plot scripts in, and what is it names. Because I
want to add China's province boundaries on the plots . The plots which
I had just had national boundaries now.

Kindly advise.

Thanks in advance.
Yuanpu




------------------------------------------------
Subject: Problem with point stat tool
From: Randy Bullock
Time: Wed Jun 03 11:12:00 2015

Hello Yuanpu -

Thanks for emailing met_help.

To answer your questions, the graphing software used to produce the
plots
is included in the MET package.  It's in the src/libcode/vx_ps
directory.
No external plotting programs or scripts are used.

The map data we use is in the data/map directory (and its
subdirectories).
If you could reformat your map data to be the same as those data
files, and
then edit the ConfigMapData file (in data/config) to use your new
files,
then you should be able to get China's provinces to show up on the
maps
that MET produces.

Let me know if you have any other questions or need any further help.

Randy

On Tue, Jun 2, 2015 at 7:08 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
>
>
> Hello met_help members,
>
> I have some questions about mode output plots. First, what graph
software
> was used to plot the PostScript plots ? Second, which directory was
the
> plot scripts in, and what is it names. Because I want to add China's
> province boundaries on the plots . The plots which I had just had
national
> boundaries now.
>
> Kindly advise.
>
> Thanks in advance.
> Yuanpu
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Thu Jun 04 02:11:07 2015

Hello Randy

Thanks for your reply.

I changed ConfigMapData (in data/config) file like below,but it also
has no China's provinces boundares. If it was necessary to recompile
the MET after change the ConfigMapData file. It also no changes in
plots when I change the line_width.

////////////////////////////////////////////////////////////////////////////////
//
// Default map data plotting information.
//
////////////////////////////////////////////////////////////////////////////////

map_data = {

   line_color = [ 25, 25, 25 ]; // rgb triple values, 0-255
   line_width = 0.5;
   line_dash  = "";

   source = [
      {
        file_name = "MET_BASE/map/legacy/world_data";
      },
      {
        file_name = "MET_BASE/map/legacy/country_data";
      },
      {
        file_name = "/data1/home/liuyp/local/met-
5.0/data/map/update_2012/admin_by_country/admin_China_data";
      }

   ];
};

////////////////////////////////////////////////////////////////////////






------------------------------------------------
Subject: Problem with point stat tool
From: Randy Bullock
Time: Thu Jun 04 14:13:18 2015

Hello again, Yuanpu -

Which MET tool are you using?  Is it MODE, or is it one of the other
tools?

Randy

On Thu, Jun 4, 2015 at 2:11 AM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
> Hello Randy
>
> Thanks for your reply.
>
> I changed ConfigMapData (in data/config) file like below,but it also
has
> no China's provinces boundares. If it was necessary to recompile the
MET
> after change the ConfigMapData file. It also no changes in plots
when I
> change the line_width.
>
>
>
////////////////////////////////////////////////////////////////////////////////
> //
> // Default map data plotting information.
> //
>
>
////////////////////////////////////////////////////////////////////////////////
>
> map_data = {
>
>    line_color = [ 25, 25, 25 ]; // rgb triple values, 0-255
>    line_width = 0.5;
>    line_dash  = "";
>
>    source = [
>       {
>         file_name = "MET_BASE/map/legacy/world_data";
>       },
>       {
>         file_name = "MET_BASE/map/legacy/country_data";
>       },
>       {
>         file_name =
> "/data1/home/liuyp/local/met-
5.0/data/map/update_2012/admin_by_country/admin_China_data";
>       }
>
>    ];
> };
>
>
////////////////////////////////////////////////////////////////////////
>
>
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Thu Jun 04 18:52:12 2015

Hello Randy

I have solved the  problems when  recompiled the MET(v5.0). By the
way, I used MODE tool.

Thanks for your reply.

I appreciate it .

Have a good day.

Yuanpu


> -----原始邮件-----
> 发件人: "Randy Bullock via RT" <met_help at ucar.edu>
> 发送时间: 2015年6月5日 星期五
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: [rt.rap.ucar.edu #71731] About mode PostScript plots
>
> Hello again, Yuanpu -
>
> Which MET tool are you using?  Is it MODE, or is it one of the other
tools?
>
> Randy
>
> On Thu, Jun 4, 2015 at 2:11 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> > Hello Randy
> >
> > Thanks for your reply.
> >
> > I changed ConfigMapData (in data/config) file like below,but it
also has
> > no China's provinces boundares. If it was necessary to recompile
the MET
> > after change the ConfigMapData file. It also no changes in plots
when I
> > change the line_width.
> >
> >
> >
////////////////////////////////////////////////////////////////////////////////
> > //
> > // Default map data plotting information.
> > //
> >
> >
////////////////////////////////////////////////////////////////////////////////
> >
> > map_data = {
> >
> >    line_color = [ 25, 25, 25 ]; // rgb triple values, 0-255
> >    line_width = 0.5;
> >    line_dash  = "";
> >
> >    source = [
> >       {
> >         file_name = "MET_BASE/map/legacy/world_data";
> >       },
> >       {
> >         file_name = "MET_BASE/map/legacy/country_data";
> >       },
> >       {
> >         file_name =
> > "/data1/home/liuyp/local/met-
5.0/data/map/update_2012/admin_by_country/admin_China_data";
> >       }
> >
> >    ];
> > };
> >
> >
////////////////////////////////////////////////////////////////////////
> >
> >
> >
> >
> >
> >
> >
>




------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Sun Aug 16 20:51:14 2015


Dear all,

I want to verify ECMWF data with point_stat tool of  MET 5.0. When I
used the original data(grib format), it has the followed error
information:

DEBUG 1: Default Config File: /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
ERROR  :
ERROR  : GribFile::read_record(GribRecord &) -> Second order packing
not implemented.
ERROR  :

 But when I  thransfered  the original data to nc format with NCL, it
has the followed error information:

DEBUG 1: Default Config File: /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
NetCDF: Attribute not found

The original ECMWF nc format data in attachment.

Kindly advise.

Thanks in advance.

Yuanpu








------------------------------------------------
Subject: Problem with point stat tool
From: John Halley Gotway
Time: Mon Aug 17 12:08:57 2015

Yuanpu,

Unfortunately, the NetCDF output of the ncl_convert2nc tool is not
CF-compliant so MET will not be able to read it.  Can you please send
me
the original ECMWF GRIB1 file that produced the error you're seeing.
It
doesn't need to be the full file, just a single record would suffice.

You can post it to our anonymous ftp site following the instructions
listed
here:
   http://www.dtcenter.org/met/users/support/met_help.php#ftp

Please write back once you've posted the data so that I can go grab
it.

Thanks,
John Halley Gotway
met_help at ucar.edu

On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
>
> Dear all,
>
> I want to verify ECMWF data with point_stat tool of  MET 5.0. When I
used
> the original data(grib format), it has the followed error
information:
>
> DEBUG 1: Default Config File:
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> ERROR  :
> ERROR  : GribFile::read_record(GribRecord &) -> Second order packing
not
> implemented.
> ERROR  :
>
>  But when I  thransfered  the original data to nc format with NCL,
it has
> the followed error information:
>
> DEBUG 1: Default Config File:
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> NetCDF: Attribute not found
>
> The original ECMWF nc format data in attachment.
>
> Kindly advise.
>
> Thanks in advance.
>
> Yuanpu
>
>
>
>
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Tue Aug 18 03:33:16 2015

Hello John

Thank you for you reply. I have put the ECMWF
data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-ACHN.bin) on
your  anonymous ftp site, it was in incoming/irap/met_help/liu_data
dirctory.but the observed data was not thransfered completed. So I put
it in attachment.



yuanpu




> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年8月18日 星期二
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
>
> Yuanpu,
>
> Unfortunately, the NetCDF output of the ncl_convert2nc tool is not
> CF-compliant so MET will not be able to read it.  Can you please
send me
> the original ECMWF GRIB1 file that produced the error you're seeing.
It
> doesn't need to be the full file, just a single record would
suffice.
>
> You can post it to our anonymous ftp site following the instructions
listed
> here:
>    http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please write back once you've posted the data so that I can go grab
it.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> >
> > Dear all,
> >
> > I want to verify ECMWF data with point_stat tool of  MET 5.0. When
I used
> > the original data(grib format), it has the followed error
information:
> >
> > DEBUG 1: Default Config File:
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > ERROR  :
> > ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing not
> > implemented.
> > ERROR  :
> >
> >  But when I  thransfered  the original data to nc format with NCL,
it has
> > the followed error information:
> >
> > DEBUG 1: Default Config File:
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > NetCDF: Attribute not found
> >
> > The original ECMWF nc format data in attachment.
> >
> > Kindly advise.
> >
> > Thanks in advance.
> >
> > Yuanpu
> >
> >
> >
> >
> >
> >
> >
> >
> >
>


从lzb.ac.cn发来的超大附件
Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
查看下载信息:https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e




------------------------------------------------
Subject: Problem with point stat tool
From: John Halley Gotway
Time: Wed Aug 19 14:31:31 2015

Yuanpu,

Thanks for sending that data.  When I run it through MET, I see the
same
error you're getting:
   ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing not
implemented.

Running it through wgrib to inventory the file, I get a similar error:
   *** Cannot decode complex packed fields n=20089***

And trying to plot the data with the Integrated Data Viewer (IDV), I
also
get a similar error:
  Couldn't get data
  Grib1BinaryDataSection: (octet 4, 1st half) not grid point data and
simple packing

Take a look at GRIB1 table 11:
   http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html

Having the second bit of BDS octet 4 set to 1 indicates that the data
uses
complex packing, which is generally not supported.

I also tried viewing the NetCDF data file you sent me, the output of
ncl_convert2nc, using the ncview utility.  While ncview is able to
open the
file, all of the data values plotted are 0.

Honestly, this is the first time I've seen this error.  My guess is
that
your GRIB1 data file is corrupted or badly created in some way.  I'd
take a
close look at how that GRIB1 file is being created.  Once you're able
to
view the data using the IDV tool, you should be able to use it in MET.

Thanks,
John

On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
> Hello John
>
> Thank you for you reply. I have put the ECMWF
> data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-ACHN.bin)
on your
> anonymous ftp site, it was in incoming/irap/met_help/liu_data
dirctory.but
> the observed data was not thransfered completed. So I put it in
attachment.
>
>
>
> yuanpu
>
>
>
>
> > -----原始邮件-----
> > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > 发送时间: 2015年8月18日 星期二
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> >
> > Yuanpu,
> >
> > Unfortunately, the NetCDF output of the ncl_convert2nc tool is not
> > CF-compliant so MET will not be able to read it.  Can you please
send me
> > the original ECMWF GRIB1 file that produced the error you're
seeing.  It
> > doesn't need to be the full file, just a single record would
suffice.
> >
> > You can post it to our anonymous ftp site following the
instructions
> listed
> > here:
> >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Please write back once you've posted the data so that I can go
grab it.
> >
> > Thanks,
> > John Halley Gotway
> > met_help at ucar.edu
> >
> > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > >
> > >
> > > Dear all,
> > >
> > > I want to verify ECMWF data with point_stat tool of  MET 5.0.
When I
> used
> > > the original data(grib format), it has the followed error
information:
> > >
> > > DEBUG 1: Default Config File:
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > ERROR  :
> > > ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing
> not
> > > implemented.
> > > ERROR  :
> > >
> > >  But when I  thransfered  the original data to nc format with
NCL, it
> has
> > > the followed error information:
> > >
> > > DEBUG 1: Default Config File:
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > NetCDF: Attribute not found
> > >
> > > The original ECMWF nc format data in attachment.
> > >
> > > Kindly advise.
> > >
> > > Thanks in advance.
> > >
> > > Yuanpu
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> >
>
>
> 从lzb.ac.cn发来的超大附件
> Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> 查看下载信息:
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Thu Aug 20 20:51:15 2015

Hello John,

Thanks for your reply. The original ECMWF data were complex packing,
so I decoded them with grib tool. But I see the new error :

DEBUG 1: Default Config File: /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
ERROR  :
ERROR  : GribFile::read_record(GribRecord &) -> found a grib record
larger than the buffer size.
ERROR  :
ERROR  :   Increase the buffer to at least 8819010 bytes.
ERROR  :
ERROR  :

I put the grib1 data in attachment.

yuanpu

> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年8月20日 星期四
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
>
> Yuanpu,
>
> Thanks for sending that data.  When I run it through MET, I see the
same
> error you're getting:
>    ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing not
> implemented.
>
> Running it through wgrib to inventory the file, I get a similar
error:
>    *** Cannot decode complex packed fields n=20089***
>
> And trying to plot the data with the Integrated Data Viewer (IDV), I
also
> get a similar error:
>   Couldn't get data
>   Grib1BinaryDataSection: (octet 4, 1st half) not grid point data
and
> simple packing
>
> Take a look at GRIB1 table 11:
>    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
>
> Having the second bit of BDS octet 4 set to 1 indicates that the
data uses
> complex packing, which is generally not supported.
>
> I also tried viewing the NetCDF data file you sent me, the output of
> ncl_convert2nc, using the ncview utility.  While ncview is able to
open the
> file, all of the data values plotted are 0.
>
> Honestly, this is the first time I've seen this error.  My guess is
that
> your GRIB1 data file is corrupted or badly created in some way.  I'd
take a
> close look at how that GRIB1 file is being created.  Once you're
able to
> view the data using the IDV tool, you should be able to use it in
MET.
>
> Thanks,
> John
>
> On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> > Hello John
> >
> > Thank you for you reply. I have put the ECMWF
> > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-ACHN.bin)
on your
> > anonymous ftp site, it was in incoming/irap/met_help/liu_data
dirctory.but
> > the observed data was not thransfered completed. So I put it in
attachment.
> >
> >
> >
> > yuanpu
> >
> >
> >
> >
> > > -----原始邮件-----
> > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > 发送时间: 2015年8月18日 星期二
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > >
> > > Yuanpu,
> > >
> > > Unfortunately, the NetCDF output of the ncl_convert2nc tool is
not
> > > CF-compliant so MET will not be able to read it.  Can you please
send me
> > > the original ECMWF GRIB1 file that produced the error you're
seeing.  It
> > > doesn't need to be the full file, just a single record would
suffice.
> > >
> > > You can post it to our anonymous ftp site following the
instructions
> > listed
> > > here:
> > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Please write back once you've posted the data so that I can go
grab it.
> > >
> > > Thanks,
> > > John Halley Gotway
> > > met_help at ucar.edu
> > >
> > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
>
> > > >
> > > >
> > > > Dear all,
> > > >
> > > > I want to verify ECMWF data with point_stat tool of  MET 5.0.
When I
> > used
> > > > the original data(grib format), it has the followed error
information:
> > > >
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > ERROR  :
> > > > ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing
> > not
> > > > implemented.
> > > > ERROR  :
> > > >
> > > >  But when I  thransfered  the original data to nc format with
NCL, it
> > has
> > > > the followed error information:
> > > >
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > NetCDF: Attribute not found
> > > >
> > > > The original ECMWF nc format data in attachment.
> > > >
> > > > Kindly advise.
> > > >
> > > > Thanks in advance.
> > > >
> > > > Yuanpu
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > >
> >
> >
> > 从lzb.ac.cn发来的超大附件
> > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > 查看下载信息:
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> >
> >
> >
> >
> >
>


从lzb.ac.cn发来的超大附件
sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
查看下载信息:https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk




------------------------------------------------
Subject: Problem with point stat tool
From: John Halley Gotway
Time: Fri Aug 21 10:11:03 2015

Yuanpu,

The attachment didn't come through.  Please post it to our anonymous
ftp
site following these instructions:
   http://www.dtcenter.org/met/users/support/met_help.php#ftp

Yes, that error you're seeing is a pesky little issue.  We do have a
development task defined to dynamically allocate the amount of space
required rather than having a fixed upper limit.

Please try editing this source code file:
   src/libcode/vx_data2d_grib/grib_classes.h

Find and set this line as follows:
   static const int gribfile_buf_size = 16777216;   //  2^24

And then from the top-level MET directory, do a "make clean" followed
by a
"make install".

Please let me know how it goes after making that change.

Thanks,
John

On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
> Hello John,
>
> Thanks for your reply. The original ECMWF data were complex packing,
so I
> decoded them with grib tool. But I see the new error :
>
> DEBUG 1: Default Config File:
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> ERROR  :
> ERROR  : GribFile::read_record(GribRecord &) -> found a grib record
larger
> than the buffer size.
> ERROR  :
> ERROR  :   Increase the buffer to at least 8819010 bytes.
> ERROR  :
> ERROR  :
>
> I put the grib1 data in attachment.
>
> yuanpu
>
> > -----原始邮件-----
> > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > 发送时间: 2015年8月20日 星期四
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> >
> > Yuanpu,
> >
> > Thanks for sending that data.  When I run it through MET, I see
the same
> > error you're getting:
> >    ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing
> not
> > implemented.
> >
> > Running it through wgrib to inventory the file, I get a similar
error:
> >    *** Cannot decode complex packed fields n=20089***
> >
> > And trying to plot the data with the Integrated Data Viewer (IDV),
I also
> > get a similar error:
> >   Couldn't get data
> >   Grib1BinaryDataSection: (octet 4, 1st half) not grid point data
and
> > simple packing
> >
> > Take a look at GRIB1 table 11:
> >    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> >
> > Having the second bit of BDS octet 4 set to 1 indicates that the
data
> uses
> > complex packing, which is generally not supported.
> >
> > I also tried viewing the NetCDF data file you sent me, the output
of
> > ncl_convert2nc, using the ncview utility.  While ncview is able to
open
> the
> > file, all of the data values plotted are 0.
> >
> > Honestly, this is the first time I've seen this error.  My guess
is that
> > your GRIB1 data file is corrupted or badly created in some way.
I'd
> take a
> > close look at how that GRIB1 file is being created.  Once you're
able to
> > view the data using the IDV tool, you should be able to use it in
MET.
> >
> > Thanks,
> > John
> >
> > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > >
> > > Hello John
> > >
> > > Thank you for you reply. I have put the ECMWF
> > > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-
ACHN.bin) on
> your
> > > anonymous ftp site, it was in incoming/irap/met_help/liu_data
> dirctory.but
> > > the observed data was not thransfered completed. So I put it in
> attachment.
> > >
> > >
> > >
> > > yuanpu
> > >
> > >
> > >
> > >
> > > > -----原始邮件-----
> > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > 发送时间: 2015年8月18日 星期二
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > > >
> > > > Yuanpu,
> > > >
> > > > Unfortunately, the NetCDF output of the ncl_convert2nc tool is
not
> > > > CF-compliant so MET will not be able to read it.  Can you
please
> send me
> > > > the original ECMWF GRIB1 file that produced the error you're
> seeing.  It
> > > > doesn't need to be the full file, just a single record would
suffice.
> > > >
> > > > You can post it to our anonymous ftp site following the
instructions
> > > listed
> > > > here:
> > > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Please write back once you've posted the data so that I can go
grab
> it.
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > > met_help at ucar.edu
> > > >
> > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > >
> > > > >
> > > > > Dear all,
> > > > >
> > > > > I want to verify ECMWF data with point_stat tool of  MET
5.0. When
> I
> > > used
> > > > > the original data(grib format), it has the followed error
> information:
> > > > >
> > > > > DEBUG 1: Default Config File:
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > ERROR  :
> > > > > ERROR  : GribFile::read_record(GribRecord &) -> Second order
> packing
> > > not
> > > > > implemented.
> > > > > ERROR  :
> > > > >
> > > > >  But when I  thransfered  the original data to nc format
with NCL,
> it
> > > has
> > > > > the followed error information:
> > > > >
> > > > > DEBUG 1: Default Config File:
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > NetCDF: Attribute not found
> > > > >
> > > > > The original ECMWF nc format data in attachment.
> > > > >
> > > > > Kindly advise.
> > > > >
> > > > > Thanks in advance.
> > > > >
> > > > > Yuanpu
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > >
> > >
> > >
> > > 从lzb.ac.cn发来的超大附件
> > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > 查看下载信息:
> > >
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > >
> > >
> > >
> > >
> > >
> >
>
>
> 从lzb.ac.cn发来的超大附件
> sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> 查看下载信息:
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Sun Aug 23 22:13:27 2015

Hello John,

Unfortunately , when I set lines as     static const int
gribfile_buf_size = 16777216;   //  2^24

and mke "make clean" and "make install". But I see the same error:

DEBUG 1: Default Config File: /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
ERROR  :
ERROR  : GribFile::read_record(GribRecord &) -> found a grib record
larger than the buffer size.
ERROR  :
ERROR  :   Increase the buffer to at least 8819010 bytes.
ERROR  :
ERROR  :

I has post the girb1 data in liu_data directory of your ftp site.

yuanpu



> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年8月22日 星期六
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
>
> Yuanpu,
>
> The attachment didn't come through.  Please post it to our anonymous
ftp
> site following these instructions:
>    http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Yes, that error you're seeing is a pesky little issue.  We do have a
> development task defined to dynamically allocate the amount of space
> required rather than having a fixed upper limit.
>
> Please try editing this source code file:
>    src/libcode/vx_data2d_grib/grib_classes.h
>
> Find and set this line as follows:
>    static const int gribfile_buf_size = 16777216;   //  2^24
>
> And then from the top-level MET directory, do a "make clean"
followed by a
> "make install".
>
> Please let me know how it goes after making that change.
>
> Thanks,
> John
>
> On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> > Hello John,
> >
> > Thanks for your reply. The original ECMWF data were complex
packing, so I
> > decoded them with grib tool. But I see the new error :
> >
> > DEBUG 1: Default Config File:
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > ERROR  :
> > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record larger
> > than the buffer size.
> > ERROR  :
> > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > ERROR  :
> > ERROR  :
> >
> > I put the grib1 data in attachment.
> >
> > yuanpu
> >
> > > -----原始邮件-----
> > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > 发送时间: 2015年8月20日 星期四
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > >
> > > Yuanpu,
> > >
> > > Thanks for sending that data.  When I run it through MET, I see
the same
> > > error you're getting:
> > >    ERROR  : GribFile::read_record(GribRecord &) -> Second order
packing
> > not
> > > implemented.
> > >
> > > Running it through wgrib to inventory the file, I get a similar
error:
> > >    *** Cannot decode complex packed fields n=20089***
> > >
> > > And trying to plot the data with the Integrated Data Viewer
(IDV), I also
> > > get a similar error:
> > >   Couldn't get data
> > >   Grib1BinaryDataSection: (octet 4, 1st half) not grid point
data and
> > > simple packing
> > >
> > > Take a look at GRIB1 table 11:
> > >    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> > >
> > > Having the second bit of BDS octet 4 set to 1 indicates that the
data
> > uses
> > > complex packing, which is generally not supported.
> > >
> > > I also tried viewing the NetCDF data file you sent me, the
output of
> > > ncl_convert2nc, using the ncview utility.  While ncview is able
to open
> > the
> > > file, all of the data values plotted are 0.
> > >
> > > Honestly, this is the first time I've seen this error.  My guess
is that
> > > your GRIB1 data file is corrupted or badly created in some way.
I'd
> > take a
> > > close look at how that GRIB1 file is being created.  Once you're
able to
> > > view the data using the IDV tool, you should be able to use it
in MET.
> > >
> > > Thanks,
> > > John
> > >
> > > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
>
> > > >
> > > > Hello John
> > > >
> > > > Thank you for you reply. I have put the ECMWF
> > > > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-
ACHN.bin) on
> > your
> > > > anonymous ftp site, it was in incoming/irap/met_help/liu_data
> > dirctory.but
> > > > the observed data was not thransfered completed. So I put it
in
> > attachment.
> > > >
> > > >
> > > >
> > > > yuanpu
> > > >
> > > >
> > > >
> > > >
> > > > > -----原始邮件-----
> > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > 发送时间: 2015年8月18日 星期二
> > > > > 收件人: lyp at lzb.ac.cn
> > > > > 抄送:
> > > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > > > >
> > > > > Yuanpu,
> > > > >
> > > > > Unfortunately, the NetCDF output of the ncl_convert2nc tool
is not
> > > > > CF-compliant so MET will not be able to read it.  Can you
please
> > send me
> > > > > the original ECMWF GRIB1 file that produced the error you're
> > seeing.  It
> > > > > doesn't need to be the full file, just a single record would
suffice.
> > > > >
> > > > > You can post it to our anonymous ftp site following the
instructions
> > > > listed
> > > > > here:
> > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > >
> > > > > Please write back once you've posted the data so that I can
go grab
> > it.
> > > > >
> > > > > Thanks,
> > > > > John Halley Gotway
> > > > > met_help at ucar.edu
> > > > >
> > > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > > >
> > > > > >
> > > > > > Dear all,
> > > > > >
> > > > > > I want to verify ECMWF data with point_stat tool of  MET
5.0. When
> > I
> > > > used
> > > > > > the original data(grib format), it has the followed error
> > information:
> > > > > >
> > > > > > DEBUG 1: Default Config File:
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > ERROR  :
> > > > > > ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> > packing
> > > > not
> > > > > > implemented.
> > > > > > ERROR  :
> > > > > >
> > > > > >  But when I  thransfered  the original data to nc format
with NCL,
> > it
> > > > has
> > > > > > the followed error information:
> > > > > >
> > > > > > DEBUG 1: Default Config File:
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > NetCDF: Attribute not found
> > > > > >
> > > > > > The original ECMWF nc format data in attachment.
> > > > > >
> > > > > > Kindly advise.
> > > > > >
> > > > > > Thanks in advance.
> > > > > >
> > > > > > Yuanpu
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > >
> > > >
> > > > 从lzb.ac.cn发来的超大附件
> > > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > > 查看下载信息:
> > > >
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > > >
> > > >
> > > >
> > > >
> > > >
> > >
> >
> >
> > 从lzb.ac.cn发来的超大附件
> > sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> > 查看下载信息:
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
> >
> >
> >
> >
> >
>




------------------------------------------------
Subject: Problem with point stat tool
From: John Halley Gotway
Time: Mon Aug 24 10:13:20 2015

Yuanpu,

Unfortunately, it looks like the GRIB file you posted to our ftp site
is
corrupted.  I get an error when trying to read it.  Please look at the
size
of the file you tried to post:

   -rw-rw-r--   1 johnhg rap 1.5M Aug 24 10:07
sfc_20150616_36_grid_simple.grib1

Is it supposed to be larger than 1.5M?

Thanks,
John


On Sun, Aug 23, 2015 at 10:13 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
> Hello John,
>
> Unfortunately , when I set lines as     static const int
gribfile_buf_size
> = 16777216;   //  2^24
>
> and mke "make clean" and "make install". But I see the same error:
>
> DEBUG 1: Default Config File:
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> ERROR  :
> ERROR  : GribFile::read_record(GribRecord &) -> found a grib record
larger
> than the buffer size.
> ERROR  :
> ERROR  :   Increase the buffer to at least 8819010 bytes.
> ERROR  :
> ERROR  :
>
> I has post the girb1 data in liu_data directory of your ftp site.
>
> yuanpu
>
>
>
> > -----原始邮件-----
> > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > 发送时间: 2015年8月22日 星期六
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> >
> > Yuanpu,
> >
> > The attachment didn't come through.  Please post it to our
anonymous ftp
> > site following these instructions:
> >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Yes, that error you're seeing is a pesky little issue.  We do have
a
> > development task defined to dynamically allocate the amount of
space
> > required rather than having a fixed upper limit.
> >
> > Please try editing this source code file:
> >    src/libcode/vx_data2d_grib/grib_classes.h
> >
> > Find and set this line as follows:
> >    static const int gribfile_buf_size = 16777216;   //  2^24
> >
> > And then from the top-level MET directory, do a "make clean"
followed by
> a
> > "make install".
> >
> > Please let me know how it goes after making that change.
> >
> > Thanks,
> > John
> >
> > On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > >
> > > Hello John,
> > >
> > > Thanks for your reply. The original ECMWF data were complex
packing,
> so I
> > > decoded them with grib tool. But I see the new error :
> > >
> > > DEBUG 1: Default Config File:
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > ERROR  :
> > > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record
> larger
> > > than the buffer size.
> > > ERROR  :
> > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > ERROR  :
> > > ERROR  :
> > >
> > > I put the grib1 data in attachment.
> > >
> > > yuanpu
> > >
> > > > -----原始邮件-----
> > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > 发送时间: 2015年8月20日 星期四
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > > >
> > > > Yuanpu,
> > > >
> > > > Thanks for sending that data.  When I run it through MET, I
see the
> same
> > > > error you're getting:
> > > >    ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> packing
> > > not
> > > > implemented.
> > > >
> > > > Running it through wgrib to inventory the file, I get a
similar
> error:
> > > >    *** Cannot decode complex packed fields n=20089***
> > > >
> > > > And trying to plot the data with the Integrated Data Viewer
(IDV), I
> also
> > > > get a similar error:
> > > >   Couldn't get data
> > > >   Grib1BinaryDataSection: (octet 4, 1st half) not grid point
data and
> > > > simple packing
> > > >
> > > > Take a look at GRIB1 table 11:
> > > >    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> > > >
> > > > Having the second bit of BDS octet 4 set to 1 indicates that
the data
> > > uses
> > > > complex packing, which is generally not supported.
> > > >
> > > > I also tried viewing the NetCDF data file you sent me, the
output of
> > > > ncl_convert2nc, using the ncview utility.  While ncview is
able to
> open
> > > the
> > > > file, all of the data values plotted are 0.
> > > >
> > > > Honestly, this is the first time I've seen this error.  My
guess is
> that
> > > > your GRIB1 data file is corrupted or badly created in some
way.  I'd
> > > take a
> > > > close look at how that GRIB1 file is being created.  Once
you're
> able to
> > > > view the data using the IDV tool, you should be able to use it
in
> MET.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > >
> > > > > Hello John
> > > > >
> > > > > Thank you for you reply. I have put the ECMWF
> > > > > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-
ACHN.bin)
> on
> > > your
> > > > > anonymous ftp site, it was in
incoming/irap/met_help/liu_data
> > > dirctory.but
> > > > > the observed data was not thransfered completed. So I put it
in
> > > attachment.
> > > > >
> > > > >
> > > > >
> > > > > yuanpu
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > > -----原始邮件-----
> > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > 发送时间: 2015年8月18日 星期二
> > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > 抄送:
> > > > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF Data
> > > > > >
> > > > > > Yuanpu,
> > > > > >
> > > > > > Unfortunately, the NetCDF output of the ncl_convert2nc
tool is
> not
> > > > > > CF-compliant so MET will not be able to read it.  Can you
please
> > > send me
> > > > > > the original ECMWF GRIB1 file that produced the error
you're
> > > seeing.  It
> > > > > > doesn't need to be the full file, just a single record
would
> suffice.
> > > > > >
> > > > > > You can post it to our anonymous ftp site following the
> instructions
> > > > > listed
> > > > > > here:
> > > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > >
> > > > > > Please write back once you've posted the data so that I
can go
> grab
> > > it.
> > > > > >
> > > > > > Thanks,
> > > > > > John Halley Gotway
> > > > > > met_help at ucar.edu
> > > > > >
> > > > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> > > wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> >
> > > > > > >
> > > > > > >
> > > > > > > Dear all,
> > > > > > >
> > > > > > > I want to verify ECMWF data with point_stat tool of  MET
5.0.
> When
> > > I
> > > > > used
> > > > > > > the original data(grib format), it has the followed
error
> > > information:
> > > > > > >
> > > > > > > DEBUG 1: Default Config File:
> > > > > > >
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > ERROR  :
> > > > > > > ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> > > packing
> > > > > not
> > > > > > > implemented.
> > > > > > > ERROR  :
> > > > > > >
> > > > > > >  But when I  thransfered  the original data to nc format
with
> NCL,
> > > it
> > > > > has
> > > > > > > the followed error information:
> > > > > > >
> > > > > > > DEBUG 1: Default Config File:
> > > > > > >
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > NetCDF: Attribute not found
> > > > > > >
> > > > > > > The original ECMWF nc format data in attachment.
> > > > > > >
> > > > > > > Kindly advise.
> > > > > > >
> > > > > > > Thanks in advance.
> > > > > > >
> > > > > > > Yuanpu
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > 从lzb.ac.cn发来的超大附件
> > > > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > > > 查看下载信息:
> > > > >
> > >
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > >
> > >
> > >
> > > 从lzb.ac.cn发来的超大附件
> > > sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> > > 查看下载信息:
> > >
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
> > >
> > >
> > >
> > >
> > >
> >
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Mon Aug 24 21:12:48 2015


Hello John

I has repost the grib data in you ftp site. It was thransfered
completed this time.

yuanpu.

> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年8月25日 星期二
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
>
> Yuanpu,
>
> Unfortunately, it looks like the GRIB file you posted to our ftp
site is
> corrupted.  I get an error when trying to read it.  Please look at
the size
> of the file you tried to post:
>
>    -rw-rw-r--   1 johnhg rap 1.5M Aug 24 10:07
> sfc_20150616_36_grid_simple.grib1
>
> Is it supposed to be larger than 1.5M?
>
> Thanks,
> John
>
>
> On Sun, Aug 23, 2015 at 10:13 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> > Hello John,
> >
> > Unfortunately , when I set lines as     static const int
gribfile_buf_size
> > = 16777216;   //  2^24
> >
> > and mke "make clean" and "make install". But I see the same error:
> >
> > DEBUG 1: Default Config File:
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > ERROR  :
> > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record larger
> > than the buffer size.
> > ERROR  :
> > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > ERROR  :
> > ERROR  :
> >
> > I has post the girb1 data in liu_data directory of your ftp site.
> >
> > yuanpu
> >
> >
> >
> > > -----原始邮件-----
> > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > 发送时间: 2015年8月22日 星期六
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > >
> > > Yuanpu,
> > >
> > > The attachment didn't come through.  Please post it to our
anonymous ftp
> > > site following these instructions:
> > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Yes, that error you're seeing is a pesky little issue.  We do
have a
> > > development task defined to dynamically allocate the amount of
space
> > > required rather than having a fixed upper limit.
> > >
> > > Please try editing this source code file:
> > >    src/libcode/vx_data2d_grib/grib_classes.h
> > >
> > > Find and set this line as follows:
> > >    static const int gribfile_buf_size = 16777216;   //  2^24
> > >
> > > And then from the top-level MET directory, do a "make clean"
followed by
> > a
> > > "make install".
> > >
> > > Please let me know how it goes after making that change.
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
>
> > > >
> > > > Hello John,
> > > >
> > > > Thanks for your reply. The original ECMWF data were complex
packing,
> > so I
> > > > decoded them with grib tool. But I see the new error :
> > > >
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > ERROR  :
> > > > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record
> > larger
> > > > than the buffer size.
> > > > ERROR  :
> > > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > > ERROR  :
> > > > ERROR  :
> > > >
> > > > I put the grib1 data in attachment.
> > > >
> > > > yuanpu
> > > >
> > > > > -----原始邮件-----
> > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > 发送时间: 2015年8月20日 星期四
> > > > > 收件人: lyp at lzb.ac.cn
> > > > > 抄送:
> > > > > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > > > >
> > > > > Yuanpu,
> > > > >
> > > > > Thanks for sending that data.  When I run it through MET, I
see the
> > same
> > > > > error you're getting:
> > > > >    ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> > packing
> > > > not
> > > > > implemented.
> > > > >
> > > > > Running it through wgrib to inventory the file, I get a
similar
> > error:
> > > > >    *** Cannot decode complex packed fields n=20089***
> > > > >
> > > > > And trying to plot the data with the Integrated Data Viewer
(IDV), I
> > also
> > > > > get a similar error:
> > > > >   Couldn't get data
> > > > >   Grib1BinaryDataSection: (octet 4, 1st half) not grid point
data and
> > > > > simple packing
> > > > >
> > > > > Take a look at GRIB1 table 11:
> > > > >    http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> > > > >
> > > > > Having the second bit of BDS octet 4 set to 1 indicates that
the data
> > > > uses
> > > > > complex packing, which is generally not supported.
> > > > >
> > > > > I also tried viewing the NetCDF data file you sent me, the
output of
> > > > > ncl_convert2nc, using the ncview utility.  While ncview is
able to
> > open
> > > > the
> > > > > file, all of the data values plotted are 0.
> > > > >
> > > > > Honestly, this is the first time I've seen this error.  My
guess is
> > that
> > > > > your GRIB1 data file is corrupted or badly created in some
way.  I'd
> > > > take a
> > > > > close look at how that GRIB1 file is being created.  Once
you're
> > able to
> > > > > view the data using the IDV tool, you should be able to use
it in
> > MET.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT
<met_help at ucar.edu>
> > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > > >
> > > > > > Hello John
> > > > > >
> > > > > > Thank you for you reply. I have put the ECMWF
> > > > > > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-
ACHN.bin)
> > on
> > > > your
> > > > > > anonymous ftp site, it was in
incoming/irap/met_help/liu_data
> > > > dirctory.but
> > > > > > the observed data was not thransfered completed. So I put
it in
> > > > attachment.
> > > > > >
> > > > > >
> > > > > >
> > > > > > yuanpu
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > > -----原始邮件-----
> > > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > > 发送时间: 2015年8月18日 星期二
> > > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > > 抄送:
> > > > > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > > > > > >
> > > > > > > Yuanpu,
> > > > > > >
> > > > > > > Unfortunately, the NetCDF output of the ncl_convert2nc
tool is
> > not
> > > > > > > CF-compliant so MET will not be able to read it.  Can
you please
> > > > send me
> > > > > > > the original ECMWF GRIB1 file that produced the error
you're
> > > > seeing.  It
> > > > > > > doesn't need to be the full file, just a single record
would
> > suffice.
> > > > > > >
> > > > > > > You can post it to our anonymous ftp site following the
> > instructions
> > > > > > listed
> > > > > > > here:
> > > > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > > >
> > > > > > > Please write back once you've posted the data so that I
can go
> > grab
> > > > it.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > John Halley Gotway
> > > > > > > met_help at ucar.edu
> > > > > > >
> > > > > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Dear all,
> > > > > > > >
> > > > > > > > I want to verify ECMWF data with point_stat tool of
MET 5.0.
> > When
> > > > I
> > > > > > used
> > > > > > > > the original data(grib format), it has the followed
error
> > > > information:
> > > > > > > >
> > > > > > > > DEBUG 1: Default Config File:
> > > > > > > >
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > > ERROR  :
> > > > > > > > ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> > > > packing
> > > > > > not
> > > > > > > > implemented.
> > > > > > > > ERROR  :
> > > > > > > >
> > > > > > > >  But when I  thransfered  the original data to nc
format with
> > NCL,
> > > > it
> > > > > > has
> > > > > > > > the followed error information:
> > > > > > > >
> > > > > > > > DEBUG 1: Default Config File:
> > > > > > > >
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > > NetCDF: Attribute not found
> > > > > > > >
> > > > > > > > The original ECMWF nc format data in attachment.
> > > > > > > >
> > > > > > > > Kindly advise.
> > > > > > > >
> > > > > > > > Thanks in advance.
> > > > > > > >
> > > > > > > > Yuanpu
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > 从lzb.ac.cn发来的超大附件
> > > > > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > > > > 查看下载信息:
> > > > > >
> > > >
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > >
> > > >
> > > > 从lzb.ac.cn发来的超大附件
> > > > sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> > > > 查看下载信息:
> > > >
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
> > > >
> > > >
> > > >
> > > >
> > > >
> > >
> >
> >
> >
> >
> >
>




------------------------------------------------
Subject: Problem with point stat tool
From: John Halley Gotway
Time: Tue Aug 25 10:30:55 2015

Yuanpu,

Thanks for sending that data file.  I was able to grab it and work
with it.

Making that change I indicated earlier did fix the problem.  Please
try the
following:

(1) Save the attached version of "grib_classes.h" to your machine.
And
place it into the met-5.0 directory:
  mv grib_classes.h met-5.0/src/libcode/vx_data2d_grib/grib_classes.h

(2) Rebuild met-5.0 (be sure to run "make clean"):
  cd met-5.0
  make clean
  make install

Once you've successfully recompiled with this change.  Try running
this
command:

met-5.0/bin/plot_data_plane sfc_20150616_36_grid_simple.grib1
2T_plot.ps
'name="var167"; level="Z0"; GRIB1_ptv=128;'

That should produce that attached image.

But why did we call it "var167"?  That's because met-5.0 does not
include
support for ECMWF GRIB tables.  Running wgrib on the file you sent, I
see
that the "2T" record is encoded as GRIB code 167 from parameter table
version number 128:

wgrib -V -d 10 sfc_20150616_36_grid_simple.grib1
rec 10:62251956:date 2015061612 2T kpds5=167 kpds6=1 kpds7=0
levels=(0,0)
grid=255 sfc 36hr fcst:
  2T=2 metre temperature [K]
  timerange 0 P1 36 P2 0 TimeU 1  nx 2880 ny 1441 GDS grid 0
num_in_ave 0
missing 0
  center 98 subcenter 0 process 145 Table 128 scan: WE:NS winds(N/S)
  latlon: lat  90.000000 to -90.000000 by 0.125000  nxny 4150080
          long 0.000000 to -0.125000 by 0.125000, (2880 x 1441) scan 0
mode
128 bdsgrid 1
  min/max data 202.203 320.359  num bits 12  BDS_Ref 202.203  DecScale
0
BinScale -5

You could edit the GRIB tables used by MET if you'd like
(met-5.0/data/table_files/nceptab_flat.txt).  Look for the line
beginning
with "167 128".  You could replace "var167" with "2T" if you'd like.

Hope that helps.

Thanks,
John


On Mon, Aug 24, 2015 at 9:12 PM, 柳媛普 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
>
>
> Hello John
>
> I has repost the grib data in you ftp site. It was thransfered
completed
> this time.
>
> yuanpu.
>
> > -----原始邮件-----
> > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > 发送时间: 2015年8月25日 星期二
> > 收件人: lyp at lzb.ac.cn
> > 抄送:
> > 主题: Re: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> >
> > Yuanpu,
> >
> > Unfortunately, it looks like the GRIB file you posted to our ftp
site is
> > corrupted.  I get an error when trying to read it.  Please look at
the
> size
> > of the file you tried to post:
> >
> >    -rw-rw-r--   1 johnhg rap 1.5M Aug 24 10:07
> > sfc_20150616_36_grid_simple.grib1
> >
> > Is it supposed to be larger than 1.5M?
> >
> > Thanks,
> > John
> >
> >
> > On Sun, Aug 23, 2015 at 10:13 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > >
> > > Hello John,
> > >
> > > Unfortunately , when I set lines as     static const int
> gribfile_buf_size
> > > = 16777216;   //  2^24
> > >
> > > and mke "make clean" and "make install". But I see the same
error:
> > >
> > > DEBUG 1: Default Config File:
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > ERROR  :
> > > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record
> larger
> > > than the buffer size.
> > > ERROR  :
> > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > ERROR  :
> > > ERROR  :
> > >
> > > I has post the girb1 data in liu_data directory of your ftp
site.
> > >
> > > yuanpu
> > >
> > >
> > >
> > > > -----原始邮件-----
> > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > 发送时间: 2015年8月22日 星期六
> > > > 收件人: lyp at lzb.ac.cn
> > > > 抄送:
> > > > 主题: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > > >
> > > > Yuanpu,
> > > >
> > > > The attachment didn't come through.  Please post it to our
anonymous
> ftp
> > > > site following these instructions:
> > > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Yes, that error you're seeing is a pesky little issue.  We do
have a
> > > > development task defined to dynamically allocate the amount of
space
> > > > required rather than having a fixed upper limit.
> > > >
> > > > Please try editing this source code file:
> > > >    src/libcode/vx_data2d_grib/grib_classes.h
> > > >
> > > > Find and set this line as follows:
> > > >    static const int gribfile_buf_size = 16777216;   //  2^24
> > > >
> > > > And then from the top-level MET directory, do a "make clean"
> followed by
> > > a
> > > > "make install".
> > > >
> > > > Please let me know how it goes after making that change.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > >
> > > > > Hello John,
> > > > >
> > > > > Thanks for your reply. The original ECMWF data were complex
> packing,
> > > so I
> > > > > decoded them with grib tool. But I see the new error :
> > > > >
> > > > > DEBUG 1: Default Config File:
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > ERROR  :
> > > > > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record
> > > larger
> > > > > than the buffer size.
> > > > > ERROR  :
> > > > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > > > ERROR  :
> > > > > ERROR  :
> > > > >
> > > > > I put the grib1 data in attachment.
> > > > >
> > > > > yuanpu
> > > > >
> > > > > > -----原始邮件-----
> > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > 发送时间: 2015年8月20日 星期四
> > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > 抄送:
> > > > > > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > > > > >
> > > > > > Yuanpu,
> > > > > >
> > > > > > Thanks for sending that data.  When I run it through MET,
I see
> the
> > > same
> > > > > > error you're getting:
> > > > > >    ERROR  : GribFile::read_record(GribRecord &) -> Second
order
> > > packing
> > > > > not
> > > > > > implemented.
> > > > > >
> > > > > > Running it through wgrib to inventory the file, I get a
similar
> > > error:
> > > > > >    *** Cannot decode complex packed fields n=20089***
> > > > > >
> > > > > > And trying to plot the data with the Integrated Data
Viewer
> (IDV), I
> > > also
> > > > > > get a similar error:
> > > > > >   Couldn't get data
> > > > > >   Grib1BinaryDataSection: (octet 4, 1st half) not grid
point
> data and
> > > > > > simple packing
> > > > > >
> > > > > > Take a look at GRIB1 table 11:
> > > > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> > > > > >
> > > > > > Having the second bit of BDS octet 4 set to 1 indicates
that the
> data
> > > > > uses
> > > > > > complex packing, which is generally not supported.
> > > > > >
> > > > > > I also tried viewing the NetCDF data file you sent me, the
> output of
> > > > > > ncl_convert2nc, using the ncview utility.  While ncview is
able
> to
> > > open
> > > > > the
> > > > > > file, all of the data values plotted are 0.
> > > > > >
> > > > > > Honestly, this is the first time I've seen this error.  My
guess
> is
> > > that
> > > > > > your GRIB1 data file is corrupted or badly created in some
way.
> I'd
> > > > > take a
> > > > > > close look at how that GRIB1 file is being created.  Once
you're
> > > able to
> > > > > > view the data using the IDV tool, you should be able to
use it in
> > > MET.
> > > > > >
> > > > > > Thanks,
> > > > > > John
> > > > > >
> > > > > > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT
<met_help at ucar.edu>
> > > wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> >
> > > > > > >
> > > > > > > Hello John
> > > > > > >
> > > > > > > Thank you for you reply. I have put the ECMWF
> > > > > > >
> data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-ACHN.bin)
> > > on
> > > > > your
> > > > > > > anonymous ftp site, it was in
incoming/irap/met_help/liu_data
> > > > > dirctory.but
> > > > > > > the observed data was not thransfered completed. So I
put it in
> > > > > attachment.
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > yuanpu
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > > -----原始邮件-----
> > > > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > > > 发送时间: 2015年8月18日 星期二
> > > > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > > > 抄送:
> > > > > > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified ECMWF
Data
> > > > > > > >
> > > > > > > > Yuanpu,
> > > > > > > >
> > > > > > > > Unfortunately, the NetCDF output of the ncl_convert2nc
tool
> is
> > > not
> > > > > > > > CF-compliant so MET will not be able to read it.  Can
you
> please
> > > > > send me
> > > > > > > > the original ECMWF GRIB1 file that produced the error
you're
> > > > > seeing.  It
> > > > > > > > doesn't need to be the full file, just a single record
would
> > > suffice.
> > > > > > > >
> > > > > > > > You can post it to our anonymous ftp site following
the
> > > instructions
> > > > > > > listed
> > > > > > > > here:
> > > > > > > >
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > > > >
> > > > > > > > Please write back once you've posted the data so that
I can
> go
> > > grab
> > > > > it.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > John Halley Gotway
> > > > > > > > met_help at ucar.edu
> > > > > > > >
> > > > > > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <
> met_help at ucar.edu>
> > > > > wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > Dear all,
> > > > > > > > >
> > > > > > > > > I want to verify ECMWF data with point_stat tool of
MET
> 5.0.
> > > When
> > > > > I
> > > > > > > used
> > > > > > > > > the original data(grib format), it has the followed
error
> > > > > information:
> > > > > > > > >
> > > > > > > > > DEBUG 1: Default Config File:
> > > > > > > > >
> > > > > > >
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > > > ERROR  :
> > > > > > > > > ERROR  : GribFile::read_record(GribRecord &) ->
Second
> order
> > > > > packing
> > > > > > > not
> > > > > > > > > implemented.
> > > > > > > > > ERROR  :
> > > > > > > > >
> > > > > > > > >  But when I  thransfered  the original data to nc
format
> with
> > > NCL,
> > > > > it
> > > > > > > has
> > > > > > > > > the followed error information:
> > > > > > > > >
> > > > > > > > > DEBUG 1: Default Config File:
> > > > > > > > >
> > > > > > >
> > > > >
> > >
> /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > > > > NetCDF: Attribute not found
> > > > > > > > >
> > > > > > > > > The original ECMWF nc format data in attachment.
> > > > > > > > >
> > > > > > > > > Kindly advise.
> > > > > > > > >
> > > > > > > > > Thanks in advance.
> > > > > > > > >
> > > > > > > > > Yuanpu
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > 从lzb.ac.cn发来的超大附件
> > > > > > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > > > > > 查看下载信息:
> > > > > > >
> > > > >
> > >
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > 从lzb.ac.cn发来的超大附件
> > > > > sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> > > > > 查看下载信息:
> > > > >
> > >
>
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > >
> > >
> > >
> > >
> > >
> > >
> >
>
>
>
>
>

------------------------------------------------
Subject: Problem with point stat tool
From: 柳媛普
Time: Wed Aug 26 01:39:04 2015

Hello John,

Thanks for your helps. The met tool read the ECMWF data correctly when
I tried to follow your indication.
I appreciated it.

yuanpu

> -----原始邮件-----
> 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> 发送时间: 2015年8月26日 星期三
> 收件人: lyp at lzb.ac.cn
> 抄送:
> 主题: Re: Re: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified
ECMWF Data
>
> Yuanpu,
>
> Thanks for sending that data file.  I was able to grab it and work
with it.
>
> Making that change I indicated earlier did fix the problem.  Please
try the
> following:
>
> (1) Save the attached version of "grib_classes.h" to your machine.
And
> place it into the met-5.0 directory:
>   mv grib_classes.h met-
5.0/src/libcode/vx_data2d_grib/grib_classes.h
>
> (2) Rebuild met-5.0 (be sure to run "make clean"):
>   cd met-5.0
>   make clean
>   make install
>
> Once you've successfully recompiled with this change.  Try running
this
> command:
>
> met-5.0/bin/plot_data_plane sfc_20150616_36_grid_simple.grib1
2T_plot.ps
> 'name="var167"; level="Z0"; GRIB1_ptv=128;'
>
> That should produce that attached image.
>
> But why did we call it "var167"?  That's because met-5.0 does not
include
> support for ECMWF GRIB tables.  Running wgrib on the file you sent,
I see
> that the "2T" record is encoded as GRIB code 167 from parameter
table
> version number 128:
>
> wgrib -V -d 10 sfc_20150616_36_grid_simple.grib1
> rec 10:62251956:date 2015061612 2T kpds5=167 kpds6=1 kpds7=0
levels=(0,0)
> grid=255 sfc 36hr fcst:
>   2T=2 metre temperature [K]
>   timerange 0 P1 36 P2 0 TimeU 1  nx 2880 ny 1441 GDS grid 0
num_in_ave 0
> missing 0
>   center 98 subcenter 0 process 145 Table 128 scan: WE:NS winds(N/S)
>   latlon: lat  90.000000 to -90.000000 by 0.125000  nxny 4150080
>           long 0.000000 to -0.125000 by 0.125000, (2880 x 1441) scan
0 mode
> 128 bdsgrid 1
>   min/max data 202.203 320.359  num bits 12  BDS_Ref 202.203
DecScale 0
> BinScale -5
>
> You could edit the GRIB tables used by MET if you'd like
> (met-5.0/data/table_files/nceptab_flat.txt).  Look for the line
beginning
> with "167 128".  You could replace "var167" with "2T" if you'd like.
>
> Hope that helps.
>
> Thanks,
> John
>
>
> On Mon, Aug 24, 2015 at 9:12 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> >
> >
> > Hello John
> >
> > I has repost the grib data in you ftp site. It was thransfered
completed
> > this time.
> >
> > yuanpu.
> >
> > > -----原始邮件-----
> > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > 发送时间: 2015年8月25日 星期二
> > > 收件人: lyp at lzb.ac.cn
> > > 抄送:
> > > 主题: Re: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified
ECMWF Data
> > >
> > > Yuanpu,
> > >
> > > Unfortunately, it looks like the GRIB file you posted to our ftp
site is
> > > corrupted.  I get an error when trying to read it.  Please look
at the
> > size
> > > of the file you tried to post:
> > >
> > >    -rw-rw-r--   1 johnhg rap 1.5M Aug 24 10:07
> > > sfc_20150616_36_grid_simple.grib1
> > >
> > > Is it supposed to be larger than 1.5M?
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Sun, Aug 23, 2015 at 10:13 PM, 柳媛普 via RT <met_help at ucar.edu>
wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
>
> > > >
> > > > Hello John,
> > > >
> > > > Unfortunately , when I set lines as     static const int
> > gribfile_buf_size
> > > > = 16777216;   //  2^24
> > > >
> > > > and mke "make clean" and "make install". But I see the same
error:
> > > >
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > ERROR  :
> > > > ERROR  : GribFile::read_record(GribRecord &) -> found a grib
record
> > larger
> > > > than the buffer size.
> > > > ERROR  :
> > > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > > ERROR  :
> > > > ERROR  :
> > > >
> > > > I has post the girb1 data in liu_data directory of your ftp
site.
> > > >
> > > > yuanpu
> > > >
> > > >
> > > >
> > > > > -----原始邮件-----
> > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > 发送时间: 2015年8月22日 星期六
> > > > > 收件人: lyp at lzb.ac.cn
> > > > > 抄送:
> > > > > 主题: Re: Re: Re: [rt.rap.ucar.edu #71731] About verified
ECMWF Data
> > > > >
> > > > > Yuanpu,
> > > > >
> > > > > The attachment didn't come through.  Please post it to our
anonymous
> > ftp
> > > > > site following these instructions:
> > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > >
> > > > > Yes, that error you're seeing is a pesky little issue.  We
do have a
> > > > > development task defined to dynamically allocate the amount
of space
> > > > > required rather than having a fixed upper limit.
> > > > >
> > > > > Please try editing this source code file:
> > > > >    src/libcode/vx_data2d_grib/grib_classes.h
> > > > >
> > > > > Find and set this line as follows:
> > > > >    static const int gribfile_buf_size = 16777216;   //  2^24
> > > > >
> > > > > And then from the top-level MET directory, do a "make clean"
> > followed by
> > > > a
> > > > > "make install".
> > > > >
> > > > > Please let me know how it goes after making that change.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Thu, Aug 20, 2015 at 8:51 PM, 柳媛普 via RT
<met_help at ucar.edu>
> > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731 >
> > > > > >
> > > > > > Hello John,
> > > > > >
> > > > > > Thanks for your reply. The original ECMWF data were
complex
> > packing,
> > > > so I
> > > > > > decoded them with grib tool. But I see the new error :
> > > > > >
> > > > > > DEBUG 1: Default Config File:
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > DEBUG 1: User Config File: ./PointStatConfig_sfc_EC
> > > > > > ERROR  :
> > > > > > ERROR  : GribFile::read_record(GribRecord &) -> found a
grib record
> > > > larger
> > > > > > than the buffer size.
> > > > > > ERROR  :
> > > > > > ERROR  :   Increase the buffer to at least 8819010 bytes.
> > > > > > ERROR  :
> > > > > > ERROR  :
> > > > > >
> > > > > > I put the grib1 data in attachment.
> > > > > >
> > > > > > yuanpu
> > > > > >
> > > > > > > -----原始邮件-----
> > > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > > 发送时间: 2015年8月20日 星期四
> > > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > > 抄送:
> > > > > > > 主题: Re: Re: [rt.rap.ucar.edu #71731] About verified
ECMWF Data
> > > > > > >
> > > > > > > Yuanpu,
> > > > > > >
> > > > > > > Thanks for sending that data.  When I run it through
MET, I see
> > the
> > > > same
> > > > > > > error you're getting:
> > > > > > >    ERROR  : GribFile::read_record(GribRecord &) ->
Second order
> > > > packing
> > > > > > not
> > > > > > > implemented.
> > > > > > >
> > > > > > > Running it through wgrib to inventory the file, I get a
similar
> > > > error:
> > > > > > >    *** Cannot decode complex packed fields n=20089***
> > > > > > >
> > > > > > > And trying to plot the data with the Integrated Data
Viewer
> > (IDV), I
> > > > also
> > > > > > > get a similar error:
> > > > > > >   Couldn't get data
> > > > > > >   Grib1BinaryDataSection: (octet 4, 1st half) not grid
point
> > data and
> > > > > > > simple packing
> > > > > > >
> > > > > > > Take a look at GRIB1 table 11:
> > > > > > >
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table11.html
> > > > > > >
> > > > > > > Having the second bit of BDS octet 4 set to 1 indicates
that the
> > data
> > > > > > uses
> > > > > > > complex packing, which is generally not supported.
> > > > > > >
> > > > > > > I also tried viewing the NetCDF data file you sent me,
the
> > output of
> > > > > > > ncl_convert2nc, using the ncview utility.  While ncview
is able
> > to
> > > > open
> > > > > > the
> > > > > > > file, all of the data values plotted are 0.
> > > > > > >
> > > > > > > Honestly, this is the first time I've seen this error.
My guess
> > is
> > > > that
> > > > > > > your GRIB1 data file is corrupted or badly created in
some way.
> > I'd
> > > > > > take a
> > > > > > > close look at how that GRIB1 file is being created.
Once you're
> > > > able to
> > > > > > > view the data using the IDV tool, you should be able to
use it in
> > > > MET.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > John
> > > > > > >
> > > > > > > On Tue, Aug 18, 2015 at 3:33 AM, 柳媛普 via RT
<met_help at ucar.edu>
> > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> > >
> > > > > > > >
> > > > > > > > Hello John
> > > > > > > >
> > > > > > > > Thank you for you reply. I have put the ECMWF
> > > > > > > >
> > data(W_NAFP_C_ECMF_20150718201205_P_C3E07181200071815001-ACHN.bin)
> > > > on
> > > > > > your
> > > > > > > > anonymous ftp site, it was in
incoming/irap/met_help/liu_data
> > > > > > dirctory.but
> > > > > > > > the observed data was not thransfered completed. So I
put it in
> > > > > > attachment.
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > yuanpu
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > > -----原始邮件-----
> > > > > > > > > 发件人: "John Halley Gotway via RT" <met_help at ucar.edu>
> > > > > > > > > 发送时间: 2015年8月18日 星期二
> > > > > > > > > 收件人: lyp at lzb.ac.cn
> > > > > > > > > 抄送:
> > > > > > > > > 主题: Re: [rt.rap.ucar.edu #71731] About verified
ECMWF Data
> > > > > > > > >
> > > > > > > > > Yuanpu,
> > > > > > > > >
> > > > > > > > > Unfortunately, the NetCDF output of the
ncl_convert2nc tool
> > is
> > > > not
> > > > > > > > > CF-compliant so MET will not be able to read it.
Can you
> > please
> > > > > > send me
> > > > > > > > > the original ECMWF GRIB1 file that produced the
error you're
> > > > > > seeing.  It
> > > > > > > > > doesn't need to be the full file, just a single
record would
> > > > suffice.
> > > > > > > > >
> > > > > > > > > You can post it to our anonymous ftp site following
the
> > > > instructions
> > > > > > > > listed
> > > > > > > > > here:
> > > > > > > > >
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > > > > >
> > > > > > > > > Please write back once you've posted the data so
that I can
> > go
> > > > grab
> > > > > > it.
> > > > > > > > >
> > > > > > > > > Thanks,
> > > > > > > > > John Halley Gotway
> > > > > > > > > met_help at ucar.edu
> > > > > > > > >
> > > > > > > > > On Sun, Aug 16, 2015 at 8:51 PM, 柳媛普 via RT <
> > met_help at ucar.edu>
> > > > > > wrote:
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71731
> > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Dear all,
> > > > > > > > > >
> > > > > > > > > > I want to verify ECMWF data with point_stat tool
of  MET
> > 5.0.
> > > > When
> > > > > > I
> > > > > > > > used
> > > > > > > > > > the original data(grib format), it has the
followed error
> > > > > > information:
> > > > > > > > > >
> > > > > > > > > > DEBUG 1: Default Config File:
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > > > DEBUG 1: User Config File:
./PointStatConfig_sfc_EC
> > > > > > > > > > ERROR  :
> > > > > > > > > > ERROR  : GribFile::read_record(GribRecord &) ->
Second
> > order
> > > > > > packing
> > > > > > > > not
> > > > > > > > > > implemented.
> > > > > > > > > > ERROR  :
> > > > > > > > > >
> > > > > > > > > >  But when I  thransfered  the original data to nc
format
> > with
> > > > NCL,
> > > > > > it
> > > > > > > > has
> > > > > > > > > > the followed error information:
> > > > > > > > > >
> > > > > > > > > > DEBUG 1: Default Config File:
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > /home/liuyp/local/met-5.0/met-
install/share/met/config/PointStatConfig_default
> > > > > > > > > > DEBUG 1: User Config File:
./PointStatConfig_sfc_EC
> > > > > > > > > > NetCDF: Attribute not found
> > > > > > > > > >
> > > > > > > > > > The original ECMWF nc format data in attachment.
> > > > > > > > > >
> > > > > > > > > > Kindly advise.
> > > > > > > > > >
> > > > > > > > > > Thanks in advance.
> > > > > > > > > >
> > > > > > > > > > Yuanpu
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > 从lzb.ac.cn发来的超大附件
> > > > > > > > Z_SURF_ascii_20150718.nc (151MB, 2015-09-02 到期)
> > > > > > > > 查看下载信息:
> > > > > > > >
> > > > > >
> > > >
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU0xU2DjjU83U4DjqUn5U3DjYU0kUfUj2U88U77jLUe7j-
7wWB96J&code=u777hw9e
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > 从lzb.ac.cn发来的超大附件
> > > > > > sfc_20150616_36_grid_simple.grib1 (166MB, 2015-09-05 到期)
> > > > > > 查看下载信息:
> > > > > >
> > > >
> >
https://mail.cstnet.cn/coremail/viewDownloadFile.jsp?key=1U31SsvkjDCdhVjUTemgo90yo9FLUu2kTDa0SstETevJon2LzsGN3gKq3XKyTgGka1AXFbGpS4APTuxE3Z0XasIX-
4lpf9fETemgo90yo9FLwujXTWCyfu2LaV7L3ZtmjeIk-
s71UUUUjU2l39EtjqanVW8XFWDZrWfur1xWr15yVZGKf8jdUZCmaUtsU18USUjgUnkU7DjBU8kUc7j0U83U4DjqUn5U3DjqUnkUSUjBU88U77jLUn8j-
7rX86ZX&code=cy7f45hk
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > >
> >
> >
> >
> >
> >
>




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


More information about the Met_help mailing list