[Met_help] [rt.rap.ucar.edu #73656] History for Re: Use of WRF-ARW files in MET
Julie Prestopnik via RT
met_help at ucar.edu
Mon Nov 30 15:59:20 MST 2015
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear all,
problem solved. Please let me follow up with another question:
when using netcdf files as input, how can I specify different levels (i.e.
the equivalent of P0-100000 using grib files)? I would like to verify all
vertical levels against all available measures for a given time.
Thanks for your help
Gianluca
On Tue, Oct 6, 2015 at 12:15 PM, Gianluca Meneghello <gianmail at gmail.com>
wrote:
> Dear all,
>
> I am trying to use MET's point_stat on WRF-ARW forecasts. I have two
> questions about how to do that.
>
> First, I understand I have to run the p_interp tool to translate WRF-ARW
> output as an input for MET. Is that still the case for MET5.0 and WRF3.7.1?
>
> Second, I cannot figure out how to write the PointStatConfig file for
> netCDF input.
>
> My current PointStatConfig follows, but point_stat seems to default to the
> default point_stat file instead of reading it.
>
> Any suggestion, or a pointer to an example PointStatConfig file for
> WRF-ARW output is more than welcome.
>
> Thanks in advance for your help and your work
>
> Gianluca
>
>
>
> ////////////////////////////////////////////////////////////////////////////////
> //
> // Point-Stat configuration file.
> //
> // For additional information, see the MET_BASE/config/README file.
> //
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Output model name to be written
> //
> model = "WRF";
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Forecast and observation fields to be verified
> //
> fcst = {
> wind_thresh = [];
> message_type = [];
> sid_exc = [];
>
> fcst_field = [ "UU(0,0,*,*)", "VV(0,0,*,*)" ];
>
> };
>
> obs = {
> wind_thresh = [];
> message_type =
> ["ADPUPA","ADPSFC","AIRCAR","AIRCFT","SATWND","PROFLR","VADWND","SATEMP","ADPSFC","SFCSHP","SFCBOG","SPSSMI","SYNDAT","ERS1DA","GOESND","QKSWND","MSONET","GPSIPW","RASSDA","WDSATR","ASCATW"];
> sid_exc = [];
>
> obs_field = [ "UGRD(0,0,*,*)", "VGRD(0,0,*,*)" ];
>
> };
>
>
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Point observation time window
> //
> obs_window = {
> beg = -5400;
> end = 5400;
> }
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Verification masking regions
> //
> mask = {
> grid = [];
> poly = ["JQN.poly"];
> sid = [];
> };
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Confidence interval settings
> //
> ci_alpha = [ 0.05 ];
>
> boot = {
> interval = PCTILE;
> rep_prop = 1.0;
> n_rep = 1000;
> rng = "mt19937";
> seed = "";
> };
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Interpolation methods
> //
> interp = {
> vld_thresh = 1.0;
>
> type = [
> {
> method = [DW_MEAN, LS_FIT, BILIN];
> width = 5;
> }
> ];
> };
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> //
> // Statistical output types
> //
> output_flag = {
> fho = NONE;
> ctc = NONE;
> cts = NONE;
> mctc = NONE;
> mcts = NONE;
> cnt = BOTH;
> sl1l2 = NONE;
> sal1l2 = NONE;
> vl1l2 = NONE;
> val1l2 = NONE;
> pct = NONE;
> pstd = NONE;
> pjc = NONE;
> prc = NONE;
> mpr = NONE;
> };
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
> obs_quality = [];
> duplicate_flag = NONE;
> rank_corr_flag = TRUE;
> tmp_dir = "/tmp";
> output_prefix = "";
> version = "V5.0";
>
>
> ////////////////////////////////////////////////////////////////////////////////
>
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: Use of WRF-ARW files in MET
From: Julie Prestopnik
Time: Wed Oct 07 10:39:43 2015
Hi Gianluca.
The p_interp utility can be used to convert WRF-ARW forecast into a
flavor
of NetCDF that MET can read, but with some limitations. WRF is run on
a
staggered grid with the wind fields staggered, while the mass points
are on
a regular grid. The p_interp output is indeed on pressure levels, but
the
wind points are still staggered, and MET is unable to read them. So
using
p_interp is not a good choice if you are interested in winds. If you
decide you are not interested in winds, the gridded p_interp output
can be
used in any place in MET that GRIB1 or GRIB2 is used.
The other WRF post-processing alternative is the Unified Post
Processor
(UPP). It's output is in GRIB1 format, which MET fully supports. If
possible, I'd suggest using UPP instead of p_interp to avoid the
staggered
grid limitations of NetCDF. Support for UPP is provided at upp-help
at
ucar.edu.
If you decide to proceed with the p_interp utility, please let us know
and
we can answer your questions.
I hope this helps.
Thanks,
Julie
On Tue, Oct 6, 2015 at 11:29 PM, Gianluca Meneghello via RT <
met_help at ucar.edu> wrote:
>
> Tue Oct 06 17:29:26 2015: Request 73656 was acted upon.
> Transaction: Ticket created by gianmail at gmail.com
> Queue: met_help
> Subject: Re: Use of WRF-ARW files in MET
> Owner: Nobody
> Requestors: gianmail at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73656 >
>
>
> Dear all,
>
> problem solved. Please let me follow up with another question:
>
> when using netcdf files as input, how can I specify different levels
(i.e.
> the equivalent of P0-100000 using grib files)? I would like to
verify all
> vertical levels against all available measures for a given time.
>
> Thanks for your help
>
> Gianluca
>
> On Tue, Oct 6, 2015 at 12:15 PM, Gianluca Meneghello
<gianmail at gmail.com>
> wrote:
>
> > Dear all,
> >
> > I am trying to use MET's point_stat on WRF-ARW forecasts. I have
two
> > questions about how to do that.
> >
> > First, I understand I have to run the p_interp tool to translate
WRF-ARW
> > output as an input for MET. Is that still the case for MET5.0 and
> WRF3.7.1?
> >
> > Second, I cannot figure out how to write the PointStatConfig file
for
> > netCDF input.
> >
> > My current PointStatConfig follows, but point_stat seems to
default to
> the
> > default point_stat file instead of reading it.
> >
> > Any suggestion, or a pointer to an example PointStatConfig file
for
> > WRF-ARW output is more than welcome.
> >
> > Thanks in advance for your help and your work
> >
> > Gianluca
> >
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > //
> > // Point-Stat configuration file.
> > //
> > // For additional information, see the MET_BASE/config/README
file.
> > //
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Output model name to be written
> > //
> > model = "WRF";
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Forecast and observation fields to be verified
> > //
> > fcst = {
> > wind_thresh = [];
> > message_type = [];
> > sid_exc = [];
> >
> > fcst_field = [ "UU(0,0,*,*)", "VV(0,0,*,*)" ];
> >
> > };
> >
> > obs = {
> > wind_thresh = [];
> > message_type =
> >
>
["ADPUPA","ADPSFC","AIRCAR","AIRCFT","SATWND","PROFLR","VADWND","SATEMP","ADPSFC","SFCSHP","SFCBOG","SPSSMI","SYNDAT","ERS1DA","GOESND","QKSWND","MSONET","GPSIPW","RASSDA","WDSATR","ASCATW"];
> > sid_exc = [];
> >
> > obs_field = [ "UGRD(0,0,*,*)", "VGRD(0,0,*,*)" ];
> >
> > };
> >
> >
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Point observation time window
> > //
> > obs_window = {
> > beg = -5400;
> > end = 5400;
> > }
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Verification masking regions
> > //
> > mask = {
> > grid = [];
> > poly = ["JQN.poly"];
> > sid = [];
> > };
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Confidence interval settings
> > //
> > ci_alpha = [ 0.05 ];
> >
> > boot = {
> > interval = PCTILE;
> > rep_prop = 1.0;
> > n_rep = 1000;
> > rng = "mt19937";
> > seed = "";
> > };
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Interpolation methods
> > //
> > interp = {
> > vld_thresh = 1.0;
> >
> > type = [
> > {
> > method = [DW_MEAN, LS_FIT, BILIN];
> > width = 5;
> > }
> > ];
> > };
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > //
> > // Statistical output types
> > //
> > output_flag = {
> > fho = NONE;
> > ctc = NONE;
> > cts = NONE;
> > mctc = NONE;
> > mcts = NONE;
> > cnt = BOTH;
> > sl1l2 = NONE;
> > sal1l2 = NONE;
> > vl1l2 = NONE;
> > val1l2 = NONE;
> > pct = NONE;
> > pstd = NONE;
> > pjc = NONE;
> > prc = NONE;
> > mpr = NONE;
> > };
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
> > obs_quality = [];
> > duplicate_flag = NONE;
> > rank_corr_flag = TRUE;
> > tmp_dir = "/tmp";
> > output_prefix = "";
> > version = "V5.0";
> >
> >
> >
>
////////////////////////////////////////////////////////////////////////////////
> >
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: Re: Use of WRF-ARW files in MET
From: Gianluca Meneghello
Time: Fri Oct 09 11:23:02 2015
Hi Julie,
thanks for your reply. I am indeed interested in winds: I will try UPP
(I
have some difficulties at compiling it right now, but I am working on
it!)
Thanks again
Gianluca
On Wed, Oct 7, 2015 at 9:39 AM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Gianluca.
>
> The p_interp utility can be used to convert WRF-ARW forecast into a
flavor
> of NetCDF that MET can read, but with some limitations. WRF is run
on a
> staggered grid with the wind fields staggered, while the mass points
are on
> a regular grid. The p_interp output is indeed on pressure levels,
but the
> wind points are still staggered, and MET is unable to read them. So
using
> p_interp is not a good choice if you are interested in winds. If
you
> decide you are not interested in winds, the gridded p_interp output
can be
> used in any place in MET that GRIB1 or GRIB2 is used.
>
> The other WRF post-processing alternative is the Unified Post
Processor
> (UPP). It's output is in GRIB1 format, which MET fully supports.
If
> possible, I'd suggest using UPP instead of p_interp to avoid the
staggered
> grid limitations of NetCDF. Support for UPP is provided at upp-help
at
> ucar.edu.
>
> If you decide to proceed with the p_interp utility, please let us
know and
> we can answer your questions.
>
> I hope this helps.
>
> Thanks,
> Julie
>
>
> On Tue, Oct 6, 2015 at 11:29 PM, Gianluca Meneghello via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Tue Oct 06 17:29:26 2015: Request 73656 was acted upon.
> > Transaction: Ticket created by gianmail at gmail.com
> > Queue: met_help
> > Subject: Re: Use of WRF-ARW files in MET
> > Owner: Nobody
> > Requestors: gianmail at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73656 >
> >
> >
> > Dear all,
> >
> > problem solved. Please let me follow up with another question:
> >
> > when using netcdf files as input, how can I specify different
levels
> (i.e.
> > the equivalent of P0-100000 using grib files)? I would like to
verify all
> > vertical levels against all available measures for a given time.
> >
> > Thanks for your help
> >
> > Gianluca
> >
> > On Tue, Oct 6, 2015 at 12:15 PM, Gianluca Meneghello
<gianmail at gmail.com
> >
> > wrote:
> >
> > > Dear all,
> > >
> > > I am trying to use MET's point_stat on WRF-ARW forecasts. I have
two
> > > questions about how to do that.
> > >
> > > First, I understand I have to run the p_interp tool to translate
> WRF-ARW
> > > output as an input for MET. Is that still the case for MET5.0
and
> > WRF3.7.1?
> > >
> > > Second, I cannot figure out how to write the PointStatConfig
file for
> > > netCDF input.
> > >
> > > My current PointStatConfig follows, but point_stat seems to
default to
> > the
> > > default point_stat file instead of reading it.
> > >
> > > Any suggestion, or a pointer to an example PointStatConfig file
for
> > > WRF-ARW output is more than welcome.
> > >
> > > Thanks in advance for your help and your work
> > >
> > > Gianluca
> > >
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > > //
> > > // Point-Stat configuration file.
> > > //
> > > // For additional information, see the MET_BASE/config/README
file.
> > > //
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Output model name to be written
> > > //
> > > model = "WRF";
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Forecast and observation fields to be verified
> > > //
> > > fcst = {
> > > wind_thresh = [];
> > > message_type = [];
> > > sid_exc = [];
> > >
> > > fcst_field = [ "UU(0,0,*,*)", "VV(0,0,*,*)" ];
> > >
> > > };
> > >
> > > obs = {
> > > wind_thresh = [];
> > > message_type =
> > >
> >
>
["ADPUPA","ADPSFC","AIRCAR","AIRCFT","SATWND","PROFLR","VADWND","SATEMP","ADPSFC","SFCSHP","SFCBOG","SPSSMI","SYNDAT","ERS1DA","GOESND","QKSWND","MSONET","GPSIPW","RASSDA","WDSATR","ASCATW"];
> > > sid_exc = [];
> > >
> > > obs_field = [ "UGRD(0,0,*,*)", "VGRD(0,0,*,*)" ];
> > >
> > > };
> > >
> > >
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Point observation time window
> > > //
> > > obs_window = {
> > > beg = -5400;
> > > end = 5400;
> > > }
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Verification masking regions
> > > //
> > > mask = {
> > > grid = [];
> > > poly = ["JQN.poly"];
> > > sid = [];
> > > };
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Confidence interval settings
> > > //
> > > ci_alpha = [ 0.05 ];
> > >
> > > boot = {
> > > interval = PCTILE;
> > > rep_prop = 1.0;
> > > n_rep = 1000;
> > > rng = "mt19937";
> > > seed = "";
> > > };
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Interpolation methods
> > > //
> > > interp = {
> > > vld_thresh = 1.0;
> > >
> > > type = [
> > > {
> > > method = [DW_MEAN, LS_FIT, BILIN];
> > > width = 5;
> > > }
> > > ];
> > > };
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > //
> > > // Statistical output types
> > > //
> > > output_flag = {
> > > fho = NONE;
> > > ctc = NONE;
> > > cts = NONE;
> > > mctc = NONE;
> > > mcts = NONE;
> > > cnt = BOTH;
> > > sl1l2 = NONE;
> > > sal1l2 = NONE;
> > > vl1l2 = NONE;
> > > val1l2 = NONE;
> > > pct = NONE;
> > > pstd = NONE;
> > > pjc = NONE;
> > > prc = NONE;
> > > mpr = NONE;
> > > };
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> > > obs_quality = [];
> > > duplicate_flag = NONE;
> > > rank_corr_flag = TRUE;
> > > tmp_dir = "/tmp";
> > > output_prefix = "";
> > > version = "V5.0";
> > >
> > >
> > >
> >
>
////////////////////////////////////////////////////////////////////////////////
> > >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
------------------------------------------------
More information about the Met_help
mailing list