[Met_help] [rt.rap.ucar.edu #99455] History for Display MET products
Minna Win via RT
met_help at ucar.edu
Thu Apr 22 14:44:35 MDT 2021
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I have several questions:
1. How could I use METView in CHEYENNE? I saw there is a Metview
/glade/p/ral/jntp/MET/METviewer
2. How to use the python plotting scripts you provided? I need to plot time
series of RMSE/BIAS etc.
cheyenne:/glade/p/ral/jntp/MET/METplus/METplus-4.0.0-beta3/ush/plotting_scripts>ls
plot_date_by_level.py plot_lead_by_level.py plot_threshold_by_lead.py
plot_lead_average.py plot_stat_by_level.py plot_time_series.py
plot_lead_by_date.py plot_threshold_average.py plot_util.py
3. Could I use Fortran code to rewrite/recompute the MET output of
point_stat_000000L_20150715_000000V.stat?
Do you have a simple subroutine to read this file?
4. What do these thresh means? The units of Obs and Fcst should be the
same, or could be different? such as obs is K, FCst can be degC.
FCST_VAR1_THRESH = <=273, >273 !for temperature
OBS_VAR1_THRESH = <=273, >273
FCST_VAR2_THRESH = >=5 !for wind
Thanks!
Zhuming
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Display MET products
From: Minna Win
Time: Wed Apr 07 15:42:54 2021
Hello,
It looks like you have a lot of different questions. I will need to
get
assistance for answering these. Thank you for your patience while I
get
answers/guidance from the appropriate sources.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Apr 7, 2021 at 2:53 PM George McCabe via RT
<met_help at ucar.edu>
wrote:
>
> Wed Apr 07 14:53:23 2021: Request 99455 was acted upon.
> Transaction: Given to minnawin (Minna Win) by mccabe
> Queue: met_help
> Subject: Display MET products
> Owner: minnawin
> Requestors: zhuming at ucar.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455 >
>
>
> This transaction appears to have no content
>
------------------------------------------------
Subject: Display MET products
From: Minna Win
Time: Wed Apr 07 16:02:15 2021
Hello Zhuming,
For question #4, Thresh is used for setting threshold criteria.
According
to George McCabe, you can indeed have one set of units for Obs, and a
different set of units for Fcst. Typically, units for both Fcst and
Obs
tend to be the same.
For question #3, I need to understand *why* you want to recompute
point-stat. The .stat file is a white-space separated text file. You
could read it into a pandas dataframe as a csv file (whitespace
separated
rather than comma-separated) and do any manipulations on the
dataframe.
The MET code for point-stat is written in C/C++.
Regards,
Minna
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Apr 7, 2021 at 11:56 AM zhuming at ucar.edu via RT
<met_help at ucar.edu>
wrote:
>
> Wed Apr 07 11:56:31 2021: Request 99455 was acted upon.
> Transaction: Ticket created by zhuming at ucar.edu
> Queue: met_help
> Subject: Display MET products
> Owner: Nobody
> Requestors: zhuming at ucar.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455 >
>
>
> Hello,
>
> I have several questions:
>
> 1. How could I use METView in CHEYENNE? I saw there is a Metview
>
> /glade/p/ral/jntp/MET/METviewer
>
> 2. How to use the python plotting scripts you provided? I need to
plot time
> series of RMSE/BIAS etc.
>
>
> cheyenne:/glade/p/ral/jntp/MET/METplus/METplus-4.0.0-
beta3/ush/plotting_scripts>ls
>
> plot_date_by_level.py plot_lead_by_level.py
plot_threshold_by_lead.py
>
> plot_lead_average.py plot_stat_by_level.py
plot_time_series.py
>
> plot_lead_by_date.py plot_threshold_average.py plot_util.py
>
> 3. Could I use Fortran code to rewrite/recompute the MET output of
> point_stat_000000L_20150715_000000V.stat?
> Do you have a simple subroutine to read this file?
>
> 4. What do these thresh means? The units of Obs and Fcst should be
the
> same, or could be different? such as obs is K, FCst can be degC.
>
> FCST_VAR1_THRESH = <=273, >273 !for temperature
>
> OBS_VAR1_THRESH = <=273, >273
>
> FCST_VAR2_THRESH = >=5 !for wind
>
>
> Thanks!
>
> Zhuming
>
>
------------------------------------------------
Subject: Display MET products
From: zhuming at ucar.edu
Time: Wed Apr 07 17:39:32 2021
< For question #4, Thresh is used for setting threshold criteria.
so "FCST_VAR2_THRESH = >=5 " means remove all the data with the value
<5 ?
On Wed, Apr 7, 2021 at 4:02 PM Minna Win via RT <met_help at ucar.edu>
wrote:
> Hello Zhuming,
>
> For question #4, Thresh is used for setting threshold criteria.
According
> to George McCabe, you can indeed have one set of units for Obs, and
a
> different set of units for Fcst. Typically, units for both Fcst and
Obs
> tend to be the same.
>
> For question #3, I need to understand *why* you want to recompute
> point-stat. The .stat file is a white-space separated text file.
You
> could read it into a pandas dataframe as a csv file (whitespace
separated
> rather than comma-separated) and do any manipulations on the
dataframe.
> The MET code for point-stat is written in C/C++.
>
> Regards,
> Minna
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
> ---------------
> Pronouns: she/her
>
>
> On Wed, Apr 7, 2021 at 11:56 AM zhuming at ucar.edu via RT
<met_help at ucar.edu
> >
> wrote:
>
> >
> > Wed Apr 07 11:56:31 2021: Request 99455 was acted upon.
> > Transaction: Ticket created by zhuming at ucar.edu
> > Queue: met_help
> > Subject: Display MET products
> > Owner: Nobody
> > Requestors: zhuming at ucar.edu
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455 >
> >
> >
> > Hello,
> >
> > I have several questions:
> >
> > 1. How could I use METView in CHEYENNE? I saw there is a Metview
> >
> > /glade/p/ral/jntp/MET/METviewer
> >
> > 2. How to use the python plotting scripts you provided? I need to
plot
> time
> > series of RMSE/BIAS etc.
> >
> >
> >
> cheyenne:/glade/p/ral/jntp/MET/METplus/METplus-4.0.0-
beta3/ush/plotting_scripts>ls
> >
> > plot_date_by_level.py plot_lead_by_level.py
> plot_threshold_by_lead.py
> >
> > plot_lead_average.py plot_stat_by_level.py
plot_time_series.py
> >
> > plot_lead_by_date.py plot_threshold_average.py plot_util.py
> >
> > 3. Could I use Fortran code to rewrite/recompute the MET output
of
> > point_stat_000000L_20150715_000000V.stat?
> > Do you have a simple subroutine to read this file?
> >
> > 4. What do these thresh means? The units of Obs and Fcst should be
the
> > same, or could be different? such as obs is K, FCst can be degC.
> >
> > FCST_VAR1_THRESH = <=273, >273 !for temperature
> >
> > OBS_VAR1_THRESH = <=273, >273
> >
> > FCST_VAR2_THRESH = >=5 !for wind
> >
> >
> > Thanks!
> >
> > Zhuming
> >
> >
>
>
------------------------------------------------
Subject: Display MET products
From: Minna Win
Time: Wed Apr 07 19:40:15 2021
Not remove, but filter out the data that doesn't fall within the
threshold
criteria. In your example, you are only interested in FCST values
greater
than or equal to 5. Does that make sense?
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Apr 7, 2021 at 5:39 PM zhuming at ucar.edu via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455 >
>
> < For question #4, Thresh is used for setting threshold criteria.
>
> so "FCST_VAR2_THRESH = >=5 " means remove all the data with the
value <5 ?
>
> On Wed, Apr 7, 2021 at 4:02 PM Minna Win via RT <met_help at ucar.edu>
wrote:
>
> > Hello Zhuming,
> >
> > For question #4, Thresh is used for setting threshold criteria.
> According
> > to George McCabe, you can indeed have one set of units for Obs,
and a
> > different set of units for Fcst. Typically, units for both Fcst
and Obs
> > tend to be the same.
> >
> > For question #3, I need to understand *why* you want to recompute
> > point-stat. The .stat file is a white-space separated text file.
You
> > could read it into a pandas dataframe as a csv file (whitespace
separated
> > rather than comma-separated) and do any manipulations on the
dataframe.
> > The MET code for point-stat is written in C/C++.
> >
> > Regards,
> > Minna
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> > ---------------
> > Pronouns: she/her
> >
> >
> > On Wed, Apr 7, 2021 at 11:56 AM zhuming at ucar.edu via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > >
> > > Wed Apr 07 11:56:31 2021: Request 99455 was acted upon.
> > > Transaction: Ticket created by zhuming at ucar.edu
> > > Queue: met_help
> > > Subject: Display MET products
> > > Owner: Nobody
> > > Requestors: zhuming at ucar.edu
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455
> >
> > >
> > >
> > > Hello,
> > >
> > > I have several questions:
> > >
> > > 1. How could I use METView in CHEYENNE? I saw there is a Metview
> > >
> > > /glade/p/ral/jntp/MET/METviewer
> > >
> > > 2. How to use the python plotting scripts you provided? I need
to plot
> > time
> > > series of RMSE/BIAS etc.
> > >
> > >
> > >
> >
> cheyenne:/glade/p/ral/jntp/MET/METplus/METplus-4.0.0-
beta3/ush/plotting_scripts>ls
> > >
> > > plot_date_by_level.py plot_lead_by_level.py
> > plot_threshold_by_lead.py
> > >
> > > plot_lead_average.py plot_stat_by_level.py
plot_time_series.py
> > >
> > > plot_lead_by_date.py plot_threshold_average.py plot_util.py
> > >
> > > 3. Could I use Fortran code to rewrite/recompute the MET output
of
> > > point_stat_000000L_20150715_000000V.stat?
> > > Do you have a simple subroutine to read this file?
> > >
> > > 4. What do these thresh means? The units of Obs and Fcst should
be the
> > > same, or could be different? such as obs is K, FCst can be degC.
> > >
> > > FCST_VAR1_THRESH = <=273, >273 !for temperature
> > >
> > > OBS_VAR1_THRESH = <=273, >273
> > >
> > > FCST_VAR2_THRESH = >=5 !for wind
> > >
> > >
> > > Thanks!
> > >
> > > Zhuming
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Display MET products
From: Minna Win
Time: Thu Apr 08 09:28:43 2021
Hi Zhuming,
Responses for questions #1 and #2:
-----------------------------------------------
I asked around, and it turns out that the METviewer instance installed
on
'cheyenne' was installed sometime last year for another project, but
can be
used by anyone, but is unsupported. The data loaded in that instance
was
used for the Skymet project.
Here is a link to the METviewer documentation, on 'Read the Docs':
https://metviewer.readthedocs.io/en/latest/Users_Guide/index.html
You can use METviewer to generate the RMSE/Bias plots you were
interested
in, here is an example output from METviewer:
[image: plot_example.png]
This example has the following details (provided by Tatiana):
METviewer can calculate/display:
BCMSE - Bias-corrected mean-squared root
BCRMSE - Bias-corrected root mean square error
The plots can be built using pre-calculated by MET values or using
sl1l2
data.
Here is an example from dakota
I hope that answers your questions.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Apr 7, 2021 at 11:56 AM zhuming at ucar.edu via RT
<met_help at ucar.edu>
wrote:
>
> Wed Apr 07 11:56:31 2021: Request 99455 was acted upon.
> Transaction: Ticket created by zhuming at ucar.edu
> Queue: met_help
> Subject: Display MET products
> Owner: Nobody
> Requestors: zhuming at ucar.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=99455 >
>
>
> Hello,
>
> I have several questions:
>
> 1. How could I use METView in CHEYENNE? I saw there is a Metview
>
> /glade/p/ral/jntp/MET/METviewer
>
> 2. How to use the python plotting scripts you provided? I need to
plot time
> series of RMSE/BIAS etc.
>
>
> cheyenne:/glade/p/ral/jntp/MET/METplus/METplus-4.0.0-
beta3/ush/plotting_scripts>ls
>
> plot_date_by_level.py plot_lead_by_level.py
plot_threshold_by_lead.py
>
> plot_lead_average.py plot_stat_by_level.py
plot_time_series.py
>
> plot_lead_by_date.py plot_threshold_average.py plot_util.py
>
> 3. Could I use Fortran code to rewrite/recompute the MET output of
> point_stat_000000L_20150715_000000V.stat?
> Do you have a simple subroutine to read this file?
>
> 4. What do these thresh means? The units of Obs and Fcst should be
the
> same, or could be different? such as obs is K, FCst can be degC.
>
> FCST_VAR1_THRESH = <=273, >273 !for temperature
>
> OBS_VAR1_THRESH = <=273, >273
>
> FCST_VAR2_THRESH = >=5 !for wind
>
>
> Thanks!
>
> Zhuming
>
>
------------------------------------------------
More information about the Met_help
mailing list