[Met_help] [rt.rap.ucar.edu #71833] History for Running MET with single radar verification

John Halley Gotway via RT met_help at ucar.edu
Mon Aug 17 12:12:41 MDT 2015


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

I am trying to do some verification of my WRF-ARW data using a single radar
as "truth" data. I have put my WRF output through UPP and have a grib
output, but I'm unsure about where to go from here (I'm actually very
unsure if I set the levels in UPP's wrf_cntrl.parm correctly). I have gone
through the MET tutorial, but I still have some doubt about how to use MET
when dealing with radar data.

I gridded the radar data using Radx and it is in a netCDF (CF-compliant)
format but am unsure if MET can handle this format or if I will have to
process this file further.

Could you please tell me if I'm headed in the right direction, if there is
anything I should watch out for, and/or if there are any additional
resources related to this type of verification?
Thank you,
Mariana


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

Subject: Running MET with single radar verification
From: John Halley Gotway
Time: Thu May 07 10:05:24 2015

Mariana,

MET "should" be able to handle CF-compliant NetCDF data, as long as it
includes the projection information for the grid on which it resides.
If
it's on a simple lat/lon grid, then MET can infer the grid by reading
the
lat/lon variables and computing their spacing.

I'd suggest trying it out by using the "plot_data_plane" utility to
visualize the data.  If plot_data_plane can read the data, MODE can
too.
Let's suppose your CF-compliant NetCDF data is stored in a variable
named
"RADAR(lat, lon)".  Just substitute in the correct variable and file
names
in the example below.  Try running:
    met-5.0/bin/plot_data_plane radar_file.nc radar_file.ps
'name="RADAR";
level="(*,*)";'

If that runs without error, you're all set!

If not, there's more work to do.  It'd probably be easiest to have you
send
us a sample data file, following the instructions posted here:
   http://www.dtcenter.org/met/users/support/met_help.php#ftp
We can take a look and let you know what changes are required to get
MET
reading your data.

Thanks,
John Halley Gotway
met_help at ucar.edu

On Thu, May 7, 2015 at 9:22 AM, Mariana Scott via RT
<met_help at ucar.edu>
wrote:

>
> Thu May 07 09:22:51 2015: Request 71833 was acted upon.
> Transaction: Ticket created by mof0001 at uah.edu
>        Queue: met_help
>      Subject: Running MET with single radar verification
>        Owner: Nobody
>   Requestors: mof0001 at uah.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
>
>
> I am trying to do some verification of my WRF-ARW data using a
single radar
> as "truth" data. I have put my WRF output through UPP and have a
grib
> output, but I'm unsure about where to go from here (I'm actually
very
> unsure if I set the levels in UPP's wrf_cntrl.parm correctly). I
have gone
> through the MET tutorial, but I still have some doubt about how to
use MET
> when dealing with radar data.
>
> I gridded the radar data using Radx and it is in a netCDF (CF-
compliant)
> format but am unsure if MET can handle this format or if I will have
to
> process this file further.
>
> Could you please tell me if I'm headed in the right direction, if
there is
> anything I should watch out for, and/or if there are any additional
> resources related to this type of verification?
> Thank you,
> Mariana
>
>

------------------------------------------------
Subject: Running MET with single radar verification
From: Mariana Scott
Time: Wed Jun 03 06:49:19 2015

Hi John,
I apologize for the delay. Here is the error I get when trying to use
plot_data_plane:


DEBUG 1: Opening data file: /home/data/radar/ncf_20140925_101248.nc
WARNING:
WARNING: NcCfFile::open() -> could not extract init time from file
name
WARNING: Using init time of 0
WARNING:
ERROR  :
ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
projection
from information in netCDF file.
ERROR  :


Does this mean that met will not work with the Radx-created netcdf?
I'll keep looking for plot_data_plane options, in case I called it
incorrectly.
Thanks,
-Mariana


On Wed, May 20, 2015 at 10:19 AM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>

------------------------------------------------
Subject: Running MET with single radar verification
From: Julie Prestopnik
Time: Wed Jun 03 09:57:39 2015

Hi Mariana.  Unfortunately, John is out of the office this week.  He
is
likely the best person to answer your question, but I'll see if I can
do
anything to help.

Could you please send us a sample data file, following the
instructions
posted here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp

along with the command that you are running for plot_data_plane.

Thanks,
Julie

On Wed, Jun 3, 2015 at 6:49 AM, Mariana Scott via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
>
> Hi John,
> I apologize for the delay. Here is the error I get when trying to
use
> plot_data_plane:
>
>
> DEBUG 1: Opening data file: /home/data/radar/ncf_20140925_101248.nc
> WARNING:
> WARNING: NcCfFile::open() -> could not extract init time from file
name
> WARNING: Using init time of 0
> WARNING:
> ERROR  :
> ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
projection
> from information in netCDF file.
> ERROR  :
>
>
> Does this mean that met will not work with the Radx-created netcdf?
> I'll keep looking for plot_data_plane options, in case I called it
> incorrectly.
> Thanks,
> -Mariana
>
>
> On Wed, May 20, 2015 at 10:19 AM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > According to our records, your request has been resolved. If you
have any
> > further questions or concerns, please respond to this message.
> >
>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

------------------------------------------------
Subject: Running MET with single radar verification
From: Mariana Scott
Time: Wed Jun 03 10:51:25 2015

Julie,
The sample data file (ncf_20140925_101248.nc) was uploaded to:
/incoming/irap/met_help/scott_data/

The command I used is a variation of (inside met-5.0 directory):
bin/plot_data_plane /path/to/file/ncf_20140925_101248.nc /path/to/met-
5.0/
trial1.ps 'name="ZZ";level="(*,*)";'

Where I attempted changing the "level" variable to level="2",
level="(1,2,*,*)", level="(0,2,*,*)", and level="z0" -- getting the
same
projection information error.
Thanks,
Mariana



On Wed, Jun 3, 2015 at 10:57 AM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Mariana.  Unfortunately, John is out of the office this week.  He
is
> likely the best person to answer your question, but I'll see if I
can do
> anything to help.
>
> Could you please send us a sample data file, following the
instructions
> posted here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> along with the command that you are running for plot_data_plane.
>
> Thanks,
> Julie
>
> On Wed, Jun 3, 2015 at 6:49 AM, Mariana Scott via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
> >
> > Hi John,
> > I apologize for the delay. Here is the error I get when trying to
use
> > plot_data_plane:
> >
> >
> > DEBUG 1: Opening data file:
/home/data/radar/ncf_20140925_101248.nc
> > WARNING:
> > WARNING: NcCfFile::open() -> could not extract init time from file
name
> > WARNING: Using init time of 0
> > WARNING:
> > ERROR  :
> > ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
projection
> > from information in netCDF file.
> > ERROR  :
> >
> >
> > Does this mean that met will not work with the Radx-created
netcdf?
> > I'll keep looking for plot_data_plane options, in case I called it
> > incorrectly.
> > Thanks,
> > -Mariana
> >
> >
> > On Wed, May 20, 2015 at 10:19 AM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > According to our records, your request has been resolved. If you
have
> any
> > > further questions or concerns, please respond to this message.
> > >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>

------------------------------------------------
Subject: Running MET with single radar verification
From: Julie Prestopnik
Time: Wed Jun 03 12:45:31 2015

Hi Mariana.  I have looked at the code and see that it checks for the
attribute "grid_mapping", which your variable "ZZ" has.  I am not
quite
sure why the code does not recognize it, however.

I believe John will be able to provide a solution and more information
to
you next week.

Julie

On Wed, Jun 3, 2015 at 10:51 AM, Mariana Scott via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
>
> Julie,
> The sample data file (ncf_20140925_101248.nc) was uploaded to:
> /incoming/irap/met_help/scott_data/
>
> The command I used is a variation of (inside met-5.0 directory):
> bin/plot_data_plane /path/to/file/ncf_20140925_101248.nc
/path/to/met-5.0/
> trial1.ps 'name="ZZ";level="(*,*)";'
>
> Where I attempted changing the "level" variable to level="2",
> level="(1,2,*,*)", level="(0,2,*,*)", and level="z0" -- getting the
same
> projection information error.
> Thanks,
> Mariana
>
>
>
> On Wed, Jun 3, 2015 at 10:57 AM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Mariana.  Unfortunately, John is out of the office this week.
He is
> > likely the best person to answer your question, but I'll see if I
can do
> > anything to help.
> >
> > Could you please send us a sample data file, following the
instructions
> > posted here:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > along with the command that you are running for plot_data_plane.
> >
> > Thanks,
> > Julie
> >
> > On Wed, Jun 3, 2015 at 6:49 AM, Mariana Scott via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
> > >
> > > Hi John,
> > > I apologize for the delay. Here is the error I get when trying
to use
> > > plot_data_plane:
> > >
> > >
> > > DEBUG 1: Opening data file:
/home/data/radar/ncf_20140925_101248.nc
> > > WARNING:
> > > WARNING: NcCfFile::open() -> could not extract init time from
file name
> > > WARNING: Using init time of 0
> > > WARNING:
> > > ERROR  :
> > > ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
projection
> > > from information in netCDF file.
> > > ERROR  :
> > >
> > >
> > > Does this mean that met will not work with the Radx-created
netcdf?
> > > I'll keep looking for plot_data_plane options, in case I called
it
> > > incorrectly.
> > > Thanks,
> > > -Mariana
> > >
> > >
> > > On Wed, May 20, 2015 at 10:19 AM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > According to our records, your request has been resolved. If
you have
> > any
> > > > further questions or concerns, please respond to this message.
> > > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

------------------------------------------------
Subject: Running MET with single radar verification
From: Mariana Scott
Time: Wed Jun 03 13:00:57 2015

Sounds good, thanks!
-Mariana

On Wed, Jun 3, 2015 at 1:45 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Mariana.  I have looked at the code and see that it checks for
the
> attribute "grid_mapping", which your variable "ZZ" has.  I am not
quite
> sure why the code does not recognize it, however.
>
> I believe John will be able to provide a solution and more
information to
> you next week.
>
> Julie
>
> On Wed, Jun 3, 2015 at 10:51 AM, Mariana Scott via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833 >
> >
> > Julie,
> > The sample data file (ncf_20140925_101248.nc) was uploaded to:
> > /incoming/irap/met_help/scott_data/
> >
> > The command I used is a variation of (inside met-5.0 directory):
> > bin/plot_data_plane /path/to/file/ncf_20140925_101248.nc
> /path/to/met-5.0/
> > trial1.ps 'name="ZZ";level="(*,*)";'
> >
> > Where I attempted changing the "level" variable to level="2",
> > level="(1,2,*,*)", level="(0,2,*,*)", and level="z0" -- getting
the same
> > projection information error.
> > Thanks,
> > Mariana
> >
> >
> >
> > On Wed, Jun 3, 2015 at 10:57 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Mariana.  Unfortunately, John is out of the office this week.
He is
> > > likely the best person to answer your question, but I'll see if
I can
> do
> > > anything to help.
> > >
> > > Could you please send us a sample data file, following the
instructions
> > > posted here:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > along with the command that you are running for plot_data_plane.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Wed, Jun 3, 2015 at 6:49 AM, Mariana Scott via RT <
> met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71833
>
> > > >
> > > > Hi John,
> > > > I apologize for the delay. Here is the error I get when trying
to use
> > > > plot_data_plane:
> > > >
> > > >
> > > > DEBUG 1: Opening data file:
/home/data/radar/ncf_20140925_101248.nc
> > > > WARNING:
> > > > WARNING: NcCfFile::open() -> could not extract init time from
file
> name
> > > > WARNING: Using init time of 0
> > > > WARNING:
> > > > ERROR  :
> > > > ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
> projection
> > > > from information in netCDF file.
> > > > ERROR  :
> > > >
> > > >
> > > > Does this mean that met will not work with the Radx-created
netcdf?
> > > > I'll keep looking for plot_data_plane options, in case I
called it
> > > > incorrectly.
> > > > Thanks,
> > > > -Mariana
> > > >
> > > >
> > > > On Wed, May 20, 2015 at 10:19 AM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > According to our records, your request has been resolved. If
you
> have
> > > any
> > > > > further questions or concerns, please respond to this
message.
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> >
> >
>
>
> --
> 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