[Met_help] [rt.rap.ucar.edu #84829] History for obtaining counts for observation threshold levels to create a heatmap?
Julie Prestopnik via RT
met_help at ucar.edu
Wed May 2 10:04:11 MDT 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi,
I'm using ASCAT and GFS surface wind speed data for verification. When I
use point-stat, I have set the cat_thresh to a variety of wind threshold
levels so I can calculate statistics within these bins.
So, I'm trying to make a heatmap of the max wind speed at a certain lat/lon
point, over a period of time, and also, how many times a certain lat/lon
point had a wind speed in a particular threshold.
So, I was wondering if there was a way to keep track of thresholds at a
certain lat/lon point, and get the count for a period of time, and, to get
the max of the wind speed at that particular lat/lon point over a period of
time?
Does MET have the capability to generate such a file that I can plot with
Python?
Thanks!
Roz
--
Rosalyn MacCracken
Support Scientist
Ocean Applications Branch
NOAA/NWS Ocean Prediction Center
NCWCP
5830 University Research Ct
College Park, MD 20740-3818
(p) 301-683-1551
rosalyn.maccracken at noaa.gov
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: Julie Prestopnik
Time: Thu Apr 19 09:25:52 2018
Hello again, Roz. My apologies for the delay in responding.
Unfortunately, John is out of the office this week, and I do not know
the
answers to your questions. Thank you for your patience in advance.
We'll
respond with an answer as soon as we can.
Thanks,
Julie
On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA Affiliate
via RT
<met_help at ucar.edu> wrote:
>
> Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> Queue: met_help
> Subject: obtaining counts for observation threshold levels to
create
> a heatmap?
> Owner: Nobody
> Requestors: rosalyn.maccracken at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
>
>
> Hi,
>
> I'm using ASCAT and GFS surface wind speed data for verification.
When I
> use point-stat, I have set the cat_thresh to a variety of wind
threshold
> levels so I can calculate statistics within these bins.
>
> So, I'm trying to make a heatmap of the max wind speed at a certain
lat/lon
> point, over a period of time, and also, how many times a certain
lat/lon
> point had a wind speed in a particular threshold.
>
> So, I was wondering if there was a way to keep track of thresholds
at a
> certain lat/lon point, and get the count for a period of time, and,
to get
> the max of the wind speed at that particular lat/lon point over a
period of
> time?
>
> Does MET have the capability to generate such a file that I can plot
with
> Python?
>
> Thanks!
>
> Roz
>
> --
> Rosalyn MacCracken
> Support Scientist
>
> Ocean Applications Branch
> NOAA/NWS Ocean Prediction Center
> NCWCP
> 5830 University Research Ct
> College Park, MD 20740-3818
>
> (p) 301-683-1551
> rosalyn.maccracken at noaa.gov
>
>
------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: Rosalyn MacCracken - NOAA Affiliate
Time: Thu Apr 19 09:35:33 2018
Hi Julie,
Thanks for getting back to me. I can wait until John gets back for
the
answer. In the meantime, I'll try and dig through the documentation
to
find the answer.
Roz
On Thu, Apr 19, 2018 at 11:25 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> wrote:
> Hello again, Roz. My apologies for the delay in responding.
>
> Unfortunately, John is out of the office this week, and I do not
know the
> answers to your questions. Thank you for your patience in advance.
We'll
> respond with an answer as soon as we can.
>
> Thanks,
> Julie
>
> On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA Affiliate
via RT
> <met_help at ucar.edu> wrote:
>
> >
> > Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> > Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> > Queue: met_help
> > Subject: obtaining counts for observation threshold levels to
create
> > a heatmap?
> > Owner: Nobody
> > Requestors: rosalyn.maccracken at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
> >
> >
> > Hi,
> >
> > I'm using ASCAT and GFS surface wind speed data for verification.
When I
> > use point-stat, I have set the cat_thresh to a variety of wind
threshold
> > levels so I can calculate statistics within these bins.
> >
> > So, I'm trying to make a heatmap of the max wind speed at a
certain
> lat/lon
> > point, over a period of time, and also, how many times a certain
lat/lon
> > point had a wind speed in a particular threshold.
> >
> > So, I was wondering if there was a way to keep track of thresholds
at a
> > certain lat/lon point, and get the count for a period of time,
and, to
> get
> > the max of the wind speed at that particular lat/lon point over a
period
> of
> > time?
> >
> > Does MET have the capability to generate such a file that I can
plot with
> > Python?
> >
> > Thanks!
> >
> > Roz
> >
> > --
> > Rosalyn MacCracken
> > Support Scientist
> >
> > Ocean Applications Branch
> > NOAA/NWS Ocean Prediction Center
> > NCWCP
> > 5830 University Research Ct
> > College Park, MD 20740-3818
> >
> > (p) 301-683-1551
> > rosalyn.maccracken at noaa.gov
> >
> >
>
>
--
Rosalyn MacCracken
Support Scientist
Ocean Applications Branch
NOAA/NWS Ocean Prediction Center
NCWCP
5830 University Research Ct
College Park, MD 20740-3818
(p) 301-683-1551
rosalyn.maccracken at noaa.gov
------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: John Halley Gotway
Time: Mon Apr 23 12:32:33 2018
Roz,
I see that you're interested in determining the maximum wind speed
value at
particular lat/lon points. It isn't crystal clear to me exactly what
you're trying to do. Do you want to see the maximum gridded
*forecast*
value or are you interested in the lat/lon locations of the ASCAT
observations?
If you're interested in the gridded forecast values, you could use the
Series-Analysis tool to process multiple gridded fields through time.
If you're interested in the lat/lon data of the ASCAT observations,
you
could use the STAT-Analysis tool to process the MPR output lines
through
time.
I'll assume you want to do the latter. I'd recommend running a
*summary*
job, like this:
stat_analysis \
-lookin directory_with_mpr_data \ # List a directory containing
many
.stat files with MPR data in them
-job aggregate_stat \ # Job type is summary
-line_type MPR \ # Input line type = MPR
-column FCST,OBS \ # Process both the FCST and
OBS
columns of data
-fcst_var WIND \ # Only process lines
where
FCST_VAR column = WIND
-fcst_init_beg 20180401 \ # Beginning of init times to
process
-fcst_init_end 20180405 \ # End of init times to
process
-by OBS_SID \ # Run summary job for
each
unique station ID name
-by MODEL,FCST_LEV,VX_MASK,INTERP_MTHD,INTERP_PNTS # Run summary
job for
each unique combo of these columns
-out WIND_by_OBS_SID.stat # Write output to an ascii file
This job will read MPR lines where FCST_VAR = WIND only including the
data
falling between FCST_INIT_BEG and FCST_INIT_END. You could also
specify
particular lead times to include using the "-fcst_lead" option. It'll
organize the input MPR lines into buckets based on the "-by" options
we
listed. Since we used "-by OBS_SID" you'll get separate output for
each
station ID. The maximum value is one of the output of the summary
job.
And it sounds like that's what you're interested in.
Your Python script could read the the output WIND_by_OBS_SID.stat and
plot
the maximum value for each station on a map.
Thanks,
John
On Thu, Apr 19, 2018 at 9:35 AM, Rosalyn MacCracken - NOAA Affiliate
via RT
<met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
>
> Hi Julie,
>
> Thanks for getting back to me. I can wait until John gets back for
the
> answer. In the meantime, I'll try and dig through the documentation
to
> find the answer.
>
> Roz
>
> On Thu, Apr 19, 2018 at 11:25 AM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hello again, Roz. My apologies for the delay in responding.
> >
> > Unfortunately, John is out of the office this week, and I do not
know the
> > answers to your questions. Thank you for your patience in
advance.
> We'll
> > respond with an answer as soon as we can.
> >
> > Thanks,
> > Julie
> >
> > On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA
Affiliate via
> RT
> > <met_help at ucar.edu> wrote:
> >
> > >
> > > Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> > > Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> > > Queue: met_help
> > > Subject: obtaining counts for observation threshold levels
to
> create
> > > a heatmap?
> > > Owner: Nobody
> > > Requestors: rosalyn.maccracken at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829
> >
> > >
> > >
> > > Hi,
> > >
> > > I'm using ASCAT and GFS surface wind speed data for
verification.
> When I
> > > use point-stat, I have set the cat_thresh to a variety of wind
> threshold
> > > levels so I can calculate statistics within these bins.
> > >
> > > So, I'm trying to make a heatmap of the max wind speed at a
certain
> > lat/lon
> > > point, over a period of time, and also, how many times a certain
> lat/lon
> > > point had a wind speed in a particular threshold.
> > >
> > > So, I was wondering if there was a way to keep track of
thresholds at a
> > > certain lat/lon point, and get the count for a period of time,
and, to
> > get
> > > the max of the wind speed at that particular lat/lon point over
a
> period
> > of
> > > time?
> > >
> > > Does MET have the capability to generate such a file that I can
plot
> with
> > > Python?
> > >
> > > Thanks!
> > >
> > > Roz
> > >
> > > --
> > > Rosalyn MacCracken
> > > Support Scientist
> > >
> > > Ocean Applications Branch
> > > NOAA/NWS Ocean Prediction Center
> > > NCWCP
> > > 5830 University Research Ct
> > > College Park, MD 20740-3818
> > >
> > > (p) 301-683-1551
> > > rosalyn.maccracken at noaa.gov
> > >
> > >
> >
> >
>
>
> --
> Rosalyn MacCracken
> Support Scientist
>
> Ocean Applications Branch
> NOAA/NWS Ocean Prediction Center
> NCWCP
> 5830 University Research Ct
> College Park, MD 20740-3818
>
> (p) 301-683-1551
> rosalyn.maccracken at noaa.gov
>
>
------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: Rosalyn MacCracken - NOAA Affiliate
Time: Mon Apr 23 13:12:48 2018
Hi John,
That sounds pretty much like what I want to do. I'm trying to make a
heatmap of :
(1) The max value at each ASCAT "station" (lat/lon point) for a
particular
period of time
and
(2) threshold counts at each station, or a period of time (which I'm
not
sure I can do with stat-analysis).
I'm in the middle of writing a program to do this, but, it's not
working
the right way yet. So, I will try this stat-analysis call and see if
it
gives me what I need.
Thanks!
Roz
On Mon, Apr 23, 2018 at 2:32 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Roz,
>
> I see that you're interested in determining the maximum wind speed
value at
> particular lat/lon points. It isn't crystal clear to me exactly
what
> you're trying to do. Do you want to see the maximum gridded
*forecast*
> value or are you interested in the lat/lon locations of the ASCAT
> observations?
>
> If you're interested in the gridded forecast values, you could use
the
> Series-Analysis tool to process multiple gridded fields through
time.
>
> If you're interested in the lat/lon data of the ASCAT observations,
you
> could use the STAT-Analysis tool to process the MPR output lines
through
> time.
>
> I'll assume you want to do the latter. I'd recommend running a
*summary*
> job, like this:
>
> stat_analysis \
> -lookin directory_with_mpr_data \ # List a directory containing
many
> .stat files with MPR data in them
> -job aggregate_stat \ # Job type is summary
> -line_type MPR \ # Input line type = MPR
> -column FCST,OBS \ # Process both the FCST and
OBS
> columns of data
> -fcst_var WIND \ # Only process lines
where
> FCST_VAR column = WIND
> -fcst_init_beg 20180401 \ # Beginning of init times
to
> process
> -fcst_init_end 20180405 \ # End of init times to
process
> -by OBS_SID \ # Run summary job for
each
> unique station ID name
> -by MODEL,FCST_LEV,VX_MASK,INTERP_MTHD,INTERP_PNTS # Run summary
job
> for
> each unique combo of these columns
> -out WIND_by_OBS_SID.stat # Write output to an ascii file
>
> This job will read MPR lines where FCST_VAR = WIND only including
the data
> falling between FCST_INIT_BEG and FCST_INIT_END. You could also
specify
> particular lead times to include using the "-fcst_lead" option.
It'll
> organize the input MPR lines into buckets based on the "-by" options
we
> listed. Since we used "-by OBS_SID" you'll get separate output for
each
> station ID. The maximum value is one of the output of the summary
job.
> And it sounds like that's what you're interested in.
>
> Your Python script could read the the output WIND_by_OBS_SID.stat
and plot
> the maximum value for each station on a map.
>
> Thanks,
> John
>
>
> On Thu, Apr 19, 2018 at 9:35 AM, Rosalyn MacCracken - NOAA Affiliate
via RT
> <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
> >
> > Hi Julie,
> >
> > Thanks for getting back to me. I can wait until John gets back
for the
> > answer. In the meantime, I'll try and dig through the
documentation to
> > find the answer.
> >
> > Roz
> >
> > On Thu, Apr 19, 2018 at 11:25 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hello again, Roz. My apologies for the delay in responding.
> > >
> > > Unfortunately, John is out of the office this week, and I do not
know
> the
> > > answers to your questions. Thank you for your patience in
advance.
> > We'll
> > > respond with an answer as soon as we can.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA
Affiliate
> via
> > RT
> > > <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> > > > Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> > > > Queue: met_help
> > > > Subject: obtaining counts for observation threshold
levels to
> > create
> > > > a heatmap?
> > > > Owner: Nobody
> > > > Requestors: rosalyn.maccracken at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=84829
> > >
> > > >
> > > >
> > > > Hi,
> > > >
> > > > I'm using ASCAT and GFS surface wind speed data for
verification.
> > When I
> > > > use point-stat, I have set the cat_thresh to a variety of wind
> > threshold
> > > > levels so I can calculate statistics within these bins.
> > > >
> > > > So, I'm trying to make a heatmap of the max wind speed at a
certain
> > > lat/lon
> > > > point, over a period of time, and also, how many times a
certain
> > lat/lon
> > > > point had a wind speed in a particular threshold.
> > > >
> > > > So, I was wondering if there was a way to keep track of
thresholds
> at a
> > > > certain lat/lon point, and get the count for a period of time,
and,
> to
> > > get
> > > > the max of the wind speed at that particular lat/lon point
over a
> > period
> > > of
> > > > time?
> > > >
> > > > Does MET have the capability to generate such a file that I
can plot
> > with
> > > > Python?
> > > >
> > > > Thanks!
> > > >
> > > > Roz
> > > >
> > > > --
> > > > Rosalyn MacCracken
> > > > Support Scientist
> > > >
> > > > Ocean Applications Branch
> > > > NOAA/NWS Ocean Prediction Center
> > > > NCWCP
> > > > 5830 University Research Ct
> > > > College Park, MD 20740-3818
> > > >
> > > > (p) 301-683-1551
> > > > rosalyn.maccracken at noaa.gov
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Rosalyn MacCracken
> > Support Scientist
> >
> > Ocean Applications Branch
> > NOAA/NWS Ocean Prediction Center
> > NCWCP
> > 5830 University Research Ct
> > College Park, MD 20740-3818
> >
> > (p) 301-683-1551
> > rosalyn.maccracken at noaa.gov
> >
> >
>
>
--
Rosalyn MacCracken
Support Scientist
Ocean Applications Branch
NOAA/NWS Ocean Prediction Center
NCWCP
5830 University Research Ct
College Park, MD 20740-3818
(p) 301-683-1551
rosalyn.maccracken at noaa.gov
------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: Rosalyn MacCracken - NOAA Affiliate
Time: Mon Apr 23 13:26:17 2018
Oh, you know, this isn't going to work. We don't use SID. So, it's
the
summary at the Lat/Lon points. And, you know, the lat/lon points
aren't
the same everyday, because of the rotation of the earth. They are
always
slightly moving, so, they aren't exactly the same. So...can I run a
summary job for both lat and lon?
Roz
On Mon, Apr 23, 2018 at 2:32 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Roz,
>
> I see that you're interested in determining the maximum wind speed
value at
> particular lat/lon points. It isn't crystal clear to me exactly
what
> you're trying to do. Do you want to see the maximum gridded
*forecast*
> value or are you interested in the lat/lon locations of the ASCAT
> observations?
>
> If you're interested in the gridded forecast values, you could use
the
> Series-Analysis tool to process multiple gridded fields through
time.
>
> If you're interested in the lat/lon data of the ASCAT observations,
you
> could use the STAT-Analysis tool to process the MPR output lines
through
> time.
>
> I'll assume you want to do the latter. I'd recommend running a
*summary*
> job, like this:
>
> stat_analysis \
> -lookin directory_with_mpr_data \ # List a directory containing
many
> .stat files with MPR data in them
> -job aggregate_stat \ # Job type is summary
> -line_type MPR \ # Input line type = MPR
> -column FCST,OBS \ # Process both the FCST and
OBS
> columns of data
> -fcst_var WIND \ # Only process lines
where
> FCST_VAR column = WIND
> -fcst_init_beg 20180401 \ # Beginning of init times
to
> process
> -fcst_init_end 20180405 \ # End of init times to
process
> -by OBS_SID \ # Run summary job for
each
> unique station ID name
> -by MODEL,FCST_LEV,VX_MASK,INTERP_MTHD,INTERP_PNTS # Run summary
job
> for
> each unique combo of these columns
> -out WIND_by_OBS_SID.stat # Write output to an ascii file
>
> This job will read MPR lines where FCST_VAR = WIND only including
the data
> falling between FCST_INIT_BEG and FCST_INIT_END. You could also
specify
> particular lead times to include using the "-fcst_lead" option.
It'll
> organize the input MPR lines into buckets based on the "-by" options
we
> listed. Since we used "-by OBS_SID" you'll get separate output for
each
> station ID. The maximum value is one of the output of the summary
job.
> And it sounds like that's what you're interested in.
>
> Your Python script could read the the output WIND_by_OBS_SID.stat
and plot
> the maximum value for each station on a map.
>
> Thanks,
> John
>
>
> On Thu, Apr 19, 2018 at 9:35 AM, Rosalyn MacCracken - NOAA Affiliate
via RT
> <met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
> >
> > Hi Julie,
> >
> > Thanks for getting back to me. I can wait until John gets back
for the
> > answer. In the meantime, I'll try and dig through the
documentation to
> > find the answer.
> >
> > Roz
> >
> > On Thu, Apr 19, 2018 at 11:25 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hello again, Roz. My apologies for the delay in responding.
> > >
> > > Unfortunately, John is out of the office this week, and I do not
know
> the
> > > answers to your questions. Thank you for your patience in
advance.
> > We'll
> > > respond with an answer as soon as we can.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA
Affiliate
> via
> > RT
> > > <met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> > > > Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> > > > Queue: met_help
> > > > Subject: obtaining counts for observation threshold
levels to
> > create
> > > > a heatmap?
> > > > Owner: Nobody
> > > > Requestors: rosalyn.maccracken at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=84829
> > >
> > > >
> > > >
> > > > Hi,
> > > >
> > > > I'm using ASCAT and GFS surface wind speed data for
verification.
> > When I
> > > > use point-stat, I have set the cat_thresh to a variety of wind
> > threshold
> > > > levels so I can calculate statistics within these bins.
> > > >
> > > > So, I'm trying to make a heatmap of the max wind speed at a
certain
> > > lat/lon
> > > > point, over a period of time, and also, how many times a
certain
> > lat/lon
> > > > point had a wind speed in a particular threshold.
> > > >
> > > > So, I was wondering if there was a way to keep track of
thresholds
> at a
> > > > certain lat/lon point, and get the count for a period of time,
and,
> to
> > > get
> > > > the max of the wind speed at that particular lat/lon point
over a
> > period
> > > of
> > > > time?
> > > >
> > > > Does MET have the capability to generate such a file that I
can plot
> > with
> > > > Python?
> > > >
> > > > Thanks!
> > > >
> > > > Roz
> > > >
> > > > --
> > > > Rosalyn MacCracken
> > > > Support Scientist
> > > >
> > > > Ocean Applications Branch
> > > > NOAA/NWS Ocean Prediction Center
> > > > NCWCP
> > > > 5830 University Research Ct
> > > > College Park, MD 20740-3818
> > > >
> > > > (p) 301-683-1551
> > > > rosalyn.maccracken at noaa.gov
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Rosalyn MacCracken
> > Support Scientist
> >
> > Ocean Applications Branch
> > NOAA/NWS Ocean Prediction Center
> > NCWCP
> > 5830 University Research Ct
> > College Park, MD 20740-3818
> >
> > (p) 301-683-1551
> > rosalyn.maccracken at noaa.gov
> >
> >
>
>
--
Rosalyn MacCracken
Support Scientist
Ocean Applications Branch
NOAA/NWS Ocean Prediction Center
NCWCP
5830 University Research Ct
College Park, MD 20740-3818
(p) 301-683-1551
rosalyn.maccracken at noaa.gov
------------------------------------------------
Subject: obtaining counts for observation threshold levels to create a heatmap?
From: John Halley Gotway
Time: Mon Apr 23 14:10:12 2018
Roz,
Take a look at the contents of the MPR line type on page 186 of the
user's
guide:
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v7.0.pdf
The OBS_LAT and OBS_LON columns list the lat/lon of the observation
point.
If you run STAT-Analysis using "-by OBS_LAT,OBS_LON" then you'll get
separate output for each unique combination of lat/lon values. Sounds
like
that may be a lot of output!
And regarding (2), yes, you should be able to do that using STAT-
Analysis.
Instead, of running "-job summary", run "-job aggregate_stat
-line_type MPR
-out_line_type CTC -out_thresh gt6 -by OBS_LAT,OBS_LON" (or substitute
whatever threshold you're interested in. That'll generate contingency
table stats for each unique combination of lat/lon value. Then you'd
choose the contingency table stat you'd like to show in your plot.
For
example, BASER is the frequency of the event in the observations and
FMEAN
is the frequency of the event in the forecast values.
Thanks,
John
On Mon, Apr 23, 2018 at 1:26 PM, Rosalyn MacCracken - NOAA Affiliate
via RT
<met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
>
> Oh, you know, this isn't going to work. We don't use SID. So, it's
the
> summary at the Lat/Lon points. And, you know, the lat/lon points
aren't
> the same everyday, because of the rotation of the earth. They are
always
> slightly moving, so, they aren't exactly the same. So...can I run a
> summary job for both lat and lon?
>
> Roz
>
> On Mon, Apr 23, 2018 at 2:32 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Roz,
> >
> > I see that you're interested in determining the maximum wind speed
value
> at
> > particular lat/lon points. It isn't crystal clear to me exactly
what
> > you're trying to do. Do you want to see the maximum gridded
*forecast*
> > value or are you interested in the lat/lon locations of the ASCAT
> > observations?
> >
> > If you're interested in the gridded forecast values, you could use
the
> > Series-Analysis tool to process multiple gridded fields through
time.
> >
> > If you're interested in the lat/lon data of the ASCAT
observations, you
> > could use the STAT-Analysis tool to process the MPR output lines
through
> > time.
> >
> > I'll assume you want to do the latter. I'd recommend running a
*summary*
> > job, like this:
> >
> > stat_analysis \
> > -lookin directory_with_mpr_data \ # List a directory
containing many
> > .stat files with MPR data in them
> > -job aggregate_stat \ # Job type is summary
> > -line_type MPR \ # Input line type =
MPR
> > -column FCST,OBS \ # Process both the FCST
and OBS
> > columns of data
> > -fcst_var WIND \ # Only process lines
where
> > FCST_VAR column = WIND
> > -fcst_init_beg 20180401 \ # Beginning of init
times to
> > process
> > -fcst_init_end 20180405 \ # End of init times to
process
> > -by OBS_SID \ # Run summary job for
each
> > unique station ID name
> > -by MODEL,FCST_LEV,VX_MASK,INTERP_MTHD,INTERP_PNTS # Run
summary job
> > for
> > each unique combo of these columns
> > -out WIND_by_OBS_SID.stat # Write output to an ascii file
> >
> > This job will read MPR lines where FCST_VAR = WIND only including
the
> data
> > falling between FCST_INIT_BEG and FCST_INIT_END. You could also
specify
> > particular lead times to include using the "-fcst_lead" option.
It'll
> > organize the input MPR lines into buckets based on the "-by"
options we
> > listed. Since we used "-by OBS_SID" you'll get separate output
for each
> > station ID. The maximum value is one of the output of the summary
job.
> > And it sounds like that's what you're interested in.
> >
> > Your Python script could read the the output WIND_by_OBS_SID.stat
and
> plot
> > the maximum value for each station on a map.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, Apr 19, 2018 at 9:35 AM, Rosalyn MacCracken - NOAA
Affiliate via
> RT
> > <met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84829 >
> > >
> > > Hi Julie,
> > >
> > > Thanks for getting back to me. I can wait until John gets back
for the
> > > answer. In the meantime, I'll try and dig through the
documentation to
> > > find the answer.
> > >
> > > Roz
> > >
> > > On Thu, Apr 19, 2018 at 11:25 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hello again, Roz. My apologies for the delay in responding.
> > > >
> > > > Unfortunately, John is out of the office this week, and I do
not know
> > the
> > > > answers to your questions. Thank you for your patience in
advance.
> > > We'll
> > > > respond with an answer as soon as we can.
> > > >
> > > > Thanks,
> > > > Julie
> > > >
> > > > On Wed, Apr 18, 2018 at 1:51 PM, Rosalyn MacCracken - NOAA
Affiliate
> > via
> > > RT
> > > > <met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Wed Apr 18 13:51:55 2018: Request 84829 was acted upon.
> > > > > Transaction: Ticket created by rosalyn.maccracken at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: obtaining counts for observation threshold
levels to
> > > create
> > > > > a heatmap?
> > > > > Owner: Nobody
> > > > > Requestors: rosalyn.maccracken at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=84829
> > > >
> > > > >
> > > > >
> > > > > Hi,
> > > > >
> > > > > I'm using ASCAT and GFS surface wind speed data for
verification.
> > > When I
> > > > > use point-stat, I have set the cat_thresh to a variety of
wind
> > > threshold
> > > > > levels so I can calculate statistics within these bins.
> > > > >
> > > > > So, I'm trying to make a heatmap of the max wind speed at a
certain
> > > > lat/lon
> > > > > point, over a period of time, and also, how many times a
certain
> > > lat/lon
> > > > > point had a wind speed in a particular threshold.
> > > > >
> > > > > So, I was wondering if there was a way to keep track of
thresholds
> > at a
> > > > > certain lat/lon point, and get the count for a period of
time, and,
> > to
> > > > get
> > > > > the max of the wind speed at that particular lat/lon point
over a
> > > period
> > > > of
> > > > > time?
> > > > >
> > > > > Does MET have the capability to generate such a file that I
can
> plot
> > > with
> > > > > Python?
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Roz
> > > > >
> > > > > --
> > > > > Rosalyn MacCracken
> > > > > Support Scientist
> > > > >
> > > > > Ocean Applications Branch
> > > > > NOAA/NWS Ocean Prediction Center
> > > > > NCWCP
> > > > > 5830 University Research Ct
> > > > > College Park, MD 20740-3818
> > > > >
> > > > > (p) 301-683-1551
> > > > > rosalyn.maccracken at noaa.gov
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Rosalyn MacCracken
> > > Support Scientist
> > >
> > > Ocean Applications Branch
> > > NOAA/NWS Ocean Prediction Center
> > > NCWCP
> > > 5830 University Research Ct
> > > College Park, MD 20740-3818
> > >
> > > (p) 301-683-1551
> > > rosalyn.maccracken at noaa.gov
> > >
> > >
> >
> >
>
>
> --
> Rosalyn MacCracken
> Support Scientist
>
> Ocean Applications Branch
> NOAA/NWS Ocean Prediction Center
> NCWCP
> 5830 University Research Ct
> College Park, MD 20740-3818
>
> (p) 301-683-1551
> rosalyn.maccracken at noaa.gov
>
>
------------------------------------------------
More information about the Met_help
mailing list