[Met_help] [rt.rap.ucar.edu #91664] History for SL1L2 output from grid_stat

John Halley Gotway via RT met_help at ucar.edu
Mon Aug 26 09:40:28 MDT 2019


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

Hi,

I'm running grid_stat in version 8.1 and trying to figure out why I get duplicate records when I output SL1L2.  Every line is written out twice, and the pair are identical for all fields.  Here is some sample output with the repeated lines:

VERSION MODEL   DESC FCST_LEAD FCST_VALID_BEG  FCST_VALID_END  OBS_LEAD OBS_VALID_BEG   OBS_VALID_END   FCST_VAR FCST_UNITS FCST_LEV OBS_VAR OBS_
UNITS OBS_LEV OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH COV_THRESH ALPHA LINE_TYPE
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 TMP      K          P400     TMP     K
      P400    ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240    241.48153    241.48947       58
488.4646        58487.39199       58489.87999  0.43695
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 TMP      K          P400     TMP     K
      P400    ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240    241.48153    241.48947       58
488.4646        58487.39199       58489.87999  0.43695
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 TMP      K          Z2       TMP     K
      Z2      ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240    277.3807     277.30876       77
472.23732       77490.92909       77455.03569  0.71251
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 TMP      K          Z2       TMP     K
      Z2      ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240    277.3807     277.30876       77
472.23732       77490.92909       77455.03569  0.71251
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 HGT      gpm        P400     HGT     gpm
      P400    ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240   7109.98602   7108.65435    50746
400.5502     50755397.41531    50737446.60268  4.27849
V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000 000000   20190401_120000 20190401_120000 HGT      gpm        P400     HGT     gpm
      P400    ANALYS FULL    NEAREST     1           NA          NA         NA         NA    SL1L2     1038240   7109.98602   7108.65435    50746
400.5502     50755397.41531    50737446.60268  4.27849

I'll post my config file below.  I don't see any obvious reason why this happens.

//
// Output model name to be written
//
model = "GALWEM1";

//
// Output description to be written
// May be set separately in each "obs.field" entry
//
desc = "NA";

//
// Output observation type to be written
//
obtype = "ANALYS";

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

//
// Verification grid
//
regrid = {
   to_grid    = NONE;
   method     = NEAREST;
   width      = 1;
   vld_thresh = 1.0;
   shape      = SQUARE;
}

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

//
// May be set separately in each "field" entry
//
censor_thresh    = [];
censor_val       = [];
cat_thresh       = [];
cnt_thresh       = [ NA ];
cnt_logic        = UNION;
wind_thresh      = [ NA ];
wind_logic       = UNION;
eclv_points      = 0.05;
nc_pairs_var_str = "";
rank_corr_flag   = FALSE;

//
// Forecast and observation fields to be verified
//
fcst = {

   file_type = GRIB2;
   field = [
      { name = "TMP"; level = "P400"; },
      { name = "TMP"; level = "Z2"; },
      { name = "HGT"; level = "P400"; },
      { name = "DPT"; level = "P850"; },
      { name = "DPT"; level = "P700"; },
      { name = "DPT"; level = "P400"; },
      { name = "DPT"; level = "Z2"; },
      { name = "WIND"; level = "Z10"; },
      { name = "WIND"; level = "P850"; },
      { name = "WIND"; level = "P400"; },
      { name = "WIND"; level = "P250"; },
      { name = "PRMSL"; level = "Z0"; }
   ];

}

obs = fcst;

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

//
// Climatology data
//
climo_mean = {

   file_name = [];
   field     = [];

   regrid = {
      method     = NEAREST;
      width      = 1;
      vld_thresh = 1.0;
      shape      = SQUARE;
   }

   time_interp_method = DW_MEAN;
   match_month        = TRUE;
   match_day          = FALSE;
   time_step          = 21600;
}

climo_stdev = climo_mean;
climo_stdev = {
   file_name = [];
}

climo_cdf_bins = 1;

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

//
// Verification masking regions
//
mask = {
   grid = [ "FULL" ];
   poly = [ "" ];
}

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

//
// Confidence interval settings
//
ci_alpha  = [ 0.05 ];

boot = {
   interval = PCTILE;
   rep_prop = 1.0;
   n_rep    = 0;
   rng      = "mt19937";
   seed     = "";
}

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

//
// Data smoothing methods
//
interp = {
   field      = NONE;
   vld_thresh = 1.0;
   shape      = SQUARE;

   type = [
      {
         method = NEAREST;
         width  = 1;
      }
   ];
}

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

//
// Neighborhood methods
//
nbrhd = {
   field      = BOTH;
   vld_thresh = 1.0;
   shape      = SQUARE;
   width      = [ 1 ];
   cov_thresh = [ >=0.5 ];
}

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

//
// Fourier decomposition
//
fourier = {
   wave_1d_beg = [];
   wave_1d_end = [];
}

//
// Gradient statistics
// May be set separately in each "obs.field" entry
//
gradient = {
   dx = [ 1 ];
   dy = [ 1 ];
}

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

//
// Statistical output types
//
output_flag = {
   fho    = NONE;
   ctc    = NONE;
   cts    = NONE;
   mctc   = NONE;
   mcts   = NONE;
   cnt    = NONE;
   sl1l2  = STAT;
   sal1l2 = NONE;
   vl1l2  = NONE;
   val1l2 = NONE;
   pct    = NONE;
   pstd   = NONE;
   pjc    = NONE;
   prc    = NONE;
   eclv   = NONE;
   nbrctc = NONE;
   nbrcts = NONE;
   nbrcnt = NONE;
   grad   = NONE;
}

//
// NetCDF matched pairs output file
//
nc_pairs_flag   = {
   latlon     = FALSE;
   raw        = FALSE;
   diff       = FALSE;
   climo      = FALSE;
   weight     = FALSE;
   nbrhd      = FALSE;
   fourier    = FALSE;
   gradient   = FALSE;
   apply_mask = FALSE;
}

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

grid_weight_flag = NONE;
tmp_dir          = "/tmp";
output_prefix    = "GALWEM1";
version          = "V8.1";

Thanks,
Matt

Matthew C. Sittel, Meteorologist
University Corporation for Atmospheric Research
Matthew.Sittel.ctr at us.af.mil



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

Subject: SL1L2 output from grid_stat
From: David Fillmore
Time: Fri Aug 23 10:23:58 2019

Hi Matt -

We will have a look.
Are you able to provide sample data so that we can replicate the issue
if
necessary?
See ftp upload instructions here:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk


thanks,
David

On Fri, Aug 23, 2019 at 9:40 AM SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
via RT <met_help at ucar.edu> wrote:

>
> Fri Aug 23 09:40:26 2019: Request 91664 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
>        Queue: met_help
>      Subject: SL1L2 output from grid_stat
>        Owner: Nobody
>   Requestors: matthew.sittel.ctr at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664 >
>
>
> Hi,
>
> I'm running grid_stat in version 8.1 and trying to figure out why I
get
> duplicate records when I output SL1L2.  Every line is written out
twice,
> and the pair are identical for all fields.  Here is some sample
output with
> the repeated lines:
>
> VERSION MODEL   DESC FCST_LEAD FCST_VALID_BEG  FCST_VALID_END
OBS_LEAD
> OBS_VALID_BEG   OBS_VALID_END   FCST_VAR FCST_UNITS FCST_LEV OBS_VAR
OBS_
> UNITS OBS_LEV OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH
> OBS_THRESH COV_THRESH ALPHA LINE_TYPE
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> 488.4646        58487.39199       58489.87999  0.43695
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> 488.4646        58487.39199       58489.87999  0.43695
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
>       Z2      ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> 472.23732       77490.92909       77455.03569  0.71251
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
>       Z2      ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> 472.23732       77490.92909       77455.03569  0.71251
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> 400.5502     50755397.41531    50737446.60268  4.27849
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> 400.5502     50755397.41531    50737446.60268  4.27849
>
> I'll post my config file below.  I don't see any obvious reason why
this
> happens.
>
> //
> // Output model name to be written
> //
> model = "GALWEM1";
>
> //
> // Output description to be written
> // May be set separately in each "obs.field" entry
> //
> desc = "NA";
>
> //
> // Output observation type to be written
> //
> obtype = "ANALYS";
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Verification grid
> //
> regrid = {
>    to_grid    = NONE;
>    method     = NEAREST;
>    width      = 1;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // May be set separately in each "field" entry
> //
> censor_thresh    = [];
> censor_val       = [];
> cat_thresh       = [];
> cnt_thresh       = [ NA ];
> cnt_logic        = UNION;
> wind_thresh      = [ NA ];
> wind_logic       = UNION;
> eclv_points      = 0.05;
> nc_pairs_var_str = "";
> rank_corr_flag   = FALSE;
>
> //
> // Forecast and observation fields to be verified
> //
> fcst = {
>
>    file_type = GRIB2;
>    field = [
>       { name = "TMP"; level = "P400"; },
>       { name = "TMP"; level = "Z2"; },
>       { name = "HGT"; level = "P400"; },
>       { name = "DPT"; level = "P850"; },
>       { name = "DPT"; level = "P700"; },
>       { name = "DPT"; level = "P400"; },
>       { name = "DPT"; level = "Z2"; },
>       { name = "WIND"; level = "Z10"; },
>       { name = "WIND"; level = "P850"; },
>       { name = "WIND"; level = "P400"; },
>       { name = "WIND"; level = "P250"; },
>       { name = "PRMSL"; level = "Z0"; }
>    ];
>
> }
>
> obs = fcst;
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Climatology data
> //
> climo_mean = {
>
>    file_name = [];
>    field     = [];
>
>    regrid = {
>       method     = NEAREST;
>       width      = 1;
>       vld_thresh = 1.0;
>       shape      = SQUARE;
>    }
>
>    time_interp_method = DW_MEAN;
>    match_month        = TRUE;
>    match_day          = FALSE;
>    time_step          = 21600;
> }
>
> climo_stdev = climo_mean;
> climo_stdev = {
>    file_name = [];
> }
>
> climo_cdf_bins = 1;
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Verification masking regions
> //
> mask = {
>    grid = [ "FULL" ];
>    poly = [ "" ];
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Confidence interval settings
> //
> ci_alpha  = [ 0.05 ];
>
> boot = {
>    interval = PCTILE;
>    rep_prop = 1.0;
>    n_rep    = 0;
>    rng      = "mt19937";
>    seed     = "";
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Data smoothing methods
> //
> interp = {
>    field      = NONE;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
>
>    type = [
>       {
>          method = NEAREST;
>          width  = 1;
>       }
>    ];
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Neighborhood methods
> //
> nbrhd = {
>    field      = BOTH;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
>    width      = [ 1 ];
>    cov_thresh = [ >=0.5 ];
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Fourier decomposition
> //
> fourier = {
>    wave_1d_beg = [];
>    wave_1d_end = [];
> }
>
> //
> // Gradient statistics
> // May be set separately in each "obs.field" entry
> //
> gradient = {
>    dx = [ 1 ];
>    dy = [ 1 ];
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> //
> // Statistical output types
> //
> output_flag = {
>    fho    = NONE;
>    ctc    = NONE;
>    cts    = NONE;
>    mctc   = NONE;
>    mcts   = NONE;
>    cnt    = NONE;
>    sl1l2  = STAT;
>    sal1l2 = NONE;
>    vl1l2  = NONE;
>    val1l2 = NONE;
>    pct    = NONE;
>    pstd   = NONE;
>    pjc    = NONE;
>    prc    = NONE;
>    eclv   = NONE;
>    nbrctc = NONE;
>    nbrcts = NONE;
>    nbrcnt = NONE;
>    grad   = NONE;
> }
>
> //
> // NetCDF matched pairs output file
> //
> nc_pairs_flag   = {
>    latlon     = FALSE;
>    raw        = FALSE;
>    diff       = FALSE;
>    climo      = FALSE;
>    weight     = FALSE;
>    nbrhd      = FALSE;
>    fourier    = FALSE;
>    gradient   = FALSE;
>    apply_mask = FALSE;
> }
>
>
>
////////////////////////////////////////////////////////////////////////////////
>
> grid_weight_flag = NONE;
> tmp_dir          = "/tmp";
> output_prefix    = "GALWEM1";
> version          = "V8.1";
>
> Thanks,
> Matt
>
> Matthew C. Sittel, Meteorologist
> University Corporation for Atmospheric Research
> Matthew.Sittel.ctr at us.af.mil
>
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2 output from grid_stat
From: SITTEL, MATTHEW C CTR USAF AFWA 16 WS/WXN
Time: Fri Aug 23 10:41:29 2019

FTP doesn't work for us; we can't get through the Air Force firewall.
We have to use AMRDEC to transfer files.  It's a bit of a hassle,
which is why I sent what we did, in hopes that the sample output
showed the nature of the problem, and the config file contains the
exact settings I used.  I'm hoping that illustrates the problem.

The SL1L2 output will be used to calculate GO Index in stat_analysis.
Duplicate records slow down the process and potentially contaminate
the results if this repetitive output isn't consistent from run to
run.

Matt

-----Original Message-----
From: David Fillmore via RT <met_help at ucar.edu>
Sent: Friday, August 23, 2019 11:24 AM
To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL
<matthew.sittel.ctr at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2 output
from grid_stat

Hi Matt -

We will have a look.
Are you able to provide sample data so that we can replicate the issue
if necessary?
See ftp upload instructions here:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk


thanks,
David

On Fri, Aug 23, 2019 at 9:40 AM SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN via RT <met_help at ucar.edu> wrote:

>
> Fri Aug 23 09:40:26 2019: Request 91664 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
>        Queue: met_help
>      Subject: SL1L2 output from grid_stat
>        Owner: Nobody
>   Requestors: matthew.sittel.ctr at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664
> >
>
>
> Hi,
>
> I'm running grid_stat in version 8.1 and trying to figure out why I
> get duplicate records when I output SL1L2.  Every line is written
out
> twice, and the pair are identical for all fields.  Here is some
sample
> output with the repeated lines:
>
> VERSION MODEL   DESC FCST_LEAD FCST_VALID_BEG  FCST_VALID_END
OBS_LEAD
> OBS_VALID_BEG   OBS_VALID_END   FCST_VAR FCST_UNITS FCST_LEV OBS_VAR
OBS_
> UNITS OBS_LEV OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH
> OBS_THRESH COV_THRESH ALPHA LINE_TYPE
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> 488.4646        58487.39199       58489.87999  0.43695
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> 488.4646        58487.39199       58489.87999  0.43695
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
>       Z2      ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> 472.23732       77490.92909       77455.03569  0.71251
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
>       Z2      ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> 472.23732       77490.92909       77455.03569  0.71251
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> 400.5502     50755397.41531    50737446.60268  4.27849
> V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
>  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
>       P400    ANALYS FULL    NEAREST     1           NA          NA
>  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> 400.5502     50755397.41531    50737446.60268  4.27849
>
> I'll post my config file below.  I don't see any obvious reason why
> this happens.
>
> //
> // Output model name to be written
> //
> model = "GALWEM1";
>
> //
> // Output description to be written
> // May be set separately in each "obs.field" entry // desc = "NA";
>
> //
> // Output observation type to be written // obtype = "ANALYS";
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Verification grid
> //
> regrid = {
>    to_grid    = NONE;
>    method     = NEAREST;
>    width      = 1;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // May be set separately in each "field" entry //
> censor_thresh    = [];
> censor_val       = [];
> cat_thresh       = [];
> cnt_thresh       = [ NA ];
> cnt_logic        = UNION;
> wind_thresh      = [ NA ];
> wind_logic       = UNION;
> eclv_points      = 0.05;
> nc_pairs_var_str = "";
> rank_corr_flag   = FALSE;
>
> //
> // Forecast and observation fields to be verified // fcst = {
>
>    file_type = GRIB2;
>    field = [
>       { name = "TMP"; level = "P400"; },
>       { name = "TMP"; level = "Z2"; },
>       { name = "HGT"; level = "P400"; },
>       { name = "DPT"; level = "P850"; },
>       { name = "DPT"; level = "P700"; },
>       { name = "DPT"; level = "P400"; },
>       { name = "DPT"; level = "Z2"; },
>       { name = "WIND"; level = "Z10"; },
>       { name = "WIND"; level = "P850"; },
>       { name = "WIND"; level = "P400"; },
>       { name = "WIND"; level = "P250"; },
>       { name = "PRMSL"; level = "Z0"; }
>    ];
>
> }
>
> obs = fcst;
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Climatology data
> //
> climo_mean = {
>
>    file_name = [];
>    field     = [];
>
>    regrid = {
>       method     = NEAREST;
>       width      = 1;
>       vld_thresh = 1.0;
>       shape      = SQUARE;
>    }
>
>    time_interp_method = DW_MEAN;
>    match_month        = TRUE;
>    match_day          = FALSE;
>    time_step          = 21600;
> }
>
> climo_stdev = climo_mean;
> climo_stdev = {
>    file_name = [];
> }
>
> climo_cdf_bins = 1;
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Verification masking regions
> //
> mask = {
>    grid = [ "FULL" ];
>    poly = [ "" ];
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Confidence interval settings
> //
> ci_alpha  = [ 0.05 ];
>
> boot = {
>    interval = PCTILE;
>    rep_prop = 1.0;
>    n_rep    = 0;
>    rng      = "mt19937";
>    seed     = "";
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Data smoothing methods
> //
> interp = {
>    field      = NONE;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
>
>    type = [
>       {
>          method = NEAREST;
>          width  = 1;
>       }
>    ];
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Neighborhood methods
> //
> nbrhd = {
>    field      = BOTH;
>    vld_thresh = 1.0;
>    shape      = SQUARE;
>    width      = [ 1 ];
>    cov_thresh = [ >=0.5 ];
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Fourier decomposition
> //
> fourier = {
>    wave_1d_beg = [];
>    wave_1d_end = [];
> }
>
> //
> // Gradient statistics
> // May be set separately in each "obs.field" entry // gradient = {
>    dx = [ 1 ];
>    dy = [ 1 ];
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> //
> // Statistical output types
> //
> output_flag = {
>    fho    = NONE;
>    ctc    = NONE;
>    cts    = NONE;
>    mctc   = NONE;
>    mcts   = NONE;
>    cnt    = NONE;
>    sl1l2  = STAT;
>    sal1l2 = NONE;
>    vl1l2  = NONE;
>    val1l2 = NONE;
>    pct    = NONE;
>    pstd   = NONE;
>    pjc    = NONE;
>    prc    = NONE;
>    eclv   = NONE;
>    nbrctc = NONE;
>    nbrcts = NONE;
>    nbrcnt = NONE;
>    grad   = NONE;
> }
>
> //
> // NetCDF matched pairs output file
> //
> nc_pairs_flag   = {
>    latlon     = FALSE;
>    raw        = FALSE;
>    diff       = FALSE;
>    climo      = FALSE;
>    weight     = FALSE;
>    nbrhd      = FALSE;
>    fourier    = FALSE;
>    gradient   = FALSE;
>    apply_mask = FALSE;
> }
>
>
>
//////////////////////////////////////////////////////////////////////
> //////////
>
> grid_weight_flag = NONE;
> tmp_dir          = "/tmp";
> output_prefix    = "GALWEM1";
> version          = "V8.1";
>
> Thanks,
> Matt
>
> Matthew C. Sittel, Meteorologist
> University Corporation for Atmospheric Research
> Matthew.Sittel.ctr at us.af.mil
>
>
>



------------------------------------------------
Subject: SL1L2 output from grid_stat
From: John Halley Gotway
Time: Fri Aug 23 11:19:23 2019

Hi Matt and David,

Take a look at this setting in the Grid-Stat config file:
mask = {
   grid = [ "FULL" ];
   poly = [ "" ];
}

Please try changing it to the following:
mask = {
   grid = [ "FULL" ];
   poly = [ ];
}

And see if that makes the duplicates go away.  If this does fix that
behavior, we'll need to run this situation through the debugger and
fix the
behavior.  Instead of producing duplicate output lines, it *should*
print
an error message stating that it can't parse an empty string as the
name of
the polyline file.

Thanks,
John

On Fri, Aug 23, 2019 at 10:42 AM SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664 >
>
> FTP doesn't work for us; we can't get through the Air Force
firewall.  We
> have to use AMRDEC to transfer files.  It's a bit of a hassle, which
is why
> I sent what we did, in hopes that the sample output showed the
nature of
> the problem, and the config file contains the exact settings I used.
I'm
> hoping that illustrates the problem.
>
> The SL1L2 output will be used to calculate GO Index in
stat_analysis.
> Duplicate records slow down the process and potentially contaminate
the
> results if this repetitive output isn't consistent from run to run.
>
> Matt
>
> -----Original Message-----
> From: David Fillmore via RT <met_help at ucar.edu>
> Sent: Friday, August 23, 2019 11:24 AM
> To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL <
> matthew.sittel.ctr at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2 output
from
> grid_stat
>
> Hi Matt -
>
> We will have a look.
> Are you able to provide sample data so that we can replicate the
issue if
> necessary?
> See ftp upload instructions here:
>
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
>
>
> thanks,
> David
>
> On Fri, Aug 23, 2019 at 9:40 AM SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN
> via RT <met_help at ucar.edu> wrote:
>
> >
> > Fri Aug 23 09:40:26 2019: Request 91664 was acted upon.
> > Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
> >        Queue: met_help
> >      Subject: SL1L2 output from grid_stat
> >        Owner: Nobody
> >   Requestors: matthew.sittel.ctr at us.af.mil
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664
> > >
> >
> >
> > Hi,
> >
> > I'm running grid_stat in version 8.1 and trying to figure out why
I
> > get duplicate records when I output SL1L2.  Every line is written
out
> > twice, and the pair are identical for all fields.  Here is some
sample
> > output with the repeated lines:
> >
> > VERSION MODEL   DESC FCST_LEAD FCST_VALID_BEG  FCST_VALID_END
OBS_LEAD
> > OBS_VALID_BEG   OBS_VALID_END   FCST_VAR FCST_UNITS FCST_LEV
OBS_VAR OBS_
> > UNITS OBS_LEV OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH
> > OBS_THRESH COV_THRESH ALPHA LINE_TYPE
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> > 488.4646        58487.39199       58489.87999  0.43695
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> > 488.4646        58487.39199       58489.87999  0.43695
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
> >       Z2      ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> > 472.23732       77490.92909       77455.03569  0.71251
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
> >       Z2      ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> > 472.23732       77490.92909       77455.03569  0.71251
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> > 400.5502     50755397.41531    50737446.60268  4.27849
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> > 400.5502     50755397.41531    50737446.60268  4.27849
> >
> > I'll post my config file below.  I don't see any obvious reason
why
> > this happens.
> >
> > //
> > // Output model name to be written
> > //
> > model = "GALWEM1";
> >
> > //
> > // Output description to be written
> > // May be set separately in each "obs.field" entry // desc = "NA";
> >
> > //
> > // Output observation type to be written // obtype = "ANALYS";
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Verification grid
> > //
> > regrid = {
> >    to_grid    = NONE;
> >    method     = NEAREST;
> >    width      = 1;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // May be set separately in each "field" entry //
> > censor_thresh    = [];
> > censor_val       = [];
> > cat_thresh       = [];
> > cnt_thresh       = [ NA ];
> > cnt_logic        = UNION;
> > wind_thresh      = [ NA ];
> > wind_logic       = UNION;
> > eclv_points      = 0.05;
> > nc_pairs_var_str = "";
> > rank_corr_flag   = FALSE;
> >
> > //
> > // Forecast and observation fields to be verified // fcst = {
> >
> >    file_type = GRIB2;
> >    field = [
> >       { name = "TMP"; level = "P400"; },
> >       { name = "TMP"; level = "Z2"; },
> >       { name = "HGT"; level = "P400"; },
> >       { name = "DPT"; level = "P850"; },
> >       { name = "DPT"; level = "P700"; },
> >       { name = "DPT"; level = "P400"; },
> >       { name = "DPT"; level = "Z2"; },
> >       { name = "WIND"; level = "Z10"; },
> >       { name = "WIND"; level = "P850"; },
> >       { name = "WIND"; level = "P400"; },
> >       { name = "WIND"; level = "P250"; },
> >       { name = "PRMSL"; level = "Z0"; }
> >    ];
> >
> > }
> >
> > obs = fcst;
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Climatology data
> > //
> > climo_mean = {
> >
> >    file_name = [];
> >    field     = [];
> >
> >    regrid = {
> >       method     = NEAREST;
> >       width      = 1;
> >       vld_thresh = 1.0;
> >       shape      = SQUARE;
> >    }
> >
> >    time_interp_method = DW_MEAN;
> >    match_month        = TRUE;
> >    match_day          = FALSE;
> >    time_step          = 21600;
> > }
> >
> > climo_stdev = climo_mean;
> > climo_stdev = {
> >    file_name = [];
> > }
> >
> > climo_cdf_bins = 1;
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Verification masking regions
> > //
> > mask = {
> >    grid = [ "FULL" ];
> >    poly = [ "" ];
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Confidence interval settings
> > //
> > ci_alpha  = [ 0.05 ];
> >
> > boot = {
> >    interval = PCTILE;
> >    rep_prop = 1.0;
> >    n_rep    = 0;
> >    rng      = "mt19937";
> >    seed     = "";
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Data smoothing methods
> > //
> > interp = {
> >    field      = NONE;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> >
> >    type = [
> >       {
> >          method = NEAREST;
> >          width  = 1;
> >       }
> >    ];
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Neighborhood methods
> > //
> > nbrhd = {
> >    field      = BOTH;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> >    width      = [ 1 ];
> >    cov_thresh = [ >=0.5 ];
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Fourier decomposition
> > //
> > fourier = {
> >    wave_1d_beg = [];
> >    wave_1d_end = [];
> > }
> >
> > //
> > // Gradient statistics
> > // May be set separately in each "obs.field" entry // gradient = {
> >    dx = [ 1 ];
> >    dy = [ 1 ];
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > //
> > // Statistical output types
> > //
> > output_flag = {
> >    fho    = NONE;
> >    ctc    = NONE;
> >    cts    = NONE;
> >    mctc   = NONE;
> >    mcts   = NONE;
> >    cnt    = NONE;
> >    sl1l2  = STAT;
> >    sal1l2 = NONE;
> >    vl1l2  = NONE;
> >    val1l2 = NONE;
> >    pct    = NONE;
> >    pstd   = NONE;
> >    pjc    = NONE;
> >    prc    = NONE;
> >    eclv   = NONE;
> >    nbrctc = NONE;
> >    nbrcts = NONE;
> >    nbrcnt = NONE;
> >    grad   = NONE;
> > }
> >
> > //
> > // NetCDF matched pairs output file
> > //
> > nc_pairs_flag   = {
> >    latlon     = FALSE;
> >    raw        = FALSE;
> >    diff       = FALSE;
> >    climo      = FALSE;
> >    weight     = FALSE;
> >    nbrhd      = FALSE;
> >    fourier    = FALSE;
> >    gradient   = FALSE;
> >    apply_mask = FALSE;
> > }
> >
> >
> >
//////////////////////////////////////////////////////////////////////
> > //////////
> >
> > grid_weight_flag = NONE;
> > tmp_dir          = "/tmp";
> > output_prefix    = "GALWEM1";
> > version          = "V8.1";
> >
> > Thanks,
> > Matt
> >
> > Matthew C. Sittel, Meteorologist
> > University Corporation for Atmospheric Research
> > Matthew.Sittel.ctr at us.af.mil
> >
> >
> >
>
>
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2 output from grid_stat
From: SITTEL, MATTHEW C CTR USAF AFWA 16 WS/WXN
Time: Mon Aug 26 05:43:53 2019

That change fixed the problem.  Thanks for your help!

Matt

-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Friday, August 23, 2019 12:19 PM
To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL
<matthew.sittel.ctr at us.af.mil>
Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2
output from grid_stat

Hi Matt and David,

Take a look at this setting in the Grid-Stat config file:
mask = {
   grid = [ "FULL" ];
   poly = [ "" ];
}

Please try changing it to the following:
mask = {
   grid = [ "FULL" ];
   poly = [ ];
}

And see if that makes the duplicates go away.  If this does fix that
behavior, we'll need to run this situation through the debugger and
fix the behavior.  Instead of producing duplicate output lines, it
*should* print an error message stating that it can't parse an empty
string as the name of the polyline file.

Thanks,
John

On Fri, Aug 23, 2019 at 10:42 AM SITTEL, MATTHEW C CTR USAF AFWA 16
WS/WXN via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664 >
>
> FTP doesn't work for us; we can't get through the Air Force
firewall.
> We have to use AMRDEC to transfer files.  It's a bit of a hassle,
> which is why I sent what we did, in hopes that the sample output
> showed the nature of the problem, and the config file contains the
> exact settings I used.  I'm hoping that illustrates the problem.
>
> The SL1L2 output will be used to calculate GO Index in
stat_analysis.
> Duplicate records slow down the process and potentially contaminate
> the results if this repetitive output isn't consistent from run to
run.
>
> Matt
>
> -----Original Message-----
> From: David Fillmore via RT <met_help at ucar.edu>
> Sent: Friday, August 23, 2019 11:24 AM
> To: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL <
> matthew.sittel.ctr at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91664] SL1L2 output
> from grid_stat
>
> Hi Matt -
>
> We will have a look.
> Are you able to provide sample data so that we can replicate the
issue
> if necessary?
> See ftp upload instructions here:
>
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
hel
> p-desk
>
>
> thanks,
> David
>
> On Fri, Aug 23, 2019 at 9:40 AM SITTEL, MATTHEW C CTR USAF AFWA 16
> WS/WXN via RT <met_help at ucar.edu> wrote:
>
> >
> > Fri Aug 23 09:40:26 2019: Request 91664 was acted upon.
> > Transaction: Ticket created by matthew.sittel.ctr at us.af.mil
> >        Queue: met_help
> >      Subject: SL1L2 output from grid_stat
> >        Owner: Nobody
> >   Requestors: matthew.sittel.ctr at us.af.mil
> >       Status: new
> >  Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91664
> > >
> >
> >
> > Hi,
> >
> > I'm running grid_stat in version 8.1 and trying to figure out why
I
> > get duplicate records when I output SL1L2.  Every line is written
> > out twice, and the pair are identical for all fields.  Here is
some
> > sample output with the repeated lines:
> >
> > VERSION MODEL   DESC FCST_LEAD FCST_VALID_BEG  FCST_VALID_END
OBS_LEAD
> > OBS_VALID_BEG   OBS_VALID_END   FCST_VAR FCST_UNITS FCST_LEV
OBS_VAR OBS_
> > UNITS OBS_LEV OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH
> > OBS_THRESH COV_THRESH ALPHA LINE_TYPE
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> > 488.4646        58487.39199       58489.87999  0.43695
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          P400     TMP
K
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    241.48153    241.48947
58
> > 488.4646        58487.39199       58489.87999  0.43695
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
> >       Z2      ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> > 472.23732       77490.92909       77455.03569  0.71251
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 TMP      K          Z2       TMP
K
> >       Z2      ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240    277.3807     277.30876
77
> > 472.23732       77490.92909       77455.03569  0.71251
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> > 400.5502     50755397.41531    50737446.60268  4.27849
> > V8.1    GALWEM1 NA   120000    20190401_120000 20190401_120000
000000
> >  20190401_120000 20190401_120000 HGT      gpm        P400     HGT
gpm
> >       P400    ANALYS FULL    NEAREST     1           NA
NA
> >  NA         NA    SL1L2     1038240   7109.98602   7108.65435
50746
> > 400.5502     50755397.41531    50737446.60268  4.27849
> >
> > I'll post my config file below.  I don't see any obvious reason
why
> > this happens.
> >
> > //
> > // Output model name to be written
> > //
> > model = "GALWEM1";
> >
> > //
> > // Output description to be written
> > // May be set separately in each "obs.field" entry // desc = "NA";
> >
> > //
> > // Output observation type to be written // obtype = "ANALYS";
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Verification grid
> > //
> > regrid = {
> >    to_grid    = NONE;
> >    method     = NEAREST;
> >    width      = 1;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // May be set separately in each "field" entry //
> > censor_thresh    = [];
> > censor_val       = [];
> > cat_thresh       = [];
> > cnt_thresh       = [ NA ];
> > cnt_logic        = UNION;
> > wind_thresh      = [ NA ];
> > wind_logic       = UNION;
> > eclv_points      = 0.05;
> > nc_pairs_var_str = "";
> > rank_corr_flag   = FALSE;
> >
> > //
> > // Forecast and observation fields to be verified // fcst = {
> >
> >    file_type = GRIB2;
> >    field = [
> >       { name = "TMP"; level = "P400"; },
> >       { name = "TMP"; level = "Z2"; },
> >       { name = "HGT"; level = "P400"; },
> >       { name = "DPT"; level = "P850"; },
> >       { name = "DPT"; level = "P700"; },
> >       { name = "DPT"; level = "P400"; },
> >       { name = "DPT"; level = "Z2"; },
> >       { name = "WIND"; level = "Z10"; },
> >       { name = "WIND"; level = "P850"; },
> >       { name = "WIND"; level = "P400"; },
> >       { name = "WIND"; level = "P250"; },
> >       { name = "PRMSL"; level = "Z0"; }
> >    ];
> >
> > }
> >
> > obs = fcst;
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Climatology data
> > //
> > climo_mean = {
> >
> >    file_name = [];
> >    field     = [];
> >
> >    regrid = {
> >       method     = NEAREST;
> >       width      = 1;
> >       vld_thresh = 1.0;
> >       shape      = SQUARE;
> >    }
> >
> >    time_interp_method = DW_MEAN;
> >    match_month        = TRUE;
> >    match_day          = FALSE;
> >    time_step          = 21600;
> > }
> >
> > climo_stdev = climo_mean;
> > climo_stdev = {
> >    file_name = [];
> > }
> >
> > climo_cdf_bins = 1;
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Verification masking regions
> > //
> > mask = {
> >    grid = [ "FULL" ];
> >    poly = [ "" ];
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Confidence interval settings
> > //
> > ci_alpha  = [ 0.05 ];
> >
> > boot = {
> >    interval = PCTILE;
> >    rep_prop = 1.0;
> >    n_rep    = 0;
> >    rng      = "mt19937";
> >    seed     = "";
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Data smoothing methods
> > //
> > interp = {
> >    field      = NONE;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> >
> >    type = [
> >       {
> >          method = NEAREST;
> >          width  = 1;
> >       }
> >    ];
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Neighborhood methods
> > //
> > nbrhd = {
> >    field      = BOTH;
> >    vld_thresh = 1.0;
> >    shape      = SQUARE;
> >    width      = [ 1 ];
> >    cov_thresh = [ >=0.5 ];
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Fourier decomposition
> > //
> > fourier = {
> >    wave_1d_beg = [];
> >    wave_1d_end = [];
> > }
> >
> > //
> > // Gradient statistics
> > // May be set separately in each "obs.field" entry // gradient = {
> >    dx = [ 1 ];
> >    dy = [ 1 ];
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > //
> > // Statistical output types
> > //
> > output_flag = {
> >    fho    = NONE;
> >    ctc    = NONE;
> >    cts    = NONE;
> >    mctc   = NONE;
> >    mcts   = NONE;
> >    cnt    = NONE;
> >    sl1l2  = STAT;
> >    sal1l2 = NONE;
> >    vl1l2  = NONE;
> >    val1l2 = NONE;
> >    pct    = NONE;
> >    pstd   = NONE;
> >    pjc    = NONE;
> >    prc    = NONE;
> >    eclv   = NONE;
> >    nbrctc = NONE;
> >    nbrcts = NONE;
> >    nbrcnt = NONE;
> >    grad   = NONE;
> > }
> >
> > //
> > // NetCDF matched pairs output file
> > //
> > nc_pairs_flag   = {
> >    latlon     = FALSE;
> >    raw        = FALSE;
> >    diff       = FALSE;
> >    climo      = FALSE;
> >    weight     = FALSE;
> >    nbrhd      = FALSE;
> >    fourier    = FALSE;
> >    gradient   = FALSE;
> >    apply_mask = FALSE;
> > }
> >
> >
> >
////////////////////////////////////////////////////////////////////
> > //
> > //////////
> >
> > grid_weight_flag = NONE;
> > tmp_dir          = "/tmp";
> > output_prefix    = "GALWEM1";
> > version          = "V8.1";
> >
> > Thanks,
> > Matt
> >
> > Matthew C. Sittel, Meteorologist
> > University Corporation for Atmospheric Research
> > Matthew.Sittel.ctr at us.af.mil
> >
> >
> >
>
>
>
>



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


More information about the Met_help mailing list