[Met_help] [rt.rap.ucar.edu #89084] History for Scripting in MET

Tue Jul 9 12:06:59 MDT 2019

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

Hi All,

I am wondering if it is possible to write some codes inside MET to get our favorite outputs.

For the beginning, I need to have the intensity of rainfall in each time-step of MTD output. MTD just give me the maximum and percentile intensity for the whole period of a storm.
Another thing I need is the vorticity or rotational velocity.

Do you have any suggestions to get these information if it is not possible to script in  MET ?

Finally, Do you have a script or any pre-written code for plotting the out put information of MTD ?

Sorry for asking too much questions,
Thanks,
Hooman

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

Subject: Scripting in MET
From: John Halley Gotway
Time: Thu Feb 28 15:47:51 2019

Hooman,

The MET package is open source and you're welcome to modify it in
whatever
way you'd like.

I see that you'd specifically like to see the range of intensity
values
within each time step of the mtd objects.  I would suggest processing
the
NetCDF output files from mtd to get this information.  They contain
the raw
forecast and observation fields, along with mask field for the
objects.
You could write a script to read data from these NetCDF files and
extract
the raw values you need.

As for vorticity/rotational velocity, I'm really not sure on that one.
I'll ask Randy Bullock, the developer of mtd to comment on that one.

We do not have support scripts for plotting mtd output, but we have
had a
user contribute some sample python scripts.  Hopefully Randy can
comment on
that as well.

Thanks,
John Halley Gotway

On Tue, Feb 26, 2019 at 1:50 AM Hooman Ayat via RT <met_help at ucar.edu>
wrote:

>
> Tue Feb 26 01:50:08 2019: Request 89084 was acted upon.
> Transaction: Ticket created by h.ayat at student.unsw.edu.au
>        Queue: met_help
>      Subject: Scripting in MET
>        Owner: Nobody
>   Requestors: h.ayat at student.unsw.edu.au
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=89084 >
>
>
> Hi All,
>
> I am wondering if it is possible to write some codes inside MET to
get our
> favorite outputs.
>
> For the beginning, I need to have the intensity of rainfall in each
> time-step of MTD output. MTD just give me the maximum and percentile
> intensity for the whole period of a storm.
> Another thing I need is the vorticity or rotational velocity.
>
> Do you have any suggestions to get these information if it is not
possible
> to script in  MET ?
>
> Finally, Do you have a script or any pre-written code for plotting
the out
> put information of MTD ?
>
> Sorry for asking too much questions,
> Thanks,
> Hooman
>
>

------------------------------------------------
Subject: Scripting in MET
From: Randy Bullock
Time: Fri Mar 01 09:08:52 2019

Hello Hooman -

Thanks for emailing us.

Regarding intensity levels at each time step:  I think that would be a
good
thing add to a future version of MTD, but, as you pointed out, it's
not
currently in the output of the tool.  I suggest you go into the netcdf
file
that MTD produces, and get intensity distributions for each time step
for
each object directly.  That's basically what the netcdf file is for
--- to
enable users to investigate or calculate things for themselves that
the
tool (in this case MTD) may not write out directly.

As for rotational velocity, that can be obtained by using the spatial
axis
angle information that MTD writes out for each time step and for each
object.  Note however that you'll have to do some anti-aliasing on the
axis
angle values.  The angle takes values in the range from -90 degrees to
+90
degrees.  So, for example, if you have a sequence of axis angle values
at
consecutive time steps that look like this: 65, 75, 85, -85, then
you'll
probably want to de-alias that -85 value and change it to 95.

Once that de-aliasing is done, you can fit a model to the values ---
wither
a constant angular velocity model, or some kind of nonlinear growth
and
decay model.   Either of these would give you an angular velocity.  In
the
first case, you would get a time-dependent angular velocity that would
apply to the whole spacetime object.  In the second case, you would
get a
time-dependent angular velocity that would vary over the lifetime of
the
object.

Finally, regarding plotting the output of MTD, we have some code that
we
use internally, but we've decided not to publicly release that.  On
the
plus side, others have used various Python utilites (like matplotlib)
to
create plots of MTD output.  You might consider looking into that.

Hope this helps.

Randy

On Thu, Feb 28, 2019 at 8:56 PM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=89084 >
>
> Hooman,
>
> The MET package is open source and you're welcome to modify it in
whatever
> way you'd like.
>
> I see that you'd specifically like to see the range of intensity
values
> within each time step of the mtd objects.  I would suggest
processing the
> NetCDF output files from mtd to get this information.  They contain
the raw
> forecast and observation fields, along with mask field for the
objects.
> You could write a script to read data from these NetCDF files and
extract
> the raw values you need.
>
> As for vorticity/rotational velocity, I'm really not sure on that
one.
> I'll ask Randy Bullock, the developer of mtd to comment on that one.
>
> We do not have support scripts for plotting mtd output, but we have
had a
> user contribute some sample python scripts.  Hopefully Randy can
comment on
> that as well.
>
> Thanks,
> John Halley Gotway
>
>
> On Tue, Feb 26, 2019 at 1:50 AM Hooman Ayat via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Tue Feb 26 01:50:08 2019: Request 89084 was acted upon.
> > Transaction: Ticket created by h.ayat at student.unsw.edu.au
> >        Queue: met_help
> >      Subject: Scripting in MET
> >        Owner: Nobody
> >   Requestors: h.ayat at student.unsw.edu.au
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=89084 >
> >
> >
> > Hi All,
> >
> > I am wondering if it is possible to write some codes inside MET to
get
> our
> > favorite outputs.
> >
> > For the beginning, I need to have the intensity of rainfall in
each
> > time-step of MTD output. MTD just give me the maximum and
percentile
> > intensity for the whole period of a storm.
> > Another thing I need is the vorticity or rotational velocity.
> >
> > Do you have any suggestions to get these information if it is not
> possible
> > to script in  MET ?
> >
> > Finally, Do you have a script or any pre-written code for plotting
the
> out
> > put information of MTD ?
> >
> > Sorry for asking too much questions,
> > Thanks,
> > Hooman
> >
> >
>
>

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