[Met_help] [rt.rap.ucar.edu #63414] History for need help to grid-stat against pressure range

John Halley Gotway via RT met_help at ucar.edu
Tue Oct 15 15:59:58 MDT 2013


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

Hi MET_HELP,

I have problem to do grad-stat against pressure level range. my grib files
include TMP from 1000 mb to 100 mb, with interval 25mb in first several
pressure levels. I would like to use grid-stat to do statistic analysis fro
pressure 1000mb to 800mb. I set level as P1000-P800, or P800-1000, it
doesn't work. I uploaded my two input grib files and configure file here.
It is at ftp.rap.ucar.edu/incoming/irap/met_help/zhu. Could you help me to
figure out the problem?

Thank you in advance,

Jiang


**********************************************************
Jiang Zhu, Ph.D.
Geographic Information Network of Alaska
University of Alaska Fairbanks
Phone: 907-474-5689, EMAIL: jiang at gina.alaska.edu
***********************************************************


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

Subject: need help to grid-stat against pressure range
From: John Halley Gotway
Time: Tue Oct 15 12:22:09 2013

Jiang,

The Grid-Stat tool can only be used to extract a single field from
your forecast file and compare it to a single field in your
observation file.  You have 9 different fields between 1000 and 800 mb
and you can't verify them all as a single level like P1000-P800.
Instead, you need to compare the forecast and observation data level-
by-level.  You can do that in a single call to Grid-Stat, and
you'll have statistics computed for each of those 9 levels in the
output.  If you'd like to aggregate those statistics over the 9
vertical levels, you could run the Stat-Analysis tool on the Grid-Stat
output.  I'll list the details about doing these steps below.

But as of METv4.1, you also have another option.  You could use the
Series-Analysis tool to define a height series.  Then for each grid
point in your domain, you could compute statistics over the 9
matched pairs in your height series.  Using this approach you should
be able to see how the model performance varies over you domain.  I'll
also list the details about this below.

Grid-Stat config file settings...

- Specify the exact levels to be verified in config file:
         name       = "TMP";
         level      = ["P1000", "P975", "P950", "P925",
                       "P900",  "P875", "P850", "P825",
                       "P800"];
         cat_thresh = [];
- Turn on the SL1L2 output line:
     sl1l2  = BOTH;

Running Stat-Analysis job to aggregate over all vertical levels:
    METv4.1/bin/stat_analysis -lookin out -job aggregate_stat
-line_type SL1L2 -out_line_type CNT -fcst_var TMP

Running Stat-Analysis job to aggregate over a subset of vertical
levels:
    METv4.1/bin/stat_analysis -lookin out -job aggregate_stat
-line_type SL1L2 -out_line_type CNT -fcst_var TMP \
                              -fcst_lev P1000 -fcst_lev P975 -fcst_lev
P950 -fcst_lev P925 \
                              -fcst_lev P900 -fcst_lev P875 -fcst_lev
P850 -fcst_lev P825  -fcst_lev P800

Running the Series-Analysis tool...

- I've attached the config file I used.

- Here's how I called it on the command line:
    METv4.1/bin/series_analysis \
    -fcst rh-polar_merged_AWIP32.2012110500.grb \
    -obs WRFPRS_CNTL_d01_2012-11-05_00_00_00.grb \
    -config SeriesAnalysisConfig \
    -out temp_height_series_P1000_P800.nc

- And here's how I created a plot of the resulting RMSE field:
    METv4.1/bin/plot_data_plane \
    temp_height_series_P1000_P800.nc \
    temp_height_series_P1000_P800.ps \
    'name="series_cnt_RMSE"; level="(*,*)";'

- I've attached the resulting image.

Hope that helps.

Thanks,
John Halley Gotway



On 10/14/2013 04:33 PM, Jiang Zhu via RT wrote:
>
> Mon Oct 14 16:33:03 2013: Request 63414 was acted upon.
> Transaction: Ticket created by jiang at gina.alaska.edu
>         Queue: met_help
>       Subject: need help to grid-stat against pressure range
>         Owner: Nobody
>    Requestors: jiang at gina.alaska.edu
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=63414 >
>
>
> Hi MET_HELP,
>
> I have problem to do grad-stat against pressure level range. my grib
files
> include TMP from 1000 mb to 100 mb, with interval 25mb in first
several
> pressure levels. I would like to use grid-stat to do statistic
analysis fro
> pressure 1000mb to 800mb. I set level as P1000-P800, or P800-1000,
it
> doesn't work. I uploaded my two input grib files and configure file
here.
> It is at ftp.rap.ucar.edu/incoming/irap/met_help/zhu. Could you help
me to
> figure out the problem?
>
> Thank you in advance,
>
> Jiang
>
>
> **********************************************************
> Jiang Zhu, Ph.D.
> Geographic Information Network of Alaska
> University of Alaska Fairbanks
> Phone: 907-474-5689, EMAIL: jiang at gina.alaska.edu
> ***********************************************************
>

------------------------------------------------
Subject: need help to grid-stat against pressure range
From: John Halley Gotway
Time: Tue Oct 15 12:22:09 2013

////////////////////////////////////////////////////////////////////////////////
//
// Series-Analysis configuration file.
//
// For additional information, see the MET_BASE/data/config/README
file.
//
////////////////////////////////////////////////////////////////////////////////

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

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

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

   cat_thresh  = [ >0.0, >=5.0 ];

   field = [
      {
        name  = "TMP";
        level = ["P1000", "P975", "P950", "P925",
                 "P900",  "P875", "P850", "P825",
                 "P800"];
      }
   ];

};
obs = fcst;

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

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

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

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

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

