[Met_help] [rt.rap.ucar.edu #79208] History for MODE error: Couldn't figure out projection from information in netCDF file.

Julie Prestopnik via RT met_help at ucar.edu
Thu Mar 2 13:38:44 MST 2017


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

Hi,

Upon running MODE, I encounter the error: Couldn't figure out projection
from information in netCDF file.

I assume because it says "in netCDF file" that it's referring to my obs
file?  In this case, that is a radar netCDF processed through radx called
ncf_20141127_010250.nc.

I'm pretty sure the projection is Mercator.  Can I just hard code that in
somewhere, since having MODE search for the projection doesn't work?

Thanks,
Ryan

-- 
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu


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

Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Jan 19 12:54:15 2017

Hi Ryan.

I see that you're having trouble with a NetCDF file using MODE.  Thank
you
for letting us know that you are using a radar NetCDF file processed
through radx.

Currently, MET has the ability to read three different flavors of
NetCDF...
 - The format used to write intermediate files in MET (i.e. the output
of
pcp_combine)
 - NetCDF files following the Climate Forecast conventions (i.e.
CF-compliant NetCDF)
 - The output of the pinterp or wrf_interp utility

I suspect that your file is not in one of those formats.

I have attached an example of a NetCDF file (out.nc) with a Mercator
projection that MET can read.  If you can get your file in that
format,
MODE should be able to read it.

I hope that helps!

Julie

On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> Transaction: Ticket created by rconne01 at gmail.com
>        Queue: met_help
>      Subject: MODE error: Couldn't figure out projection from
information
> in netCDF file.
>        Owner: Nobody
>   Requestors: rconne01 at gmail.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
>
> Hi,
>
> Upon running MODE, I encounter the error: Couldn't figure out
projection
> from information in netCDF file.
>
> I assume because it says "in netCDF file" that it's referring to my
obs
> file?  In this case, that is a radar netCDF processed through radx
called
> ncf_20141127_010250.nc.
>
> I'm pretty sure the projection is Mercator.  Can I just hard code
that in
> somewhere, since having MODE search for the projection doesn't work?
>
> Thanks,
> Ryan
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Thu Jan 19 16:25:46 2017

Hi Julie,

Is there a way to make radar netCDF files that MODE will read out of
the
native .gz files?  I like what I've been doing so far with radx since
I
want to do my analysis on constant height surfaces, rather than worry
about
elevation angle, and radx allows me to do that.

Unrelated to that issue, when I run with your out.nc file, it tells me
that
the GRIB1 field abbreviation isn't recognized.  It's for radar
reflectivity
AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc., and none
of
those is right, apparently.  Any idea what it is or how I can find
out?  (I
know for sure that that field was computed when running the UPP, as I
can
see the field REFD in my gdinfo in GEMPAK.)

Thanks,
Ryan

On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Ryan.
>
> I see that you're having trouble with a NetCDF file using MODE.
Thank you
> for letting us know that you are using a radar NetCDF file processed
> through radx.
>
> Currently, MET has the ability to read three different flavors of
NetCDF...
>  - The format used to write intermediate files in MET (i.e. the
output of
> pcp_combine)
>  - NetCDF files following the Climate Forecast conventions (i.e.
> CF-compliant NetCDF)
>  - The output of the pinterp or wrf_interp utility
>
> I suspect that your file is not in one of those formats.
>
> I have attached an example of a NetCDF file (out.nc) with a Mercator
> projection that MET can read.  If you can get your file in that
format,
> MODE should be able to read it.
>
> I hope that helps!
>
> Julie
>
> On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > Transaction: Ticket created by rconne01 at gmail.com
> >        Queue: met_help
> >      Subject: MODE error: Couldn't figure out projection from
information
> > in netCDF file.
> >        Owner: Nobody
> >   Requestors: rconne01 at gmail.com
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> >
> > Hi,
> >
> > Upon running MODE, I encounter the error: Couldn't figure out
projection
> > from information in netCDF file.
> >
> > I assume because it says "in netCDF file" that it's referring to
my obs
> > file?  In this case, that is a radar netCDF processed through radx
called
> > ncf_20141127_010250.nc.
> >
> > I'm pretty sure the projection is Mercator.  Can I just hard code
that in
> > somewhere, since having MODE search for the projection doesn't
work?
> >
> > Thanks,
> > Ryan
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>


--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Fri Jan 20 09:24:59 2017

Hi Ryan.

Unfortunately, it is the user's responsibility to get the files into a
format that MET can read.  I am not familiar with the radar native .gz
files and do not know of any existing conversion utilities.

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.
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!

The out.nc file that I provided to you was only meant as an example of
the
format that MET can read.  It actually contains temperature values and
not
any radar related values.

Julie

On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> Is there a way to make radar netCDF files that MODE will read out of
the
> native .gz files?  I like what I've been doing so far with radx
since I
> want to do my analysis on constant height surfaces, rather than
worry about
> elevation angle, and radx allows me to do that.
>
> Unrelated to that issue, when I run with your out.nc file, it tells
me
> that
> the GRIB1 field abbreviation isn't recognized.  It's for radar
reflectivity
> AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc., and none
of
> those is right, apparently.  Any idea what it is or how I can find
out?  (I
> know for sure that that field was computed when running the UPP, as
I can
> see the field REFD in my gdinfo in GEMPAK.)
>
> Thanks,
> Ryan
>
> On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Ryan.
> >
> > I see that you're having trouble with a NetCDF file using MODE.
Thank
> you
> > for letting us know that you are using a radar NetCDF file
processed
> > through radx.
> >
> > Currently, MET has the ability to read three different flavors of
> NetCDF...
> >  - The format used to write intermediate files in MET (i.e. the
output of
> > pcp_combine)
> >  - NetCDF files following the Climate Forecast conventions (i.e.
> > CF-compliant NetCDF)
> >  - The output of the pinterp or wrf_interp utility
> >
> > I suspect that your file is not in one of those formats.
> >
> > I have attached an example of a NetCDF file (out.nc) with a
Mercator
> > projection that MET can read.  If you can get your file in that
format,
> > MODE should be able to read it.
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > Transaction: Ticket created by rconne01 at gmail.com
> > >        Queue: met_help
> > >      Subject: MODE error: Couldn't figure out projection from
> information
> > > in netCDF file.
> > >        Owner: Nobody
> > >   Requestors: rconne01 at gmail.com
> > >       Status: new
> > >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
> >
> > >
> > >
> > > Hi,
> > >
> > > Upon running MODE, I encounter the error: Couldn't figure out
> projection
> > > from information in netCDF file.
> > >
> > > I assume because it says "in netCDF file" that it's referring to
my obs
> > > file?  In this case, that is a radar netCDF processed through
radx
> called
> > > ncf_20141127_010250.nc.
> > >
> > > I'm pretty sure the projection is Mercator.  Can I just hard
code that
> in
> > > somewhere, since having MODE search for the projection doesn't
work?
> > >
> > > Thanks,
> > > Ryan
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Fri Jan 20 14:21:01 2017

Hi Julie,

The CfRadial files I have don't work either, which is perplexing,
because
they should be cf-compliant, right?

Anyway, right now I'm trying to run pcp_combine to use the output of
it as
the input for MODE as you suggested in your first email.  But even
that is
returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure out
projection from information in netCDF file., on both the cartesian
netCDF,
and the (supposedly) cf-compliant CfRadial files.

An ncdump of the CfRadial file revelas the following projection
information:

int grid_mapping ;
                grid_mapping:grid_mapping_name =
"radar_lidar_radial_scan" ;
                grid_mapping:longitude_of_projection_origin =
-70.2563629150391 ;
                grid_mapping:latitude_of_projection_origin =
43.8913040161133 ;
                grid_mapping:altitude_of_projection_origin = 144. ;
                grid_mapping:false_northing = 0. ;
                grid_mapping:false_easting = 0. ;
        double latitude ;
                latitude:long_name = "latitude" ;
                latitude:units = "degrees_north" ;
                latitude:_FillValue = -9999. ;
        double longitude ;
                longitude:long_name = "longitude" ;
                longitude:units = "degrees_east" ;
                longitude:_FillValue = -9999. ;

and an ncdump of the cartesian ncf contains:

int grid_mapping_0 ;
                grid_mapping_0:grid_mapping_name =
"azimuthal_equidistant" ;
                grid_mapping_0:longitude_of_projection_origin =
-70.25636f ;
                grid_mapping_0:latitude_of_projection_origin =
43.8913f ;
                grid_mapping_0:false_easting = 0.f ;
                grid_mapping_0:false_northing = 0.f ;

So why is even pcp_combine not able to figure out projection info?

Thanks,
Ryan

On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> wrote:

> Hi Ryan.
>
> Unfortunately, it is the user's responsibility to get the files into
a
> format that MET can read.  I am not familiar with the radar native
.gz
> files and do not know of any existing conversion utilities.
>
> 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. 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!
>
> The out.nc file that I provided to you was only meant as an example
of the
> format that MET can read.  It actually contains temperature values
and not
> any radar related values.
>
> Julie
>
> On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> > Hi Julie,
> >
> > Is there a way to make radar netCDF files that MODE will read out
of the
> > native .gz files?  I like what I've been doing so far with radx
since I
> > want to do my analysis on constant height surfaces, rather than
worry
> about
> > elevation angle, and radx allows me to do that.
> >
> > Unrelated to that issue, when I run with your out.nc file, it
tells me
> > that
> > the GRIB1 field abbreviation isn't recognized.  It's for radar
> reflectivity
> > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc., and
none of
> > those is right, apparently.  Any idea what it is or how I can find
out?
> (I
> > know for sure that that field was computed when running the UPP,
as I can
> > see the field REFD in my gdinfo in GEMPAK.)
> >
> > Thanks,
> > Ryan
> >
> > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Ryan.
> > >
> > > I see that you're having trouble with a NetCDF file using MODE.
Thank
> > you
> > > for letting us know that you are using a radar NetCDF file
processed
> > > through radx.
> > >
> > > Currently, MET has the ability to read three different flavors
of
> > NetCDF...
> > >  - The format used to write intermediate files in MET (i.e. the
output
> of
> > > pcp_combine)
> > >  - NetCDF files following the Climate Forecast conventions (i.e.
> > > CF-compliant NetCDF)
> > >  - The output of the pinterp or wrf_interp utility
> > >
> > > I suspect that your file is not in one of those formats.
> > >
> > > I have attached an example of a NetCDF file (out.nc) with a
Mercator
> > > projection that MET can read.  If you can get your file in that
format,
> > > MODE should be able to read it.
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > > Transaction: Ticket created by rconne01 at gmail.com
> > > >        Queue: met_help
> > > >      Subject: MODE error: Couldn't figure out projection from
> > information
> > > > in netCDF file.
> > > >        Owner: Nobody
> > > >   Requestors: rconne01 at gmail.com
> > > >       Status: new
> > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=79208
> > >
> > > >
> > > >
> > > > Hi,
> > > >
> > > > Upon running MODE, I encounter the error: Couldn't figure out
> > projection
> > > > from information in netCDF file.
> > > >
> > > > I assume because it says "in netCDF file" that it's referring
to my
> obs
> > > > file?  In this case, that is a radar netCDF processed through
radx
> > called
> > > > ncf_20141127_010250.nc.
> > > >
> > > > I'm pretty sure the projection is Mercator.  Can I just hard
code
> that
> > in
> > > > somewhere, since having MODE search for the projection doesn't
work?
> > > >
> > > > Thanks,
> > > > Ryan
> > > >
> > > > --
> > > > Ryan Connelly
> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > University
> > > > rconne01 at gmail.com
> > > > ryan.connelly at stonybrook.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>


--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Fri Jan 20 15:28:18 2017

Hi Ryan.

It would be helpful for us to take a look at your data.  Could you
please
upload your data used in the pcp_combine command to our anonymous ftp
site
as follows?

   http://www.dtcenter.org/met/users/support/met_help.php#ftp

Could you please also reply with the command that you used to run
pcp_combine?

Thanks,
Julie


On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> The CfRadial files I have don't work either, which is perplexing,
because
> they should be cf-compliant, right?
>
> Anyway, right now I'm trying to run pcp_combine to use the output of
it as
> the input for MODE as you suggested in your first email.  But even
that is
> returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure
out
> projection from information in netCDF file., on both the cartesian
netCDF,
> and the (supposedly) cf-compliant CfRadial files.
>
> An ncdump of the CfRadial file revelas the following projection
> information:
>
> int grid_mapping ;
>                 grid_mapping:grid_mapping_name =
"radar_lidar_radial_scan"
> ;
>                 grid_mapping:longitude_of_projection_origin =
> -70.2563629150391 ;
>                 grid_mapping:latitude_of_projection_origin =
> 43.8913040161133 ;
>                 grid_mapping:altitude_of_projection_origin = 144. ;
>                 grid_mapping:false_northing = 0. ;
>                 grid_mapping:false_easting = 0. ;
>         double latitude ;
>                 latitude:long_name = "latitude" ;
>                 latitude:units = "degrees_north" ;
>                 latitude:_FillValue = -9999. ;
>         double longitude ;
>                 longitude:long_name = "longitude" ;
>                 longitude:units = "degrees_east" ;
>                 longitude:_FillValue = -9999. ;
>
> and an ncdump of the cartesian ncf contains:
>
> int grid_mapping_0 ;
>                 grid_mapping_0:grid_mapping_name =
> "azimuthal_equidistant" ;
>                 grid_mapping_0:longitude_of_projection_origin =
> -70.25636f ;
>                 grid_mapping_0:latitude_of_projection_origin =
43.8913f ;
>                 grid_mapping_0:false_easting = 0.f ;
>                 grid_mapping_0:false_northing = 0.f ;
>
> So why is even pcp_combine not able to figure out projection info?
>
> Thanks,
> Ryan
>
> On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Ryan.
> >
> > Unfortunately, it is the user's responsibility to get the files
into a
> > format that MET can read.  I am not familiar with the radar native
.gz
> > files and do not know of any existing conversion utilities.
> >
> > 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. 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!
> >
> > The out.nc file that I provided to you was only meant as an
example of
> the
> > format that MET can read.  It actually contains temperature values
and
> not
> > any radar related values.
> >
> > Julie
> >
> > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT
<met_help at ucar.edu
> >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > Is there a way to make radar netCDF files that MODE will read
out of
> the
> > > native .gz files?  I like what I've been doing so far with radx
since I
> > > want to do my analysis on constant height surfaces, rather than
worry
> > about
> > > elevation angle, and radx allows me to do that.
> > >
> > > Unrelated to that issue, when I run with your out.nc file, it
tells me
> > > that
> > > the GRIB1 field abbreviation isn't recognized.  It's for radar
> > reflectivity
> > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc., and
none of
> > > those is right, apparently.  Any idea what it is or how I can
find out?
> > (I
> > > know for sure that that field was computed when running the UPP,
as I
> can
> > > see the field REFD in my gdinfo in GEMPAK.)
> > >
> > > Thanks,
> > > Ryan
> > >
> > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Ryan.
> > > >
> > > > I see that you're having trouble with a NetCDF file using
MODE.
> Thank
> > > you
> > > > for letting us know that you are using a radar NetCDF file
processed
> > > > through radx.
> > > >
> > > > Currently, MET has the ability to read three different flavors
of
> > > NetCDF...
> > > >  - The format used to write intermediate files in MET (i.e.
the
> output
> > of
> > > > pcp_combine)
> > > >  - NetCDF files following the Climate Forecast conventions
(i.e.
> > > > CF-compliant NetCDF)
> > > >  - The output of the pinterp or wrf_interp utility
> > > >
> > > > I suspect that your file is not in one of those formats.
> > > >
> > > > I have attached an example of a NetCDF file (out.nc) with a
Mercator
> > > > projection that MET can read.  If you can get your file in
that
> format,
> > > > MODE should be able to read it.
> > > >
> > > > I hope that helps!
> > > >
> > > > Julie
> > > >
> > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > >        Queue: met_help
> > > > >      Subject: MODE error: Couldn't figure out projection
from
> > > information
> > > > > in netCDF file.
> > > > >        Owner: Nobody
> > > > >   Requestors: rconne01 at gmail.com
> > > > >       Status: new
> > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=79208
> > > >
> > > > >
> > > > >
> > > > > Hi,
> > > > >
> > > > > Upon running MODE, I encounter the error: Couldn't figure
out
> > > projection
> > > > > from information in netCDF file.
> > > > >
> > > > > I assume because it says "in netCDF file" that it's
referring to my
> > obs
> > > > > file?  In this case, that is a radar netCDF processed
through radx
> > > called
> > > > > ncf_20141127_010250.nc.
> > > > >
> > > > > I'm pretty sure the projection is Mercator.  Can I just hard
code
> > that
> > > in
> > > > > somewhere, since having MODE search for the projection
doesn't
> work?
> > > > >
> > > > > Thanks,
> > > > > Ryan
> > > > >
> > > > > --
> > > > > Ryan Connelly
> > > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > > University
> > > > > rconne01 at gmail.com
> > > > > ryan.connelly at stonybrook.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Mon Jan 30 11:36:09 2017

Hi Julie,

I apologize for the delay.  I was away at the general meeting.

The commands I used to run pcp_combine are:

bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
/D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc -pcpdir
/D2/ryan/Research/20141126_Case/temp

bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
/D2/ryan/Research/20141126_Case/temp/
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
-pcpdir /D2/ryan/Research/20141126_Case/temp


The files - a polar coordinate cfrad and a Cartesian gridded ncf, are
located in Connelly_data on the ftp server.

Thanks,
Ryan

On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Ryan.
>
> It would be helpful for us to take a look at your data.  Could you
please
> upload your data used in the pcp_combine command to our anonymous
ftp site
> as follows?
>
>    http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Could you please also reply with the command that you used to run
> pcp_combine?
>
> Thanks,
> Julie
>
>
> On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> > Hi Julie,
> >
> > The CfRadial files I have don't work either, which is perplexing,
because
> > they should be cf-compliant, right?
> >
> > Anyway, right now I'm trying to run pcp_combine to use the output
of it
> as
> > the input for MODE as you suggested in your first email.  But even
that
> is
> > returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't figure
out
> > projection from information in netCDF file., on both the cartesian
> netCDF,
> > and the (supposedly) cf-compliant CfRadial files.
> >
> > An ncdump of the CfRadial file revelas the following projection
> > information:
> >
> > int grid_mapping ;
> >                 grid_mapping:grid_mapping_name =
> "radar_lidar_radial_scan"
> > ;
> >                 grid_mapping:longitude_of_projection_origin =
> > -70.2563629150391 ;
> >                 grid_mapping:latitude_of_projection_origin =
> > 43.8913040161133 ;
> >                 grid_mapping:altitude_of_projection_origin = 144.
;
> >                 grid_mapping:false_northing = 0. ;
> >                 grid_mapping:false_easting = 0. ;
> >         double latitude ;
> >                 latitude:long_name = "latitude" ;
> >                 latitude:units = "degrees_north" ;
> >                 latitude:_FillValue = -9999. ;
> >         double longitude ;
> >                 longitude:long_name = "longitude" ;
> >                 longitude:units = "degrees_east" ;
> >                 longitude:_FillValue = -9999. ;
> >
> > and an ncdump of the cartesian ncf contains:
> >
> > int grid_mapping_0 ;
> >                 grid_mapping_0:grid_mapping_name =
> > "azimuthal_equidistant" ;
> >                 grid_mapping_0:longitude_of_projection_origin =
> > -70.25636f ;
> >                 grid_mapping_0:latitude_of_projection_origin =
43.8913f
> ;
> >                 grid_mapping_0:false_easting = 0.f ;
> >                 grid_mapping_0:false_northing = 0.f ;
> >
> > So why is even pcp_combine not able to figure out projection info?
> >
> > Thanks,
> > Ryan
> >
> > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Ryan.
> > >
> > > Unfortunately, it is the user's responsibility to get the files
into a
> > > format that MET can read.  I am not familiar with the radar
native .gz
> > > files and do not know of any existing conversion utilities.
> > >
> > > 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.
> 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!
> > >
> > > The out.nc file that I provided to you was only meant as an
example of
> > the
> > > format that MET can read.  It actually contains temperature
values and
> > not
> > > any radar related values.
> > >
> > > Julie
> > >
> > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> met_help at ucar.edu
> > >
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
>
> > > >
> > > > Hi Julie,
> > > >
> > > > Is there a way to make radar netCDF files that MODE will read
out of
> > the
> > > > native .gz files?  I like what I've been doing so far with
radx
> since I
> > > > want to do my analysis on constant height surfaces, rather
than worry
> > > about
> > > > elevation angle, and radx allows me to do that.
> > > >
> > > > Unrelated to that issue, when I run with your out.nc file, it
tells
> me
> > > > that
> > > > the GRIB1 field abbreviation isn't recognized.  It's for radar
> > > reflectivity
> > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc.,
and none
> of
> > > > those is right, apparently.  Any idea what it is or how I can
find
> out?
> > > (I
> > > > know for sure that that field was computed when running the
UPP, as I
> > can
> > > > see the field REFD in my gdinfo in GEMPAK.)
> > > >
> > > > Thanks,
> > > > Ryan
> > > >
> > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > > Hi Ryan.
> > > > >
> > > > > I see that you're having trouble with a NetCDF file using
MODE.
> > Thank
> > > > you
> > > > > for letting us know that you are using a radar NetCDF file
> processed
> > > > > through radx.
> > > > >
> > > > > Currently, MET has the ability to read three different
flavors of
> > > > NetCDF...
> > > > >  - The format used to write intermediate files in MET (i.e.
the
> > output
> > > of
> > > > > pcp_combine)
> > > > >  - NetCDF files following the Climate Forecast conventions
(i.e.
> > > > > CF-compliant NetCDF)
> > > > >  - The output of the pinterp or wrf_interp utility
> > > > >
> > > > > I suspect that your file is not in one of those formats.
> > > > >
> > > > > I have attached an example of a NetCDF file (out.nc) with a
> Mercator
> > > > > projection that MET can read.  If you can get your file in
that
> > format,
> > > > > MODE should be able to read it.
> > > > >
> > > > > I hope that helps!
> > > > >
> > > > > Julie
> > > > >
> > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > >        Queue: met_help
> > > > > >      Subject: MODE error: Couldn't figure out projection
from
> > > > information
> > > > > > in netCDF file.
> > > > > >        Owner: Nobody
> > > > > >   Requestors: rconne01 at gmail.com
> > > > > >       Status: new
> > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > Ticket/Display.html?id=79208
> > > > >
> > > > > >
> > > > > >
> > > > > > Hi,
> > > > > >
> > > > > > Upon running MODE, I encounter the error: Couldn't figure
out
> > > > projection
> > > > > > from information in netCDF file.
> > > > > >
> > > > > > I assume because it says "in netCDF file" that it's
referring to
> my
> > > obs
> > > > > > file?  In this case, that is a radar netCDF processed
through
> radx
> > > > called
> > > > > > ncf_20141127_010250.nc.
> > > > > >
> > > > > > I'm pretty sure the projection is Mercator.  Can I just
hard code
> > > that
> > > > in
> > > > > > somewhere, since having MODE search for the projection
doesn't
> > work?
> > > > > >
> > > > > > Thanks,
> > > > > > Ryan
> > > > > >
> > > > > > --
> > > > > > Ryan Connelly
> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> > > > > > University
> > > > > > rconne01 at gmail.com
> > > > > > ryan.connelly at stonybrook.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Ryan Connelly
> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > University
> > > > rconne01 at gmail.com
> > > > ryan.connelly at stonybrook.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>


--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Tue Jan 31 12:23:37 2017

Hi Ryan.

Thank you for uploading your files to our ftp server.  It was helpful
to
have those.

Regarding the file:
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
It looks like this is CF compliant point data.  Unfortunately, MET
only
supports CF compliant gridded data.   Our suggestion would be to get
the
format in the point observation ascii format as described on this page
of
the tutorial:

http://www.dtcenter.org/met/users/support/online_tutorial/METv5.2/tutorial.php?name=ascii2nc&category=index

and then use ascii2nc to convert the data to a NetCDF format that MET
can
read.

Regarding the file: ncf_20141127_010250.nc
METv5.2 had some trouble parsing the file, however, our developmental
version of METv6.0 was able to parse it.  Unfortunately, this file is
using
the azimuthal equidistant projection, which is currently not supported
by
MET.  The projections that MET currently supports are Lat Lon (aka
cylindrical equidistant), Lambert Conformal, Mercator, and Polar
Stereographic.  So, your file would need to use one of those
projections
before MET could work with it.  If we had some funding from SBU, we
could
add support for the azimuthal equidistant projection.

Julie


On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> I apologize for the delay.  I was away at the general meeting.
>
> The commands I used to run pcp_combine are:
>
> bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc -pcpdir
> /D2/ryan/Research/20141126_Case/temp
>
> bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> /D2/ryan/Research/20141126_Case/temp/
>
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
> -pcpdir /D2/ryan/Research/20141126_Case/temp
>
>
> The files - a polar coordinate cfrad and a Cartesian gridded ncf,
are
> located in Connelly_data on the ftp server.
>
> Thanks,
> Ryan
>
> On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Ryan.
> >
> > It would be helpful for us to take a look at your data.  Could you
please
> > upload your data used in the pcp_combine command to our anonymous
ftp
> site
> > as follows?
> >
> >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Could you please also reply with the command that you used to run
> > pcp_combine?
> >
> > Thanks,
> > Julie
> >
> >
> > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT
<met_help at ucar.edu
> >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > The CfRadial files I have don't work either, which is
perplexing,
> because
> > > they should be cf-compliant, right?
> > >
> > > Anyway, right now I'm trying to run pcp_combine to use the
output of it
> > as
> > > the input for MODE as you suggested in your first email.  But
even that
> > is
> > > returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't
figure out
> > > projection from information in netCDF file., on both the
cartesian
> > netCDF,
> > > and the (supposedly) cf-compliant CfRadial files.
> > >
> > > An ncdump of the CfRadial file revelas the following projection
> > > information:
> > >
> > > int grid_mapping ;
> > >                 grid_mapping:grid_mapping_name =
> > "radar_lidar_radial_scan"
> > > ;
> > >                 grid_mapping:longitude_of_projection_origin =
> > > -70.2563629150391 ;
> > >                 grid_mapping:latitude_of_projection_origin =
> > > 43.8913040161133 ;
> > >                 grid_mapping:altitude_of_projection_origin =
144. ;
> > >                 grid_mapping:false_northing = 0. ;
> > >                 grid_mapping:false_easting = 0. ;
> > >         double latitude ;
> > >                 latitude:long_name = "latitude" ;
> > >                 latitude:units = "degrees_north" ;
> > >                 latitude:_FillValue = -9999. ;
> > >         double longitude ;
> > >                 longitude:long_name = "longitude" ;
> > >                 longitude:units = "degrees_east" ;
> > >                 longitude:_FillValue = -9999. ;
> > >
> > > and an ncdump of the cartesian ncf contains:
> > >
> > > int grid_mapping_0 ;
> > >                 grid_mapping_0:grid_mapping_name =
> > > "azimuthal_equidistant" ;
> > >                 grid_mapping_0:longitude_of_projection_origin =
> > > -70.25636f ;
> > >                 grid_mapping_0:latitude_of_projection_origin =
> 43.8913f
> > ;
> > >                 grid_mapping_0:false_easting = 0.f ;
> > >                 grid_mapping_0:false_northing = 0.f ;
> > >
> > > So why is even pcp_combine not able to figure out projection
info?
> > >
> > > Thanks,
> > > Ryan
> > >
> > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Ryan.
> > > >
> > > > Unfortunately, it is the user's responsibility to get the
files into
> a
> > > > format that MET can read.  I am not familiar with the radar
native
> .gz
> > > > files and do not know of any existing conversion utilities.
> > > >
> > > > 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.
> > 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!
> > > >
> > > > The out.nc file that I provided to you was only meant as an
example
> of
> > > the
> > > > format that MET can read.  It actually contains temperature
values
> and
> > > not
> > > > any radar related values.
> > > >
> > > > Julie
> > > >
> > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> > met_help at ucar.edu
> > > >
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > > >
> > > > > Hi Julie,
> > > > >
> > > > > Is there a way to make radar netCDF files that MODE will
read out
> of
> > > the
> > > > > native .gz files?  I like what I've been doing so far with
radx
> > since I
> > > > > want to do my analysis on constant height surfaces, rather
than
> worry
> > > > about
> > > > > elevation angle, and radx allows me to do that.
> > > > >
> > > > > Unrelated to that issue, when I run with your out.nc file,
it
> tells
> > me
> > > > > that
> > > > > the GRIB1 field abbreviation isn't recognized.  It's for
radar
> > > > reflectivity
> > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc.,
and
> none
> > of
> > > > > those is right, apparently.  Any idea what it is or how I
can find
> > out?
> > > > (I
> > > > > know for sure that that field was computed when running the
UPP,
> as I
> > > can
> > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > > >
> > > > > Thanks,
> > > > > Ryan
> > > > >
> > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > > Hi Ryan.
> > > > > >
> > > > > > I see that you're having trouble with a NetCDF file using
MODE.
> > > Thank
> > > > > you
> > > > > > for letting us know that you are using a radar NetCDF file
> > processed
> > > > > > through radx.
> > > > > >
> > > > > > Currently, MET has the ability to read three different
flavors of
> > > > > NetCDF...
> > > > > >  - The format used to write intermediate files in MET
(i.e. the
> > > output
> > > > of
> > > > > > pcp_combine)
> > > > > >  - NetCDF files following the Climate Forecast conventions
(i.e.
> > > > > > CF-compliant NetCDF)
> > > > > >  - The output of the pinterp or wrf_interp utility
> > > > > >
> > > > > > I suspect that your file is not in one of those formats.
> > > > > >
> > > > > > I have attached an example of a NetCDF file (out.nc) with
a
> > Mercator
> > > > > > projection that MET can read.  If you can get your file in
that
> > > format,
> > > > > > MODE should be able to read it.
> > > > > >
> > > > > > I hope that helps!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > > > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > > >        Queue: met_help
> > > > > > >      Subject: MODE error: Couldn't figure out projection
from
> > > > > information
> > > > > > > in netCDF file.
> > > > > > >        Owner: Nobody
> > > > > > >   Requestors: rconne01 at gmail.com
> > > > > > >       Status: new
> > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=79208
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hi,
> > > > > > >
> > > > > > > Upon running MODE, I encounter the error: Couldn't
figure out
> > > > > projection
> > > > > > > from information in netCDF file.
> > > > > > >
> > > > > > > I assume because it says "in netCDF file" that it's
referring
> to
> > my
> > > > obs
> > > > > > > file?  In this case, that is a radar netCDF processed
through
> > radx
> > > > > called
> > > > > > > ncf_20141127_010250.nc.
> > > > > > >
> > > > > > > I'm pretty sure the projection is Mercator.  Can I just
hard
> code
> > > > that
> > > > > in
> > > > > > > somewhere, since having MODE search for the projection
doesn't
> > > work?
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Ryan
> > > > > > >
> > > > > > > --
> > > > > > > Ryan Connelly
> > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> > Valparaiso
> > > > > > > University
> > > > > > > rconne01 at gmail.com
> > > > > > > ryan.connelly at stonybrook.edu
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik
> > > > > > National Center for Atmospheric Research
> > > > > > Research Applications Laboratory
> > > > > > Phone: 303.497.8399
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Ryan Connelly
> > > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > > University
> > > > > rconne01 at gmail.com
> > > > > ryan.connelly at stonybrook.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Wed Feb 01 12:25:24 2017

Hi Julie,

I remade a netcdf file with Lambert Conformal as the projection, and
added
the required global attributes as described in this help list email:
http://mailman.ucar.edu/pipermail/met_help/2014-September/002166.html.
However, pcp_combine still throws the error that projection info can't
be
found in the file.

That file is now in Connelly_data on the ftp site.  Note that it is a
different time: ncf_20141127_001023.nc

Thanks,
Ryan


On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

Hi Ryan.

Thank you for uploading your files to our ftp server.  It was helpful
to
have those.

Regarding the file:
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
It looks like this is CF compliant point data.  Unfortunately, MET
only
supports CF compliant gridded data.   Our suggestion would be to get
the
format in the point observation ascii format as described on this page
of
the tutorial:

http://www.dtcenter.org/met/users/support/online_tutorial/METv5.2/tutorial.php?name=ascii2nc&category=index

and then use ascii2nc to convert the data to a NetCDF format that MET
can
read.

Regarding the file: ncf_20141127_010250.nc
METv5.2 had some trouble parsing the file, however, our developmental
version of METv6.0 was able to parse it.  Unfortunately, this file is
using
the azimuthal equidistant projection, which is currently not supported
by
MET.  The projections that MET currently supports are Lat Lon (aka
cylindrical equidistant), Lambert Conformal, Mercator, and Polar
Stereographic.  So, your file would need to use one of those
projections
before MET could work with it.  If we had some funding from SBU, we
could
add support for the azimuthal equidistant projection.

Julie


On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> I apologize for the delay.  I was away at the general meeting.
>
> The commands I used to run pcp_combine are:
>
> bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc -pcpdir
> /D2/ryan/Research/20141126_Case/temp
>
> bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> /D2/ryan/Research/20141126_Case/temp/
>
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
> -pcpdir /D2/ryan/Research/20141126_Case/temp
>
>
> The files - a polar coordinate cfrad and a Cartesian gridded ncf,
are
> located in Connelly_data on the ftp server.
>
> Thanks,
> Ryan
>
> On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Ryan.
> >
> > It would be helpful for us to take a look at your data.  Could you
please
> > upload your data used in the pcp_combine command to our anonymous
ftp
> site
> > as follows?
> >
> >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Could you please also reply with the command that you used to run
> > pcp_combine?
> >
> > Thanks,
> > Julie
> >
> >
> > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT
<met_help at ucar.edu
> >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > The CfRadial files I have don't work either, which is
perplexing,
> because
> > > they should be cf-compliant, right?
> > >
> > > Anyway, right now I'm trying to run pcp_combine to use the
output of
it
> > as
> > > the input for MODE as you suggested in your first email.  But
even
that
> > is
> > > returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't
figure out
> > > projection from information in netCDF file., on both the
cartesian
> > netCDF,
> > > and the (supposedly) cf-compliant CfRadial files.
> > >
> > > An ncdump of the CfRadial file revelas the following projection
> > > information:
> > >
> > > int grid_mapping ;
> > >                 grid_mapping:grid_mapping_name =
> > "radar_lidar_radial_scan"
> > > ;
> > >                 grid_mapping:longitude_of_projection_origin =
> > > -70.2563629150391 ;
> > >                 grid_mapping:latitude_of_projection_origin =
> > > 43.8913040161133 ;
> > >                 grid_mapping:altitude_of_projection_origin =
144. ;
> > >                 grid_mapping:false_northing = 0. ;
> > >                 grid_mapping:false_easting = 0. ;
> > >         double latitude ;
> > >                 latitude:long_name = "latitude" ;
> > >                 latitude:units = "degrees_north" ;
> > >                 latitude:_FillValue = -9999. ;
> > >         double longitude ;
> > >                 longitude:long_name = "longitude" ;
> > >                 longitude:units = "degrees_east" ;
> > >                 longitude:_FillValue = -9999. ;
> > >
> > > and an ncdump of the cartesian ncf contains:
> > >
> > > int grid_mapping_0 ;
> > >                 grid_mapping_0:grid_mapping_name =
> > > "azimuthal_equidistant" ;
> > >                 grid_mapping_0:longitude_of_projection_origin =
> > > -70.25636f ;
> > >                 grid_mapping_0:latitude_of_projection_origin =
> 43.8913f
> > ;
> > >                 grid_mapping_0:false_easting = 0.f ;
> > >                 grid_mapping_0:false_northing = 0.f ;
> > >
> > > So why is even pcp_combine not able to figure out projection
info?
> > >
> > > Thanks,
> > > Ryan
> > >
> > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Ryan.
> > > >
> > > > Unfortunately, it is the user's responsibility to get the
files into
> a
> > > > format that MET can read.  I am not familiar with the radar
native
> .gz
> > > > files and do not know of any existing conversion utilities.
> > > >
> > > > 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.
> > 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!
> > > >
> > > > The out.nc file that I provided to you was only meant as an
example
> of
> > > the
> > > > format that MET can read.  It actually contains temperature
values
> and
> > > not
> > > > any radar related values.
> > > >
> > > > Julie
> > > >
> > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> > met_help at ucar.edu
> > > >
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > > >
> > > > > Hi Julie,
> > > > >
> > > > > Is there a way to make radar netCDF files that MODE will
read out
> of
> > > the
> > > > > native .gz files?  I like what I've been doing so far with
radx
> > since I
> > > > > want to do my analysis on constant height surfaces, rather
than
> worry
> > > > about
> > > > > elevation angle, and radx allows me to do that.
> > > > >
> > > > > Unrelated to that issue, when I run with your out.nc file,
it
> tells
> > me
> > > > > that
> > > > > the GRIB1 field abbreviation isn't recognized.  It's for
radar
> > > > reflectivity
> > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM, etc.,
and
> none
> > of
> > > > > those is right, apparently.  Any idea what it is or how I
can find
> > out?
> > > > (I
> > > > > know for sure that that field was computed when running the
UPP,
> as I
> > > can
> > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > > >
> > > > > Thanks,
> > > > > Ryan
> > > > >
> > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > > Hi Ryan.
> > > > > >
> > > > > > I see that you're having trouble with a NetCDF file using
MODE.
> > > Thank
> > > > > you
> > > > > > for letting us know that you are using a radar NetCDF file
> > processed
> > > > > > through radx.
> > > > > >
> > > > > > Currently, MET has the ability to read three different
flavors
of
> > > > > NetCDF...
> > > > > >  - The format used to write intermediate files in MET
(i.e. the
> > > output
> > > > of
> > > > > > pcp_combine)
> > > > > >  - NetCDF files following the Climate Forecast conventions
(i.e.
> > > > > > CF-compliant NetCDF)
> > > > > >  - The output of the pinterp or wrf_interp utility
> > > > > >
> > > > > > I suspect that your file is not in one of those formats.
> > > > > >
> > > > > > I have attached an example of a NetCDF file (out.nc) with
a
> > Mercator
> > > > > > projection that MET can read.  If you can get your file in
that
> > > format,
> > > > > > MODE should be able to read it.
> > > > > >
> > > > > > I hope that helps!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > > > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted upon.
> > > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > > >        Queue: met_help
> > > > > > >      Subject: MODE error: Couldn't figure out projection
from
> > > > > information
> > > > > > > in netCDF file.
> > > > > > >        Owner: Nobody
> > > > > > >   Requestors: rconne01 at gmail.com
> > > > > > >       Status: new
> > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=79208
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hi,
> > > > > > >
> > > > > > > Upon running MODE, I encounter the error: Couldn't
figure out
> > > > > projection
> > > > > > > from information in netCDF file.
> > > > > > >
> > > > > > > I assume because it says "in netCDF file" that it's
referring
> to
> > my
> > > > obs
> > > > > > > file?  In this case, that is a radar netCDF processed
through
> > radx
> > > > > called
> > > > > > > ncf_20141127_010250.nc.
> > > > > > >
> > > > > > > I'm pretty sure the projection is Mercator.  Can I just
hard
> code
> > > > that
> > > > > in
> > > > > > > somewhere, since having MODE search for the projection
doesn't
> > > work?
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Ryan
> > > > > > >
> > > > > > > --
> > > > > > > Ryan Connelly
> > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> > Valparaiso
> > > > > > > University
> > > > > > > rconne01 at gmail.com
> > > > > > > ryan.connelly at stonybrook.edu
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik
> > > > > > National Center for Atmospheric Research
> > > > > > Research Applications Laboratory
> > > > > > Phone: 303.497.8399
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Ryan Connelly
> > > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > > University
> > > > > rconne01 at gmail.com
> > > > > ryan.connelly at stonybrook.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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




--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu
--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Feb 02 11:34:27 2017

Hi Ryan.  Thank you for uploading the file. Unfortunately, METv5.2 had
some
trouble parsing the file, however, our developmental version of
METv6.0 was
able to parse it.

We will be working on a bug fix for METv5.2 to fix this issue so that
MET
will be able to parse the file.  We will let you know as soon as the
fix is
ready (likely sometime next week).

Thanks,
Julie

On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> I remade a netcdf file with Lambert Conformal as the projection, and
added
> the required global attributes as described in this help list email:
> http://mailman.ucar.edu/pipermail/met_help/2014-
September/002166.html.
> However, pcp_combine still throws the error that projection info
can't be
> found in the file.
>
> That file is now in Connelly_data on the ftp site.  Note that it is
a
> different time: ncf_20141127_001023.nc
>
> Thanks,
> Ryan
>
>
> On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> Hi Ryan.
>
> Thank you for uploading your files to our ftp server.  It was
helpful to
> have those.
>
> Regarding the file:
>
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_Surveillance_SUR.nc
> It looks like this is CF compliant point data.  Unfortunately, MET
only
> supports CF compliant gridded data.   Our suggestion would be to get
the
> format in the point observation ascii format as described on this
page of
> the tutorial:
>
> http://www.dtcenter.org/met/users/support/online_tutorial/
> METv5.2/tutorial.php?name=ascii2nc&category=index
>
> and then use ascii2nc to convert the data to a NetCDF format that
MET can
> read.
>
> Regarding the file: ncf_20141127_010250.nc
> METv5.2 had some trouble parsing the file, however, our
developmental
> version of METv6.0 was able to parse it.  Unfortunately, this file
is using
> the azimuthal equidistant projection, which is currently not
supported by
> MET.  The projections that MET currently supports are Lat Lon (aka
> cylindrical equidistant), Lambert Conformal, Mercator, and Polar
> Stereographic.  So, your file would need to use one of those
projections
> before MET could work with it.  If we had some funding from SBU, we
could
> add support for the azimuthal equidistant projection.
>
> Julie
>
>
> On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> > Hi Julie,
> >
> > I apologize for the delay.  I was away at the general meeting.
> >
> > The commands I used to run pcp_combine are:
> >
> > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
-pcpdir
> > /D2/ryan/Research/20141126_Case/temp
> >
> > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > /D2/ryan/Research/20141126_Case/temp/
> > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> Surveillance_SUR.nc
> > -pcpdir /D2/ryan/Research/20141126_Case/temp
> >
> >
> > The files - a polar coordinate cfrad and a Cartesian gridded ncf,
are
> > located in Connelly_data on the ftp server.
> >
> > Thanks,
> > Ryan
> >
> > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Ryan.
> > >
> > > It would be helpful for us to take a look at your data.  Could
you
> please
> > > upload your data used in the pcp_combine command to our
anonymous ftp
> > site
> > > as follows?
> > >
> > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Could you please also reply with the command that you used to
run
> > > pcp_combine?
> > >
> > > Thanks,
> > > Julie
> > >
> > >
> > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT <
> met_help at ucar.edu
> > >
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
>
> > > >
> > > > Hi Julie,
> > > >
> > > > The CfRadial files I have don't work either, which is
perplexing,
> > because
> > > > they should be cf-compliant, right?
> > > >
> > > > Anyway, right now I'm trying to run pcp_combine to use the
output of
> it
> > > as
> > > > the input for MODE as you suggested in your first email.  But
even
> that
> > > is
> > > > returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't
figure
> out
> > > > projection from information in netCDF file., on both the
cartesian
> > > netCDF,
> > > > and the (supposedly) cf-compliant CfRadial files.
> > > >
> > > > An ncdump of the CfRadial file revelas the following
projection
> > > > information:
> > > >
> > > > int grid_mapping ;
> > > >                 grid_mapping:grid_mapping_name =
> > > "radar_lidar_radial_scan"
> > > > ;
> > > >                 grid_mapping:longitude_of_projection_origin =
> > > > -70.2563629150391 ;
> > > >                 grid_mapping:latitude_of_projection_origin =
> > > > 43.8913040161133 ;
> > > >                 grid_mapping:altitude_of_projection_origin =
144. ;
> > > >                 grid_mapping:false_northing = 0. ;
> > > >                 grid_mapping:false_easting = 0. ;
> > > >         double latitude ;
> > > >                 latitude:long_name = "latitude" ;
> > > >                 latitude:units = "degrees_north" ;
> > > >                 latitude:_FillValue = -9999. ;
> > > >         double longitude ;
> > > >                 longitude:long_name = "longitude" ;
> > > >                 longitude:units = "degrees_east" ;
> > > >                 longitude:_FillValue = -9999. ;
> > > >
> > > > and an ncdump of the cartesian ncf contains:
> > > >
> > > > int grid_mapping_0 ;
> > > >                 grid_mapping_0:grid_mapping_name =
> > > > "azimuthal_equidistant" ;
> > > >                 grid_mapping_0:longitude_of_projection_origin
=
> > > > -70.25636f ;
> > > >                 grid_mapping_0:latitude_of_projection_origin =
> > 43.8913f
> > > ;
> > > >                 grid_mapping_0:false_easting = 0.f ;
> > > >                 grid_mapping_0:false_northing = 0.f ;
> > > >
> > > > So why is even pcp_combine not able to figure out projection
info?
> > > >
> > > > Thanks,
> > > > Ryan
> > > >
> > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > > Hi Ryan.
> > > > >
> > > > > Unfortunately, it is the user's responsibility to get the
files
> into
> > a
> > > > > format that MET can read.  I am not familiar with the radar
native
> > .gz
> > > > > files and do not know of any existing conversion utilities.
> > > > >
> > > > > 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.
> > > 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!
> > > > >
> > > > > The out.nc file that I provided to you was only meant as an
> example
> > of
> > > > the
> > > > > format that MET can read.  It actually contains temperature
values
> > and
> > > > not
> > > > > any radar related values.
> > > > >
> > > > > Julie
> > > > >
> > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> > > met_help at ucar.edu
> > > > >
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > > > >
> > > > > > Hi Julie,
> > > > > >
> > > > > > Is there a way to make radar netCDF files that MODE will
read out
> > of
> > > > the
> > > > > > native .gz files?  I like what I've been doing so far with
radx
> > > since I
> > > > > > want to do my analysis on constant height surfaces, rather
than
> > worry
> > > > > about
> > > > > > elevation angle, and radx allows me to do that.
> > > > > >
> > > > > > Unrelated to that issue, when I run with your out.nc file,
it
> > tells
> > > me
> > > > > > that
> > > > > > the GRIB1 field abbreviation isn't recognized.  It's for
radar
> > > > > reflectivity
> > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM,
etc., and
> > none
> > > of
> > > > > > those is right, apparently.  Any idea what it is or how I
can
> find
> > > out?
> > > > > (I
> > > > > > know for sure that that field was computed when running
the UPP,
> > as I
> > > > can
> > > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > > > >
> > > > > > Thanks,
> > > > > > Ryan
> > > > > >
> > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT <
> > > > > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Ryan.
> > > > > > >
> > > > > > > I see that you're having trouble with a NetCDF file
using MODE.
> > > > Thank
> > > > > > you
> > > > > > > for letting us know that you are using a radar NetCDF
file
> > > processed
> > > > > > > through radx.
> > > > > > >
> > > > > > > Currently, MET has the ability to read three different
flavors
> of
> > > > > > NetCDF...
> > > > > > >  - The format used to write intermediate files in MET
(i.e. the
> > > > output
> > > > > of
> > > > > > > pcp_combine)
> > > > > > >  - NetCDF files following the Climate Forecast
conventions
> (i.e.
> > > > > > > CF-compliant NetCDF)
> > > > > > >  - The output of the pinterp or wrf_interp utility
> > > > > > >
> > > > > > > I suspect that your file is not in one of those formats.
> > > > > > >
> > > > > > > I have attached an example of a NetCDF file (out.nc)
with a
> > > Mercator
> > > > > > > projection that MET can read.  If you can get your file
in that
> > > > format,
> > > > > > > MODE should be able to read it.
> > > > > > >
> > > > > > > I hope that helps!
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT <
> > > > > > met_help at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted
upon.
> > > > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > > > >        Queue: met_help
> > > > > > > >      Subject: MODE error: Couldn't figure out
projection from
> > > > > > information
> > > > > > > > in netCDF file.
> > > > > > > >        Owner: Nobody
> > > > > > > >   Requestors: rconne01 at gmail.com
> > > > > > > >       Status: new
> > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > Ticket/Display.html?id=79208
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Hi,
> > > > > > > >
> > > > > > > > Upon running MODE, I encounter the error: Couldn't
figure out
> > > > > > projection
> > > > > > > > from information in netCDF file.
> > > > > > > >
> > > > > > > > I assume because it says "in netCDF file" that it's
referring
> > to
> > > my
> > > > > obs
> > > > > > > > file?  In this case, that is a radar netCDF processed
through
> > > radx
> > > > > > called
> > > > > > > > ncf_20141127_010250.nc.
> > > > > > > >
> > > > > > > > I'm pretty sure the projection is Mercator.  Can I
just hard
> > code
> > > > > that
> > > > > > in
> > > > > > > > somewhere, since having MODE search for the projection
> doesn't
> > > > work?
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > Ryan
> > > > > > > >
> > > > > > > > --
> > > > > > > > Ryan Connelly
> > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > Valparaiso
> > > > > > > > University
> > > > > > > > rconne01 at gmail.com
> > > > > > > > ryan.connelly at stonybrook.edu
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Julie Prestopnik
> > > > > > > National Center for Atmospheric Research
> > > > > > > Research Applications Laboratory
> > > > > > > Phone: 303.497.8399
> > > > > > > Email: jpresto at ucar.edu
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Ryan Connelly
> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> > > > > > University
> > > > > > rconne01 at gmail.com
> > > > > > ryan.connelly at stonybrook.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Ryan Connelly
> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > University
> > > > rconne01 at gmail.com
> > > > ryan.connelly at stonybrook.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Thu Feb 02 11:40:39 2017

Hi Julie,

Thanks for trying it out.  It is just a big - not something I can fix
by
choosing a different map projection or editing/creating more
attributes or
something like that?  If so, I'm perfectly content to wait until next
week;
I just wanted to see if I could move ahead now already.

Thanks for your help!
Ryan

On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Ryan.  Thank you for uploading the file. Unfortunately, METv5.2
had some
> trouble parsing the file, however, our developmental version of
METv6.0 was
> able to parse it.
>
> We will be working on a bug fix for METv5.2 to fix this issue so
that MET
> will be able to parse the file.  We will let you know as soon as the
fix is
> ready (likely sometime next week).
>
> Thanks,
> Julie
>
> On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> > Hi Julie,
> >
> > I remade a netcdf file with Lambert Conformal as the projection,
and
> added
> > the required global attributes as described in this help list
email:
> > http://mailman.ucar.edu/pipermail/met_help/2014-
September/002166.html.
> > However, pcp_combine still throws the error that projection info
can't be
> > found in the file.
> >
> > That file is now in Connelly_data on the ftp site.  Note that it
is a
> > different time: ncf_20141127_001023.nc
> >
> > Thanks,
> > Ryan
> >
> >
> > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > Hi Ryan.
> >
> > Thank you for uploading your files to our ftp server.  It was
helpful to
> > have those.
> >
> > Regarding the file:
> > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> Surveillance_SUR.nc
> > It looks like this is CF compliant point data.  Unfortunately, MET
only
> > supports CF compliant gridded data.   Our suggestion would be to
get the
> > format in the point observation ascii format as described on this
page of
> > the tutorial:
> >
> > http://www.dtcenter.org/met/users/support/online_tutorial/
> > METv5.2/tutorial.php?name=ascii2nc&category=index
> >
> > and then use ascii2nc to convert the data to a NetCDF format that
MET can
> > read.
> >
> > Regarding the file: ncf_20141127_010250.nc
> > METv5.2 had some trouble parsing the file, however, our
developmental
> > version of METv6.0 was able to parse it.  Unfortunately, this file
is
> using
> > the azimuthal equidistant projection, which is currently not
supported by
> > MET.  The projections that MET currently supports are Lat Lon (aka
> > cylindrical equidistant), Lambert Conformal, Mercator, and Polar
> > Stereographic.  So, your file would need to use one of those
projections
> > before MET could work with it.  If we had some funding from SBU,
we could
> > add support for the azimuthal equidistant projection.
> >
> > Julie
> >
> >
> > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
> met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > I apologize for the delay.  I was away at the general meeting.
> > >
> > > The commands I used to run pcp_combine are:
> > >
> > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > > /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
-pcpdir
> > > /D2/ryan/Research/20141126_Case/temp
> > >
> > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > > /D2/ryan/Research/20141126_Case/temp/
> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > Surveillance_SUR.nc
> > > -pcpdir /D2/ryan/Research/20141126_Case/temp
> > >
> > >
> > > The files - a polar coordinate cfrad and a Cartesian gridded
ncf, are
> > > located in Connelly_data on the ftp server.
> > >
> > > Thanks,
> > > Ryan
> > >
> > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Ryan.
> > > >
> > > > It would be helpful for us to take a look at your data.  Could
you
> > please
> > > > upload your data used in the pcp_combine command to our
anonymous ftp
> > > site
> > > > as follows?
> > > >
> > > >    http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Could you please also reply with the command that you used to
run
> > > > pcp_combine?
> > > >
> > > > Thanks,
> > > > Julie
> > > >
> > > >
> > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT <
> > met_help at ucar.edu
> > > >
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > > >
> > > > > Hi Julie,
> > > > >
> > > > > The CfRadial files I have don't work either, which is
perplexing,
> > > because
> > > > > they should be cf-compliant, right?
> > > > >
> > > > > Anyway, right now I'm trying to run pcp_combine to use the
output
> of
> > it
> > > > as
> > > > > the input for MODE as you suggested in your first email.
But even
> > that
> > > > is
> > > > > returning ERROR  : NcCfFile::read_netcdf_grid() -> Couldn't
figure
> > out
> > > > > projection from information in netCDF file., on both the
cartesian
> > > > netCDF,
> > > > > and the (supposedly) cf-compliant CfRadial files.
> > > > >
> > > > > An ncdump of the CfRadial file revelas the following
projection
> > > > > information:
> > > > >
> > > > > int grid_mapping ;
> > > > >                 grid_mapping:grid_mapping_name =
> > > > "radar_lidar_radial_scan"
> > > > > ;
> > > > >                 grid_mapping:longitude_of_projection_origin
=
> > > > > -70.2563629150391 ;
> > > > >                 grid_mapping:latitude_of_projection_origin =
> > > > > 43.8913040161133 ;
> > > > >                 grid_mapping:altitude_of_projection_origin =
144.
> ;
> > > > >                 grid_mapping:false_northing = 0. ;
> > > > >                 grid_mapping:false_easting = 0. ;
> > > > >         double latitude ;
> > > > >                 latitude:long_name = "latitude" ;
> > > > >                 latitude:units = "degrees_north" ;
> > > > >                 latitude:_FillValue = -9999. ;
> > > > >         double longitude ;
> > > > >                 longitude:long_name = "longitude" ;
> > > > >                 longitude:units = "degrees_east" ;
> > > > >                 longitude:_FillValue = -9999. ;
> > > > >
> > > > > and an ncdump of the cartesian ncf contains:
> > > > >
> > > > > int grid_mapping_0 ;
> > > > >                 grid_mapping_0:grid_mapping_name =
> > > > > "azimuthal_equidistant" ;
> > > > >
grid_mapping_0:longitude_of_projection_origin =
> > > > > -70.25636f ;
> > > > >                 grid_mapping_0:latitude_of_projection_origin
=
> > > 43.8913f
> > > > ;
> > > > >                 grid_mapping_0:false_easting = 0.f ;
> > > > >                 grid_mapping_0:false_northing = 0.f ;
> > > > >
> > > > > So why is even pcp_combine not able to figure out projection
info?
> > > > >
> > > > > Thanks,
> > > > > Ryan
> > > > >
> > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > > Hi Ryan.
> > > > > >
> > > > > > Unfortunately, it is the user's responsibility to get the
files
> > into
> > > a
> > > > > > format that MET can read.  I am not familiar with the
radar
> native
> > > .gz
> > > > > > files and do not know of any existing conversion
utilities.
> > > > > >
> > > > > > 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.
> > > > 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!
> > > > > >
> > > > > > The out.nc file that I provided to you was only meant as
an
> > example
> > > of
> > > > > the
> > > > > > format that MET can read.  It actually contains
temperature
> values
> > > and
> > > > > not
> > > > > > any radar related values.
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> > > > met_help at ucar.edu
> > > > > >
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
> >
> > > > > > >
> > > > > > > Hi Julie,
> > > > > > >
> > > > > > > Is there a way to make radar netCDF files that MODE will
read
> out
> > > of
> > > > > the
> > > > > > > native .gz files?  I like what I've been doing so far
with radx
> > > > since I
> > > > > > > want to do my analysis on constant height surfaces,
rather than
> > > worry
> > > > > > about
> > > > > > > elevation angle, and radx allows me to do that.
> > > > > > >
> > > > > > > Unrelated to that issue, when I run with your out.nc
file, it
> > > tells
> > > > me
> > > > > > > that
> > > > > > > the GRIB1 field abbreviation isn't recognized.  It's for
radar
> > > > > > reflectivity
> > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM,
etc., and
> > > none
> > > > of
> > > > > > > those is right, apparently.  Any idea what it is or how
I can
> > find
> > > > out?
> > > > > > (I
> > > > > > > know for sure that that field was computed when running
the
> UPP,
> > > as I
> > > > > can
> > > > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Ryan
> > > > > > >
> > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via RT
<
> > > > > > > met_help at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > > Hi Ryan.
> > > > > > > >
> > > > > > > > I see that you're having trouble with a NetCDF file
using
> MODE.
> > > > > Thank
> > > > > > > you
> > > > > > > > for letting us know that you are using a radar NetCDF
file
> > > > processed
> > > > > > > > through radx.
> > > > > > > >
> > > > > > > > Currently, MET has the ability to read three different
> flavors
> > of
> > > > > > > NetCDF...
> > > > > > > >  - The format used to write intermediate files in MET
(i.e.
> the
> > > > > output
> > > > > > of
> > > > > > > > pcp_combine)
> > > > > > > >  - NetCDF files following the Climate Forecast
conventions
> > (i.e.
> > > > > > > > CF-compliant NetCDF)
> > > > > > > >  - The output of the pinterp or wrf_interp utility
> > > > > > > >
> > > > > > > > I suspect that your file is not in one of those
formats.
> > > > > > > >
> > > > > > > > I have attached an example of a NetCDF file (out.nc)
with a
> > > > Mercator
> > > > > > > > projection that MET can read.  If you can get your
file in
> that
> > > > > format,
> > > > > > > > MODE should be able to read it.
> > > > > > > >
> > > > > > > > I hope that helps!
> > > > > > > >
> > > > > > > > Julie
> > > > > > > >
> > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via RT
<
> > > > > > > met_help at ucar.edu>
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted
upon.
> > > > > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > > > > >        Queue: met_help
> > > > > > > > >      Subject: MODE error: Couldn't figure out
projection
> from
> > > > > > > information
> > > > > > > > > in netCDF file.
> > > > > > > > >        Owner: Nobody
> > > > > > > > >   Requestors: rconne01 at gmail.com
> > > > > > > > >       Status: new
> > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > Ticket/Display.html?id=79208
> > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > Hi,
> > > > > > > > >
> > > > > > > > > Upon running MODE, I encounter the error: Couldn't
figure
> out
> > > > > > > projection
> > > > > > > > > from information in netCDF file.
> > > > > > > > >
> > > > > > > > > I assume because it says "in netCDF file" that it's
> referring
> > > to
> > > > my
> > > > > > obs
> > > > > > > > > file?  In this case, that is a radar netCDF
processed
> through
> > > > radx
> > > > > > > called
> > > > > > > > > ncf_20141127_010250.nc.
> > > > > > > > >
> > > > > > > > > I'm pretty sure the projection is Mercator.  Can I
just
> hard
> > > code
> > > > > > that
> > > > > > > in
> > > > > > > > > somewhere, since having MODE search for the
projection
> > doesn't
> > > > > work?
> > > > > > > > >
> > > > > > > > > Thanks,
> > > > > > > > > Ryan
> > > > > > > > >
> > > > > > > > > --
> > > > > > > > > Ryan Connelly
> > > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
> University
> > > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > > Valparaiso
> > > > > > > > > University
> > > > > > > > > rconne01 at gmail.com
> > > > > > > > > ryan.connelly at stonybrook.edu
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > --
> > > > > > > > Julie Prestopnik
> > > > > > > > National Center for Atmospheric Research
> > > > > > > > Research Applications Laboratory
> > > > > > > > Phone: 303.497.8399
> > > > > > > > Email: jpresto at ucar.edu
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Ryan Connelly
> > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> > Valparaiso
> > > > > > > University
> > > > > > > rconne01 at gmail.com
> > > > > > > ryan.connelly at stonybrook.edu
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik
> > > > > > National Center for Atmospheric Research
> > > > > > Research Applications Laboratory
> > > > > > Phone: 303.497.8399
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Ryan Connelly
> > > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > > University
> > > > > rconne01 at gmail.com
> > > > > ryan.connelly at stonybrook.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>


--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Feb 02 14:53:15 2017

Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.

I'll let you know once it's ready.

Thanks,
Julie

On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> Thanks for trying it out.  It is just a big - not something I can
fix by
> choosing a different map projection or editing/creating more
attributes or
> something like that?  If so, I'm perfectly content to wait until
next week;
> I just wanted to see if I could move ahead now already.
>
> Thanks for your help!
> Ryan
>
> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Ryan.  Thank you for uploading the file. Unfortunately, METv5.2
had
> some
> > trouble parsing the file, however, our developmental version of
METv6.0
> was
> > able to parse it.
> >
> > We will be working on a bug fix for METv5.2 to fix this issue so
that MET
> > will be able to parse the file.  We will let you know as soon as
the fix
> is
> > ready (likely sometime next week).
> >
> > Thanks,
> > Julie
> >
> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT
<met_help at ucar.edu
> >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > I remade a netcdf file with Lambert Conformal as the projection,
and
> > added
> > > the required global attributes as described in this help list
email:
> > > http://mailman.ucar.edu/pipermail/met_help/2014-
September/002166.html.
> > > However, pcp_combine still throws the error that projection info
can't
> be
> > > found in the file.
> > >
> > > That file is now in Connelly_data on the ftp site.  Note that it
is a
> > > different time: ncf_20141127_001023.nc
> > >
> > > Thanks,
> > > Ryan
> > >
> > >
> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > Hi Ryan.
> > >
> > > Thank you for uploading your files to our ftp server.  It was
helpful
> to
> > > have those.
> > >
> > > Regarding the file:
> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > Surveillance_SUR.nc
> > > It looks like this is CF compliant point data.  Unfortunately,
MET only
> > > supports CF compliant gridded data.   Our suggestion would be to
get
> the
> > > format in the point observation ascii format as described on
this page
> of
> > > the tutorial:
> > >
> > > http://www.dtcenter.org/met/users/support/online_tutorial/
> > > METv5.2/tutorial.php?name=ascii2nc&category=index
> > >
> > > and then use ascii2nc to convert the data to a NetCDF format
that MET
> can
> > > read.
> > >
> > > Regarding the file: ncf_20141127_010250.nc
> > > METv5.2 had some trouble parsing the file, however, our
developmental
> > > version of METv6.0 was able to parse it.  Unfortunately, this
file is
> > using
> > > the azimuthal equidistant projection, which is currently not
supported
> by
> > > MET.  The projections that MET currently supports are Lat Lon
(aka
> > > cylindrical equidistant), Lambert Conformal, Mercator, and Polar
> > > Stereographic.  So, your file would need to use one of those
> projections
> > > before MET could work with it.  If we had some funding from SBU,
we
> could
> > > add support for the azimuthal equidistant projection.
> > >
> > > Julie
> > >
> > >
> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
> > met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
>
> > > >
> > > > Hi Julie,
> > > >
> > > > I apologize for the delay.  I was away at the general meeting.
> > > >
> > > > The commands I used to run pcp_combine are:
> > > >
> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > > > /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
-pcpdir
> > > > /D2/ryan/Research/20141126_Case/temp
> > > >
> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> > > > /D2/ryan/Research/20141126_Case/temp/
> > > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > > Surveillance_SUR.nc
> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
> > > >
> > > >
> > > > The files - a polar coordinate cfrad and a Cartesian gridded
ncf, are
> > > > located in Connelly_data on the ftp server.
> > > >
> > > > Thanks,
> > > > Ryan
> > > >
> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > > Hi Ryan.
> > > > >
> > > > > It would be helpful for us to take a look at your data.
Could you
> > > please
> > > > > upload your data used in the pcp_combine command to our
anonymous
> ftp
> > > > site
> > > > > as follows?
> > > > >
> > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > >
> > > > > Could you please also reply with the command that you used
to run
> > > > > pcp_combine?
> > > > >
> > > > > Thanks,
> > > > > Julie
> > > > >
> > > > >
> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT <
> > > met_help at ucar.edu
> > > > >
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > > > >
> > > > > > Hi Julie,
> > > > > >
> > > > > > The CfRadial files I have don't work either, which is
perplexing,
> > > > because
> > > > > > they should be cf-compliant, right?
> > > > > >
> > > > > > Anyway, right now I'm trying to run pcp_combine to use the
output
> > of
> > > it
> > > > > as
> > > > > > the input for MODE as you suggested in your first email.
But
> even
> > > that
> > > > > is
> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid() ->
Couldn't
> figure
> > > out
> > > > > > projection from information in netCDF file., on both the
> cartesian
> > > > > netCDF,
> > > > > > and the (supposedly) cf-compliant CfRadial files.
> > > > > >
> > > > > > An ncdump of the CfRadial file revelas the following
projection
> > > > > > information:
> > > > > >
> > > > > > int grid_mapping ;
> > > > > >                 grid_mapping:grid_mapping_name =
> > > > > "radar_lidar_radial_scan"
> > > > > > ;
> > > > > >
grid_mapping:longitude_of_projection_origin =
> > > > > > -70.2563629150391 ;
> > > > > >                 grid_mapping:latitude_of_projection_origin
=
> > > > > > 43.8913040161133 ;
> > > > > >                 grid_mapping:altitude_of_projection_origin
=
> 144.
> > ;
> > > > > >                 grid_mapping:false_northing = 0. ;
> > > > > >                 grid_mapping:false_easting = 0. ;
> > > > > >         double latitude ;
> > > > > >                 latitude:long_name = "latitude" ;
> > > > > >                 latitude:units = "degrees_north" ;
> > > > > >                 latitude:_FillValue = -9999. ;
> > > > > >         double longitude ;
> > > > > >                 longitude:long_name = "longitude" ;
> > > > > >                 longitude:units = "degrees_east" ;
> > > > > >                 longitude:_FillValue = -9999. ;
> > > > > >
> > > > > > and an ncdump of the cartesian ncf contains:
> > > > > >
> > > > > > int grid_mapping_0 ;
> > > > > >                 grid_mapping_0:grid_mapping_name =
> > > > > > "azimuthal_equidistant" ;
> > > > > >
grid_mapping_0:longitude_of_projection_origin =
> > > > > > -70.25636f ;
> > > > > >
grid_mapping_0:latitude_of_projection_origin =
> > > > 43.8913f
> > > > > ;
> > > > > >                 grid_mapping_0:false_easting = 0.f ;
> > > > > >                 grid_mapping_0:false_northing = 0.f ;
> > > > > >
> > > > > > So why is even pcp_combine not able to figure out
projection
> info?
> > > > > >
> > > > > > Thanks,
> > > > > > Ryan
> > > > > >
> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT
<
> > > > > > met_help at ucar.edu
> > > > > > > wrote:
> > > > > >
> > > > > > > Hi Ryan.
> > > > > > >
> > > > > > > Unfortunately, it is the user's responsibility to get
the files
> > > into
> > > > a
> > > > > > > format that MET can read.  I am not familiar with the
radar
> > native
> > > > .gz
> > > > > > > files and do not know of any existing conversion
utilities.
> > > > > > >
> > > > > > > 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.
> > > > > 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!
> > > > > > >
> > > > > > > The out.nc file that I provided to you was only meant as
an
> > > example
> > > > of
> > > > > > the
> > > > > > > format that MET can read.  It actually contains
temperature
> > values
> > > > and
> > > > > > not
> > > > > > > any radar related values.
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
> > > > > met_help at ucar.edu
> > > > > > >
> > > > > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=79208
> > >
> > > > > > > >
> > > > > > > > Hi Julie,
> > > > > > > >
> > > > > > > > Is there a way to make radar netCDF files that MODE
will read
> > out
> > > > of
> > > > > > the
> > > > > > > > native .gz files?  I like what I've been doing so far
with
> radx
> > > > > since I
> > > > > > > > want to do my analysis on constant height surfaces,
rather
> than
> > > > worry
> > > > > > > about
> > > > > > > > elevation angle, and radx allows me to do that.
> > > > > > > >
> > > > > > > > Unrelated to that issue, when I run with your out.nc
file,
> it
> > > > tells
> > > > > me
> > > > > > > > that
> > > > > > > > the GRIB1 field abbreviation isn't recognized.  It's
for
> radar
> > > > > > > reflectivity
> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM,
etc.,
> and
> > > > none
> > > > > of
> > > > > > > > those is right, apparently.  Any idea what it is or
how I can
> > > find
> > > > > out?
> > > > > > > (I
> > > > > > > > know for sure that that field was computed when
running the
> > UPP,
> > > > as I
> > > > > > can
> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > Ryan
> > > > > > > >
> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via
RT <
> > > > > > > > met_help at ucar.edu>
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > > Hi Ryan.
> > > > > > > > >
> > > > > > > > > I see that you're having trouble with a NetCDF file
using
> > MODE.
> > > > > > Thank
> > > > > > > > you
> > > > > > > > > for letting us know that you are using a radar
NetCDF file
> > > > > processed
> > > > > > > > > through radx.
> > > > > > > > >
> > > > > > > > > Currently, MET has the ability to read three
different
> > flavors
> > > of
> > > > > > > > NetCDF...
> > > > > > > > >  - The format used to write intermediate files in
MET (i.e.
> > the
> > > > > > output
> > > > > > > of
> > > > > > > > > pcp_combine)
> > > > > > > > >  - NetCDF files following the Climate Forecast
conventions
> > > (i.e.
> > > > > > > > > CF-compliant NetCDF)
> > > > > > > > >  - The output of the pinterp or wrf_interp utility
> > > > > > > > >
> > > > > > > > > I suspect that your file is not in one of those
formats.
> > > > > > > > >
> > > > > > > > > I have attached an example of a NetCDF file (out.nc)
with
> a
> > > > > Mercator
> > > > > > > > > projection that MET can read.  If you can get your
file in
> > that
> > > > > > format,
> > > > > > > > > MODE should be able to read it.
> > > > > > > > >
> > > > > > > > > I hope that helps!
> > > > > > > > >
> > > > > > > > > Julie
> > > > > > > > >
> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via
RT <
> > > > > > > > met_help at ucar.edu>
> > > > > > > > > wrote:
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted
upon.
> > > > > > > > > > Transaction: Ticket created by rconne01 at gmail.com
> > > > > > > > > >        Queue: met_help
> > > > > > > > > >      Subject: MODE error: Couldn't figure out
projection
> > from
> > > > > > > > information
> > > > > > > > > > in netCDF file.
> > > > > > > > > >        Owner: Nobody
> > > > > > > > > >   Requestors: rconne01 at gmail.com
> > > > > > > > > >       Status: new
> > > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > > Ticket/Display.html?id=79208
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Hi,
> > > > > > > > > >
> > > > > > > > > > Upon running MODE, I encounter the error: Couldn't
figure
> > out
> > > > > > > > projection
> > > > > > > > > > from information in netCDF file.
> > > > > > > > > >
> > > > > > > > > > I assume because it says "in netCDF file" that
it's
> > referring
> > > > to
> > > > > my
> > > > > > > obs
> > > > > > > > > > file?  In this case, that is a radar netCDF
processed
> > through
> > > > > radx
> > > > > > > > called
> > > > > > > > > > ncf_20141127_010250.nc.
> > > > > > > > > >
> > > > > > > > > > I'm pretty sure the projection is Mercator.  Can I
just
> > hard
> > > > code
> > > > > > > that
> > > > > > > > in
> > > > > > > > > > somewhere, since having MODE search for the
projection
> > > doesn't
> > > > > > work?
> > > > > > > > > >
> > > > > > > > > > Thanks,
> > > > > > > > > > Ryan
> > > > > > > > > >
> > > > > > > > > > --
> > > > > > > > > > Ryan Connelly
> > > > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
> > University
> > > > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > > > Valparaiso
> > > > > > > > > > University
> > > > > > > > > > rconne01 at gmail.com
> > > > > > > > > > ryan.connelly at stonybrook.edu
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > --
> > > > > > > > > Julie Prestopnik
> > > > > > > > > National Center for Atmospheric Research
> > > > > > > > > Research Applications Laboratory
> > > > > > > > > Phone: 303.497.8399
> > > > > > > > > Email: jpresto at ucar.edu
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > --
> > > > > > > > Ryan Connelly
> > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > Valparaiso
> > > > > > > > University
> > > > > > > > rconne01 at gmail.com
> > > > > > > > ryan.connelly at stonybrook.edu
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Julie Prestopnik
> > > > > > > National Center for Atmospheric Research
> > > > > > > Research Applications Laboratory
> > > > > > > Phone: 303.497.8399
> > > > > > > Email: jpresto at ucar.edu
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Ryan Connelly
> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> > > > > > University
> > > > > > rconne01 at gmail.com
> > > > > > ryan.connelly at stonybrook.edu
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Ryan Connelly
> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > University
> > > > rconne01 at gmail.com
> > > > ryan.connelly at stonybrook.edu
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>


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

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Feb 09 09:52:16 2017

Hi Ryan.  We have fixed the code and made a patch available.  Please
take a
look at the following page and follow the directions for installing
the
patch:

http://www.dtcenter.org/met/users/support/known_issues/METv5.2/index.php

Please let us know if you have any questions.  Thank you for your
patience
on this issue.

Julie



On Thu, Feb 2, 2017 at 2:52 PM, Julie Prestopnik <jpresto at ucar.edu>
wrote:

> Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.
>
> I'll let you know once it's ready.
>
> Thanks,
> Julie
>
> On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>>
>> Hi Julie,
>>
>> Thanks for trying it out.  It is just a big - not something I can
fix by
>> choosing a different map projection or editing/creating more
attributes or
>> something like that?  If so, I'm perfectly content to wait until
next
>> week;
>> I just wanted to see if I could move ahead now already.
>>
>> Thanks for your help!
>> Ryan
>>
>> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT <
>> met_help at ucar.edu>
>> wrote:
>>
>> > Hi Ryan.  Thank you for uploading the file. Unfortunately,
METv5.2 had
>> some
>> > trouble parsing the file, however, our developmental version of
METv6.0
>> was
>> > able to parse it.
>> >
>> > We will be working on a bug fix for METv5.2 to fix this issue so
that
>> MET
>> > will be able to parse the file.  We will let you know as soon as
the
>> fix is
>> > ready (likely sometime next week).
>> >
>> > Thanks,
>> > Julie
>> >
>> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT <
>> met_help at ucar.edu>
>> > wrote:
>> >
>> > >
>> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>> > >
>> > > Hi Julie,
>> > >
>> > > I remade a netcdf file with Lambert Conformal as the
projection, and
>> > added
>> > > the required global attributes as described in this help list
email:
>> > > http://mailman.ucar.edu/pipermail/met_help/2014-
September/002166.html
>> .
>> > > However, pcp_combine still throws the error that projection
info
>> can't be
>> > > found in the file.
>> > >
>> > > That file is now in Connelly_data on the ftp site.  Note that
it is a
>> > > different time: ncf_20141127_001023.nc
>> > >
>> > > Thanks,
>> > > Ryan
>> > >
>> > >
>> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
>> > > met_help at ucar.edu>
>> > > wrote:
>> > >
>> > > Hi Ryan.
>> > >
>> > > Thank you for uploading your files to our ftp server.  It was
helpful
>> to
>> > > have those.
>> > >
>> > > Regarding the file:
>> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
>> > Surveillance_SUR.nc
>> > > It looks like this is CF compliant point data.  Unfortunately,
MET
>> only
>> > > supports CF compliant gridded data.   Our suggestion would be
to get
>> the
>> > > format in the point observation ascii format as described on
this
>> page of
>> > > the tutorial:
>> > >
>> > > http://www.dtcenter.org/met/users/support/online_tutorial/
>> > > METv5.2/tutorial.php?name=ascii2nc&category=index
>> > >
>> > > and then use ascii2nc to convert the data to a NetCDF format
that MET
>> can
>> > > read.
>> > >
>> > > Regarding the file: ncf_20141127_010250.nc
>> > > METv5.2 had some trouble parsing the file, however, our
developmental
>> > > version of METv6.0 was able to parse it.  Unfortunately, this
file is
>> > using
>> > > the azimuthal equidistant projection, which is currently not
>> supported by
>> > > MET.  The projections that MET currently supports are Lat Lon
(aka
>> > > cylindrical equidistant), Lambert Conformal, Mercator, and
Polar
>> > > Stereographic.  So, your file would need to use one of those
>> projections
>> > > before MET could work with it.  If we had some funding from
SBU, we
>> could
>> > > add support for the azimuthal equidistant projection.
>> > >
>> > > Julie
>> > >
>> > >
>> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
>> > met_help at ucar.edu>
>> > > wrote:
>> > >
>> > > >
>> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
>
>> > > >
>> > > > Hi Julie,
>> > > >
>> > > > I apologize for the delay.  I was away at the general
meeting.
>> > > >
>> > > > The commands I used to run pcp_combine are:
>> > > >
>> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
>> > > > /D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
-pcpdir
>> > > > /D2/ryan/Research/20141126_Case/temp
>> > > >
>> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
>> > > > /D2/ryan/Research/20141126_Case/temp/
>> > > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
>> > > Surveillance_SUR.nc
>> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
>> > > >
>> > > >
>> > > > The files - a polar coordinate cfrad and a Cartesian gridded
ncf,
>> are
>> > > > located in Connelly_data on the ftp server.
>> > > >
>> > > > Thanks,
>> > > > Ryan
>> > > >
>> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
>> > > > met_help at ucar.edu>
>> > > > wrote:
>> > > >
>> > > > > Hi Ryan.
>> > > > >
>> > > > > It would be helpful for us to take a look at your data.
Could you
>> > > please
>> > > > > upload your data used in the pcp_combine command to our
anonymous
>> ftp
>> > > > site
>> > > > > as follows?
>> > > > >
>> > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> > > > >
>> > > > > Could you please also reply with the command that you used
to run
>> > > > > pcp_combine?
>> > > > >
>> > > > > Thanks,
>> > > > > Julie
>> > > > >
>> > > > >
>> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT <
>> > > met_help at ucar.edu
>> > > > >
>> > > > > wrote:
>> > > > >
>> > > > > >
>> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>> > > > > >
>> > > > > > Hi Julie,
>> > > > > >
>> > > > > > The CfRadial files I have don't work either, which is
>> perplexing,
>> > > > because
>> > > > > > they should be cf-compliant, right?
>> > > > > >
>> > > > > > Anyway, right now I'm trying to run pcp_combine to use
the
>> output
>> > of
>> > > it
>> > > > > as
>> > > > > > the input for MODE as you suggested in your first email.
But
>> even
>> > > that
>> > > > > is
>> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid() ->
Couldn't
>> figure
>> > > out
>> > > > > > projection from information in netCDF file., on both the
>> cartesian
>> > > > > netCDF,
>> > > > > > and the (supposedly) cf-compliant CfRadial files.
>> > > > > >
>> > > > > > An ncdump of the CfRadial file revelas the following
projection
>> > > > > > information:
>> > > > > >
>> > > > > > int grid_mapping ;
>> > > > > >                 grid_mapping:grid_mapping_name =
>> > > > > "radar_lidar_radial_scan"
>> > > > > > ;
>> > > > > >
grid_mapping:longitude_of_projection_origin =
>> > > > > > -70.2563629150391 ;
>> > > > > >
grid_mapping:latitude_of_projection_origin =
>> > > > > > 43.8913040161133 ;
>> > > > > >
grid_mapping:altitude_of_projection_origin =
>> 144.
>> > ;
>> > > > > >                 grid_mapping:false_northing = 0. ;
>> > > > > >                 grid_mapping:false_easting = 0. ;
>> > > > > >         double latitude ;
>> > > > > >                 latitude:long_name = "latitude" ;
>> > > > > >                 latitude:units = "degrees_north" ;
>> > > > > >                 latitude:_FillValue = -9999. ;
>> > > > > >         double longitude ;
>> > > > > >                 longitude:long_name = "longitude" ;
>> > > > > >                 longitude:units = "degrees_east" ;
>> > > > > >                 longitude:_FillValue = -9999. ;
>> > > > > >
>> > > > > > and an ncdump of the cartesian ncf contains:
>> > > > > >
>> > > > > > int grid_mapping_0 ;
>> > > > > >                 grid_mapping_0:grid_mapping_name =
>> > > > > > "azimuthal_equidistant" ;
>> > > > > >
grid_mapping_0:longitude_of_projection_origin =
>> > > > > > -70.25636f ;
>> > > > > >
grid_mapping_0:latitude_of_projection_origin =
>> > > > 43.8913f
>> > > > > ;
>> > > > > >                 grid_mapping_0:false_easting = 0.f ;
>> > > > > >                 grid_mapping_0:false_northing = 0.f ;
>> > > > > >
>> > > > > > So why is even pcp_combine not able to figure out
projection
>> info?
>> > > > > >
>> > > > > > Thanks,
>> > > > > > Ryan
>> > > > > >
>> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via RT
<
>> > > > > > met_help at ucar.edu
>> > > > > > > wrote:
>> > > > > >
>> > > > > > > Hi Ryan.
>> > > > > > >
>> > > > > > > Unfortunately, it is the user's responsibility to get
the
>> files
>> > > into
>> > > > a
>> > > > > > > format that MET can read.  I am not familiar with the
radar
>> > native
>> > > > .gz
>> > > > > > > files and do not know of any existing conversion
utilities.
>> > > > > > >
>> > > > > > > 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.
>> > > > > 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!
>> > > > > > >
>> > > > > > > The out.nc file that I provided to you was only meant
as an
>> > > example
>> > > > of
>> > > > > > the
>> > > > > > > format that MET can read.  It actually contains
temperature
>> > values
>> > > > and
>> > > > > > not
>> > > > > > > any radar related values.
>> > > > > > >
>> > > > > > > Julie
>> > > > > > >
>> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via RT <
>> > > > > met_help at ucar.edu
>> > > > > > >
>> > > > > > > wrote:
>> > > > > > >
>> > > > > > > >
>> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
>> ket/Display.html?id=79208
>> > >
>> > > > > > > >
>> > > > > > > > Hi Julie,
>> > > > > > > >
>> > > > > > > > Is there a way to make radar netCDF files that MODE
will
>> read
>> > out
>> > > > of
>> > > > > > the
>> > > > > > > > native .gz files?  I like what I've been doing so far
with
>> radx
>> > > > > since I
>> > > > > > > > want to do my analysis on constant height surfaces,
rather
>> than
>> > > > worry
>> > > > > > > about
>> > > > > > > > elevation angle, and radx allows me to do that.
>> > > > > > > >
>> > > > > > > > Unrelated to that issue, when I run with your out.nc
file,
>> it
>> > > > tells
>> > > > > me
>> > > > > > > > that
>> > > > > > > > the GRIB1 field abbreviation isn't recognized.  It's
for
>> radar
>> > > > > > > reflectivity
>> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL, REFL_10CM,
etc.,
>> and
>> > > > none
>> > > > > of
>> > > > > > > > those is right, apparently.  Any idea what it is or
how I
>> can
>> > > find
>> > > > > out?
>> > > > > > > (I
>> > > > > > > > know for sure that that field was computed when
running the
>> > UPP,
>> > > > as I
>> > > > > > can
>> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
>> > > > > > > >
>> > > > > > > > Thanks,
>> > > > > > > > Ryan
>> > > > > > > >
>> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik via
RT <
>> > > > > > > > met_help at ucar.edu>
>> > > > > > > > wrote:
>> > > > > > > >
>> > > > > > > > > Hi Ryan.
>> > > > > > > > >
>> > > > > > > > > I see that you're having trouble with a NetCDF file
using
>> > MODE.
>> > > > > > Thank
>> > > > > > > > you
>> > > > > > > > > for letting us know that you are using a radar
NetCDF file
>> > > > > processed
>> > > > > > > > > through radx.
>> > > > > > > > >
>> > > > > > > > > Currently, MET has the ability to read three
different
>> > flavors
>> > > of
>> > > > > > > > NetCDF...
>> > > > > > > > >  - The format used to write intermediate files in
MET
>> (i.e.
>> > the
>> > > > > > output
>> > > > > > > of
>> > > > > > > > > pcp_combine)
>> > > > > > > > >  - NetCDF files following the Climate Forecast
conventions
>> > > (i.e.
>> > > > > > > > > CF-compliant NetCDF)
>> > > > > > > > >  - The output of the pinterp or wrf_interp utility
>> > > > > > > > >
>> > > > > > > > > I suspect that your file is not in one of those
formats.
>> > > > > > > > >
>> > > > > > > > > I have attached an example of a NetCDF file
(out.nc)
>> with a
>> > > > > Mercator
>> > > > > > > > > projection that MET can read.  If you can get your
file in
>> > that
>> > > > > > format,
>> > > > > > > > > MODE should be able to read it.
>> > > > > > > > >
>> > > > > > > > > I hope that helps!
>> > > > > > > > >
>> > > > > > > > > Julie
>> > > > > > > > >
>> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly via
RT <
>> > > > > > > > met_help at ucar.edu>
>> > > > > > > > > wrote:
>> > > > > > > > >
>> > > > > > > > > >
>> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was acted
upon.
>> > > > > > > > > > Transaction: Ticket created by rconne01 at gmail.com
>> > > > > > > > > >        Queue: met_help
>> > > > > > > > > >      Subject: MODE error: Couldn't figure out
projection
>> > from
>> > > > > > > > information
>> > > > > > > > > > in netCDF file.
>> > > > > > > > > >        Owner: Nobody
>> > > > > > > > > >   Requestors: rconne01 at gmail.com
>> > > > > > > > > >       Status: new
>> > > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
>> > > > > > > Ticket/Display.html?id=79208
>> > > > > > > > >
>> > > > > > > > > >
>> > > > > > > > > >
>> > > > > > > > > > Hi,
>> > > > > > > > > >
>> > > > > > > > > > Upon running MODE, I encounter the error:
Couldn't
>> figure
>> > out
>> > > > > > > > projection
>> > > > > > > > > > from information in netCDF file.
>> > > > > > > > > >
>> > > > > > > > > > I assume because it says "in netCDF file" that
it's
>> > referring
>> > > > to
>> > > > > my
>> > > > > > > obs
>> > > > > > > > > > file?  In this case, that is a radar netCDF
processed
>> > through
>> > > > > radx
>> > > > > > > > called
>> > > > > > > > > > ncf_20141127_010250.nc.
>> > > > > > > > > >
>> > > > > > > > > > I'm pretty sure the projection is Mercator.  Can
I just
>> > hard
>> > > > code
>> > > > > > > that
>> > > > > > > > in
>> > > > > > > > > > somewhere, since having MODE search for the
projection
>> > > doesn't
>> > > > > > work?
>> > > > > > > > > >
>> > > > > > > > > > Thanks,
>> > > > > > > > > > Ryan
>> > > > > > > > > >
>> > > > > > > > > > --
>> > > > > > > > > > Ryan Connelly
>> > > > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
>> > University
>> > > > > > > > > > B.S. in Meteorology with Minors in Mathematics
and GIS,
>> > > > > Valparaiso
>> > > > > > > > > > University
>> > > > > > > > > > rconne01 at gmail.com
>> > > > > > > > > > ryan.connelly at stonybrook.edu
>> > > > > > > > > >
>> > > > > > > > > >
>> > > > > > > > >
>> > > > > > > > >
>> > > > > > > > > --
>> > > > > > > > > Julie Prestopnik
>> > > > > > > > > National Center for Atmospheric Research
>> > > > > > > > > Research Applications Laboratory
>> > > > > > > > > Phone: 303.497.8399
>> > > > > > > > > Email: jpresto at ucar.edu
>> > > > > > > > >
>> > > > > > > > >
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > --
>> > > > > > > > Ryan Connelly
>> > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
>> > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
>> > > Valparaiso
>> > > > > > > > University
>> > > > > > > > rconne01 at gmail.com
>> > > > > > > > ryan.connelly at stonybrook.edu
>> > > > > > > >
>> > > > > > > >
>> > > > > > >
>> > > > > > >
>> > > > > > > --
>> > > > > > > Julie Prestopnik
>> > > > > > > National Center for Atmospheric Research
>> > > > > > > Research Applications Laboratory
>> > > > > > > Phone: 303.497.8399
>> > > > > > > Email: jpresto at ucar.edu
>> > > > > > >
>> > > > > > >
>> > > > > >
>> > > > > >
>> > > > > > --
>> > > > > > Ryan Connelly
>> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
>> > > > > > B.S. in Meteorology with Minors in Mathematics and GIS,
>> Valparaiso
>> > > > > > University
>> > > > > > rconne01 at gmail.com
>> > > > > > ryan.connelly at stonybrook.edu
>> > > > > >
>> > > > > >
>> > > > >
>> > > > >
>> > > > > --
>> > > > > Julie Prestopnik
>> > > > > National Center for Atmospheric Research
>> > > > > Research Applications Laboratory
>> > > > > Phone: 303.497.8399
>> > > > > Email: jpresto at ucar.edu
>> > > > >
>> > > > >
>> > > >
>> > > >
>> > > > --
>> > > > Ryan Connelly
>> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
>> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> > > > University
>> > > > rconne01 at gmail.com
>> > > > ryan.connelly at stonybrook.edu
>> > > >
>> > > >
>> > >
>> > >
>> > > --
>> > > Julie Prestopnik
>> > > National Center for Atmospheric Research
>> > > Research Applications Laboratory
>> > > Phone: 303.497.8399
>> > > Email: jpresto at ucar.edu
>> > >
>> > >
>> > >
>> > >
>> > > --
>> > > Ryan Connelly
>> > > M.S. Student in Atmospheric Sciences, Stony Brook University
>> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> > > University
>> > > rconne01 at gmail.com
>> > > ryan.connelly at stonybrook.edu
>> > > --
>> > > Ryan Connelly
>> > > M.S. Student in Atmospheric Sciences, Stony Brook University
>> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> > > University
>> > > rconne01 at gmail.com
>> > > ryan.connelly at stonybrook.edu
>> > >
>> > >
>> >
>> >
>> > --
>> > Julie Prestopnik
>> > National Center for Atmospheric Research
>> > Research Applications Laboratory
>> > Phone: 303.497.8399
>> > Email: jpresto at ucar.edu
>> >
>> >
>>
>>
>> --
>> Ryan Connelly
>> M.S. Student in Atmospheric Sciences, Stony Brook University
>> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
>> University
>> rconne01 at gmail.com
>> ryan.connelly at stonybrook.edu
>>
>>
>
>
> --
> 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: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Wed Mar 01 12:16:05 2017

Hi Julie,

I have installed the patch.  However, now unpon running pcp_combine, I
encounter a new error:

ERROR  : NcCfFile::get_grid_mapping_azimuthal_equidistant() ->
Azimuthal
equidistant grid not handled in MET.


This error is repeated for radar netcdf files whose map projection I
have
set through radx as (1) cylindrical equidistant, (2) lambert
conformal, and
(3) mercator, which doesn't seem to make sense to me.  I've verified
through ncdump that these ncf files have the projection:

int grid_mapping_0 ;
                grid_mapping_0:grid_mapping_name = "mercator" ;

(for example)

so I don't really know why it's throwing an error about an Azimuthal
equidistant grid.

Thanks,
Ryan

On Fri, Feb 10, 2017 at 2:41 PM, Ryan Connelly <rconne01 at gmail.com>
wrote:

> Hi Julie,
>
> When I gunzip, and then tar -xvf the patches tar file, it recreates
the
> /data, /doc, and /src directories, but not /bin, so that all of the
tools
> are not updated.  What am I doing wrong?
>
> Ryan
>
> On Thu, Feb 9, 2017 at 11:52 AM, Julie Prestopnik via RT <
> met_help at ucar.edu> wrote:
>
>> Hi Ryan.  We have fixed the code and made a patch available.
Please take
>> a
>>
>> look at the following page and follow the directions for installing
the
>> patch:
>>
>>
http://www.dtcenter.org/met/users/support/known_issues/METv5.2/index.php
>>
>> Please let us know if you have any questions.  Thank you for your
patience
>> on this issue.
>>
>> Julie
>>
>>
>>
>> On Thu, Feb 2, 2017 at 2:52 PM, Julie Prestopnik <jpresto at ucar.edu>
>> wrote:
>>
>> > Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.
>> >
>> > I'll let you know once it's ready.
>> >
>> > Thanks,
>> > Julie
>> >
>> > On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT <
>> met_help at ucar.edu>
>> > wrote:
>> >
>> >>
>> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>> >>
>> >> Hi Julie,
>> >>
>> >> Thanks for trying it out.  It is just a big - not something I
can fix
>> by
>> >> choosing a different map projection or editing/creating more
>> attributes or
>> >> something like that?  If so, I'm perfectly content to wait until
next
>> >> week;
>> >> I just wanted to see if I could move ahead now already.
>> >>
>> >> Thanks for your help!
>> >> Ryan
>> >>
>> >> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT <
>> >> met_help at ucar.edu>
>> >> wrote:
>> >>
>> >> > Hi Ryan.  Thank you for uploading the file. Unfortunately,
METv5.2
>> had
>> >> some
>> >> > trouble parsing the file, however, our developmental version
of
>> METv6.0
>> >> was
>> >> > able to parse it.
>> >> >
>> >> > We will be working on a bug fix for METv5.2 to fix this issue
so that
>> >> MET
>> >> > will be able to parse the file.  We will let you know as soon
as the
>> >> fix is
>> >> > ready (likely sometime next week).
>> >> >
>> >> > Thanks,
>> >> > Julie
>> >> >
>> >> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT <
>> >> met_help at ucar.edu>
>> >> > wrote:
>> >> >
>> >> > >
>> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>> >> > >
>> >> > > Hi Julie,
>> >> > >
>> >> > > I remade a netcdf file with Lambert Conformal as the
projection,
>> and
>> >> > added
>> >> > > the required global attributes as described in this help
list
>> email:
>> >> > > http://mailman.ucar.edu/pipermail/met_help/2014-September/
>> 002166.html
>> >> .
>> >> > > However, pcp_combine still throws the error that projection
info
>> >> can't be
>> >> > > found in the file.
>> >> > >
>> >> > > That file is now in Connelly_data on the ftp site.  Note
that it
>> is a
>> >> > > different time: ncf_20141127_001023.nc
>> >> > >
>> >> > > Thanks,
>> >> > > Ryan
>> >> > >
>> >> > >
>> >> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
>> >> > > met_help at ucar.edu>
>> >> > > wrote:
>> >> > >
>> >> > > Hi Ryan.
>> >> > >
>> >> > > Thank you for uploading your files to our ftp server.  It
was
>> helpful
>> >> to
>> >> > > have those.
>> >> > >
>> >> > > Regarding the file:
>> >> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
>> >> > Surveillance_SUR.nc
>> >> > > It looks like this is CF compliant point data.
Unfortunately, MET
>> >> only
>> >> > > supports CF compliant gridded data.   Our suggestion would
be to
>> get
>> >> the
>> >> > > format in the point observation ascii format as described on
this
>> >> page of
>> >> > > the tutorial:
>> >> > >
>> >> > > http://www.dtcenter.org/met/users/support/online_tutorial/
>> >> > > METv5.2/tutorial.php?name=ascii2nc&category=index
>> >> > >
>> >> > > and then use ascii2nc to convert the data to a NetCDF format
that
>> MET
>> >> can
>> >> > > read.
>> >> > >
>> >> > > Regarding the file: ncf_20141127_010250.nc
>> >> > > METv5.2 had some trouble parsing the file, however, our
>> developmental
>> >> > > version of METv6.0 was able to parse it.  Unfortunately,
this file
>> is
>> >> > using
>> >> > > the azimuthal equidistant projection, which is currently not
>> >> supported by
>> >> > > MET.  The projections that MET currently supports are Lat
Lon (aka
>> >> > > cylindrical equidistant), Lambert Conformal, Mercator, and
Polar
>> >> > > Stereographic.  So, your file would need to use one of those
>> >> projections
>> >> > > before MET could work with it.  If we had some funding from
SBU, we
>> >> could
>> >> > > add support for the azimuthal equidistant projection.
>> >> > >
>> >> > > Julie
>> >> > >
>> >> > >
>> >> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
>> >> > met_help at ucar.edu>
>> >> > > wrote:
>> >> > >
>> >> > > >
>> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>> >> > > >
>> >> > > > Hi Julie,
>> >> > > >
>> >> > > > I apologize for the delay.  I was away at the general
meeting.
>> >> > > >
>> >> > > > The commands I used to run pcp_combine are:
>> >> > > >
>> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
>> >> > > >
/D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
>> -pcpdir
>> >> > > > /D2/ryan/Research/20141126_Case/temp
>> >> > > >
>> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
>> >> > > > /D2/ryan/Research/20141126_Case/temp/
>> >> > > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
>> >> > > Surveillance_SUR.nc
>> >> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
>> >> > > >
>> >> > > >
>> >> > > > The files - a polar coordinate cfrad and a Cartesian
gridded ncf,
>> >> are
>> >> > > > located in Connelly_data on the ftp server.
>> >> > > >
>> >> > > > Thanks,
>> >> > > > Ryan
>> >> > > >
>> >> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT <
>> >> > > > met_help at ucar.edu>
>> >> > > > wrote:
>> >> > > >
>> >> > > > > Hi Ryan.
>> >> > > > >
>> >> > > > > It would be helpful for us to take a look at your data.
Could
>> you
>> >> > > please
>> >> > > > > upload your data used in the pcp_combine command to our
>> anonymous
>> >> ftp
>> >> > > > site
>> >> > > > > as follows?
>> >> > > > >
>> >> > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> >> > > > >
>> >> > > > > Could you please also reply with the command that you
used to
>> run
>> >> > > > > pcp_combine?
>> >> > > > >
>> >> > > > > Thanks,
>> >> > > > > Julie
>> >> > > > >
>> >> > > > >
>> >> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT <
>> >> > > met_help at ucar.edu
>> >> > > > >
>> >> > > > > wrote:
>> >> > > > >
>> >> > > > > >
>> >> > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
>> ket/Display.html?id=79208 >
>> >> > > > > >
>> >> > > > > > Hi Julie,
>> >> > > > > >
>> >> > > > > > The CfRadial files I have don't work either, which is
>> >> perplexing,
>> >> > > > because
>> >> > > > > > they should be cf-compliant, right?
>> >> > > > > >
>> >> > > > > > Anyway, right now I'm trying to run pcp_combine to use
the
>> >> output
>> >> > of
>> >> > > it
>> >> > > > > as
>> >> > > > > > the input for MODE as you suggested in your first
email.  But
>> >> even
>> >> > > that
>> >> > > > > is
>> >> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid() ->
Couldn't
>> >> figure
>> >> > > out
>> >> > > > > > projection from information in netCDF file., on both
the
>> >> cartesian
>> >> > > > > netCDF,
>> >> > > > > > and the (supposedly) cf-compliant CfRadial files.
>> >> > > > > >
>> >> > > > > > An ncdump of the CfRadial file revelas the following
>> projection
>> >> > > > > > information:
>> >> > > > > >
>> >> > > > > > int grid_mapping ;
>> >> > > > > >                 grid_mapping:grid_mapping_name =
>> >> > > > > "radar_lidar_radial_scan"
>> >> > > > > > ;
>> >> > > > > >
grid_mapping:longitude_of_projection_origin
>> =
>> >> > > > > > -70.2563629150391 ;
>> >> > > > > >
grid_mapping:latitude_of_projection_origin =
>> >> > > > > > 43.8913040161133 ;
>> >> > > > > >
grid_mapping:altitude_of_projection_origin =
>> >> 144.
>> >> > ;
>> >> > > > > >                 grid_mapping:false_northing = 0. ;
>> >> > > > > >                 grid_mapping:false_easting = 0. ;
>> >> > > > > >         double latitude ;
>> >> > > > > >                 latitude:long_name = "latitude" ;
>> >> > > > > >                 latitude:units = "degrees_north" ;
>> >> > > > > >                 latitude:_FillValue = -9999. ;
>> >> > > > > >         double longitude ;
>> >> > > > > >                 longitude:long_name = "longitude" ;
>> >> > > > > >                 longitude:units = "degrees_east" ;
>> >> > > > > >                 longitude:_FillValue = -9999. ;
>> >> > > > > >
>> >> > > > > > and an ncdump of the cartesian ncf contains:
>> >> > > > > >
>> >> > > > > > int grid_mapping_0 ;
>> >> > > > > >                 grid_mapping_0:grid_mapping_name =
>> >> > > > > > "azimuthal_equidistant" ;
>> >> > > > > >
grid_mapping_0:longitude_of_projection_origin
>> =
>> >> > > > > > -70.25636f ;
>> >> > > > > >
grid_mapping_0:latitude_of_projection_origin
>> =
>> >> > > > 43.8913f
>> >> > > > > ;
>> >> > > > > >                 grid_mapping_0:false_easting = 0.f ;
>> >> > > > > >                 grid_mapping_0:false_northing = 0.f ;
>> >> > > > > >
>> >> > > > > > So why is even pcp_combine not able to figure out
projection
>> >> info?
>> >> > > > > >
>> >> > > > > > Thanks,
>> >> > > > > > Ryan
>> >> > > > > >
>> >> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik via
RT <
>> >> > > > > > met_help at ucar.edu
>> >> > > > > > > wrote:
>> >> > > > > >
>> >> > > > > > > Hi Ryan.
>> >> > > > > > >
>> >> > > > > > > Unfortunately, it is the user's responsibility to
get the
>> >> files
>> >> > > into
>> >> > > > a
>> >> > > > > > > format that MET can read.  I am not familiar with
the radar
>> >> > native
>> >> > > > .gz
>> >> > > > > > > files and do not know of any existing conversion
utilities.
>> >> > > > > > >
>> >> > > > > > > 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.
>> >> > > > > 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!
>> >> > > > > > >
>> >> > > > > > > The out.nc file that I provided to you was only
meant as
>> an
>> >> > > example
>> >> > > > of
>> >> > > > > > the
>> >> > > > > > > format that MET can read.  It actually contains
temperature
>> >> > values
>> >> > > > and
>> >> > > > > > not
>> >> > > > > > > any radar related values.
>> >> > > > > > >
>> >> > > > > > > Julie
>> >> > > > > > >
>> >> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via
RT <
>> >> > > > > met_help at ucar.edu
>> >> > > > > > >
>> >> > > > > > > wrote:
>> >> > > > > > >
>> >> > > > > > > >
>> >> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
>> >> ket/Display.html?id=79208
>> >> > >
>> >> > > > > > > >
>> >> > > > > > > > Hi Julie,
>> >> > > > > > > >
>> >> > > > > > > > Is there a way to make radar netCDF files that
MODE will
>> >> read
>> >> > out
>> >> > > > of
>> >> > > > > > the
>> >> > > > > > > > native .gz files?  I like what I've been doing so
far
>> with
>> >> radx
>> >> > > > > since I
>> >> > > > > > > > want to do my analysis on constant height
surfaces,
>> rather
>> >> than
>> >> > > > worry
>> >> > > > > > > about
>> >> > > > > > > > elevation angle, and radx allows me to do that.
>> >> > > > > > > >
>> >> > > > > > > > Unrelated to that issue, when I run with your
out.nc
>> file,
>> >> it
>> >> > > > tells
>> >> > > > > me
>> >> > > > > > > > that
>> >> > > > > > > > the GRIB1 field abbreviation isn't recognized.
It's for
>> >> radar
>> >> > > > > > > reflectivity
>> >> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL,
REFL_10CM,
>> etc.,
>> >> and
>> >> > > > none
>> >> > > > > of
>> >> > > > > > > > those is right, apparently.  Any idea what it is
or how I
>> >> can
>> >> > > find
>> >> > > > > out?
>> >> > > > > > > (I
>> >> > > > > > > > know for sure that that field was computed when
running
>> the
>> >> > UPP,
>> >> > > > as I
>> >> > > > > > can
>> >> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
>> >> > > > > > > >
>> >> > > > > > > > Thanks,
>> >> > > > > > > > Ryan
>> >> > > > > > > >
>> >> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie Prestopnik
via RT
>> <
>> >> > > > > > > > met_help at ucar.edu>
>> >> > > > > > > > wrote:
>> >> > > > > > > >
>> >> > > > > > > > > Hi Ryan.
>> >> > > > > > > > >
>> >> > > > > > > > > I see that you're having trouble with a NetCDF
file
>> using
>> >> > MODE.
>> >> > > > > > Thank
>> >> > > > > > > > you
>> >> > > > > > > > > for letting us know that you are using a radar
NetCDF
>> file
>> >> > > > > processed
>> >> > > > > > > > > through radx.
>> >> > > > > > > > >
>> >> > > > > > > > > Currently, MET has the ability to read three
different
>> >> > flavors
>> >> > > of
>> >> > > > > > > > NetCDF...
>> >> > > > > > > > >  - The format used to write intermediate files
in MET
>> >> (i.e.
>> >> > the
>> >> > > > > > output
>> >> > > > > > > of
>> >> > > > > > > > > pcp_combine)
>> >> > > > > > > > >  - NetCDF files following the Climate Forecast
>> conventions
>> >> > > (i.e.
>> >> > > > > > > > > CF-compliant NetCDF)
>> >> > > > > > > > >  - The output of the pinterp or wrf_interp
utility
>> >> > > > > > > > >
>> >> > > > > > > > > I suspect that your file is not in one of those
>> formats.
>> >> > > > > > > > >
>> >> > > > > > > > > I have attached an example of a NetCDF file
(out.nc)
>> >> with a
>> >> > > > > Mercator
>> >> > > > > > > > > projection that MET can read.  If you can get
your
>> file in
>> >> > that
>> >> > > > > > format,
>> >> > > > > > > > > MODE should be able to read it.
>> >> > > > > > > > >
>> >> > > > > > > > > I hope that helps!
>> >> > > > > > > > >
>> >> > > > > > > > > Julie
>> >> > > > > > > > >
>> >> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan Connelly
via RT
>> <
>> >> > > > > > > > met_help at ucar.edu>
>> >> > > > > > > > > wrote:
>> >> > > > > > > > >
>> >> > > > > > > > > >
>> >> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was
acted
>> upon.
>> >> > > > > > > > > > Transaction: Ticket created by
rconne01 at gmail.com
>> >> > > > > > > > > >        Queue: met_help
>> >> > > > > > > > > >      Subject: MODE error: Couldn't figure out
>> projection
>> >> > from
>> >> > > > > > > > information
>> >> > > > > > > > > > in netCDF file.
>> >> > > > > > > > > >        Owner: Nobody
>> >> > > > > > > > > >   Requestors: rconne01 at gmail.com
>> >> > > > > > > > > >       Status: new
>> >> > > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
>> >> > > > > > > Ticket/Display.html?id=79208
>> >> > > > > > > > >
>> >> > > > > > > > > >
>> >> > > > > > > > > >
>> >> > > > > > > > > > Hi,
>> >> > > > > > > > > >
>> >> > > > > > > > > > Upon running MODE, I encounter the error:
Couldn't
>> >> figure
>> >> > out
>> >> > > > > > > > projection
>> >> > > > > > > > > > from information in netCDF file.
>> >> > > > > > > > > >
>> >> > > > > > > > > > I assume because it says "in netCDF file" that
it's
>> >> > referring
>> >> > > > to
>> >> > > > > my
>> >> > > > > > > obs
>> >> > > > > > > > > > file?  In this case, that is a radar netCDF
processed
>> >> > through
>> >> > > > > radx
>> >> > > > > > > > called
>> >> > > > > > > > > > ncf_20141127_010250.nc.
>> >> > > > > > > > > >
>> >> > > > > > > > > > I'm pretty sure the projection is Mercator.
Can I
>> just
>> >> > hard
>> >> > > > code
>> >> > > > > > > that
>> >> > > > > > > > in
>> >> > > > > > > > > > somewhere, since having MODE search for the
>> projection
>> >> > > doesn't
>> >> > > > > > work?
>> >> > > > > > > > > >
>> >> > > > > > > > > > Thanks,
>> >> > > > > > > > > > Ryan
>> >> > > > > > > > > >
>> >> > > > > > > > > > --
>> >> > > > > > > > > > Ryan Connelly
>> >> > > > > > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
>> >> > University
>> >> > > > > > > > > > B.S. in Meteorology with Minors in Mathematics
and
>> GIS,
>> >> > > > > Valparaiso
>> >> > > > > > > > > > University
>> >> > > > > > > > > > rconne01 at gmail.com
>> >> > > > > > > > > > ryan.connelly at stonybrook.edu
>> >> > > > > > > > > >
>> >> > > > > > > > > >
>> >> > > > > > > > >
>> >> > > > > > > > >
>> >> > > > > > > > > --
>> >> > > > > > > > > Julie Prestopnik
>> >> > > > > > > > > National Center for Atmospheric Research
>> >> > > > > > > > > Research Applications Laboratory
>> >> > > > > > > > > Phone: 303.497.8399
>> >> > > > > > > > > Email: jpresto at ucar.edu
>> >> > > > > > > > >
>> >> > > > > > > > >
>> >> > > > > > > >
>> >> > > > > > > >
>> >> > > > > > > > --
>> >> > > > > > > > Ryan Connelly
>> >> > > > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
>> University
>> >> > > > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
>> >> > > Valparaiso
>> >> > > > > > > > University
>> >> > > > > > > > rconne01 at gmail.com
>> >> > > > > > > > ryan.connelly at stonybrook.edu
>> >> > > > > > > >
>> >> > > > > > > >
>> >> > > > > > >
>> >> > > > > > >
>> >> > > > > > > --
>> >> > > > > > > Julie Prestopnik
>> >> > > > > > > National Center for Atmospheric Research
>> >> > > > > > > Research Applications Laboratory
>> >> > > > > > > Phone: 303.497.8399
>> >> > > > > > > Email: jpresto at ucar.edu
>> >> > > > > > >
>> >> > > > > > >
>> >> > > > > >
>> >> > > > > >
>> >> > > > > > --
>> >> > > > > > Ryan Connelly
>> >> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
>> >> > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
>> >> Valparaiso
>> >> > > > > > University
>> >> > > > > > rconne01 at gmail.com
>> >> > > > > > ryan.connelly at stonybrook.edu
>> >> > > > > >
>> >> > > > > >
>> >> > > > >
>> >> > > > >
>> >> > > > > --
>> >> > > > > Julie Prestopnik
>> >> > > > > National Center for Atmospheric Research
>> >> > > > > Research Applications Laboratory
>> >> > > > > Phone: 303.497.8399
>> >> > > > > Email: jpresto at ucar.edu
>> >> > > > >
>> >> > > > >
>> >> > > >
>> >> > > >
>> >> > > > --
>> >> > > > Ryan Connelly
>> >> > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
>> >> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
>> Valparaiso
>> >> > > > University
>> >> > > > rconne01 at gmail.com
>> >> > > > ryan.connelly at stonybrook.edu
>> >> > > >
>> >> > > >
>> >> > >
>> >> > >
>> >> > > --
>> >> > > Julie Prestopnik
>> >> > > National Center for Atmospheric Research
>> >> > > Research Applications Laboratory
>> >> > > Phone: 303.497.8399
>> >> > > Email: jpresto at ucar.edu
>> >> > >
>> >> > >
>> >> > >
>> >> > >
>> >> > > --
>> >> > > Ryan Connelly
>> >> > > M.S. Student in Atmospheric Sciences, Stony Brook University
>> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> >> > > University
>> >> > > rconne01 at gmail.com
>> >> > > ryan.connelly at stonybrook.edu
>> >> > > --
>> >> > > Ryan Connelly
>> >> > > M.S. Student in Atmospheric Sciences, Stony Brook University
>> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> >> > > University
>> >> > > rconne01 at gmail.com
>> >> > > ryan.connelly at stonybrook.edu
>> >> > >
>> >> > >
>> >> >
>> >> >
>> >> > --
>> >> > Julie Prestopnik
>> >> > National Center for Atmospheric Research
>> >> > Research Applications Laboratory
>> >> > Phone: 303.497.8399
>> >> > Email: jpresto at ucar.edu
>> >> >
>> >> >
>> >>
>> >>
>> >> --
>> >> Ryan Connelly
>> >> M.S. Student in Atmospheric Sciences, Stony Brook University
>> >> B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
>> >> University
>> >> rconne01 at gmail.com
>> >> ryan.connelly at stonybrook.edu
>> >>
>> >>
>> >
>> >
>> > --
>> > 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
>>
>>
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>



--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Mar 02 10:07:30 2017

Hi Ryan.  I apologize for the delay in responding. I was out of the
office.

Could you please upload the files in question to our ftp server and
also
send me the full command that you are running so that I can try to
reproduce what you are seeing?

Thanks!

Julie

On Wed, Mar 1, 2017 at 12:16 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie,
>
> I have installed the patch.  However, now unpon running pcp_combine,
I
> encounter a new error:
>
> ERROR  : NcCfFile::get_grid_mapping_azimuthal_equidistant() ->
Azimuthal
> equidistant grid not handled in MET.
>
>
> This error is repeated for radar netcdf files whose map projection I
have
> set through radx as (1) cylindrical equidistant, (2) lambert
conformal, and
> (3) mercator, which doesn't seem to make sense to me.  I've verified
> through ncdump that these ncf files have the projection:
>
> int grid_mapping_0 ;
>                 grid_mapping_0:grid_mapping_name = "mercator" ;
>
> (for example)
>
> so I don't really know why it's throwing an error about an Azimuthal
> equidistant grid.
>
> Thanks,
> Ryan
>
> On Fri, Feb 10, 2017 at 2:41 PM, Ryan Connelly <rconne01 at gmail.com>
wrote:
>
> > Hi Julie,
> >
> > When I gunzip, and then tar -xvf the patches tar file, it
recreates the
> > /data, /doc, and /src directories, but not /bin, so that all of
the tools
> > are not updated.  What am I doing wrong?
> >
> > Ryan
> >
> > On Thu, Feb 9, 2017 at 11:52 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Hi Ryan.  We have fixed the code and made a patch available.
Please
> take
> >> a
> >>
> >> look at the following page and follow the directions for
installing the
> >> patch:
> >>
> >> http://www.dtcenter.org/met/users/support/known_issues/
> METv5.2/index.php
> >>
> >> Please let us know if you have any questions.  Thank you for your
> patience
> >> on this issue.
> >>
> >> Julie
> >>
> >>
> >>
> >> On Thu, Feb 2, 2017 at 2:52 PM, Julie Prestopnik
<jpresto at ucar.edu>
> >> wrote:
> >>
> >> > Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.
> >> >
> >> > I'll let you know once it's ready.
> >> >
> >> > Thanks,
> >> > Julie
> >> >
> >> > On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT <
> >> met_help at ucar.edu>
> >> > wrote:
> >> >
> >> >>
> >> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
>
> >> >>
> >> >> Hi Julie,
> >> >>
> >> >> Thanks for trying it out.  It is just a big - not something I
can fix
> >> by
> >> >> choosing a different map projection or editing/creating more
> >> attributes or
> >> >> something like that?  If so, I'm perfectly content to wait
until next
> >> >> week;
> >> >> I just wanted to see if I could move ahead now already.
> >> >>
> >> >> Thanks for your help!
> >> >> Ryan
> >> >>
> >> >> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT <
> >> >> met_help at ucar.edu>
> >> >> wrote:
> >> >>
> >> >> > Hi Ryan.  Thank you for uploading the file. Unfortunately,
METv5.2
> >> had
> >> >> some
> >> >> > trouble parsing the file, however, our developmental version
of
> >> METv6.0
> >> >> was
> >> >> > able to parse it.
> >> >> >
> >> >> > We will be working on a bug fix for METv5.2 to fix this
issue so
> that
> >> >> MET
> >> >> > will be able to parse the file.  We will let you know as
soon as
> the
> >> >> fix is
> >> >> > ready (likely sometime next week).
> >> >> >
> >> >> > Thanks,
> >> >> > Julie
> >> >> >
> >> >> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT <
> >> >> met_help at ucar.edu>
> >> >> > wrote:
> >> >> >
> >> >> > >
> >> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >> >> > >
> >> >> > > Hi Julie,
> >> >> > >
> >> >> > > I remade a netcdf file with Lambert Conformal as the
projection,
> >> and
> >> >> > added
> >> >> > > the required global attributes as described in this help
list
> >> email:
> >> >> > > http://mailman.ucar.edu/pipermail/met_help/2014-September/
> >> 002166.html
> >> >> .
> >> >> > > However, pcp_combine still throws the error that
projection info
> >> >> can't be
> >> >> > > found in the file.
> >> >> > >
> >> >> > > That file is now in Connelly_data on the ftp site.  Note
that it
> >> is a
> >> >> > > different time: ncf_20141127_001023.nc
> >> >> > >
> >> >> > > Thanks,
> >> >> > > Ryan
> >> >> > >
> >> >> > >
> >> >> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT <
> >> >> > > met_help at ucar.edu>
> >> >> > > wrote:
> >> >> > >
> >> >> > > Hi Ryan.
> >> >> > >
> >> >> > > Thank you for uploading your files to our ftp server.  It
was
> >> helpful
> >> >> to
> >> >> > > have those.
> >> >> > >
> >> >> > > Regarding the file:
> >> >> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> >> >> > Surveillance_SUR.nc
> >> >> > > It looks like this is CF compliant point data.
Unfortunately,
> MET
> >> >> only
> >> >> > > supports CF compliant gridded data.   Our suggestion would
be to
> >> get
> >> >> the
> >> >> > > format in the point observation ascii format as described
on this
> >> >> page of
> >> >> > > the tutorial:
> >> >> > >
> >> >> > > http://www.dtcenter.org/met/users/support/online_tutorial/
> >> >> > > METv5.2/tutorial.php?name=ascii2nc&category=index
> >> >> > >
> >> >> > > and then use ascii2nc to convert the data to a NetCDF
format that
> >> MET
> >> >> can
> >> >> > > read.
> >> >> > >
> >> >> > > Regarding the file: ncf_20141127_010250.nc
> >> >> > > METv5.2 had some trouble parsing the file, however, our
> >> developmental
> >> >> > > version of METv6.0 was able to parse it.  Unfortunately,
this
> file
> >> is
> >> >> > using
> >> >> > > the azimuthal equidistant projection, which is currently
not
> >> >> supported by
> >> >> > > MET.  The projections that MET currently supports are Lat
Lon
> (aka
> >> >> > > cylindrical equidistant), Lambert Conformal, Mercator, and
Polar
> >> >> > > Stereographic.  So, your file would need to use one of
those
> >> >> projections
> >> >> > > before MET could work with it.  If we had some funding
from SBU,
> we
> >> >> could
> >> >> > > add support for the azimuthal equidistant projection.
> >> >> > >
> >> >> > > Julie
> >> >> > >
> >> >> > >
> >> >> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
> >> >> > met_help at ucar.edu>
> >> >> > > wrote:
> >> >> > >
> >> >> > > >
> >> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
> >
> >> >> > > >
> >> >> > > > Hi Julie,
> >> >> > > >
> >> >> > > > I apologize for the delay.  I was away at the general
meeting.
> >> >> > > >
> >> >> > > > The commands I used to run pcp_combine are:
> >> >> > > >
> >> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> >> >> > > >
/D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
> >> -pcpdir
> >> >> > > > /D2/ryan/Research/20141126_Case/temp
> >> >> > > >
> >> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250 1
> >> >> > > > /D2/ryan/Research/20141126_Case/temp/
> >> >> > > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> >> >> > > Surveillance_SUR.nc
> >> >> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
> >> >> > > >
> >> >> > > >
> >> >> > > > The files - a polar coordinate cfrad and a Cartesian
gridded
> ncf,
> >> >> are
> >> >> > > > located in Connelly_data on the ftp server.
> >> >> > > >
> >> >> > > > Thanks,
> >> >> > > > Ryan
> >> >> > > >
> >> >> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via RT
<
> >> >> > > > met_help at ucar.edu>
> >> >> > > > wrote:
> >> >> > > >
> >> >> > > > > Hi Ryan.
> >> >> > > > >
> >> >> > > > > It would be helpful for us to take a look at your
data.
> Could
> >> you
> >> >> > > please
> >> >> > > > > upload your data used in the pcp_combine command to
our
> >> anonymous
> >> >> ftp
> >> >> > > > site
> >> >> > > > > as follows?
> >> >> > > > >
> >> >> > > > >    http://www.dtcenter.org/met/
> users/support/met_help.php#ftp
> >> >> > > > >
> >> >> > > > > Could you please also reply with the command that you
used to
> >> run
> >> >> > > > > pcp_combine?
> >> >> > > > >
> >> >> > > > > Thanks,
> >> >> > > > > Julie
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via RT
<
> >> >> > > met_help at ucar.edu
> >> >> > > > >
> >> >> > > > > wrote:
> >> >> > > > >
> >> >> > > > > >
> >> >> > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> >> ket/Display.html?id=79208 >
> >> >> > > > > >
> >> >> > > > > > Hi Julie,
> >> >> > > > > >
> >> >> > > > > > The CfRadial files I have don't work either, which
is
> >> >> perplexing,
> >> >> > > > because
> >> >> > > > > > they should be cf-compliant, right?
> >> >> > > > > >
> >> >> > > > > > Anyway, right now I'm trying to run pcp_combine to
use the
> >> >> output
> >> >> > of
> >> >> > > it
> >> >> > > > > as
> >> >> > > > > > the input for MODE as you suggested in your first
email.
> But
> >> >> even
> >> >> > > that
> >> >> > > > > is
> >> >> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid() ->
Couldn't
> >> >> figure
> >> >> > > out
> >> >> > > > > > projection from information in netCDF file., on both
the
> >> >> cartesian
> >> >> > > > > netCDF,
> >> >> > > > > > and the (supposedly) cf-compliant CfRadial files.
> >> >> > > > > >
> >> >> > > > > > An ncdump of the CfRadial file revelas the following
> >> projection
> >> >> > > > > > information:
> >> >> > > > > >
> >> >> > > > > > int grid_mapping ;
> >> >> > > > > >                 grid_mapping:grid_mapping_name =
> >> >> > > > > "radar_lidar_radial_scan"
> >> >> > > > > > ;
> >> >> > > > > >                 grid_mapping:longitude_of_
> projection_origin
> >> =
> >> >> > > > > > -70.2563629150391 ;
> >> >> > > > > >
grid_mapping:latitude_of_projection_origin
> =
> >> >> > > > > > 43.8913040161133 ;
> >> >> > > > > >
grid_mapping:altitude_of_projection_origin
> =
> >> >> 144.
> >> >> > ;
> >> >> > > > > >                 grid_mapping:false_northing = 0. ;
> >> >> > > > > >                 grid_mapping:false_easting = 0. ;
> >> >> > > > > >         double latitude ;
> >> >> > > > > >                 latitude:long_name = "latitude" ;
> >> >> > > > > >                 latitude:units = "degrees_north" ;
> >> >> > > > > >                 latitude:_FillValue = -9999. ;
> >> >> > > > > >         double longitude ;
> >> >> > > > > >                 longitude:long_name = "longitude" ;
> >> >> > > > > >                 longitude:units = "degrees_east" ;
> >> >> > > > > >                 longitude:_FillValue = -9999. ;
> >> >> > > > > >
> >> >> > > > > > and an ncdump of the cartesian ncf contains:
> >> >> > > > > >
> >> >> > > > > > int grid_mapping_0 ;
> >> >> > > > > >                 grid_mapping_0:grid_mapping_name =
> >> >> > > > > > "azimuthal_equidistant" ;
> >> >> > > > > >                 grid_mapping_0:longitude_of_
> projection_origin
> >> =
> >> >> > > > > > -70.25636f ;
> >> >> > > > > >                 grid_mapping_0:latitude_of_
> projection_origin
> >> =
> >> >> > > > 43.8913f
> >> >> > > > > ;
> >> >> > > > > >                 grid_mapping_0:false_easting = 0.f ;
> >> >> > > > > >                 grid_mapping_0:false_northing = 0.f
;
> >> >> > > > > >
> >> >> > > > > > So why is even pcp_combine not able to figure out
> projection
> >> >> info?
> >> >> > > > > >
> >> >> > > > > > Thanks,
> >> >> > > > > > Ryan
> >> >> > > > > >
> >> >> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik
via RT <
> >> >> > > > > > met_help at ucar.edu
> >> >> > > > > > > wrote:
> >> >> > > > > >
> >> >> > > > > > > Hi Ryan.
> >> >> > > > > > >
> >> >> > > > > > > Unfortunately, it is the user's responsibility to
get the
> >> >> files
> >> >> > > into
> >> >> > > > a
> >> >> > > > > > > format that MET can read.  I am not familiar with
the
> radar
> >> >> > native
> >> >> > > > .gz
> >> >> > > > > > > files and do not know of any existing conversion
> utilities.
> >> >> > > > > > >
> >> >> > > > > > > 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.
> >> >> > > > > 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!
> >> >> > > > > > >
> >> >> > > > > > > The out.nc file that I provided to you was only
meant as
> >> an
> >> >> > > example
> >> >> > > > of
> >> >> > > > > > the
> >> >> > > > > > > format that MET can read.  It actually contains
> temperature
> >> >> > values
> >> >> > > > and
> >> >> > > > > > not
> >> >> > > > > > > any radar related values.
> >> >> > > > > > >
> >> >> > > > > > > Julie
> >> >> > > > > > >
> >> >> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly via
RT <
> >> >> > > > > met_help at ucar.edu
> >> >> > > > > > >
> >> >> > > > > > > wrote:
> >> >> > > > > > >
> >> >> > > > > > > >
> >> >> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> >> >> ket/Display.html?id=79208
> >> >> > >
> >> >> > > > > > > >
> >> >> > > > > > > > Hi Julie,
> >> >> > > > > > > >
> >> >> > > > > > > > Is there a way to make radar netCDF files that
MODE
> will
> >> >> read
> >> >> > out
> >> >> > > > of
> >> >> > > > > > the
> >> >> > > > > > > > native .gz files?  I like what I've been doing
so far
> >> with
> >> >> radx
> >> >> > > > > since I
> >> >> > > > > > > > want to do my analysis on constant height
surfaces,
> >> rather
> >> >> than
> >> >> > > > worry
> >> >> > > > > > > about
> >> >> > > > > > > > elevation angle, and radx allows me to do that.
> >> >> > > > > > > >
> >> >> > > > > > > > Unrelated to that issue, when I run with your
out.nc
> >> file,
> >> >> it
> >> >> > > > tells
> >> >> > > > > me
> >> >> > > > > > > > that
> >> >> > > > > > > > the GRIB1 field abbreviation isn't recognized.
It's
> for
> >> >> radar
> >> >> > > > > > > reflectivity
> >> >> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL,
REFL_10CM,
> >> etc.,
> >> >> and
> >> >> > > > none
> >> >> > > > > of
> >> >> > > > > > > > those is right, apparently.  Any idea what it is
or
> how I
> >> >> can
> >> >> > > find
> >> >> > > > > out?
> >> >> > > > > > > (I
> >> >> > > > > > > > know for sure that that field was computed when
running
> >> the
> >> >> > UPP,
> >> >> > > > as I
> >> >> > > > > > can
> >> >> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
> >> >> > > > > > > >
> >> >> > > > > > > > Thanks,
> >> >> > > > > > > > Ryan
> >> >> > > > > > > >
> >> >> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie
Prestopnik via
> RT
> >> <
> >> >> > > > > > > > met_help at ucar.edu>
> >> >> > > > > > > > wrote:
> >> >> > > > > > > >
> >> >> > > > > > > > > Hi Ryan.
> >> >> > > > > > > > >
> >> >> > > > > > > > > I see that you're having trouble with a NetCDF
file
> >> using
> >> >> > MODE.
> >> >> > > > > > Thank
> >> >> > > > > > > > you
> >> >> > > > > > > > > for letting us know that you are using a radar
NetCDF
> >> file
> >> >> > > > > processed
> >> >> > > > > > > > > through radx.
> >> >> > > > > > > > >
> >> >> > > > > > > > > Currently, MET has the ability to read three
> different
> >> >> > flavors
> >> >> > > of
> >> >> > > > > > > > NetCDF...
> >> >> > > > > > > > >  - The format used to write intermediate files
in MET
> >> >> (i.e.
> >> >> > the
> >> >> > > > > > output
> >> >> > > > > > > of
> >> >> > > > > > > > > pcp_combine)
> >> >> > > > > > > > >  - NetCDF files following the Climate Forecast
> >> conventions
> >> >> > > (i.e.
> >> >> > > > > > > > > CF-compliant NetCDF)
> >> >> > > > > > > > >  - The output of the pinterp or wrf_interp
utility
> >> >> > > > > > > > >
> >> >> > > > > > > > > I suspect that your file is not in one of
those
> >> formats.
> >> >> > > > > > > > >
> >> >> > > > > > > > > I have attached an example of a NetCDF file
(out.nc)
> >> >> with a
> >> >> > > > > Mercator
> >> >> > > > > > > > > projection that MET can read.  If you can get
your
> >> file in
> >> >> > that
> >> >> > > > > > format,
> >> >> > > > > > > > > MODE should be able to read it.
> >> >> > > > > > > > >
> >> >> > > > > > > > > I hope that helps!
> >> >> > > > > > > > >
> >> >> > > > > > > > > Julie
> >> >> > > > > > > > >
> >> >> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan
Connelly via
> RT
> >> <
> >> >> > > > > > > > met_help at ucar.edu>
> >> >> > > > > > > > > wrote:
> >> >> > > > > > > > >
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208 was
acted
> >> upon.
> >> >> > > > > > > > > > Transaction: Ticket created by
rconne01 at gmail.com
> >> >> > > > > > > > > >        Queue: met_help
> >> >> > > > > > > > > >      Subject: MODE error: Couldn't figure
out
> >> projection
> >> >> > from
> >> >> > > > > > > > information
> >> >> > > > > > > > > > in netCDF file.
> >> >> > > > > > > > > >        Owner: Nobody
> >> >> > > > > > > > > >   Requestors: rconne01 at gmail.com
> >> >> > > > > > > > > >       Status: new
> >> >> > > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> >> >> > > > > > > Ticket/Display.html?id=79208
> >> >> > > > > > > > >
> >> >> > > > > > > > > >
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > Hi,
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > Upon running MODE, I encounter the error:
Couldn't
> >> >> figure
> >> >> > out
> >> >> > > > > > > > projection
> >> >> > > > > > > > > > from information in netCDF file.
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > I assume because it says "in netCDF file"
that it's
> >> >> > referring
> >> >> > > > to
> >> >> > > > > my
> >> >> > > > > > > obs
> >> >> > > > > > > > > > file?  In this case, that is a radar netCDF
> processed
> >> >> > through
> >> >> > > > > radx
> >> >> > > > > > > > called
> >> >> > > > > > > > > > ncf_20141127_010250.nc.
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > I'm pretty sure the projection is Mercator.
Can I
> >> just
> >> >> > hard
> >> >> > > > code
> >> >> > > > > > > that
> >> >> > > > > > > > in
> >> >> > > > > > > > > > somewhere, since having MODE search for the
> >> projection
> >> >> > > doesn't
> >> >> > > > > > work?
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > Thanks,
> >> >> > > > > > > > > > Ryan
> >> >> > > > > > > > > >
> >> >> > > > > > > > > > --
> >> >> > > > > > > > > > Ryan Connelly
> >> >> > > > > > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
> >> >> > University
> >> >> > > > > > > > > > B.S. in Meteorology with Minors in
Mathematics and
> >> GIS,
> >> >> > > > > Valparaiso
> >> >> > > > > > > > > > University
> >> >> > > > > > > > > > rconne01 at gmail.com
> >> >> > > > > > > > > > ryan.connelly at stonybrook.edu
> >> >> > > > > > > > > >
> >> >> > > > > > > > > >
> >> >> > > > > > > > >
> >> >> > > > > > > > >
> >> >> > > > > > > > > --
> >> >> > > > > > > > > Julie Prestopnik
> >> >> > > > > > > > > National Center for Atmospheric Research
> >> >> > > > > > > > > Research Applications Laboratory
> >> >> > > > > > > > > Phone: 303.497.8399
> >> >> > > > > > > > > Email: jpresto at ucar.edu
> >> >> > > > > > > > >
> >> >> > > > > > > > >
> >> >> > > > > > > >
> >> >> > > > > > > >
> >> >> > > > > > > > --
> >> >> > > > > > > > Ryan Connelly
> >> >> > > > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
> >> University
> >> >> > > > > > > > B.S. in Meteorology with Minors in Mathematics
and GIS,
> >> >> > > Valparaiso
> >> >> > > > > > > > University
> >> >> > > > > > > > rconne01 at gmail.com
> >> >> > > > > > > > ryan.connelly at stonybrook.edu
> >> >> > > > > > > >
> >> >> > > > > > > >
> >> >> > > > > > >
> >> >> > > > > > >
> >> >> > > > > > > --
> >> >> > > > > > > Julie Prestopnik
> >> >> > > > > > > National Center for Atmospheric Research
> >> >> > > > > > > Research Applications Laboratory
> >> >> > > > > > > Phone: 303.497.8399
> >> >> > > > > > > Email: jpresto at ucar.edu
> >> >> > > > > > >
> >> >> > > > > > >
> >> >> > > > > >
> >> >> > > > > >
> >> >> > > > > > --
> >> >> > > > > > Ryan Connelly
> >> >> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
> University
> >> >> > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> >> >> Valparaiso
> >> >> > > > > > University
> >> >> > > > > > rconne01 at gmail.com
> >> >> > > > > > ryan.connelly at stonybrook.edu
> >> >> > > > > >
> >> >> > > > > >
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > --
> >> >> > > > > Julie Prestopnik
> >> >> > > > > National Center for Atmospheric Research
> >> >> > > > > Research Applications Laboratory
> >> >> > > > > Phone: 303.497.8399
> >> >> > > > > Email: jpresto at ucar.edu
> >> >> > > > >
> >> >> > > > >
> >> >> > > >
> >> >> > > >
> >> >> > > > --
> >> >> > > > Ryan Connelly
> >> >> > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> >> >> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
> >> Valparaiso
> >> >> > > > University
> >> >> > > > rconne01 at gmail.com
> >> >> > > > ryan.connelly at stonybrook.edu
> >> >> > > >
> >> >> > > >
> >> >> > >
> >> >> > >
> >> >> > > --
> >> >> > > Julie Prestopnik
> >> >> > > National Center for Atmospheric Research
> >> >> > > Research Applications Laboratory
> >> >> > > Phone: 303.497.8399
> >> >> > > Email: jpresto at ucar.edu
> >> >> > >
> >> >> > >
> >> >> > >
> >> >> > >
> >> >> > > --
> >> >> > > Ryan Connelly
> >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> >> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> >> >> > > University
> >> >> > > rconne01 at gmail.com
> >> >> > > ryan.connelly at stonybrook.edu
> >> >> > > --
> >> >> > > Ryan Connelly
> >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> >> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> >> >> > > University
> >> >> > > rconne01 at gmail.com
> >> >> > > ryan.connelly at stonybrook.edu
> >> >> > >
> >> >> > >
> >> >> >
> >> >> >
> >> >> > --
> >> >> > Julie Prestopnik
> >> >> > National Center for Atmospheric Research
> >> >> > Research Applications Laboratory
> >> >> > Phone: 303.497.8399
> >> >> > Email: jpresto at ucar.edu
> >> >> >
> >> >> >
> >> >>
> >> >>
> >> >> --
> >> >> Ryan Connelly
> >> >> M.S. Student in Atmospheric Sciences, Stony Brook University
> >> >> B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> >> >> University
> >> >> rconne01 at gmail.com
> >> >> ryan.connelly at stonybrook.edu
> >> >>
> >> >>
> >> >
> >> >
> >> > --
> >> > 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
> >>
> >>
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
>
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Ryan Connelly
Time: Thu Mar 02 13:30:17 2017

Hi Julie.  No worries.

My problems have been in running pcp_combine.  However, I don't even
remember why I'm doing that anymore.  Instead, I checked with
plot_data_plane, and a Lambert Conformal projection works!  AND I was
able
to succesfully run regrid_data_plane, and get believable output!

So I think we can call the 5.2 patch you guys released successful, and
close this ticket.

I'm sure I'll come up with more questions when running MODE, though.

Thanks,
Ryan

On Thu, Mar 2, 2017 at 12:07 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:

> Hi Ryan.  I apologize for the delay in responding. I was out of the
office.
>
> Could you please upload the files in question to our ftp server and
also
> send me the full command that you are running so that I can try to
> reproduce what you are seeing?
>
> Thanks!
>
> Julie
>
> On Wed, Mar 1, 2017 at 12:16 PM, Ryan Connelly via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> >
> > Hi Julie,
> >
> > I have installed the patch.  However, now unpon running
pcp_combine, I
> > encounter a new error:
> >
> > ERROR  : NcCfFile::get_grid_mapping_azimuthal_equidistant() ->
Azimuthal
> > equidistant grid not handled in MET.
> >
> >
> > This error is repeated for radar netcdf files whose map projection
I have
> > set through radx as (1) cylindrical equidistant, (2) lambert
conformal,
> and
> > (3) mercator, which doesn't seem to make sense to me.  I've
verified
> > through ncdump that these ncf files have the projection:
> >
> > int grid_mapping_0 ;
> >                 grid_mapping_0:grid_mapping_name = "mercator" ;
> >
> > (for example)
> >
> > so I don't really know why it's throwing an error about an
Azimuthal
> > equidistant grid.
> >
> > Thanks,
> > Ryan
> >
> > On Fri, Feb 10, 2017 at 2:41 PM, Ryan Connelly
<rconne01 at gmail.com>
> wrote:
> >
> > > Hi Julie,
> > >
> > > When I gunzip, and then tar -xvf the patches tar file, it
recreates the
> > > /data, /doc, and /src directories, but not /bin, so that all of
the
> tools
> > > are not updated.  What am I doing wrong?
> > >
> > > Ryan
> > >
> > > On Thu, Feb 9, 2017 at 11:52 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Hi Ryan.  We have fixed the code and made a patch available.
Please
> > take
> > >> a
> > >>
> > >> look at the following page and follow the directions for
installing
> the
> > >> patch:
> > >>
> > >> http://www.dtcenter.org/met/users/support/known_issues/
> > METv5.2/index.php
> > >>
> > >> Please let us know if you have any questions.  Thank you for
your
> > patience
> > >> on this issue.
> > >>
> > >> Julie
> > >>
> > >>
> > >>
> > >> On Thu, Feb 2, 2017 at 2:52 PM, Julie Prestopnik
<jpresto at ucar.edu>
> > >> wrote:
> > >>
> > >> > Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.
> > >> >
> > >> > I'll let you know once it's ready.
> > >> >
> > >> > Thanks,
> > >> > Julie
> > >> >
> > >> > On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT <
> > >> met_help at ucar.edu>
> > >> > wrote:
> > >> >
> > >> >>
> > >> >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >> >>
> > >> >> Hi Julie,
> > >> >>
> > >> >> Thanks for trying it out.  It is just a big - not something
I can
> fix
> > >> by
> > >> >> choosing a different map projection or editing/creating more
> > >> attributes or
> > >> >> something like that?  If so, I'm perfectly content to wait
until
> next
> > >> >> week;
> > >> >> I just wanted to see if I could move ahead now already.
> > >> >>
> > >> >> Thanks for your help!
> > >> >> Ryan
> > >> >>
> > >> >> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT <
> > >> >> met_help at ucar.edu>
> > >> >> wrote:
> > >> >>
> > >> >> > Hi Ryan.  Thank you for uploading the file. Unfortunately,
> METv5.2
> > >> had
> > >> >> some
> > >> >> > trouble parsing the file, however, our developmental
version of
> > >> METv6.0
> > >> >> was
> > >> >> > able to parse it.
> > >> >> >
> > >> >> > We will be working on a bug fix for METv5.2 to fix this
issue so
> > that
> > >> >> MET
> > >> >> > will be able to parse the file.  We will let you know as
soon as
> > the
> > >> >> fix is
> > >> >> > ready (likely sometime next week).
> > >> >> >
> > >> >> > Thanks,
> > >> >> > Julie
> > >> >> >
> > >> >> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT <
> > >> >> met_help at ucar.edu>
> > >> >> > wrote:
> > >> >> >
> > >> >> > >
> > >> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208
> >
> > >> >> > >
> > >> >> > > Hi Julie,
> > >> >> > >
> > >> >> > > I remade a netcdf file with Lambert Conformal as the
> projection,
> > >> and
> > >> >> > added
> > >> >> > > the required global attributes as described in this help
list
> > >> email:
> > >> >> > > http://mailman.ucar.edu/pipermail/met_help/2014-
September/
> > >> 002166.html
> > >> >> .
> > >> >> > > However, pcp_combine still throws the error that
projection
> info
> > >> >> can't be
> > >> >> > > found in the file.
> > >> >> > >
> > >> >> > > That file is now in Connelly_data on the ftp site.  Note
that
> it
> > >> is a
> > >> >> > > different time: ncf_20141127_001023.nc
> > >> >> > >
> > >> >> > > Thanks,
> > >> >> > > Ryan
> > >> >> > >
> > >> >> > >
> > >> >> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via RT
<
> > >> >> > > met_help at ucar.edu>
> > >> >> > > wrote:
> > >> >> > >
> > >> >> > > Hi Ryan.
> > >> >> > >
> > >> >> > > Thank you for uploading your files to our ftp server.
It was
> > >> helpful
> > >> >> to
> > >> >> > > have those.
> > >> >> > >
> > >> >> > > Regarding the file:
> > >> >> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > >> >> > Surveillance_SUR.nc
> > >> >> > > It looks like this is CF compliant point data.
Unfortunately,
> > MET
> > >> >> only
> > >> >> > > supports CF compliant gridded data.   Our suggestion
would be
> to
> > >> get
> > >> >> the
> > >> >> > > format in the point observation ascii format as
described on
> this
> > >> >> page of
> > >> >> > > the tutorial:
> > >> >> > >
> > >> >> > >
http://www.dtcenter.org/met/users/support/online_tutorial/
> > >> >> > > METv5.2/tutorial.php?name=ascii2nc&category=index
> > >> >> > >
> > >> >> > > and then use ascii2nc to convert the data to a NetCDF
format
> that
> > >> MET
> > >> >> can
> > >> >> > > read.
> > >> >> > >
> > >> >> > > Regarding the file: ncf_20141127_010250.nc
> > >> >> > > METv5.2 had some trouble parsing the file, however, our
> > >> developmental
> > >> >> > > version of METv6.0 was able to parse it.  Unfortunately,
this
> > file
> > >> is
> > >> >> > using
> > >> >> > > the azimuthal equidistant projection, which is currently
not
> > >> >> supported by
> > >> >> > > MET.  The projections that MET currently supports are
Lat Lon
> > (aka
> > >> >> > > cylindrical equidistant), Lambert Conformal, Mercator,
and
> Polar
> > >> >> > > Stereographic.  So, your file would need to use one of
those
> > >> >> projections
> > >> >> > > before MET could work with it.  If we had some funding
from
> SBU,
> > we
> > >> >> could
> > >> >> > > add support for the azimuthal equidistant projection.
> > >> >> > >
> > >> >> > > Julie
> > >> >> > >
> > >> >> > >
> > >> >> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT <
> > >> >> > met_help at ucar.edu>
> > >> >> > > wrote:
> > >> >> > >
> > >> >> > > >
> > >> >> > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=79208
> > >
> > >> >> > > >
> > >> >> > > > Hi Julie,
> > >> >> > > >
> > >> >> > > > I apologize for the delay.  I was away at the general
> meeting.
> > >> >> > > >
> > >> >> > > > The commands I used to run pcp_combine are:
> > >> >> > > >
> > >> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250
1
> > >> >> > > >
/D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.nc
> > >> -pcpdir
> > >> >> > > > /D2/ryan/Research/20141126_Case/temp
> > >> >> > > >
> > >> >> > > > bin/pcp_combine -sum 00000000_000000 1 20141127_010250
1
> > >> >> > > > /D2/ryan/Research/20141126_Case/temp/
> > >> >> > > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > >> >> > > Surveillance_SUR.nc
> > >> >> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
> > >> >> > > >
> > >> >> > > >
> > >> >> > > > The files - a polar coordinate cfrad and a Cartesian
gridded
> > ncf,
> > >> >> are
> > >> >> > > > located in Connelly_data on the ftp server.
> > >> >> > > >
> > >> >> > > > Thanks,
> > >> >> > > > Ryan
> > >> >> > > >
> > >> >> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik via
RT <
> > >> >> > > > met_help at ucar.edu>
> > >> >> > > > wrote:
> > >> >> > > >
> > >> >> > > > > Hi Ryan.
> > >> >> > > > >
> > >> >> > > > > It would be helpful for us to take a look at your
data.
> > Could
> > >> you
> > >> >> > > please
> > >> >> > > > > upload your data used in the pcp_combine command to
our
> > >> anonymous
> > >> >> ftp
> > >> >> > > > site
> > >> >> > > > > as follows?
> > >> >> > > > >
> > >> >> > > > >    http://www.dtcenter.org/met/
> > users/support/met_help.php#ftp
> > >> >> > > > >
> > >> >> > > > > Could you please also reply with the command that
you used
> to
> > >> run
> > >> >> > > > > pcp_combine?
> > >> >> > > > >
> > >> >> > > > > Thanks,
> > >> >> > > > > Julie
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via
RT <
> > >> >> > > met_help at ucar.edu
> > >> >> > > > >
> > >> >> > > > > wrote:
> > >> >> > > > >
> > >> >> > > > > >
> > >> >> > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> > >> ket/Display.html?id=79208 >
> > >> >> > > > > >
> > >> >> > > > > > Hi Julie,
> > >> >> > > > > >
> > >> >> > > > > > The CfRadial files I have don't work either, which
is
> > >> >> perplexing,
> > >> >> > > > because
> > >> >> > > > > > they should be cf-compliant, right?
> > >> >> > > > > >
> > >> >> > > > > > Anyway, right now I'm trying to run pcp_combine to
use
> the
> > >> >> output
> > >> >> > of
> > >> >> > > it
> > >> >> > > > > as
> > >> >> > > > > > the input for MODE as you suggested in your first
email.
> > But
> > >> >> even
> > >> >> > > that
> > >> >> > > > > is
> > >> >> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid() ->
> Couldn't
> > >> >> figure
> > >> >> > > out
> > >> >> > > > > > projection from information in netCDF file., on
both the
> > >> >> cartesian
> > >> >> > > > > netCDF,
> > >> >> > > > > > and the (supposedly) cf-compliant CfRadial files.
> > >> >> > > > > >
> > >> >> > > > > > An ncdump of the CfRadial file revelas the
following
> > >> projection
> > >> >> > > > > > information:
> > >> >> > > > > >
> > >> >> > > > > > int grid_mapping ;
> > >> >> > > > > >                 grid_mapping:grid_mapping_name =
> > >> >> > > > > "radar_lidar_radial_scan"
> > >> >> > > > > > ;
> > >> >> > > > > >                 grid_mapping:longitude_of_
> > projection_origin
> > >> =
> > >> >> > > > > > -70.2563629150391 ;
> > >> >> > > > > >                 grid_mapping:latitude_of_
> projection_origin
> > =
> > >> >> > > > > > 43.8913040161133 ;
> > >> >> > > > > >                 grid_mapping:altitude_of_
> projection_origin
> > =
> > >> >> 144.
> > >> >> > ;
> > >> >> > > > > >                 grid_mapping:false_northing = 0. ;
> > >> >> > > > > >                 grid_mapping:false_easting = 0. ;
> > >> >> > > > > >         double latitude ;
> > >> >> > > > > >                 latitude:long_name = "latitude" ;
> > >> >> > > > > >                 latitude:units = "degrees_north" ;
> > >> >> > > > > >                 latitude:_FillValue = -9999. ;
> > >> >> > > > > >         double longitude ;
> > >> >> > > > > >                 longitude:long_name = "longitude"
;
> > >> >> > > > > >                 longitude:units = "degrees_east" ;
> > >> >> > > > > >                 longitude:_FillValue = -9999. ;
> > >> >> > > > > >
> > >> >> > > > > > and an ncdump of the cartesian ncf contains:
> > >> >> > > > > >
> > >> >> > > > > > int grid_mapping_0 ;
> > >> >> > > > > >                 grid_mapping_0:grid_mapping_name =
> > >> >> > > > > > "azimuthal_equidistant" ;
> > >> >> > > > > >                 grid_mapping_0:longitude_of_
> > projection_origin
> > >> =
> > >> >> > > > > > -70.25636f ;
> > >> >> > > > > >                 grid_mapping_0:latitude_of_
> > projection_origin
> > >> =
> > >> >> > > > 43.8913f
> > >> >> > > > > ;
> > >> >> > > > > >                 grid_mapping_0:false_easting = 0.f
;
> > >> >> > > > > >                 grid_mapping_0:false_northing =
0.f ;
> > >> >> > > > > >
> > >> >> > > > > > So why is even pcp_combine not able to figure out
> > projection
> > >> >> info?
> > >> >> > > > > >
> > >> >> > > > > > Thanks,
> > >> >> > > > > > Ryan
> > >> >> > > > > >
> > >> >> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie Prestopnik
via
> RT <
> > >> >> > > > > > met_help at ucar.edu
> > >> >> > > > > > > wrote:
> > >> >> > > > > >
> > >> >> > > > > > > Hi Ryan.
> > >> >> > > > > > >
> > >> >> > > > > > > Unfortunately, it is the user's responsibility
to get
> the
> > >> >> files
> > >> >> > > into
> > >> >> > > > a
> > >> >> > > > > > > format that MET can read.  I am not familiar
with the
> > radar
> > >> >> > native
> > >> >> > > > .gz
> > >> >> > > > > > > files and do not know of any existing conversion
> > utilities.
> > >> >> > > > > > >
> > >> >> > > > > > > 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.
> > >> >> > > > > 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!
> > >> >> > > > > > >
> > >> >> > > > > > > The out.nc file that I provided to you was only
meant
> as
> > >> an
> > >> >> > > example
> > >> >> > > > of
> > >> >> > > > > > the
> > >> >> > > > > > > format that MET can read.  It actually contains
> > temperature
> > >> >> > values
> > >> >> > > > and
> > >> >> > > > > > not
> > >> >> > > > > > > any radar related values.
> > >> >> > > > > > >
> > >> >> > > > > > > Julie
> > >> >> > > > > > >
> > >> >> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly
via RT <
> > >> >> > > > > met_help at ucar.edu
> > >> >> > > > > > >
> > >> >> > > > > > > wrote:
> > >> >> > > > > > >
> > >> >> > > > > > > >
> > >> >> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> > >> >> ket/Display.html?id=79208
> > >> >> > >
> > >> >> > > > > > > >
> > >> >> > > > > > > > Hi Julie,
> > >> >> > > > > > > >
> > >> >> > > > > > > > Is there a way to make radar netCDF files that
MODE
> > will
> > >> >> read
> > >> >> > out
> > >> >> > > > of
> > >> >> > > > > > the
> > >> >> > > > > > > > native .gz files?  I like what I've been doing
so far
> > >> with
> > >> >> radx
> > >> >> > > > > since I
> > >> >> > > > > > > > want to do my analysis on constant height
surfaces,
> > >> rather
> > >> >> than
> > >> >> > > > worry
> > >> >> > > > > > > about
> > >> >> > > > > > > > elevation angle, and radx allows me to do
that.
> > >> >> > > > > > > >
> > >> >> > > > > > > > Unrelated to that issue, when I run with your
out.nc
> > >> file,
> > >> >> it
> > >> >> > > > tells
> > >> >> > > > > me
> > >> >> > > > > > > > that
> > >> >> > > > > > > > the GRIB1 field abbreviation isn't recognized.
It's
> > for
> > >> >> radar
> > >> >> > > > > > > reflectivity
> > >> >> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL,
REFL_10CM,
> > >> etc.,
> > >> >> and
> > >> >> > > > none
> > >> >> > > > > of
> > >> >> > > > > > > > those is right, apparently.  Any idea what it
is or
> > how I
> > >> >> can
> > >> >> > > find
> > >> >> > > > > out?
> > >> >> > > > > > > (I
> > >> >> > > > > > > > know for sure that that field was computed
when
> running
> > >> the
> > >> >> > UPP,
> > >> >> > > > as I
> > >> >> > > > > > can
> > >> >> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
> > >> >> > > > > > > >
> > >> >> > > > > > > > Thanks,
> > >> >> > > > > > > > Ryan
> > >> >> > > > > > > >
> > >> >> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie
Prestopnik via
> > RT
> > >> <
> > >> >> > > > > > > > met_help at ucar.edu>
> > >> >> > > > > > > > wrote:
> > >> >> > > > > > > >
> > >> >> > > > > > > > > Hi Ryan.
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > I see that you're having trouble with a
NetCDF file
> > >> using
> > >> >> > MODE.
> > >> >> > > > > > Thank
> > >> >> > > > > > > > you
> > >> >> > > > > > > > > for letting us know that you are using a
radar
> NetCDF
> > >> file
> > >> >> > > > > processed
> > >> >> > > > > > > > > through radx.
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > Currently, MET has the ability to read three
> > different
> > >> >> > flavors
> > >> >> > > of
> > >> >> > > > > > > > NetCDF...
> > >> >> > > > > > > > >  - The format used to write intermediate
files in
> MET
> > >> >> (i.e.
> > >> >> > the
> > >> >> > > > > > output
> > >> >> > > > > > > of
> > >> >> > > > > > > > > pcp_combine)
> > >> >> > > > > > > > >  - NetCDF files following the Climate
Forecast
> > >> conventions
> > >> >> > > (i.e.
> > >> >> > > > > > > > > CF-compliant NetCDF)
> > >> >> > > > > > > > >  - The output of the pinterp or wrf_interp
utility
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > I suspect that your file is not in one of
those
> > >> formats.
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > I have attached an example of a NetCDF file
(
> out.nc)
> > >> >> with a
> > >> >> > > > > Mercator
> > >> >> > > > > > > > > projection that MET can read.  If you can
get your
> > >> file in
> > >> >> > that
> > >> >> > > > > > format,
> > >> >> > > > > > > > > MODE should be able to read it.
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > I hope that helps!
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > Julie
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan
Connelly via
> > RT
> > >> <
> > >> >> > > > > > > > met_help at ucar.edu>
> > >> >> > > > > > > > > wrote:
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208
was acted
> > >> upon.
> > >> >> > > > > > > > > > Transaction: Ticket created by
> rconne01 at gmail.com
> > >> >> > > > > > > > > >        Queue: met_help
> > >> >> > > > > > > > > >      Subject: MODE error: Couldn't figure
out
> > >> projection
> > >> >> > from
> > >> >> > > > > > > > information
> > >> >> > > > > > > > > > in netCDF file.
> > >> >> > > > > > > > > >        Owner: Nobody
> > >> >> > > > > > > > > >   Requestors: rconne01 at gmail.com
> > >> >> > > > > > > > > >       Status: new
> > >> >> > > > > > > > > >  Ticket <URL: https://rt.rap.ucar.edu/rt/
> > >> >> > > > > > > Ticket/Display.html?id=79208
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > Hi,
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > Upon running MODE, I encounter the error:
> Couldn't
> > >> >> figure
> > >> >> > out
> > >> >> > > > > > > > projection
> > >> >> > > > > > > > > > from information in netCDF file.
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > I assume because it says "in netCDF file"
that
> it's
> > >> >> > referring
> > >> >> > > > to
> > >> >> > > > > my
> > >> >> > > > > > > obs
> > >> >> > > > > > > > > > file?  In this case, that is a radar
netCDF
> > processed
> > >> >> > through
> > >> >> > > > > radx
> > >> >> > > > > > > > called
> > >> >> > > > > > > > > > ncf_20141127_010250.nc.
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > I'm pretty sure the projection is
Mercator.  Can
> I
> > >> just
> > >> >> > hard
> > >> >> > > > code
> > >> >> > > > > > > that
> > >> >> > > > > > > > in
> > >> >> > > > > > > > > > somewhere, since having MODE search for
the
> > >> projection
> > >> >> > > doesn't
> > >> >> > > > > > work?
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > Thanks,
> > >> >> > > > > > > > > > Ryan
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > > --
> > >> >> > > > > > > > > > Ryan Connelly
> > >> >> > > > > > > > > > M.S. Student in Atmospheric Sciences,
Stony Brook
> > >> >> > University
> > >> >> > > > > > > > > > B.S. in Meteorology with Minors in
Mathematics
> and
> > >> GIS,
> > >> >> > > > > Valparaiso
> > >> >> > > > > > > > > > University
> > >> >> > > > > > > > > > rconne01 at gmail.com
> > >> >> > > > > > > > > > ryan.connelly at stonybrook.edu
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > > >
> > >> >> > > > > > > > >
> > >> >> > > > > > > > >
> > >> >> > > > > > > > > --
> > >> >> > > > > > > > > Julie Prestopnik
> > >> >> > > > > > > > > National Center for Atmospheric Research
> > >> >> > > > > > > > > Research Applications Laboratory
> > >> >> > > > > > > > > Phone: 303.497.8399
> > >> >> > > > > > > > > Email: jpresto at ucar.edu
> > >> >> > > > > > > > >
> > >> >> > > > > > > > >
> > >> >> > > > > > > >
> > >> >> > > > > > > >
> > >> >> > > > > > > > --
> > >> >> > > > > > > > Ryan Connelly
> > >> >> > > > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
> > >> University
> > >> >> > > > > > > > B.S. in Meteorology with Minors in Mathematics
and
> GIS,
> > >> >> > > Valparaiso
> > >> >> > > > > > > > University
> > >> >> > > > > > > > rconne01 at gmail.com
> > >> >> > > > > > > > ryan.connelly at stonybrook.edu
> > >> >> > > > > > > >
> > >> >> > > > > > > >
> > >> >> > > > > > >
> > >> >> > > > > > >
> > >> >> > > > > > > --
> > >> >> > > > > > > Julie Prestopnik
> > >> >> > > > > > > National Center for Atmospheric Research
> > >> >> > > > > > > Research Applications Laboratory
> > >> >> > > > > > > Phone: 303.497.8399
> > >> >> > > > > > > Email: jpresto at ucar.edu
> > >> >> > > > > > >
> > >> >> > > > > > >
> > >> >> > > > > >
> > >> >> > > > > >
> > >> >> > > > > > --
> > >> >> > > > > > Ryan Connelly
> > >> >> > > > > > M.S. Student in Atmospheric Sciences, Stony Brook
> > University
> > >> >> > > > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > >> >> Valparaiso
> > >> >> > > > > > University
> > >> >> > > > > > rconne01 at gmail.com
> > >> >> > > > > > ryan.connelly at stonybrook.edu
> > >> >> > > > > >
> > >> >> > > > > >
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > --
> > >> >> > > > > Julie Prestopnik
> > >> >> > > > > National Center for Atmospheric Research
> > >> >> > > > > Research Applications Laboratory
> > >> >> > > > > Phone: 303.497.8399
> > >> >> > > > > Email: jpresto at ucar.edu
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > >
> > >> >> > > >
> > >> >> > > > --
> > >> >> > > > Ryan Connelly
> > >> >> > > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > >> >> > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > >> Valparaiso
> > >> >> > > > University
> > >> >> > > > rconne01 at gmail.com
> > >> >> > > > ryan.connelly at stonybrook.edu
> > >> >> > > >
> > >> >> > > >
> > >> >> > >
> > >> >> > >
> > >> >> > > --
> > >> >> > > Julie Prestopnik
> > >> >> > > National Center for Atmospheric Research
> > >> >> > > Research Applications Laboratory
> > >> >> > > Phone: 303.497.8399
> > >> >> > > Email: jpresto at ucar.edu
> > >> >> > >
> > >> >> > >
> > >> >> > >
> > >> >> > >
> > >> >> > > --
> > >> >> > > Ryan Connelly
> > >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > >> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
> > Valparaiso
> > >> >> > > University
> > >> >> > > rconne01 at gmail.com
> > >> >> > > ryan.connelly at stonybrook.edu
> > >> >> > > --
> > >> >> > > Ryan Connelly
> > >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > >> >> > > B.S. in Meteorology with Minors in Mathematics and GIS,
> > Valparaiso
> > >> >> > > University
> > >> >> > > rconne01 at gmail.com
> > >> >> > > ryan.connelly at stonybrook.edu
> > >> >> > >
> > >> >> > >
> > >> >> >
> > >> >> >
> > >> >> > --
> > >> >> > Julie Prestopnik
> > >> >> > National Center for Atmospheric Research
> > >> >> > Research Applications Laboratory
> > >> >> > Phone: 303.497.8399
> > >> >> > Email: jpresto at ucar.edu
> > >> >> >
> > >> >> >
> > >> >>
> > >> >>
> > >> >> --
> > >> >> Ryan Connelly
> > >> >> M.S. Student in Atmospheric Sciences, Stony Brook University
> > >> >> B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > >> >> University
> > >> >> rconne01 at gmail.com
> > >> >> ryan.connelly at stonybrook.edu
> > >> >>
> > >> >>
> > >> >
> > >> >
> > >> > --
> > >> > 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
> > >>
> > >>
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> >
>
>


--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu

------------------------------------------------
Subject: MODE error: Couldn't figure out projection from information in netCDF file.
From: Julie Prestopnik
Time: Thu Mar 02 13:37:02 2017

Thanks, Ryan.  I'll go ahead and close this ticket.   I just got your
follow up "thank you", and you're very welcome!  Please let us know if
you
encounter any problems or have any questions.

Julie

On Thu, Mar 2, 2017 at 1:30 PM, Ryan Connelly via RT
<met_help at ucar.edu>
wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
>
> Hi Julie.  No worries.
>
> My problems have been in running pcp_combine.  However, I don't even
> remember why I'm doing that anymore.  Instead, I checked with
> plot_data_plane, and a Lambert Conformal projection works!  AND I
was able
> to succesfully run regrid_data_plane, and get believable output!
>
> So I think we can call the 5.2 patch you guys released successful,
and
> close this ticket.
>
> I'm sure I'll come up with more questions when running MODE, though.
>
> Thanks,
> Ryan
>
> On Thu, Mar 2, 2017 at 12:07 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Ryan.  I apologize for the delay in responding. I was out of
the
> office.
> >
> > Could you please upload the files in question to our ftp server
and also
> > send me the full command that you are running so that I can try to
> > reproduce what you are seeing?
> >
> > Thanks!
> >
> > Julie
> >
> > On Wed, Mar 1, 2017 at 12:16 PM, Ryan Connelly via RT
<met_help at ucar.edu
> >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > >
> > > Hi Julie,
> > >
> > > I have installed the patch.  However, now unpon running
pcp_combine, I
> > > encounter a new error:
> > >
> > > ERROR  : NcCfFile::get_grid_mapping_azimuthal_equidistant() ->
> Azimuthal
> > > equidistant grid not handled in MET.
> > >
> > >
> > > This error is repeated for radar netcdf files whose map
projection I
> have
> > > set through radx as (1) cylindrical equidistant, (2) lambert
conformal,
> > and
> > > (3) mercator, which doesn't seem to make sense to me.  I've
verified
> > > through ncdump that these ncf files have the projection:
> > >
> > > int grid_mapping_0 ;
> > >                 grid_mapping_0:grid_mapping_name = "mercator" ;
> > >
> > > (for example)
> > >
> > > so I don't really know why it's throwing an error about an
Azimuthal
> > > equidistant grid.
> > >
> > > Thanks,
> > > Ryan
> > >
> > > On Fri, Feb 10, 2017 at 2:41 PM, Ryan Connelly
<rconne01 at gmail.com>
> > wrote:
> > >
> > > > Hi Julie,
> > > >
> > > > When I gunzip, and then tar -xvf the patches tar file, it
recreates
> the
> > > > /data, /doc, and /src directories, but not /bin, so that all
of the
> > tools
> > > > are not updated.  What am I doing wrong?
> > > >
> > > > Ryan
> > > >
> > > > On Thu, Feb 9, 2017 at 11:52 AM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> Hi Ryan.  We have fixed the code and made a patch available.
Please
> > > take
> > > >> a
> > > >>
> > > >> look at the following page and follow the directions for
installing
> > the
> > > >> patch:
> > > >>
> > > >> http://www.dtcenter.org/met/users/support/known_issues/
> > > METv5.2/index.php
> > > >>
> > > >> Please let us know if you have any questions.  Thank you for
your
> > > patience
> > > >> on this issue.
> > > >>
> > > >> Julie
> > > >>
> > > >>
> > > >>
> > > >> On Thu, Feb 2, 2017 at 2:52 PM, Julie Prestopnik
<jpresto at ucar.edu>
> > > >> wrote:
> > > >>
> > > >> > Hi Ryan.  Yes, it is a bug we need to fix, unfortunately.
> > > >> >
> > > >> > I'll let you know once it's ready.
> > > >> >
> > > >> > Thanks,
> > > >> > Julie
> > > >> >
> > > >> > On Thu, Feb 2, 2017 at 11:40 AM, Ryan Connelly via RT <
> > > >> met_help at ucar.edu>
> > > >> > wrote:
> > > >> >
> > > >> >>
> > > >> >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=79208 >
> > > >> >>
> > > >> >> Hi Julie,
> > > >> >>
> > > >> >> Thanks for trying it out.  It is just a big - not
something I can
> > fix
> > > >> by
> > > >> >> choosing a different map projection or editing/creating
more
> > > >> attributes or
> > > >> >> something like that?  If so, I'm perfectly content to wait
until
> > next
> > > >> >> week;
> > > >> >> I just wanted to see if I could move ahead now already.
> > > >> >>
> > > >> >> Thanks for your help!
> > > >> >> Ryan
> > > >> >>
> > > >> >> On Thu, Feb 2, 2017 at 1:34 PM, Julie Prestopnik via RT <
> > > >> >> met_help at ucar.edu>
> > > >> >> wrote:
> > > >> >>
> > > >> >> > Hi Ryan.  Thank you for uploading the file.
Unfortunately,
> > METv5.2
> > > >> had
> > > >> >> some
> > > >> >> > trouble parsing the file, however, our developmental
version of
> > > >> METv6.0
> > > >> >> was
> > > >> >> > able to parse it.
> > > >> >> >
> > > >> >> > We will be working on a bug fix for METv5.2 to fix this
issue
> so
> > > that
> > > >> >> MET
> > > >> >> > will be able to parse the file.  We will let you know as
soon
> as
> > > the
> > > >> >> fix is
> > > >> >> > ready (likely sometime next week).
> > > >> >> >
> > > >> >> > Thanks,
> > > >> >> > Julie
> > > >> >> >
> > > >> >> > On Wed, Feb 1, 2017 at 12:25 PM, Ryan Connelly via RT <
> > > >> >> met_help at ucar.edu>
> > > >> >> > wrote:
> > > >> >> >
> > > >> >> > >
> > > >> >> > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=79208
> > >
> > > >> >> > >
> > > >> >> > > Hi Julie,
> > > >> >> > >
> > > >> >> > > I remade a netcdf file with Lambert Conformal as the
> > projection,
> > > >> and
> > > >> >> > added
> > > >> >> > > the required global attributes as described in this
help list
> > > >> email:
> > > >> >> > > http://mailman.ucar.edu/pipermail/met_help/2014-
September/
> > > >> 002166.html
> > > >> >> .
> > > >> >> > > However, pcp_combine still throws the error that
projection
> > info
> > > >> >> can't be
> > > >> >> > > found in the file.
> > > >> >> > >
> > > >> >> > > That file is now in Connelly_data on the ftp site.
Note that
> > it
> > > >> is a
> > > >> >> > > different time: ncf_20141127_001023.nc
> > > >> >> > >
> > > >> >> > > Thanks,
> > > >> >> > > Ryan
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > On Tue, Jan 31, 2017 at 2:23 PM, Julie Prestopnik via
RT <
> > > >> >> > > met_help at ucar.edu>
> > > >> >> > > wrote:
> > > >> >> > >
> > > >> >> > > Hi Ryan.
> > > >> >> > >
> > > >> >> > > Thank you for uploading your files to our ftp server.
It was
> > > >> helpful
> > > >> >> to
> > > >> >> > > have those.
> > > >> >> > >
> > > >> >> > > Regarding the file:
> > > >> >> > > cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > > >> >> > Surveillance_SUR.nc
> > > >> >> > > It looks like this is CF compliant point data.
> Unfortunately,
> > > MET
> > > >> >> only
> > > >> >> > > supports CF compliant gridded data.   Our suggestion
would be
> > to
> > > >> get
> > > >> >> the
> > > >> >> > > format in the point observation ascii format as
described on
> > this
> > > >> >> page of
> > > >> >> > > the tutorial:
> > > >> >> > >
> > > >> >> > >
http://www.dtcenter.org/met/users/support/online_tutorial/
> > > >> >> > > METv5.2/tutorial.php?name=ascii2nc&category=index
> > > >> >> > >
> > > >> >> > > and then use ascii2nc to convert the data to a NetCDF
format
> > that
> > > >> MET
> > > >> >> can
> > > >> >> > > read.
> > > >> >> > >
> > > >> >> > > Regarding the file: ncf_20141127_010250.nc
> > > >> >> > > METv5.2 had some trouble parsing the file, however,
our
> > > >> developmental
> > > >> >> > > version of METv6.0 was able to parse it.
Unfortunately, this
> > > file
> > > >> is
> > > >> >> > using
> > > >> >> > > the azimuthal equidistant projection, which is
currently not
> > > >> >> supported by
> > > >> >> > > MET.  The projections that MET currently supports are
Lat Lon
> > > (aka
> > > >> >> > > cylindrical equidistant), Lambert Conformal, Mercator,
and
> > Polar
> > > >> >> > > Stereographic.  So, your file would need to use one of
those
> > > >> >> projections
> > > >> >> > > before MET could work with it.  If we had some funding
from
> > SBU,
> > > we
> > > >> >> could
> > > >> >> > > add support for the azimuthal equidistant projection.
> > > >> >> > >
> > > >> >> > > Julie
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > On Mon, Jan 30, 2017 at 11:36 AM, Ryan Connelly via RT
<
> > > >> >> > met_help at ucar.edu>
> > > >> >> > > wrote:
> > > >> >> > >
> > > >> >> > > >
> > > >> >> > > > <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=79208
> > > >
> > > >> >> > > >
> > > >> >> > > > Hi Julie,
> > > >> >> > > >
> > > >> >> > > > I apologize for the delay.  I was away at the
general
> > meeting.
> > > >> >> > > >
> > > >> >> > > > The commands I used to run pcp_combine are:
> > > >> >> > > >
> > > >> >> > > > bin/pcp_combine -sum 00000000_000000 1
20141127_010250 1
> > > >> >> > > >
/D2/ryan/Research/20141126_Case/temp/ncf_20141127_010250.
> nc
> > > >> -pcpdir
> > > >> >> > > > /D2/ryan/Research/20141126_Case/temp
> > > >> >> > > >
> > > >> >> > > > bin/pcp_combine -sum 00000000_000000 1
20141127_010250 1
> > > >> >> > > > /D2/ryan/Research/20141126_Case/temp/
> > > >> >> > > >
cfrad.20141127_005742.193_to_20141127_010250.191_KGYX_
> > > >> >> > > Surveillance_SUR.nc
> > > >> >> > > > -pcpdir /D2/ryan/Research/20141126_Case/temp
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > > The files - a polar coordinate cfrad and a Cartesian
> gridded
> > > ncf,
> > > >> >> are
> > > >> >> > > > located in Connelly_data on the ftp server.
> > > >> >> > > >
> > > >> >> > > > Thanks,
> > > >> >> > > > Ryan
> > > >> >> > > >
> > > >> >> > > > On Fri, Jan 20, 2017 at 5:28 PM, Julie Prestopnik
via RT <
> > > >> >> > > > met_help at ucar.edu>
> > > >> >> > > > wrote:
> > > >> >> > > >
> > > >> >> > > > > Hi Ryan.
> > > >> >> > > > >
> > > >> >> > > > > It would be helpful for us to take a look at your
data.
> > > Could
> > > >> you
> > > >> >> > > please
> > > >> >> > > > > upload your data used in the pcp_combine command
to our
> > > >> anonymous
> > > >> >> ftp
> > > >> >> > > > site
> > > >> >> > > > > as follows?
> > > >> >> > > > >
> > > >> >> > > > >    http://www.dtcenter.org/met/
> > > users/support/met_help.php#ftp
> > > >> >> > > > >
> > > >> >> > > > > Could you please also reply with the command that
you
> used
> > to
> > > >> run
> > > >> >> > > > > pcp_combine?
> > > >> >> > > > >
> > > >> >> > > > > Thanks,
> > > >> >> > > > > Julie
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > On Fri, Jan 20, 2017 at 2:21 PM, Ryan Connelly via
RT <
> > > >> >> > > met_help at ucar.edu
> > > >> >> > > > >
> > > >> >> > > > > wrote:
> > > >> >> > > > >
> > > >> >> > > > > >
> > > >> >> > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> > > >> ket/Display.html?id=79208 >
> > > >> >> > > > > >
> > > >> >> > > > > > Hi Julie,
> > > >> >> > > > > >
> > > >> >> > > > > > The CfRadial files I have don't work either,
which is
> > > >> >> perplexing,
> > > >> >> > > > because
> > > >> >> > > > > > they should be cf-compliant, right?
> > > >> >> > > > > >
> > > >> >> > > > > > Anyway, right now I'm trying to run pcp_combine
to use
> > the
> > > >> >> output
> > > >> >> > of
> > > >> >> > > it
> > > >> >> > > > > as
> > > >> >> > > > > > the input for MODE as you suggested in your
first
> email.
> > > But
> > > >> >> even
> > > >> >> > > that
> > > >> >> > > > > is
> > > >> >> > > > > > returning ERROR  : NcCfFile::read_netcdf_grid()
->
> > Couldn't
> > > >> >> figure
> > > >> >> > > out
> > > >> >> > > > > > projection from information in netCDF file., on
both
> the
> > > >> >> cartesian
> > > >> >> > > > > netCDF,
> > > >> >> > > > > > and the (supposedly) cf-compliant CfRadial
files.
> > > >> >> > > > > >
> > > >> >> > > > > > An ncdump of the CfRadial file revelas the
following
> > > >> projection
> > > >> >> > > > > > information:
> > > >> >> > > > > >
> > > >> >> > > > > > int grid_mapping ;
> > > >> >> > > > > >                 grid_mapping:grid_mapping_name =
> > > >> >> > > > > "radar_lidar_radial_scan"
> > > >> >> > > > > > ;
> > > >> >> > > > > >                 grid_mapping:longitude_of_
> > > projection_origin
> > > >> =
> > > >> >> > > > > > -70.2563629150391 ;
> > > >> >> > > > > >                 grid_mapping:latitude_of_
> > projection_origin
> > > =
> > > >> >> > > > > > 43.8913040161133 ;
> > > >> >> > > > > >                 grid_mapping:altitude_of_
> > projection_origin
> > > =
> > > >> >> 144.
> > > >> >> > ;
> > > >> >> > > > > >                 grid_mapping:false_northing = 0.
;
> > > >> >> > > > > >                 grid_mapping:false_easting = 0.
;
> > > >> >> > > > > >         double latitude ;
> > > >> >> > > > > >                 latitude:long_name = "latitude"
;
> > > >> >> > > > > >                 latitude:units = "degrees_north"
;
> > > >> >> > > > > >                 latitude:_FillValue = -9999. ;
> > > >> >> > > > > >         double longitude ;
> > > >> >> > > > > >                 longitude:long_name =
"longitude" ;
> > > >> >> > > > > >                 longitude:units = "degrees_east"
;
> > > >> >> > > > > >                 longitude:_FillValue = -9999. ;
> > > >> >> > > > > >
> > > >> >> > > > > > and an ncdump of the cartesian ncf contains:
> > > >> >> > > > > >
> > > >> >> > > > > > int grid_mapping_0 ;
> > > >> >> > > > > >                 grid_mapping_0:grid_mapping_name
=
> > > >> >> > > > > > "azimuthal_equidistant" ;
> > > >> >> > > > > >                 grid_mapping_0:longitude_of_
> > > projection_origin
> > > >> =
> > > >> >> > > > > > -70.25636f ;
> > > >> >> > > > > >                 grid_mapping_0:latitude_of_
> > > projection_origin
> > > >> =
> > > >> >> > > > 43.8913f
> > > >> >> > > > > ;
> > > >> >> > > > > >                 grid_mapping_0:false_easting =
0.f ;
> > > >> >> > > > > >                 grid_mapping_0:false_northing =
0.f ;
> > > >> >> > > > > >
> > > >> >> > > > > > So why is even pcp_combine not able to figure
out
> > > projection
> > > >> >> info?
> > > >> >> > > > > >
> > > >> >> > > > > > Thanks,
> > > >> >> > > > > > Ryan
> > > >> >> > > > > >
> > > >> >> > > > > > On Fri, Jan 20, 2017 at 11:24 AM, Julie
Prestopnik via
> > RT <
> > > >> >> > > > > > met_help at ucar.edu
> > > >> >> > > > > > > wrote:
> > > >> >> > > > > >
> > > >> >> > > > > > > Hi Ryan.
> > > >> >> > > > > > >
> > > >> >> > > > > > > Unfortunately, it is the user's responsibility
to get
> > the
> > > >> >> files
> > > >> >> > > into
> > > >> >> > > > a
> > > >> >> > > > > > > format that MET can read.  I am not familiar
with the
> > > radar
> > > >> >> > native
> > > >> >> > > > .gz
> > > >> >> > > > > > > files and do not know of any existing
conversion
> > > utilities.
> > > >> >> > > > > > >
> > > >> >> > > > > > > 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.
> > > >> >> > > > > 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!
> > > >> >> > > > > > >
> > > >> >> > > > > > > The out.nc file that I provided to you was
only
> meant
> > as
> > > >> an
> > > >> >> > > example
> > > >> >> > > > of
> > > >> >> > > > > > the
> > > >> >> > > > > > > format that MET can read.  It actually
contains
> > > temperature
> > > >> >> > values
> > > >> >> > > > and
> > > >> >> > > > > > not
> > > >> >> > > > > > > any radar related values.
> > > >> >> > > > > > >
> > > >> >> > > > > > > Julie
> > > >> >> > > > > > >
> > > >> >> > > > > > > On Thu, Jan 19, 2017 at 4:25 PM, Ryan Connelly
via
> RT <
> > > >> >> > > > > met_help at ucar.edu
> > > >> >> > > > > > >
> > > >> >> > > > > > > wrote:
> > > >> >> > > > > > >
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/Tic
> > > >> >> ket/Display.html?id=79208
> > > >> >> > >
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > Hi Julie,
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > Is there a way to make radar netCDF files
that MODE
> > > will
> > > >> >> read
> > > >> >> > out
> > > >> >> > > > of
> > > >> >> > > > > > the
> > > >> >> > > > > > > > native .gz files?  I like what I've been
doing so
> far
> > > >> with
> > > >> >> radx
> > > >> >> > > > > since I
> > > >> >> > > > > > > > want to do my analysis on constant height
surfaces,
> > > >> rather
> > > >> >> than
> > > >> >> > > > worry
> > > >> >> > > > > > > about
> > > >> >> > > > > > > > elevation angle, and radx allows me to do
that.
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > Unrelated to that issue, when I run with
your
> out.nc
> > > >> file,
> > > >> >> it
> > > >> >> > > > tells
> > > >> >> > > > > me
> > > >> >> > > > > > > > that
> > > >> >> > > > > > > > the GRIB1 field abbreviation isn't
recognized.
> It's
> > > for
> > > >> >> radar
> > > >> >> > > > > > > reflectivity
> > > >> >> > > > > > > > AGL; I've tried REFD, RADAR REF AGL, REFL,
> REFL_10CM,
> > > >> etc.,
> > > >> >> and
> > > >> >> > > > none
> > > >> >> > > > > of
> > > >> >> > > > > > > > those is right, apparently.  Any idea what
it is or
> > > how I
> > > >> >> can
> > > >> >> > > find
> > > >> >> > > > > out?
> > > >> >> > > > > > > (I
> > > >> >> > > > > > > > know for sure that that field was computed
when
> > running
> > > >> the
> > > >> >> > UPP,
> > > >> >> > > > as I
> > > >> >> > > > > > can
> > > >> >> > > > > > > > see the field REFD in my gdinfo in GEMPAK.)
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > Thanks,
> > > >> >> > > > > > > > Ryan
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > On Thu, Jan 19, 2017 at 2:54 PM, Julie
Prestopnik
> via
> > > RT
> > > >> <
> > > >> >> > > > > > > > met_help at ucar.edu>
> > > >> >> > > > > > > > wrote:
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > > Hi Ryan.
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > I see that you're having trouble with a
NetCDF
> file
> > > >> using
> > > >> >> > MODE.
> > > >> >> > > > > > Thank
> > > >> >> > > > > > > > you
> > > >> >> > > > > > > > > for letting us know that you are using a
radar
> > NetCDF
> > > >> file
> > > >> >> > > > > processed
> > > >> >> > > > > > > > > through radx.
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > Currently, MET has the ability to read
three
> > > different
> > > >> >> > flavors
> > > >> >> > > of
> > > >> >> > > > > > > > NetCDF...
> > > >> >> > > > > > > > >  - The format used to write intermediate
files in
> > MET
> > > >> >> (i.e.
> > > >> >> > the
> > > >> >> > > > > > output
> > > >> >> > > > > > > of
> > > >> >> > > > > > > > > pcp_combine)
> > > >> >> > > > > > > > >  - NetCDF files following the Climate
Forecast
> > > >> conventions
> > > >> >> > > (i.e.
> > > >> >> > > > > > > > > CF-compliant NetCDF)
> > > >> >> > > > > > > > >  - The output of the pinterp or wrf_interp
> utility
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > I suspect that your file is not in one of
those
> > > >> formats.
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > I have attached an example of a NetCDF
file (
> > out.nc)
> > > >> >> with a
> > > >> >> > > > > Mercator
> > > >> >> > > > > > > > > projection that MET can read.  If you can
get
> your
> > > >> file in
> > > >> >> > that
> > > >> >> > > > > > format,
> > > >> >> > > > > > > > > MODE should be able to read it.
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > I hope that helps!
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > Julie
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > On Thu, Jan 19, 2017 at 12:15 PM, Ryan
Connelly
> via
> > > RT
> > > >> <
> > > >> >> > > > > > > > met_help at ucar.edu>
> > > >> >> > > > > > > > > wrote:
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > Thu Jan 19 12:15:19 2017: Request 79208
was
> acted
> > > >> upon.
> > > >> >> > > > > > > > > > Transaction: Ticket created by
> > rconne01 at gmail.com
> > > >> >> > > > > > > > > >        Queue: met_help
> > > >> >> > > > > > > > > >      Subject: MODE error: Couldn't
figure out
> > > >> projection
> > > >> >> > from
> > > >> >> > > > > > > > information
> > > >> >> > > > > > > > > > in netCDF file.
> > > >> >> > > > > > > > > >        Owner: Nobody
> > > >> >> > > > > > > > > >   Requestors: rconne01 at gmail.com
> > > >> >> > > > > > > > > >       Status: new
> > > >> >> > > > > > > > > >  Ticket <URL:
https://rt.rap.ucar.edu/rt/
> > > >> >> > > > > > > Ticket/Display.html?id=79208
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > Hi,
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > Upon running MODE, I encounter the
error:
> > Couldn't
> > > >> >> figure
> > > >> >> > out
> > > >> >> > > > > > > > projection
> > > >> >> > > > > > > > > > from information in netCDF file.
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > I assume because it says "in netCDF
file" that
> > it's
> > > >> >> > referring
> > > >> >> > > > to
> > > >> >> > > > > my
> > > >> >> > > > > > > obs
> > > >> >> > > > > > > > > > file?  In this case, that is a radar
netCDF
> > > processed
> > > >> >> > through
> > > >> >> > > > > radx
> > > >> >> > > > > > > > called
> > > >> >> > > > > > > > > > ncf_20141127_010250.nc.
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > I'm pretty sure the projection is
Mercator.
> Can
> > I
> > > >> just
> > > >> >> > hard
> > > >> >> > > > code
> > > >> >> > > > > > > that
> > > >> >> > > > > > > > in
> > > >> >> > > > > > > > > > somewhere, since having MODE search for
the
> > > >> projection
> > > >> >> > > doesn't
> > > >> >> > > > > > work?
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > Thanks,
> > > >> >> > > > > > > > > > Ryan
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > > --
> > > >> >> > > > > > > > > > Ryan Connelly
> > > >> >> > > > > > > > > > M.S. Student in Atmospheric Sciences,
Stony
> Brook
> > > >> >> > University
> > > >> >> > > > > > > > > > B.S. in Meteorology with Minors in
Mathematics
> > and
> > > >> GIS,
> > > >> >> > > > > Valparaiso
> > > >> >> > > > > > > > > > University
> > > >> >> > > > > > > > > > rconne01 at gmail.com
> > > >> >> > > > > > > > > > ryan.connelly at stonybrook.edu
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > > >
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > > --
> > > >> >> > > > > > > > > Julie Prestopnik
> > > >> >> > > > > > > > > National Center for Atmospheric Research
> > > >> >> > > > > > > > > Research Applications Laboratory
> > > >> >> > > > > > > > > Phone: 303.497.8399
> > > >> >> > > > > > > > > Email: jpresto at ucar.edu
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > > >
> > > >> >> > > > > > > >
> > > >> >> > > > > > > >
> > > >> >> > > > > > > > --
> > > >> >> > > > > > > > Ryan Connelly
> > > >> >> > > > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
> > > >> University
> > > >> >> > > > > > > > B.S. in Meteorology with Minors in
Mathematics and
> > GIS,
> > > >> >> > > Valparaiso
> > > >> >> > > > > > > > University
> > > >> >> > > > > > > > rconne01 at gmail.com
> > > >> >> > > > > > > > ryan.connelly at stonybrook.edu
> > > >> >> > > > > > > >
> > > >> >> > > > > > > >
> > > >> >> > > > > > >
> > > >> >> > > > > > >
> > > >> >> > > > > > > --
> > > >> >> > > > > > > Julie Prestopnik
> > > >> >> > > > > > > National Center for Atmospheric Research
> > > >> >> > > > > > > Research Applications Laboratory
> > > >> >> > > > > > > Phone: 303.497.8399
> > > >> >> > > > > > > Email: jpresto at ucar.edu
> > > >> >> > > > > > >
> > > >> >> > > > > > >
> > > >> >> > > > > >
> > > >> >> > > > > >
> > > >> >> > > > > > --
> > > >> >> > > > > > Ryan Connelly
> > > >> >> > > > > > M.S. Student in Atmospheric Sciences, Stony
Brook
> > > University
> > > >> >> > > > > > B.S. in Meteorology with Minors in Mathematics
and GIS,
> > > >> >> Valparaiso
> > > >> >> > > > > > University
> > > >> >> > > > > > rconne01 at gmail.com
> > > >> >> > > > > > ryan.connelly at stonybrook.edu
> > > >> >> > > > > >
> > > >> >> > > > > >
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > --
> > > >> >> > > > > Julie Prestopnik
> > > >> >> > > > > National Center for Atmospheric Research
> > > >> >> > > > > Research Applications Laboratory
> > > >> >> > > > > Phone: 303.497.8399
> > > >> >> > > > > Email: jpresto at ucar.edu
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > > --
> > > >> >> > > > Ryan Connelly
> > > >> >> > > > M.S. Student in Atmospheric Sciences, Stony Brook
> University
> > > >> >> > > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > >> Valparaiso
> > > >> >> > > > University
> > > >> >> > > > rconne01 at gmail.com
> > > >> >> > > > ryan.connelly at stonybrook.edu
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > --
> > > >> >> > > Julie Prestopnik
> > > >> >> > > National Center for Atmospheric Research
> > > >> >> > > Research Applications Laboratory
> > > >> >> > > Phone: 303.497.8399
> > > >> >> > > Email: jpresto at ucar.edu
> > > >> >> > >
> > > >> >> > >
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > --
> > > >> >> > > Ryan Connelly
> > > >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > >> >> > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > Valparaiso
> > > >> >> > > University
> > > >> >> > > rconne01 at gmail.com
> > > >> >> > > ryan.connelly at stonybrook.edu
> > > >> >> > > --
> > > >> >> > > Ryan Connelly
> > > >> >> > > M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > >> >> > > B.S. in Meteorology with Minors in Mathematics and
GIS,
> > > Valparaiso
> > > >> >> > > University
> > > >> >> > > rconne01 at gmail.com
> > > >> >> > > ryan.connelly at stonybrook.edu
> > > >> >> > >
> > > >> >> > >
> > > >> >> >
> > > >> >> >
> > > >> >> > --
> > > >> >> > Julie Prestopnik
> > > >> >> > National Center for Atmospheric Research
> > > >> >> > Research Applications Laboratory
> > > >> >> > Phone: 303.497.8399
> > > >> >> > Email: jpresto at ucar.edu
> > > >> >> >
> > > >> >> >
> > > >> >>
> > > >> >>
> > > >> >> --
> > > >> >> Ryan Connelly
> > > >> >> M.S. Student in Atmospheric Sciences, Stony Brook
University
> > > >> >> B.S. in Meteorology with Minors in Mathematics and GIS,
> Valparaiso
> > > >> >> University
> > > >> >> rconne01 at gmail.com
> > > >> >> ryan.connelly at stonybrook.edu
> > > >> >>
> > > >> >>
> > > >> >
> > > >> >
> > > >> > --
> > > >> > 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
> > > >>
> > > >>
> > > >
> > > >
> > > > --
> > > > Ryan Connelly
> > > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > > University
> > > > rconne01 at gmail.com
> > > > ryan.connelly at stonybrook.edu
> > > >
> > >
> > >
> > >
> > > --
> > > Ryan Connelly
> > > M.S. Student in Atmospheric Sciences, Stony Brook University
> > > B.S. in Meteorology with Minors in Mathematics and GIS,
Valparaiso
> > > University
> > > rconne01 at gmail.com
> > > ryan.connelly at stonybrook.edu
> > >
> > >
> >
> >
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
>

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


More information about the Met_help mailing list