[Met_help] [rt.rap.ucar.edu #98227] History for plot_data_plane and NetCDF input file
George McCabe via RT
met_help at ucar.edu
Wed Jan 20 08:43:32 MST 2021
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I am having trouble using plot_data_plane with a NetCDF file with the header information below: plot_data_plane 2020080100_ENAAPS_total_aod_001.nc total_aod_001.ps 'name="total_aod"; level="(*,*)"; start_time="20200801_00"; end_time="20200801_06";' -log pdp.log -v 4. The error is:
DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new Met2dDataFile object of type "FileType_NcCF".
DEBUG 4: NcCfFile::open() -> parsing units for the time variable "seconds since 2020-08-01 00:00:00"
DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time unit string "seconds since 2020-08-01 00:00:00" as a reference time of 20200801_000000 and 1 second(s) per time step.
DEBUG 4: NcCfFile::open() -> could not extract init time from the "forecast_reference_time" variable.
DEBUG 4: NcCfFile::open() -> could not extract init time from file name.
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of type "FileType_NcCF".
WARNING:
WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) -> returns the first available time for "total_aod" variable).
WARNING:
ERROR :
ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &) const -> needed 4 arguments for variable total_aod, got 2
ERROR :
netcdf \2020080100_ENAAPS_total_aod_001 {
dimensions:
copy = 1 ;
lat = 180 ;
lon = 360 ;
time = 21 ;
variables:
int copy(copy) ;
copy:long_name = "[Ensemble Mean, Ensemble StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
float lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
float lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
int time(time) ;
time:units = "seconds since 2020-08-01 00:00:00" ;
time:standard_name = "time" ;
float total_aod(time, copy, lat, lon) ;
// global attributes:
:Conventions = "CF-1.0" ;
:history = "Mon Jan 11 18:49:58 2021: ncks -v total_aod,copy -d copy,7 /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc" ;
:NCO = "netCDF Operators version 4.7.7 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)" ;
}
Efren A. Serra (Contractor)
Physicist
DeVine Consulting, Inc.
Naval Research Laboratory
Marine Meteorology Division
7 Grace Hopper Ave., STOP 2
Monterey, CA 93943
Code 7542
Mobile: 408-425-5027
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: plot_data_plane and NetCDF input file
From: George McCabe
Time: Wed Jan 13 13:16:13 2021
Hi Efren,
I see you are having trouble reading in NetCDF data with Plot-Data-
Plane.
It looks like there are 4 dimensions for the total_aod field. The
level you
specified only contains 2 dimensions.
level="(*,*)"
You will need to specify dimensions that match the data. The first
dimension is time (size 21) and the 2nd dimension is copy (size 1).
Also,
you are only able to read a single time into Plot-Data-Plane, so you
should
remove the start_time and end_time values.
If you want to read the first time in the file, you could set
level="(0, 0, *,*)"
If you want MET to find the appropriate time for you, you can specify
the
time value in the time dimension. I believe the syntax would be:
level="(20200801_00, 0, *,*)"
Let me know if you continue to have issues and I can take a closer
look.
Thanks,
George
On Wed, Jan 13, 2021 at 12:57 PM efren.serra.ctr at nrlmry.navy.mil via
RT <
met_help at ucar.edu> wrote:
>
> Wed Jan 13 12:56:53 2021: Request 98227 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> Queue: met_help
> Subject: plot_data_plane and NetCDF input file
> Owner: Nobody
> Requestors: efren.serra.ctr at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227 >
>
>
> I am having trouble using plot_data_plane with a NetCDF file with
the
> header information below: plot_data_plane
2020080100_ENAAPS_total_aod_
> 001.nc total_aod_001.ps 'name="total_aod"; level="(*,*)";
> start_time="20200801_00"; end_time="20200801_06";' -log pdp.log -v
4. The
> error is:
>
> DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcCF".
> DEBUG 4: NcCfFile::open() -> parsing units for the time variable
"seconds
> since 2020-08-01 00:00:00"
> DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time
unit
> string "seconds since 2020-08-01 00:00:00" as a reference time of
> 20200801_000000 and 1 second(s) per time step.
> DEBUG 4: NcCfFile::open() -> could not extract init time from the
> "forecast_reference_time" variable.
> DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcCF".
> WARNING:
> WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) ->
returns
> the first available time for "total_aod" variable).
> WARNING:
> ERROR :
> ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &)
const
> -> needed 4 arguments for variable total_aod, got 2
> ERROR :
>
> netcdf \2020080100_ENAAPS_total_aod_001 {
> dimensions:
> copy = 1 ;
> lat = 180 ;
> lon = 360 ;
> time = 21 ;
> variables:
> int copy(copy) ;
> copy:long_name = "[Ensemble Mean, Ensemble
> StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
> float lat(lat) ;
> lat:units = "degrees_north" ;
> lat:long_name = "latitude" ;
> float lon(lon) ;
> lon:units = "degrees_east" ;
> lon:long_name = "longitude" ;
> int time(time) ;
> time:units = "seconds since 2020-08-01 00:00:00" ;
> time:standard_name = "time" ;
> float total_aod(time, copy, lat, lon) ;
>
> // global attributes:
> :Conventions = "CF-1.0" ;
> :history = "Mon Jan 11 18:49:58 2021: ncks -v
> total_aod,copy -d copy,7
> /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc
> /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc" ;
> :NCO = "netCDF Operators version 4.7.7 (Homepage =
> http://nco.sf.net, Code = http://github.com/nco/nco)" ;
> }
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>
--
George McCabe - Software Engineer III
National Center for Atmospheric Research
Research Applications Laboratory
303-497-2768
---
My working day may not be your working day. Please do not feel obliged
to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF input file
From: efren.serra.ctr at nrlmry.navy.mil
Time: Wed Jan 13 13:29:44 2021
Hi George - This is nice! Thanks for the response below. Also, I take
it that I could specify in the level ="()" a list of levels such as
level="[(20200801_00, 0, *,*), (20200801_06, 0, *,*), ...];" for each
field element.
-----Original Message-----
From: George McCabe via RT <met_help at ucar.edu>
Sent: Wednesday, January 13, 2021 12:16 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF input
file
Hi Efren,
I see you are having trouble reading in NetCDF data with Plot-Data-
Plane.
It looks like there are 4 dimensions for the total_aod field. The
level you specified only contains 2 dimensions.
level="(*,*)"
You will need to specify dimensions that match the data. The first
dimension is time (size 21) and the 2nd dimension is copy (size 1).
Also, you are only able to read a single time into Plot-Data-Plane, so
you should remove the start_time and end_time values.
If you want to read the first time in the file, you could set
level="(0, 0, *,*)"
If you want MET to find the appropriate time for you, you can specify
the time value in the time dimension. I believe the syntax would be:
level="(20200801_00, 0, *,*)"
Let me know if you continue to have issues and I can take a closer
look.
Thanks,
George
On Wed, Jan 13, 2021 at 12:57 PM efren.serra.ctr at nrlmry.navy.mil via
RT < met_help at ucar.edu> wrote:
>
> Wed Jan 13 12:56:53 2021: Request 98227 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> Queue: met_help
> Subject: plot_data_plane and NetCDF input file
> Owner: Nobody
> Requestors: efren.serra.ctr at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227
> >
>
>
> I am having trouble using plot_data_plane with a NetCDF file with
the
> header information below: plot_data_plane
2020080100_ENAAPS_total_aod_
> 001.nc total_aod_001.ps 'name="total_aod"; level="(*,*)";
> start_time="20200801_00"; end_time="20200801_06";' -log pdp.log -v
4.
> The error is:
>
> DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcCF".
> DEBUG 4: NcCfFile::open() -> parsing units for the time variable
> "seconds since 2020-08-01 00:00:00"
> DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention time
> unit string "seconds since 2020-08-01 00:00:00" as a reference time
of
> 20200801_000000 and 1 second(s) per time step.
> DEBUG 4: NcCfFile::open() -> could not extract init time from the
> "forecast_reference_time" variable.
> DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of type "FileType_NcCF".
> WARNING:
> WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) ->
> returns the first available time for "total_aod" variable).
> WARNING:
> ERROR :
> ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &)
> const
> -> needed 4 arguments for variable total_aod, got 2
> ERROR :
>
> netcdf \2020080100_ENAAPS_total_aod_001 {
> dimensions:
> copy = 1 ;
> lat = 180 ;
> lon = 360 ;
> time = 21 ;
> variables:
> int copy(copy) ;
> copy:long_name = "[Ensemble Mean, Ensemble
> StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
> float lat(lat) ;
> lat:units = "degrees_north" ;
> lat:long_name = "latitude" ;
> float lon(lon) ;
> lon:units = "degrees_east" ;
> lon:long_name = "longitude" ;
> int time(time) ;
> time:units = "seconds since 2020-08-01 00:00:00" ;
> time:standard_name = "time" ;
> float total_aod(time, copy, lat, lon) ;
>
> // global attributes:
> :Conventions = "CF-1.0" ;
> :history = "Mon Jan 11 18:49:58 2021: ncks -v
> total_aod,copy -d copy,7
> /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc
> /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc" ;
> :NCO = "netCDF Operators version 4.7.7 (Homepage =
> http://nco.sf.net, Code = http://github.com/nco/nco)" ; }
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>
--
George McCabe - Software Engineer III
National Center for Atmospheric Research Research Applications
Laboratory
303-497-2768
---
My working day may not be your working day. Please do not feel obliged
to reply to this email outside of your normal working hours.
------------------------------------------------
Subject: plot_data_plane and NetCDF input file
From: George McCabe
Time: Wed Jan 13 13:36:24 2021
Hi Efren,
For some MET tools, yes, you can specify a list of levels the way you
described. However, Plot-Data-Plane only takes a single name/level to
plot,
so it will not work for that tool.
- George
On Wed, Jan 13, 2021 at 1:30 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227 >
>
> Hi George - This is nice! Thanks for the response below. Also, I
take it
> that I could specify in the level ="()" a list of levels such as
> level="[(20200801_00, 0, *,*), (20200801_06, 0, *,*), ...];" for
each field
> element.
>
> -----Original Message-----
> From: George McCabe via RT <met_help at ucar.edu>
> Sent: Wednesday, January 13, 2021 12:16 PM
> To: Serra, Mr. Efren, Contractor, Code 7531 <
> efren.serra.ctr at nrlmry.navy.mil>
> Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF
input
> file
>
> Hi Efren,
>
> I see you are having trouble reading in NetCDF data with Plot-Data-
Plane.
> It looks like there are 4 dimensions for the total_aod field. The
level
> you specified only contains 2 dimensions.
>
> level="(*,*)"
>
> You will need to specify dimensions that match the data. The first
> dimension is time (size 21) and the 2nd dimension is copy (size 1).
Also,
> you are only able to read a single time into Plot-Data-Plane, so you
should
> remove the start_time and end_time values.
>
> If you want to read the first time in the file, you could set
>
> level="(0, 0, *,*)"
>
> If you want MET to find the appropriate time for you, you can
specify the
> time value in the time dimension. I believe the syntax would be:
>
> level="(20200801_00, 0, *,*)"
>
> Let me know if you continue to have issues and I can take a closer
look.
>
> Thanks,
> George
>
> On Wed, Jan 13, 2021 at 12:57 PM efren.serra.ctr at nrlmry.navy.mil via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > Wed Jan 13 12:56:53 2021: Request 98227 was acted upon.
> > Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> > Queue: met_help
> > Subject: plot_data_plane and NetCDF input file
> > Owner: Nobody
> > Requestors: efren.serra.ctr at nrlmry.navy.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227
> > >
> >
> >
> > I am having trouble using plot_data_plane with a NetCDF file with
the
> > header information below: plot_data_plane
2020080100_ENAAPS_total_aod_
> > 001.nc total_aod_001.ps 'name="total_aod"; level="(*,*)";
> > start_time="20200801_00"; end_time="20200801_06";' -log pdp.log -v
4.
> > The error is:
> >
> > DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_NcCF".
> > DEBUG 4: NcCfFile::open() -> parsing units for the time variable
> > "seconds since 2020-08-01 00:00:00"
> > DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention
time
> > unit string "seconds since 2020-08-01 00:00:00" as a reference
time of
> > 20200801_000000 and 1 second(s) per time step.
> > DEBUG 4: NcCfFile::open() -> could not extract init time from the
> > "forecast_reference_time" variable.
> > DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> > of type "FileType_NcCF".
> > WARNING:
> > WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) ->
> > returns the first available time for "total_aod" variable).
> > WARNING:
> > ERROR :
> > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&)
> > const
> > -> needed 4 arguments for variable total_aod, got 2
> > ERROR :
> >
> > netcdf \2020080100_ENAAPS_total_aod_001 {
> > dimensions:
> > copy = 1 ;
> > lat = 180 ;
> > lon = 360 ;
> > time = 21 ;
> > variables:
> > int copy(copy) ;
> > copy:long_name = "[Ensemble Mean, Ensemble
> > StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
> > float lat(lat) ;
> > lat:units = "degrees_north" ;
> > lat:long_name = "latitude" ;
> > float lon(lon) ;
> > lon:units = "degrees_east" ;
> > lon:long_name = "longitude" ;
> > int time(time) ;
> > time:units = "seconds since 2020-08-01 00:00:00" ;
> > time:standard_name = "time" ;
> > float total_aod(time, copy, lat, lon) ;
> >
> > // global attributes:
> > :Conventions = "CF-1.0" ;
> > :history = "Mon Jan 11 18:49:58 2021: ncks -v
> > total_aod,copy -d copy,7
> > /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc
> > /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc" ;
> > :NCO = "netCDF Operators version 4.7.7 (Homepage =
> > http://nco.sf.net, Code = http://github.com/nco/nco)" ; }
> >
> > Efren A. Serra (Contractor)
> > Physicist
> >
> > DeVine Consulting, Inc.
> > Naval Research Laboratory
> > Marine Meteorology Division
> > 7 Grace Hopper Ave., STOP 2
> > Monterey, CA 93943
> > Code 7542
> > Mobile: 408-425-5027
> >
> >
> >
>
> --
> George McCabe - Software Engineer III
> National Center for Atmospheric Research Research Applications
Laboratory
> 303-497-2768
> ---
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
>
>
--
George McCabe - Software Engineer III
National Center for Atmospheric Research
Research Applications Laboratory
303-497-2768
---
My working day may not be your working day. Please do not feel obliged
to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF input file
From: efren.serra.ctr at nrlmry.navy.mil
Time: Wed Jan 13 13:39:35 2021
How about for ensemble_stat or grid_stat? Those are the tools I'm
thinking. Feel free to close this ticket after your reply.
-----Original Message-----
From: George McCabe via RT <met_help at ucar.edu>
Sent: Wednesday, January 13, 2021 12:36 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF input
file
Hi Efren,
For some MET tools, yes, you can specify a list of levels the way you
described. However, Plot-Data-Plane only takes a single name/level to
plot, so it will not work for that tool.
- George
On Wed, Jan 13, 2021 at 1:30 PM efren.serra.ctr at nrlmry.navy.mil via RT
< met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227 >
>
> Hi George - This is nice! Thanks for the response below. Also, I
take
> it that I could specify in the level ="()" a list of levels such as
> level="[(20200801_00, 0, *,*), (20200801_06, 0, *,*), ...];" for
each
> field element.
>
> -----Original Message-----
> From: George McCabe via RT <met_help at ucar.edu>
> Sent: Wednesday, January 13, 2021 12:16 PM
> To: Serra, Mr. Efren, Contractor, Code 7531 <
> efren.serra.ctr at nrlmry.navy.mil>
> Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF
input
> file
>
> Hi Efren,
>
> I see you are having trouble reading in NetCDF data with Plot-Data-
Plane.
> It looks like there are 4 dimensions for the total_aod field. The
> level you specified only contains 2 dimensions.
>
> level="(*,*)"
>
> You will need to specify dimensions that match the data. The first
> dimension is time (size 21) and the 2nd dimension is copy (size 1).
> Also, you are only able to read a single time into Plot-Data-Plane,
so
> you should remove the start_time and end_time values.
>
> If you want to read the first time in the file, you could set
>
> level="(0, 0, *,*)"
>
> If you want MET to find the appropriate time for you, you can
specify
> the time value in the time dimension. I believe the syntax would be:
>
> level="(20200801_00, 0, *,*)"
>
> Let me know if you continue to have issues and I can take a closer
look.
>
> Thanks,
> George
>
> On Wed, Jan 13, 2021 at 12:57 PM efren.serra.ctr at nrlmry.navy.mil via
> RT < met_help at ucar.edu> wrote:
>
> >
> > Wed Jan 13 12:56:53 2021: Request 98227 was acted upon.
> > Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> > Queue: met_help
> > Subject: plot_data_plane and NetCDF input file
> > Owner: Nobody
> > Requestors: efren.serra.ctr at nrlmry.navy.mil
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227
> > >
> >
> >
> > I am having trouble using plot_data_plane with a NetCDF file with
> > the header information below: plot_data_plane
> > 2020080100_ENAAPS_total_aod_ 001.nc total_aod_001.ps
> > 'name="total_aod"; level="(*,*)"; start_time="20200801_00";
end_time="20200801_06";' -log pdp.log -v 4.
> > The error is:
> >
> > DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_NcCF".
> > DEBUG 4: NcCfFile::open() -> parsing units for the time variable
> > "seconds since 2020-08-01 00:00:00"
> > DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention
time
> > unit string "seconds since 2020-08-01 00:00:00" as a reference
time
> > of
> > 20200801_000000 and 1 second(s) per time step.
> > DEBUG 4: NcCfFile::open() -> could not extract init time from the
> > "forecast_reference_time" variable.
> > DEBUG 4: NcCfFile::open() -> could not extract init time from file
name.
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
> > object of type "FileType_NcCF".
> > WARNING:
> > WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) ->
> > returns the first available time for "total_aod" variable).
> > WARNING:
> > ERROR :
> > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&)
> > const
> > -> needed 4 arguments for variable total_aod, got 2
> > ERROR :
> >
> > netcdf \2020080100_ENAAPS_total_aod_001 {
> > dimensions:
> > copy = 1 ;
> > lat = 180 ;
> > lon = 360 ;
> > time = 21 ;
> > variables:
> > int copy(copy) ;
> > copy:long_name = "[Ensemble Mean, Ensemble
> > StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
> > float lat(lat) ;
> > lat:units = "degrees_north" ;
> > lat:long_name = "latitude" ;
> > float lon(lon) ;
> > lon:units = "degrees_east" ;
> > lon:long_name = "longitude" ;
> > int time(time) ;
> > time:units = "seconds since 2020-08-01 00:00:00" ;
> > time:standard_name = "time" ;
> > float total_aod(time, copy, lat, lon) ;
> >
> > // global attributes:
> > :Conventions = "CF-1.0" ;
> > :history = "Mon Jan 11 18:49:58 2021: ncks -v
> > total_aod,copy -d copy,7
> > /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc
> > /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc" ;
> > :NCO = "netCDF Operators version 4.7.7 (Homepage =
> > http://nco.sf.net, Code = http://github.com/nco/nco)" ; }
> >
> > Efren A. Serra (Contractor)
> > Physicist
> >
> > DeVine Consulting, Inc.
> > Naval Research Laboratory
> > Marine Meteorology Division
> > 7 Grace Hopper Ave., STOP 2
> > Monterey, CA 93943
> > Code 7542
> > Mobile: 408-425-5027
> >
> >
> >
>
> --
> George McCabe - Software Engineer III
> National Center for Atmospheric Research Research Applications
> Laboratory
> 303-497-2768
> ---
> My working day may not be your working day. Please do not feel
obliged
> to reply to this email outside of your normal working hours.
>
>
>
>
--
George McCabe - Software Engineer III
National Center for Atmospheric Research Research Applications
Laboratory
303-497-2768
---
My working day may not be your working day. Please do not feel obliged
to reply to this email outside of your normal working hours.
------------------------------------------------
Subject: plot_data_plane and NetCDF input file
From: George McCabe
Time: Wed Jan 13 13:51:10 2021
Yes, those tools do support a list of field levels.
On Wed, Jan 13, 2021 at 1:39 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227 >
>
> How about for ensemble_stat or grid_stat? Those are the tools I'm
> thinking. Feel free to close this ticket after your reply.
>
> -----Original Message-----
> From: George McCabe via RT <met_help at ucar.edu>
> Sent: Wednesday, January 13, 2021 12:36 PM
> To: Serra, Mr. Efren, Contractor, Code 7531 <
> efren.serra.ctr at nrlmry.navy.mil>
> Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF
input
> file
>
> Hi Efren,
>
> For some MET tools, yes, you can specify a list of levels the way
you
> described. However, Plot-Data-Plane only takes a single name/level
to plot,
> so it will not work for that tool.
>
> - George
>
> On Wed, Jan 13, 2021 at 1:30 PM efren.serra.ctr at nrlmry.navy.mil via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227 >
> >
> > Hi George - This is nice! Thanks for the response below. Also, I
take
> > it that I could specify in the level ="()" a list of levels such
as
> > level="[(20200801_00, 0, *,*), (20200801_06, 0, *,*), ...];" for
each
> > field element.
> >
> > -----Original Message-----
> > From: George McCabe via RT <met_help at ucar.edu>
> > Sent: Wednesday, January 13, 2021 12:16 PM
> > To: Serra, Mr. Efren, Contractor, Code 7531 <
> > efren.serra.ctr at nrlmry.navy.mil>
> > Subject: Re: [rt.rap.ucar.edu #98227] plot_data_plane and NetCDF
input
> > file
> >
> > Hi Efren,
> >
> > I see you are having trouble reading in NetCDF data with Plot-
Data-Plane.
> > It looks like there are 4 dimensions for the total_aod field. The
> > level you specified only contains 2 dimensions.
> >
> > level="(*,*)"
> >
> > You will need to specify dimensions that match the data. The first
> > dimension is time (size 21) and the 2nd dimension is copy (size
1).
> > Also, you are only able to read a single time into Plot-Data-
Plane, so
> > you should remove the start_time and end_time values.
> >
> > If you want to read the first time in the file, you could set
> >
> > level="(0, 0, *,*)"
> >
> > If you want MET to find the appropriate time for you, you can
specify
> > the time value in the time dimension. I believe the syntax would
be:
> >
> > level="(20200801_00, 0, *,*)"
> >
> > Let me know if you continue to have issues and I can take a closer
look.
> >
> > Thanks,
> > George
> >
> > On Wed, Jan 13, 2021 at 12:57 PM efren.serra.ctr at nrlmry.navy.mil
via
> > RT < met_help at ucar.edu> wrote:
> >
> > >
> > > Wed Jan 13 12:56:53 2021: Request 98227 was acted upon.
> > > Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> > > Queue: met_help
> > > Subject: plot_data_plane and NetCDF input file
> > > Owner: Nobody
> > > Requestors: efren.serra.ctr at nrlmry.navy.mil
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98227
> > > >
> > >
> > >
> > > I am having trouble using plot_data_plane with a NetCDF file
with
> > > the header information below: plot_data_plane
> > > 2020080100_ENAAPS_total_aod_ 001.nc total_aod_001.ps
> > > 'name="total_aod"; level="(*,*)"; start_time="20200801_00";
> end_time="20200801_06";' -log pdp.log -v 4.
> > > The error is:
> > >
> > > DEBUG 1: Opening data file: 2020080100_ENAAPS_total_aod_001.nc
> > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > > Met2dDataFile object of type "FileType_NcCF".
> > > DEBUG 4: NcCfFile::open() -> parsing units for the time variable
> > > "seconds since 2020-08-01 00:00:00"
> > > DEBUG 4: parse_cf_time_string() -> parsed NetCDF CF convention
time
> > > unit string "seconds since 2020-08-01 00:00:00" as a reference
time
> > > of
> > > 20200801_000000 and 1 second(s) per time step.
> > > DEBUG 4: NcCfFile::open() -> could not extract init time from
the
> > > "forecast_reference_time" variable.
> > > DEBUG 4: NcCfFile::open() -> could not extract init time from
file
> name.
> > > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
> > > object of type "FileType_NcCF".
> > > WARNING:
> > > WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) ->
> > > returns the first available time for "total_aod" variable).
> > > WARNING:
> > > ERROR :
> > > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&)
> > > const
> > > -> needed 4 arguments for variable total_aod, got 2
> > > ERROR :
> > >
> > > netcdf \2020080100_ENAAPS_total_aod_001 {
> > > dimensions:
> > > copy = 1 ;
> > > lat = 180 ;
> > > lon = 360 ;
> > > time = 21 ;
> > > variables:
> > > int copy(copy) ;
> > > copy:long_name = "[Ensemble Mean, Ensemble
> > > StDev,10th,25th,50th,75th,90th Percentile, Ensemble i....n]" ;
> > > float lat(lat) ;
> > > lat:units = "degrees_north" ;
> > > lat:long_name = "latitude" ;
> > > float lon(lon) ;
> > > lon:units = "degrees_east" ;
> > > lon:long_name = "longitude" ;
> > > int time(time) ;
> > > time:units = "seconds since 2020-08-01 00:00:00"
;
> > > time:standard_name = "time" ;
> > > float total_aod(time, copy, lat, lon) ;
> > >
> > > // global attributes:
> > > :Conventions = "CF-1.0" ;
> > > :history = "Mon Jan 11 18:49:58 2021: ncks -v
> > > total_aod,copy -d copy,7
> > > /p/home/serrae/verification/Forecast/2020080100_ENAAPS.nc
> > > /p/home/serrae/verification/2020080100_ENAAPS_total_aod_001.nc"
;
> > > :NCO = "netCDF Operators version 4.7.7 (Homepage
=
> > > http://nco.sf.net, Code = http://github.com/nco/nco)" ; }
> > >
> > > Efren A. Serra (Contractor)
> > > Physicist
> > >
> > > DeVine Consulting, Inc.
> > > Naval Research Laboratory
> > > Marine Meteorology Division
> > > 7 Grace Hopper Ave., STOP 2
> > > Monterey, CA 93943
> > > Code 7542
> > > Mobile: 408-425-5027
> > >
> > >
> > >
> >
> > --
> > George McCabe - Software Engineer III
> > National Center for Atmospheric Research Research Applications
> > Laboratory
> > 303-497-2768
> > ---
> > My working day may not be your working day. Please do not feel
obliged
> > to reply to this email outside of your normal working hours.
> >
> >
> >
> >
>
> --
> George McCabe - Software Engineer III
> National Center for Atmospheric Research Research Applications
Laboratory
> 303-497-2768
> ---
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
>
>
--
George McCabe - Software Engineer III
National Center for Atmospheric Research
Research Applications Laboratory
303-497-2768
---
My working day may not be your working day. Please do not feel obliged
to
reply to this email outside of your normal working hours.
------------------------------------------------
More information about the Met_help
mailing list