//
// Number of grid points to be processed concurrently.  Set smaller to
use
// less memory but increase the number of passes through the data.
//
block_size = 1024;

//
// Ratio of valid matched pairs to compute statistics for a grid point
//
vld_thresh = 1.0;

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

//
// Statistical output types
//
output_stats = {
   fho    = [];
   ctc    = [];
   cts    = [];
   mctc   = [];
   mcts   = [];
   cnt    = [ "RMSE", "FBAR", "OBAR" ];
   sl1l2  = [];
   pct    = [];
   pstd   = [];
   pjc    = [];
   prc    = [];
};

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

rank_corr_flag = TRUE;
tmp_dir        = "/tmp";
version        = "V4.1";

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

------------------------------------------------
Subject: need help to grid-stat against pressure range
From: Jiang Zhu
Time: Tue Oct 15 15:56:30 2013

Hi John,

What you provided are very helpful for me. I will follow your guidance
to
do aggregation of multiple levels analysis later.

Thank you very much!

Jiang


On Tue, Oct 15, 2013 at 10:22 AM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Jiang,
>
> The Grid-Stat tool can only be used to extract a single field from
your
> forecast file and compare it to a single field in your observation
file.
>  You have 9 different fields between 1000 and 800 mb
> and you can't verify them all as a single level like P1000-P800.
Instead,
> you need to compare the forecast and observation data level-by-
level.  You
> can do that in a single call to Grid-Stat, and
> you'll have statistics computed for each of those 9 levels in the
output.
>  If you'd like to aggregate those statistics over the 9 vertical
levels,
> you could run the Stat-Analysis tool on the Grid-Stat
> output.  I'll list the details about doing these steps below.
>
> But as of METv4.1, you also have another option.  You could use the
> Series-Analysis tool to define a height series.  Then for each grid
point
> in your domain, you could compute statistics over the 9
> matched pairs in your height series.  Using this approach you should
be
> able to see how the model performance varies over you domain.  I'll
also
> list the details about this below.
>
> Grid-Stat config file settings...
>
> - Specify the exact levels to be verified in config file:
>          name       = "TMP";
>          level      = ["P1000", "P975", "P950", "P925",
>                        "P900",  "P875", "P850", "P825",
>                        "P800"];
>          cat_thresh = [];
> - Turn on the SL1L2 output line:
>      sl1l2  = BOTH;
>
> Running Stat-Analysis job to aggregate over all vertical levels:
>     METv4.1/bin/stat_analysis -lookin out -job aggregate_stat
-line_type
> SL1L2 -out_line_type CNT -fcst_var TMP
>
> Running Stat-Analysis job to aggregate over a subset of vertical
levels:
>     METv4.1/bin/stat_analysis -lookin out -job aggregate_stat
-line_type
> SL1L2 -out_line_type CNT -fcst_var TMP \
>                               -fcst_lev P1000 -fcst_lev P975
-fcst_lev
> P950 -fcst_lev P925 \
>                               -fcst_lev P900 -fcst_lev P875
-fcst_lev P850
> -fcst_lev P825  -fcst_lev P800
>
> Running the Series-Analysis tool...
>
> - I've attached the config file I used.
>
> - Here's how I called it on the command line:
>     METv4.1/bin/series_analysis \
>     -fcst rh-polar_merged_AWIP32.2012110500.grb \
>     -obs WRFPRS_CNTL_d01_2012-11-05_00_00_00.grb \
>     -config SeriesAnalysisConfig \
>     -out temp_height_series_P1000_P800.nc
>
> - And here's how I created a plot of the resulting RMSE field:
>     METv4.1/bin/plot_data_plane \
>     temp_height_series_P1000_P800.nc \
>     temp_height_series_P1000_P800.ps \
>     'name="series_cnt_RMSE"; level="(*,*)";'
>
> - I've attached the resulting image.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
>
>
>
> On 10/14/2013 04:33 PM, Jiang Zhu via RT wrote:
> >
> > Mon Oct 14 16:33:03 2013: Request 63414 was acted upon.
> > Transaction: Ticket created by jiang at gina.alaska.edu
> >         Queue: met_help
> >       Subject: need help to grid-stat against pressure range
> >         Owner: Nobody
> >    Requestors: jiang at gina.alaska.edu
> >        Status: new
> >   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=63414 >
> >
> >
> > Hi MET_HELP,
> >
> > I have problem to do grad-stat against pressure level range. my
grib
> files
> > include TMP from 1000 mb to 100 mb, with interval 25mb in first
several
> > pressure levels. I would like to use grid-stat to do statistic
analysis
> fro
> > pressure 1000mb to 800mb. I set level as P1000-P800, or P800-1000,
it
> > doesn't work. I uploaded my two input grib files and configure
file here.
> > It is at ftp.rap.ucar.edu/incoming/irap/met_help/zhu. Could you
help me
> to
> > figure out the problem?
> >
> > Thank you in advance,
> >
> > Jiang
> >
> >
> > **********************************************************
> > Jiang Zhu, Ph.D.
> > Geographic Information Network of Alaska
> > University of Alaska Fairbanks
> > Phone: 907-474-5689, EMAIL: jiang at gina.alaska.edu
> > ***********************************************************
> >
>
>


--
**********************************************************
Jiang Zhu, Ph.D.
Geographic Information Network of Alaska
University of Alaska Fairbanks
Phone: 907-474-5689, EMAIL: jiang at gina.alaska.edu
***********************************************************

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


More information about the Met_help mailing list