[Met_help] [rt.rap.ucar.edu #81369] History for Forecast Regrid
John Halley Gotway via RT
met_help at ucar.edu
Wed Aug 9 10:14:45 MDT 2017
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I am a new user to the MET program, and attempting to regrid a 1051x813
HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km using the budget
regridding method. I am wondering if you have a simple, straightforward
example of how to do this in MET. So far, I have managed to figure out this
is needed in the config file:
regrid = {
to_grid = NONE;
method = BUDGET;
width = 1;
vld_thresh = 0.5;
}
Any other help would be much appreciated.
Thanks in advance!
Dan
--
Daniel W. Plumb, Jr.
Scientific Programmer/Meteorologist
Ace Info Solutions, Inc. (contractor)
National Weather Service (NWS)- Meteorological Development Lab (MDL)
Location:
SSMC2 – Room 10124
1325 East West Highway
Silver Spring, MD 20910
Tel: (301) 427-9453
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Forecast Regrid
From: John Halley Gotway
Time: Thu Jul 27 12:50:34 2017
Hi Dan,
Tara gave me a heads up that you've started using MET at MDL. We're
happy
to answer any questions you have or issues that come up.
I see that you'd like to use the regridding feature within the MET
config
files. I'd actually suggest starting by running the regrid_data_plane
tool
first, and once you understand what it's doing, move those settings
into
the config file. Both the regrid_data_plane tool and the regrid
section of
the MET config files call the same underlying library code anyway.
I was hoping to grab some sample MDL data to demonstrate some
commands, but
all the links on this website aren't working:
http://www.nws.noaa.gov/mdl/hrqpf/download.php
So instead of using HRAP and NDFD data, I'll use sample GFS and NAM
data.
# Get sample GFS and NAM files
wget
http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.2017072700/gfs.t00z.pgrb2.0p25.f012
wget
http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
# Run regrid_data_plane to regrid 6-hour precip from GFS to NAM nest
domain
met-6.0/bin/regrid_data_plane \
gfs.t00z.pgrb2.0p25.f012 nam.t12z.conusnest.hiresf12.tm00.grib2
gfs_APCP_06_nam_nest_domain.nc \
-field 'name="APCP"; level="A06";' -interp_mthd BUDGET -interp_pnts
2
# Run plot_data_plane to plot the 6-hour precip both before/after the
interpolation
met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP"; level="A06";' \
-color_table met-6.0/share/met/colortables/wrf_precip2.ctable
met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06"; level="(*,*)";' \
-color_table met-6.0/share/met/colortables/wrf_precip2.ctable
The resulting plots (converted to png) are attached to this email.
In the regridding step, I defined the "to_grid" by providing the path
to a
file that's already on the grid I'd like to use. MET reads that file,
grabs the grid definition, and uses that in the regridding step.
There are
multiple ways of defining "to_grid" (as described in
met-6.0/data/config/README and the User's Guide), but I find providing
a
path to a sample file is pretty convenient.
Or, even more convenient, if I'm running grid_stat for example and
pass in
the GFS as the forecast and NAM as the observation, I can just say
"to_grid
= OBS;" and grid_stat will automatically regrid the forecast data
(GFS) to
the observation grid (NAM).
So, long story short, the following would work:
regrid = {
to_grid = "/path/to/sample/NDFD/data/file"; // Or set this to
FCST or
OBS if that's easier
method = BUDGET;
width = 2; // Budget should always have a width of 2 to use
the 4
closest points
vld_thresh = 0.5;
}
Hope that helps get you going.
Thanks,
John
On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> Transaction: Ticket created by daniel.plumb at noaa.gov
> Queue: met_help
> Subject: Forecast Regrid
> Owner: Nobody
> Requestors: daniel.plumb at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
>
> Hello,
>
> I am a new user to the MET program, and attempting to regrid a
1051x813
> HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km using the
budget
> regridding method. I am wondering if you have a simple,
straightforward
> example of how to do this in MET. So far, I have managed to figure
out this
> is needed in the config file:
>
> regrid = {
> to_grid = NONE;
> method = BUDGET;
> width = 1;
> vld_thresh = 0.5;
> }
>
> Any other help would be much appreciated.
>
> Thanks in advance!
>
> Dan
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
>
------------------------------------------------
Subject: Forecast Regrid
From: Daniel Plumb - NOAA Affiliate
Time: Thu Jul 27 14:53:25 2017
Hi John,
Thank you very much for the reply and examples.
I am trying to replicate you example for regrid_data_plane (which is
exactly the example I needed), but there seems to be something wrong
with
MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3 udunits
ncview
met/6.0 grib_util/v1.0.3", it is giving me an error saying:
"ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with the
currently
loaded module(s) 'ics/12.1'
ics/15.0.3(16):ERROR:102: Tcl command execution failed: conflict ics"
I'm looking at the instructions under
http://www.dtcenter.org/met/users/support/online_tutorial/ncwcp_201704.php.
Thanks,
Dan
On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Dan,
>
> Tara gave me a heads up that you've started using MET at MDL. We're
happy
> to answer any questions you have or issues that come up.
>
> I see that you'd like to use the regridding feature within the MET
config
> files. I'd actually suggest starting by running the
regrid_data_plane tool
> first, and once you understand what it's doing, move those settings
into
> the config file. Both the regrid_data_plane tool and the regrid
section of
> the MET config files call the same underlying library code anyway.
>
> I was hoping to grab some sample MDL data to demonstrate some
commands, but
> all the links on this website aren't working:
> http://www.nws.noaa.gov/mdl/hrqpf/download.php
>
> So instead of using HRAP and NDFD data, I'll use sample GFS and NAM
data.
>
> # Get sample GFS and NAM files
> wget
> http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> 2017072700/gfs.t00z.pgrb2.0p25.f012
> wget
> http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
>
> # Run regrid_data_plane to regrid 6-hour precip from GFS to NAM nest
domain
> met-6.0/bin/regrid_data_plane \
> gfs.t00z.pgrb2.0p25.f012 nam.t12z.conusnest.hiresf12.tm00.grib2
> gfs_APCP_06_nam_nest_domain.nc \
> -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
-interp_pnts 2
>
> # Run plot_data_plane to plot the 6-hour precip both before/after
the
> interpolation
> met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP"; level="A06";' \
> -color_table met-6.0/share/met/colortables/wrf_precip2.ctable
>
> met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
> gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06"; level="(*,*)";' \
> -color_table met-6.0/share/met/colortables/wrf_precip2.ctable
>
> The resulting plots (converted to png) are attached to this email.
>
> In the regridding step, I defined the "to_grid" by providing the
path to a
> file that's already on the grid I'd like to use. MET reads that
file,
> grabs the grid definition, and uses that in the regridding step.
There are
> multiple ways of defining "to_grid" (as described in
> met-6.0/data/config/README and the User's Guide), but I find
providing a
> path to a sample file is pretty convenient.
>
> Or, even more convenient, if I'm running grid_stat for example and
pass in
> the GFS as the forecast and NAM as the observation, I can just say
"to_grid
> = OBS;" and grid_stat will automatically regrid the forecast data
(GFS) to
> the observation grid (NAM).
>
> So, long story short, the following would work:
>
> regrid = {
> to_grid = "/path/to/sample/NDFD/data/file"; // Or set this to
FCST
> or
> OBS if that's easier
> method = BUDGET;
> width = 2; // Budget should always have a width of 2 to
use the 4
> closest points
> vld_thresh = 0.5;
> }
>
> Hope that helps get you going.
>
> Thanks,
> John
>
>
> On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> > Transaction: Ticket created by daniel.plumb at noaa.gov
> > Queue: met_help
> > Subject: Forecast Regrid
> > Owner: Nobody
> > Requestors: daniel.plumb at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> >
> >
> > Hello,
> >
> > I am a new user to the MET program, and attempting to regrid a
1051x813
> > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km using the
budget
> > regridding method. I am wondering if you have a simple,
straightforward
> > example of how to do this in MET. So far, I have managed to figure
out
> this
> > is needed in the config file:
> >
> > regrid = {
> > to_grid = NONE;
> > method = BUDGET;
> > width = 1;
> > vld_thresh = 0.5;
> > }
> >
> > Any other help would be much appreciated.
> >
> > Thanks in advance!
> >
> > Dan
> >
> > --
> > Daniel W. Plumb, Jr.
> > Scientific Programmer/Meteorologist
> > Ace Info Solutions, Inc. (contractor)
> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> >
> > Location:
> > SSMC2 – Room 10124
> > 1325 East West Highway
> > Silver Spring, MD 20910
> > Tel: (301) 427-9453
> >
> >
>
>
--
Daniel W. Plumb, Jr.
Scientific Programmer/Meteorologist
Ace Info Solutions, Inc. (contractor)
National Weather Service (NWS)- Meteorological Development Lab (MDL)
Location:
SSMC2 – Room 10124
1325 East West Highway
Silver Spring, MD 20910
Tel: (301) 427-9453
------------------------------------------------
Subject: Forecast Regrid
From: John Halley Gotway
Time: Thu Jul 27 15:02:44 2017
Dan,
I don't have access to WCOSS to replicate this problem, but by
colleague
Julie does.
Julie, could you please take a look at this?
Thanks,
John
On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
> Hi John,
>
> Thank you very much for the reply and examples.
>
> I am trying to replicate you example for regrid_data_plane (which is
> exactly the example I needed), but there seems to be something wrong
with
> MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3 udunits
ncview
> met/6.0 grib_util/v1.0.3", it is giving me an error saying:
>
> "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with the
currently
> loaded module(s) 'ics/12.1'
> ics/15.0.3(16):ERROR:102: Tcl command execution failed: conflict
ics"
>
> I'm looking at the instructions under
>
http://www.dtcenter.org/met/users/support/online_tutorial/ncwcp_201704.php
> .
>
> Thanks,
> Dan
>
> On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hi Dan,
> >
> > Tara gave me a heads up that you've started using MET at MDL.
We're
> happy
> > to answer any questions you have or issues that come up.
> >
> > I see that you'd like to use the regridding feature within the MET
config
> > files. I'd actually suggest starting by running the
regrid_data_plane
> tool
> > first, and once you understand what it's doing, move those
settings into
> > the config file. Both the regrid_data_plane tool and the regrid
section
> of
> > the MET config files call the same underlying library code anyway.
> >
> > I was hoping to grab some sample MDL data to demonstrate some
commands,
> but
> > all the links on this website aren't working:
> > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> >
> > So instead of using HRAP and NDFD data, I'll use sample GFS and
NAM data.
> >
> > # Get sample GFS and NAM files
> > wget
> > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > wget
> > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> >
> > # Run regrid_data_plane to regrid 6-hour precip from GFS to NAM
nest
> domain
> > met-6.0/bin/regrid_data_plane \
> > gfs.t00z.pgrb2.0p25.f012 nam.t12z.conusnest.hiresf12.tm00.grib2
> > gfs_APCP_06_nam_nest_domain.nc \
> > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
-interp_pnts 2
> >
> > # Run plot_data_plane to plot the 6-hour precip both before/after
the
> > interpolation
> > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP"; level="A06";' \
> > -color_table met-6.0/share/met/colortables/wrf_precip2.ctable
> >
> > met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
> > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06"; level="(*,*)";' \
> > -color_table met-6.0/share/met/colortables/wrf_precip2.ctable
> >
> > The resulting plots (converted to png) are attached to this email.
> >
> > In the regridding step, I defined the "to_grid" by providing the
path to
> a
> > file that's already on the grid I'd like to use. MET reads that
file,
> > grabs the grid definition, and uses that in the regridding step.
There
> are
> > multiple ways of defining "to_grid" (as described in
> > met-6.0/data/config/README and the User's Guide), but I find
providing a
> > path to a sample file is pretty convenient.
> >
> > Or, even more convenient, if I'm running grid_stat for example and
pass
> in
> > the GFS as the forecast and NAM as the observation, I can just say
> "to_grid
> > = OBS;" and grid_stat will automatically regrid the forecast data
(GFS)
> to
> > the observation grid (NAM).
> >
> > So, long story short, the following would work:
> >
> > regrid = {
> > to_grid = "/path/to/sample/NDFD/data/file"; // Or set this
to FCST
> > or
> > OBS if that's easier
> > method = BUDGET;
> > width = 2; // Budget should always have a width of 2 to
use
> the 4
> > closest points
> > vld_thresh = 0.5;
> > }
> >
> > Hope that helps get you going.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> > > Transaction: Ticket created by daniel.plumb at noaa.gov
> > > Queue: met_help
> > > Subject: Forecast Regrid
> > > Owner: Nobody
> > > Requestors: daniel.plumb at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
> >
> > >
> > >
> > > Hello,
> > >
> > > I am a new user to the MET program, and attempting to regrid a
1051x813
> > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km using
the
> budget
> > > regridding method. I am wondering if you have a simple,
straightforward
> > > example of how to do this in MET. So far, I have managed to
figure out
> > this
> > > is needed in the config file:
> > >
> > > regrid = {
> > > to_grid = NONE;
> > > method = BUDGET;
> > > width = 1;
> > > vld_thresh = 0.5;
> > > }
> > >
> > > Any other help would be much appreciated.
> > >
> > > Thanks in advance!
> > >
> > > Dan
> > >
> > > --
> > > Daniel W. Plumb, Jr.
> > > Scientific Programmer/Meteorologist
> > > Ace Info Solutions, Inc. (contractor)
> > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > >
> > > Location:
> > > SSMC2 – Room 10124
> > > 1325 East West Highway
> > > Silver Spring, MD 20910
> > > Tel: (301) 427-9453
> > >
> > >
> >
> >
>
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
>
------------------------------------------------
Subject: Forecast Regrid
From: Julie Prestopnik
Time: Thu Jul 27 15:16:05 2017
I only have developer access to WCOSS, but I don't have access to
WCOSS
(gyre) at this time, as I believe it is still down for developers.
Please try running "module purge" and then run:
- module use /usrx/local/dev/modulefiles
- module load met
and see if that resolves the problem. I believe there is a conflict
with
the intel module you already have loaded. I know purging isn't the
best
method to solve this problem, however, I do believe that should work.
Please let me know if you have any questions.
Thanks,
Julie
On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
> Dan,
>
> I don't have access to WCOSS to replicate this problem, but by
colleague
> Julie does.
>
> Julie, could you please take a look at this?
>
> Thanks,
> John
>
>
> On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> >
> > Hi John,
> >
> > Thank you very much for the reply and examples.
> >
> > I am trying to replicate you example for regrid_data_plane (which
is
> > exactly the example I needed), but there seems to be something
wrong with
> > MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3
udunits
> ncview
> > met/6.0 grib_util/v1.0.3", it is giving me an error saying:
> >
> > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with the
> currently
> > loaded module(s) 'ics/12.1'
> > ics/15.0.3(16):ERROR:102: Tcl command execution failed: conflict
ics"
> >
> > I'm looking at the instructions under
> > http://www.dtcenter.org/met/users/support/online_tutorial/
> ncwcp_201704.php
> > .
> >
> > Thanks,
> > Dan
> >
> > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Hi Dan,
> > >
> > > Tara gave me a heads up that you've started using MET at MDL.
We're
> > happy
> > > to answer any questions you have or issues that come up.
> > >
> > > I see that you'd like to use the regridding feature within the
MET
> config
> > > files. I'd actually suggest starting by running the
regrid_data_plane
> > tool
> > > first, and once you understand what it's doing, move those
settings
> into
> > > the config file. Both the regrid_data_plane tool and the regrid
> section
> > of
> > > the MET config files call the same underlying library code
anyway.
> > >
> > > I was hoping to grab some sample MDL data to demonstrate some
commands,
> > but
> > > all the links on this website aren't working:
> > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> > >
> > > So instead of using HRAP and NDFD data, I'll use sample GFS and
NAM
> data.
> > >
> > > # Get sample GFS and NAM files
> > > wget
> > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > > wget
> > > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> > >
> > > # Run regrid_data_plane to regrid 6-hour precip from GFS to NAM
nest
> > domain
> > > met-6.0/bin/regrid_data_plane \
> > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.tm00.grib2
> > > gfs_APCP_06_nam_nest_domain.nc \
> > > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
-interp_pnts
> 2
> > >
> > > # Run plot_data_plane to plot the 6-hour precip both
before/after the
> > > interpolation
> > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP"; level="A06";'
\
> > > -color_table met-6.0/share/met/colortables/wrf_precip2.ctable
> > >
> > > met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
> > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06"; level="(*,*)";'
\
> > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
> > >
> > > The resulting plots (converted to png) are attached to this
email.
> > >
> > > In the regridding step, I defined the "to_grid" by providing the
path
> to
> > a
> > > file that's already on the grid I'd like to use. MET reads that
file,
> > > grabs the grid definition, and uses that in the regridding step.
There
> > are
> > > multiple ways of defining "to_grid" (as described in
> > > met-6.0/data/config/README and the User's Guide), but I find
providing
> a
> > > path to a sample file is pretty convenient.
> > >
> > > Or, even more convenient, if I'm running grid_stat for example
and pass
> > in
> > > the GFS as the forecast and NAM as the observation, I can just
say
> > "to_grid
> > > = OBS;" and grid_stat will automatically regrid the forecast
data (GFS)
> > to
> > > the observation grid (NAM).
> > >
> > > So, long story short, the following would work:
> > >
> > > regrid = {
> > > to_grid = "/path/to/sample/NDFD/data/file"; // Or set this
to
> FCST
> > > or
> > > OBS if that's easier
> > > method = BUDGET;
> > > width = 2; // Budget should always have a width of 2
to use
> > the 4
> > > closest points
> > > vld_thresh = 0.5;
> > > }
> > >
> > > Hope that helps get you going.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA Affiliate
via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> > > > Transaction: Ticket created by daniel.plumb at noaa.gov
> > > > Queue: met_help
> > > > Subject: Forecast Regrid
> > > > Owner: Nobody
> > > > Requestors: daniel.plumb at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81369
> > >
> > > >
> > > >
> > > > Hello,
> > > >
> > > > I am a new user to the MET program, and attempting to regrid a
> 1051x813
> > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km using
the
> > budget
> > > > regridding method. I am wondering if you have a simple,
> straightforward
> > > > example of how to do this in MET. So far, I have managed to
figure
> out
> > > this
> > > > is needed in the config file:
> > > >
> > > > regrid = {
> > > > to_grid = NONE;
> > > > method = BUDGET;
> > > > width = 1;
> > > > vld_thresh = 0.5;
> > > > }
> > > >
> > > > Any other help would be much appreciated.
> > > >
> > > > Thanks in advance!
> > > >
> > > > Dan
> > > >
> > > > --
> > > > Daniel W. Plumb, Jr.
> > > > Scientific Programmer/Meteorologist
> > > > Ace Info Solutions, Inc. (contractor)
> > > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > > >
> > > > Location:
> > > > SSMC2 – Room 10124
> > > > 1325 East West Highway
> > > > Silver Spring, MD 20910
> > > > Tel: (301) 427-9453
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Daniel W. Plumb, Jr.
> > Scientific Programmer/Meteorologist
> > Ace Info Solutions, Inc. (contractor)
> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> >
> > Location:
> > SSMC2 – Room 10124
> > 1325 East West Highway
> > Silver Spring, MD 20910
> > Tel: (301) 427-9453
> >
> >
>
>
------------------------------------------------
Subject: Forecast Regrid
From: Julie Prestopnik
Time: Thu Jul 27 15:18:31 2017
I meant to point out, before sending that last email, that we are
pointing
to a new location for METv6.0 on the WCOSS machines. You can see the
new
locations here:
http://www.dtcenter.org/met/users/downloads/index.php#met_builds
Julie
On Thu, Jul 27, 2017 at 3:15 PM, Julie Prestopnik <jpresto at ucar.edu>
wrote:
> I only have developer access to WCOSS, but I don't have access to
WCOSS
> (gyre) at this time, as I believe it is still down for developers.
>
> Please try running "module purge" and then run:
>
> - module use /usrx/local/dev/modulefiles
> - module load met
>
> and see if that resolves the problem. I believe there is a conflict
with
> the intel module you already have loaded. I know purging isn't the
best
> method to solve this problem, however, I do believe that should
work.
>
> Please let me know if you have any questions.
>
> Thanks,
>
> Julie
>
>
>
>
>
> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>>
>> Dan,
>>
>> I don't have access to WCOSS to replicate this problem, but by
colleague
>> Julie does.
>>
>> Julie, could you please take a look at this?
>>
>> Thanks,
>> John
>>
>>
>> On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA Affiliate via
RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>> >
>> > Hi John,
>> >
>> > Thank you very much for the reply and examples.
>> >
>> > I am trying to replicate you example for regrid_data_plane (which
is
>> > exactly the example I needed), but there seems to be something
wrong
>> with
>> > MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3
udunits
>> ncview
>> > met/6.0 grib_util/v1.0.3", it is giving me an error saying:
>> >
>> > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with the
>> currently
>> > loaded module(s) 'ics/12.1'
>> > ics/15.0.3(16):ERROR:102: Tcl command execution failed: conflict
ics"
>> >
>> > I'm looking at the instructions under
>> > http://www.dtcenter.org/met/users/support/online_tutorial/nc
>> wcp_201704.php
>> > .
>> >
>> > Thanks,
>> > Dan
>> >
>> > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
>> > met_help at ucar.edu> wrote:
>> >
>> > > Hi Dan,
>> > >
>> > > Tara gave me a heads up that you've started using MET at MDL.
We're
>> > happy
>> > > to answer any questions you have or issues that come up.
>> > >
>> > > I see that you'd like to use the regridding feature within the
MET
>> config
>> > > files. I'd actually suggest starting by running the
regrid_data_plane
>> > tool
>> > > first, and once you understand what it's doing, move those
settings
>> into
>> > > the config file. Both the regrid_data_plane tool and the
regrid
>> section
>> > of
>> > > the MET config files call the same underlying library code
anyway.
>> > >
>> > > I was hoping to grab some sample MDL data to demonstrate some
>> commands,
>> > but
>> > > all the links on this website aren't working:
>> > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
>> > >
>> > > So instead of using HRAP and NDFD data, I'll use sample GFS and
NAM
>> data.
>> > >
>> > > # Get sample GFS and NAM files
>> > > wget
>> > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
>> > > 2017072700/gfs.t00z.pgrb2.0p25.f012
>> > > wget
>> > > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
>> > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
>> > >
>> > > # Run regrid_data_plane to regrid 6-hour precip from GFS to NAM
nest
>> > domain
>> > > met-6.0/bin/regrid_data_plane \
>> > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.tm00.grib2
>> > > gfs_APCP_06_nam_nest_domain.nc \
>> > > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
>> -interp_pnts 2
>> > >
>> > > # Run plot_data_plane to plot the 6-hour precip both
before/after the
>> > > interpolation
>> > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
>> > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP"; level="A06";'
\
>> > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
>> > >
>> > > met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
>> > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06"; level="(*,*)";'
\
>> > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
>> > >
>> > > The resulting plots (converted to png) are attached to this
email.
>> > >
>> > > In the regridding step, I defined the "to_grid" by providing
the path
>> to
>> > a
>> > > file that's already on the grid I'd like to use. MET reads
that file,
>> > > grabs the grid definition, and uses that in the regridding
step.
>> There
>> > are
>> > > multiple ways of defining "to_grid" (as described in
>> > > met-6.0/data/config/README and the User's Guide), but I find
>> providing a
>> > > path to a sample file is pretty convenient.
>> > >
>> > > Or, even more convenient, if I'm running grid_stat for example
and
>> pass
>> > in
>> > > the GFS as the forecast and NAM as the observation, I can just
say
>> > "to_grid
>> > > = OBS;" and grid_stat will automatically regrid the forecast
data
>> (GFS)
>> > to
>> > > the observation grid (NAM).
>> > >
>> > > So, long story short, the following would work:
>> > >
>> > > regrid = {
>> > > to_grid = "/path/to/sample/NDFD/data/file"; // Or set
this to
>> FCST
>> > > or
>> > > OBS if that's easier
>> > > method = BUDGET;
>> > > width = 2; // Budget should always have a width of 2
to use
>> > the 4
>> > > closest points
>> > > vld_thresh = 0.5;
>> > > }
>> > >
>> > > Hope that helps get you going.
>> > >
>> > > Thanks,
>> > > John
>> > >
>> > >
>> > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA Affiliate
via RT
>> <
>> > > met_help at ucar.edu> wrote:
>> > >
>> > > >
>> > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
>> > > > Transaction: Ticket created by daniel.plumb at noaa.gov
>> > > > Queue: met_help
>> > > > Subject: Forecast Regrid
>> > > > Owner: Nobody
>> > > > Requestors: daniel.plumb at noaa.gov
>> > > > Status: new
>> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
>> ket/Display.html?id=81369
>> > >
>> > > >
>> > > >
>> > > > Hello,
>> > > >
>> > > > I am a new user to the MET program, and attempting to regrid
a
>> 1051x813
>> > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km
using the
>> > budget
>> > > > regridding method. I am wondering if you have a simple,
>> straightforward
>> > > > example of how to do this in MET. So far, I have managed to
figure
>> out
>> > > this
>> > > > is needed in the config file:
>> > > >
>> > > > regrid = {
>> > > > to_grid = NONE;
>> > > > method = BUDGET;
>> > > > width = 1;
>> > > > vld_thresh = 0.5;
>> > > > }
>> > > >
>> > > > Any other help would be much appreciated.
>> > > >
>> > > > Thanks in advance!
>> > > >
>> > > > Dan
>> > > >
>> > > > --
>> > > > Daniel W. Plumb, Jr.
>> > > > Scientific Programmer/Meteorologist
>> > > > Ace Info Solutions, Inc. (contractor)
>> > > > National Weather Service (NWS)- Meteorological Development
Lab (MDL)
>> > > >
>> > > > Location:
>> > > > SSMC2 – Room 10124
>> > > > 1325 East West Highway
>> > > > Silver Spring, MD 20910
>> > > > Tel: (301) 427-9453
>> > > >
>> > > >
>> > >
>> > >
>> >
>> >
>> > --
>> > Daniel W. Plumb, Jr.
>> > Scientific Programmer/Meteorologist
>> > Ace Info Solutions, Inc. (contractor)
>> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
>> >
>> > Location:
>> > SSMC2 – Room 10124
>> > 1325 East West Highway
>> > Silver Spring, MD 20910
>> > Tel: (301) 427-9453
>> >
>> >
>>
>>
>
------------------------------------------------
Subject: Forecast Regrid
From: Daniel Plumb - NOAA Affiliate
Time: Fri Jul 28 14:05:35 2017
Please let me know if you were able to get all of my attachments. I'm
not
sure if they all went through.
Thanks,
Dan
On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA Affiliate <
daniel.plumb at noaa.gov> wrote:
> Hi Julie and John,
>
> Thank you for the help so far. I am able to run regrid_data_plane
running
> now. However, it is not giving me the desired output I would like
right
> now. I have my command below and have attached the input dataset
> (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the dataset
with
> the desired grid specs I want (ndfdcogb2p517072700qpfd13_006
> <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2 format*):
>
> regrid_data_plane qpf_06hr_2017072512f006Z.nc
> ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field 'name="APCP";
> level="A06";' -method BUDGET -width 2
>
> DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> ERROR :
> ERROR : NcCfFile::get_grid_mapping_polar_stereographic() -> Polar
> stereographic grid not handled in MET.
> ERROR :
>
> Is MET not able to handle polar stereographic grids? Or maybe the
NetCDF
> library is too old on WCOSS (dated: Dec 5 2012)?
>
> Thanks,
> Dan
>
>
>
>
> On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> met_help at ucar.edu> wrote:
>
>> I only have developer access to WCOSS, but I don't have access to
WCOSS
>> (gyre) at this time, as I believe it is still down for developers.
>>
>> Please try running "module purge" and then run:
>>
>> - module use /usrx/local/dev/modulefiles
>> - module load met
>>
>> and see if that resolves the problem. I believe there is a
conflict with
>> the intel module you already have loaded. I know purging isn't the
best
>> method to solve this problem, however, I do believe that should
work.
>>
>> Please let me know if you have any questions.
>>
>> Thanks,
>>
>> Julie
>>
>>
>>
>>
>>
>> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>> >
>> > Dan,
>> >
>> > I don't have access to WCOSS to replicate this problem, but by
colleague
>> > Julie does.
>> >
>> > Julie, could you please take a look at this?
>> >
>> > Thanks,
>> > John
>> >
>> >
>> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA Affiliate
via RT <
>> > met_help at ucar.edu> wrote:
>> >
>> > >
>> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>> > >
>> > > Hi John,
>> > >
>> > > Thank you very much for the reply and examples.
>> > >
>> > > I am trying to replicate you example for regrid_data_plane
(which is
>> > > exactly the example I needed), but there seems to be something
wrong
>> with
>> > > MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3
udunits
>> > ncview
>> > > met/6.0 grib_util/v1.0.3", it is giving me an error saying:
>> > >
>> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with
the
>> > currently
>> > > loaded module(s) 'ics/12.1'
>> > > ics/15.0.3(16):ERROR:102: Tcl command execution failed:
conflict ics"
>> > >
>> > > I'm looking at the instructions under
>> > > http://www.dtcenter.org/met/users/support/online_tutorial/
>> > ncwcp_201704.php
>> > > .
>> > >
>> > > Thanks,
>> > > Dan
>> > >
>> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
>> > > met_help at ucar.edu> wrote:
>> > >
>> > > > Hi Dan,
>> > > >
>> > > > Tara gave me a heads up that you've started using MET at MDL.
We're
>> > > happy
>> > > > to answer any questions you have or issues that come up.
>> > > >
>> > > > I see that you'd like to use the regridding feature within
the MET
>> > config
>> > > > files. I'd actually suggest starting by running the
>> regrid_data_plane
>> > > tool
>> > > > first, and once you understand what it's doing, move those
settings
>> > into
>> > > > the config file. Both the regrid_data_plane tool and the
regrid
>> > section
>> > > of
>> > > > the MET config files call the same underlying library code
anyway.
>> > > >
>> > > > I was hoping to grab some sample MDL data to demonstrate some
>> commands,
>> > > but
>> > > > all the links on this website aren't working:
>> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
>> > > >
>> > > > So instead of using HRAP and NDFD data, I'll use sample GFS
and NAM
>> > data.
>> > > >
>> > > > # Get sample GFS and NAM files
>> > > > wget
>> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
>> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
>> > > > wget
>> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
>> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
>> > > >
>> > > > # Run regrid_data_plane to regrid 6-hour precip from GFS to
NAM nest
>> > > domain
>> > > > met-6.0/bin/regrid_data_plane \
>> > > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.tm00.grib2
>> > > > gfs_APCP_06_nam_nest_domain.nc \
>> > > > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
>> -interp_pnts
>> > 2
>> > > >
>> > > > # Run plot_data_plane to plot the 6-hour precip both
before/after
>> the
>> > > > interpolation
>> > > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
>> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
level="A06";' \
>> > > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
>> > > >
>> > > > met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
>> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
level="(*,*)";' \
>> > > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
>> > > >
>> > > > The resulting plots (converted to png) are attached to this
email.
>> > > >
>> > > > In the regridding step, I defined the "to_grid" by providing
the
>> path
>> > to
>> > > a
>> > > > file that's already on the grid I'd like to use. MET reads
that
>> file,
>> > > > grabs the grid definition, and uses that in the regridding
step.
>> There
>> > > are
>> > > > multiple ways of defining "to_grid" (as described in
>> > > > met-6.0/data/config/README and the User's Guide), but I find
>> providing
>> > a
>> > > > path to a sample file is pretty convenient.
>> > > >
>> > > > Or, even more convenient, if I'm running grid_stat for
example and
>> pass
>> > > in
>> > > > the GFS as the forecast and NAM as the observation, I can
just say
>> > > "to_grid
>> > > > = OBS;" and grid_stat will automatically regrid the forecast
data
>> (GFS)
>> > > to
>> > > > the observation grid (NAM).
>> > > >
>> > > > So, long story short, the following would work:
>> > > >
>> > > > regrid = {
>> > > > to_grid = "/path/to/sample/NDFD/data/file"; // Or set
this to
>> > FCST
>> > > > or
>> > > > OBS if that's easier
>> > > > method = BUDGET;
>> > > > width = 2; // Budget should always have a width of
2 to
>> use
>> > > the 4
>> > > > closest points
>> > > > vld_thresh = 0.5;
>> > > > }
>> > > >
>> > > > Hope that helps get you going.
>> > > >
>> > > > Thanks,
>> > > > John
>> > > >
>> > > >
>> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA
Affiliate via
>> RT <
>> > > > met_help at ucar.edu> wrote:
>> > > >
>> > > > >
>> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
>> > > > > Transaction: Ticket created by daniel.plumb at noaa.gov
>> > > > > Queue: met_help
>> > > > > Subject: Forecast Regrid
>> > > > > Owner: Nobody
>> > > > > Requestors: daniel.plumb at noaa.gov
>> > > > > Status: new
>> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
>> > Ticket/Display.html?id=81369
>> > > >
>> > > > >
>> > > > >
>> > > > > Hello,
>> > > > >
>> > > > > I am a new user to the MET program, and attempting to
regrid a
>> > 1051x813
>> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km
using the
>> > > budget
>> > > > > regridding method. I am wondering if you have a simple,
>> > straightforward
>> > > > > example of how to do this in MET. So far, I have managed to
figure
>> > out
>> > > > this
>> > > > > is needed in the config file:
>> > > > >
>> > > > > regrid = {
>> > > > > to_grid = NONE;
>> > > > > method = BUDGET;
>> > > > > width = 1;
>> > > > > vld_thresh = 0.5;
>> > > > > }
>> > > > >
>> > > > > Any other help would be much appreciated.
>> > > > >
>> > > > > Thanks in advance!
>> > > > >
>> > > > > Dan
>> > > > >
>> > > > > --
>> > > > > Daniel W. Plumb, Jr.
>> > > > > Scientific Programmer/Meteorologist
>> > > > > Ace Info Solutions, Inc. (contractor)
>> > > > > National Weather Service (NWS)- Meteorological Development
Lab
>> (MDL)
>> > > > >
>> > > > > Location:
>> > > > > SSMC2 – Room 10124
>> > > > > 1325 East West Highway
>> > > > > Silver Spring, MD 20910
>> > > > > Tel: (301) 427-9453
>> > > > >
>> > > > >
>> > > >
>> > > >
>> > >
>> > >
>> > > --
>> > > Daniel W. Plumb, Jr.
>> > > Scientific Programmer/Meteorologist
>> > > Ace Info Solutions, Inc. (contractor)
>> > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
>> > >
>> > > Location:
>> > > SSMC2 – Room 10124
>> > > 1325 East West Highway
>> > > Silver Spring, MD 20910
>> > > Tel: (301) 427-9453
>> > >
>> > >
>> >
>> >
>>
>>
>
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
--
Daniel W. Plumb, Jr.
Scientific Programmer/Meteorologist
Ace Info Solutions, Inc. (contractor)
National Weather Service (NWS)- Meteorological Development Lab (MDL)
Location:
SSMC2 – Room 10124
1325 East West Highway
Silver Spring, MD 20910
Tel: (301) 427-9453
------------------------------------------------
Subject: Forecast Regrid
From: Julie Prestopnik
Time: Fri Jul 28 15:34:44 2017
Hi Dan.
Thank you for sending your second email. We never received your first
one,
possibly because the attachments were too large.
Unfortunately, MET does not currently support polar stereographic
grids for
CF compliant files, but it does support polar stereographic grids for
GRIB1, GRIB2, and MET-NetCDF format files.
We are aware that this functionality is desired and are planning on
adding
support for it in METv6.1. When we have a development version of 6.1
ready, you could try it out if you'd like.
I wish this were more helpful to you right now...
Julie
On Fri, Jul 28, 2017 at 2:05 PM, Daniel Plumb - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
> Please let me know if you were able to get all of my attachments.
I'm not
> sure if they all went through.
>
> Thanks,
> Dan
>
> On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA Affiliate <
> daniel.plumb at noaa.gov> wrote:
>
> > Hi Julie and John,
> >
> > Thank you for the help so far. I am able to run regrid_data_plane
running
> > now. However, it is not giving me the desired output I would like
right
> > now. I have my command below and have attached the input dataset
> > (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the dataset
with
> > the desired grid specs I want (ndfdcogb2p517072700qpfd13_006
> > <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2 format*):
> >
> > regrid_data_plane qpf_06hr_2017072512f006Z.nc
> > ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field
'name="APCP";
> > level="A06";' -method BUDGET -width 2
> >
> > DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> > ERROR :
> > ERROR : NcCfFile::get_grid_mapping_polar_stereographic() -> Polar
> > stereographic grid not handled in MET.
> > ERROR :
> >
> > Is MET not able to handle polar stereographic grids? Or maybe the
NetCDF
> > library is too old on WCOSS (dated: Dec 5 2012)?
> >
> > Thanks,
> > Dan
> >
> >
> >
> >
> > On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu> wrote:
> >
> >> I only have developer access to WCOSS, but I don't have access to
WCOSS
> >> (gyre) at this time, as I believe it is still down for
developers.
> >>
> >> Please try running "module purge" and then run:
> >>
> >> - module use /usrx/local/dev/modulefiles
> >> - module load met
> >>
> >> and see if that resolves the problem. I believe there is a
conflict
> with
> >> the intel module you already have loaded. I know purging isn't
the best
> >> method to solve this problem, however, I do believe that should
work.
> >>
> >> Please let me know if you have any questions.
> >>
> >> Thanks,
> >>
> >> Julie
> >>
> >>
> >>
> >>
> >>
> >> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >> >
> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> >> >
> >> > Dan,
> >> >
> >> > I don't have access to WCOSS to replicate this problem, but by
> colleague
> >> > Julie does.
> >> >
> >> > Julie, could you please take a look at this?
> >> >
> >> > Thanks,
> >> > John
> >> >
> >> >
> >> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA Affiliate
via RT
> <
> >> > met_help at ucar.edu> wrote:
> >> >
> >> > >
> >> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
>
> >> > >
> >> > > Hi John,
> >> > >
> >> > > Thank you very much for the reply and examples.
> >> > >
> >> > > I am trying to replicate you example for regrid_data_plane
(which is
> >> > > exactly the example I needed), but there seems to be
something wrong
> >> with
> >> > > MET on WCOSS tide. When I run "module load NetCDF ics/15.0.3
udunits
> >> > ncview
> >> > > met/6.0 grib_util/v1.0.3", it is giving me an error saying:
> >> > >
> >> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts with
the
> >> > currently
> >> > > loaded module(s) 'ics/12.1'
> >> > > ics/15.0.3(16):ERROR:102: Tcl command execution failed:
conflict
> ics"
> >> > >
> >> > > I'm looking at the instructions under
> >> > > http://www.dtcenter.org/met/users/support/online_tutorial/
> >> > ncwcp_201704.php
> >> > > .
> >> > >
> >> > > Thanks,
> >> > > Dan
> >> > >
> >> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT <
> >> > > met_help at ucar.edu> wrote:
> >> > >
> >> > > > Hi Dan,
> >> > > >
> >> > > > Tara gave me a heads up that you've started using MET at
MDL.
> We're
> >> > > happy
> >> > > > to answer any questions you have or issues that come up.
> >> > > >
> >> > > > I see that you'd like to use the regridding feature within
the MET
> >> > config
> >> > > > files. I'd actually suggest starting by running the
> >> regrid_data_plane
> >> > > tool
> >> > > > first, and once you understand what it's doing, move those
> settings
> >> > into
> >> > > > the config file. Both the regrid_data_plane tool and the
regrid
> >> > section
> >> > > of
> >> > > > the MET config files call the same underlying library code
anyway.
> >> > > >
> >> > > > I was hoping to grab some sample MDL data to demonstrate
some
> >> commands,
> >> > > but
> >> > > > all the links on this website aren't working:
> >> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> >> > > >
> >> > > > So instead of using HRAP and NDFD data, I'll use sample GFS
and
> NAM
> >> > data.
> >> > > >
> >> > > > # Get sample GFS and NAM files
> >> > > > wget
> >> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> >> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> >> > > > wget
> >> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> >> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> >> > > >
> >> > > > # Run regrid_data_plane to regrid 6-hour precip from GFS to
NAM
> nest
> >> > > domain
> >> > > > met-6.0/bin/regrid_data_plane \
> >> > > > gfs.t00z.pgrb2.0p25.f012 nam.t12z.conusnest.hiresf12.
> tm00.grib2
> >> > > > gfs_APCP_06_nam_nest_domain.nc \
> >> > > > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
> >> -interp_pnts
> >> > 2
> >> > > >
> >> > > > # Run plot_data_plane to plot the 6-hour precip both
before/after
> >> the
> >> > > > interpolation
> >> > > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> >> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
level="A06";' \
> >> > > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
> >> > > >
> >> > > > met-6.0/bin/plot_data_plane gfs_APCP_06_nam_nest_domain.nc
> >> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
level="(*,*)";' \
> >> > > > -color_table met-
6.0/share/met/colortables/wrf_precip2.ctable
> >> > > >
> >> > > > The resulting plots (converted to png) are attached to this
email.
> >> > > >
> >> > > > In the regridding step, I defined the "to_grid" by
providing the
> >> path
> >> > to
> >> > > a
> >> > > > file that's already on the grid I'd like to use. MET reads
that
> >> file,
> >> > > > grabs the grid definition, and uses that in the regridding
step.
> >> There
> >> > > are
> >> > > > multiple ways of defining "to_grid" (as described in
> >> > > > met-6.0/data/config/README and the User's Guide), but I
find
> >> providing
> >> > a
> >> > > > path to a sample file is pretty convenient.
> >> > > >
> >> > > > Or, even more convenient, if I'm running grid_stat for
example and
> >> pass
> >> > > in
> >> > > > the GFS as the forecast and NAM as the observation, I can
just say
> >> > > "to_grid
> >> > > > = OBS;" and grid_stat will automatically regrid the
forecast data
> >> (GFS)
> >> > > to
> >> > > > the observation grid (NAM).
> >> > > >
> >> > > > So, long story short, the following would work:
> >> > > >
> >> > > > regrid = {
> >> > > > to_grid = "/path/to/sample/NDFD/data/file"; // Or set
this
> to
> >> > FCST
> >> > > > or
> >> > > > OBS if that's easier
> >> > > > method = BUDGET;
> >> > > > width = 2; // Budget should always have a width
of 2 to
> >> use
> >> > > the 4
> >> > > > closest points
> >> > > > vld_thresh = 0.5;
> >> > > > }
> >> > > >
> >> > > > Hope that helps get you going.
> >> > > >
> >> > > > Thanks,
> >> > > > John
> >> > > >
> >> > > >
> >> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA
Affiliate via
> >> RT <
> >> > > > met_help at ucar.edu> wrote:
> >> > > >
> >> > > > >
> >> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> >> > > > > Transaction: Ticket created by daniel.plumb at noaa.gov
> >> > > > > Queue: met_help
> >> > > > > Subject: Forecast Regrid
> >> > > > > Owner: Nobody
> >> > > > > Requestors: daniel.plumb at noaa.gov
> >> > > > > Status: new
> >> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> >> > Ticket/Display.html?id=81369
> >> > > >
> >> > > > >
> >> > > > >
> >> > > > > Hello,
> >> > > > >
> >> > > > > I am a new user to the MET program, and attempting to
regrid a
> >> > 1051x813
> >> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54 km
using
> the
> >> > > budget
> >> > > > > regridding method. I am wondering if you have a simple,
> >> > straightforward
> >> > > > > example of how to do this in MET. So far, I have managed
to
> figure
> >> > out
> >> > > > this
> >> > > > > is needed in the config file:
> >> > > > >
> >> > > > > regrid = {
> >> > > > > to_grid = NONE;
> >> > > > > method = BUDGET;
> >> > > > > width = 1;
> >> > > > > vld_thresh = 0.5;
> >> > > > > }
> >> > > > >
> >> > > > > Any other help would be much appreciated.
> >> > > > >
> >> > > > > Thanks in advance!
> >> > > > >
> >> > > > > Dan
> >> > > > >
> >> > > > > --
> >> > > > > Daniel W. Plumb, Jr.
> >> > > > > Scientific Programmer/Meteorologist
> >> > > > > Ace Info Solutions, Inc. (contractor)
> >> > > > > National Weather Service (NWS)- Meteorological
Development Lab
> >> (MDL)
> >> > > > >
> >> > > > > Location:
> >> > > > > SSMC2 – Room 10124
> >> > > > > 1325 East West Highway
> >> > > > > Silver Spring, MD 20910
> >> > > > > Tel: (301) 427-9453
> >> > > > >
> >> > > > >
> >> > > >
> >> > > >
> >> > >
> >> > >
> >> > > --
> >> > > Daniel W. Plumb, Jr.
> >> > > Scientific Programmer/Meteorologist
> >> > > Ace Info Solutions, Inc. (contractor)
> >> > > National Weather Service (NWS)- Meteorological Development
Lab (MDL)
> >> > >
> >> > > Location:
> >> > > SSMC2 – Room 10124
> >> > > 1325 East West Highway
> >> > > Silver Spring, MD 20910
> >> > > Tel: (301) 427-9453
> >> > >
> >> > >
> >> >
> >> >
> >>
> >>
> >
> >
> > --
> > Daniel W. Plumb, Jr.
> > Scientific Programmer/Meteorologist
> > Ace Info Solutions, Inc. (contractor)
> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> >
> > Location:
> > SSMC2 – Room 10124
> > 1325 East West Highway
> > Silver Spring, MD 20910
> > Tel: (301) 427-9453
> >
>
>
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
>
------------------------------------------------
Subject: Forecast Regrid
From: Daniel Plumb - NOAA Affiliate
Time: Fri Jul 28 15:40:37 2017
Hi Julie,
That would be very helpful if I could assess MET v6.1. Is this
currently
installed on WCOSS?
Thanks,
Dan
On Fri, Jul 28, 2017 at 5:34 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Dan.
>
> Thank you for sending your second email. We never received your
first one,
> possibly because the attachments were too large.
>
> Unfortunately, MET does not currently support polar stereographic
grids for
> CF compliant files, but it does support polar stereographic grids
for
> GRIB1, GRIB2, and MET-NetCDF format files.
>
> We are aware that this functionality is desired and are planning on
adding
> support for it in METv6.1. When we have a development version of
6.1
> ready, you could try it out if you'd like.
>
> I wish this were more helpful to you right now...
>
> Julie
>
>
>
> On Fri, Jul 28, 2017 at 2:05 PM, Daniel Plumb - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> >
> > Please let me know if you were able to get all of my attachments.
I'm not
> > sure if they all went through.
> >
> > Thanks,
> > Dan
> >
> > On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA Affiliate <
> > daniel.plumb at noaa.gov> wrote:
> >
> > > Hi Julie and John,
> > >
> > > Thank you for the help so far. I am able to run
regrid_data_plane
> running
> > > now. However, it is not giving me the desired output I would
like right
> > > now. I have my command below and have attached the input dataset
> > > (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the
dataset with
> > > the desired grid specs I want (ndfdcogb2p517072700qpfd13_006
> > > <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2 format*):
> > >
> > > regrid_data_plane qpf_06hr_2017072512f006Z.nc
> > > ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field
'name="APCP";
> > > level="A06";' -method BUDGET -width 2
> > >
> > > DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> > > ERROR :
> > > ERROR : NcCfFile::get_grid_mapping_polar_stereographic() ->
Polar
> > > stereographic grid not handled in MET.
> > > ERROR :
> > >
> > > Is MET not able to handle polar stereographic grids? Or maybe
the
> NetCDF
> > > library is too old on WCOSS (dated: Dec 5 2012)?
> > >
> > > Thanks,
> > > Dan
> > >
> > >
> > >
> > >
> > > On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> I only have developer access to WCOSS, but I don't have access
to
> WCOSS
> > >> (gyre) at this time, as I believe it is still down for
developers.
> > >>
> > >> Please try running "module purge" and then run:
> > >>
> > >> - module use /usrx/local/dev/modulefiles
> > >> - module load met
> > >>
> > >> and see if that resolves the problem. I believe there is a
conflict
> > with
> > >> the intel module you already have loaded. I know purging isn't
the
> best
> > >> method to solve this problem, however, I do believe that should
work.
> > >>
> > >> Please let me know if you have any questions.
> > >>
> > >> Thanks,
> > >>
> > >> Julie
> > >>
> > >>
> > >>
> > >>
> > >>
> > >> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
> > >> met_help at ucar.edu> wrote:
> > >>
> > >> >
> > >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
>
> > >> >
> > >> > Dan,
> > >> >
> > >> > I don't have access to WCOSS to replicate this problem, but
by
> > colleague
> > >> > Julie does.
> > >> >
> > >> > Julie, could you please take a look at this?
> > >> >
> > >> > Thanks,
> > >> > John
> > >> >
> > >> >
> > >> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA
Affiliate via
> RT
> > <
> > >> > met_help at ucar.edu> wrote:
> > >> >
> > >> > >
> > >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > >> > >
> > >> > > Hi John,
> > >> > >
> > >> > > Thank you very much for the reply and examples.
> > >> > >
> > >> > > I am trying to replicate you example for regrid_data_plane
(which
> is
> > >> > > exactly the example I needed), but there seems to be
something
> wrong
> > >> with
> > >> > > MET on WCOSS tide. When I run "module load NetCDF
ics/15.0.3
> udunits
> > >> > ncview
> > >> > > met/6.0 grib_util/v1.0.3", it is giving me an error saying:
> > >> > >
> > >> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts
with the
> > >> > currently
> > >> > > loaded module(s) 'ics/12.1'
> > >> > > ics/15.0.3(16):ERROR:102: Tcl command execution failed:
conflict
> > ics"
> > >> > >
> > >> > > I'm looking at the instructions under
> > >> > > http://www.dtcenter.org/met/users/support/online_tutorial/
> > >> > ncwcp_201704.php
> > >> > > .
> > >> > >
> > >> > > Thanks,
> > >> > > Dan
> > >> > >
> > >> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via RT
<
> > >> > > met_help at ucar.edu> wrote:
> > >> > >
> > >> > > > Hi Dan,
> > >> > > >
> > >> > > > Tara gave me a heads up that you've started using MET at
MDL.
> > We're
> > >> > > happy
> > >> > > > to answer any questions you have or issues that come up.
> > >> > > >
> > >> > > > I see that you'd like to use the regridding feature
within the
> MET
> > >> > config
> > >> > > > files. I'd actually suggest starting by running the
> > >> regrid_data_plane
> > >> > > tool
> > >> > > > first, and once you understand what it's doing, move
those
> > settings
> > >> > into
> > >> > > > the config file. Both the regrid_data_plane tool and the
regrid
> > >> > section
> > >> > > of
> > >> > > > the MET config files call the same underlying library
code
> anyway.
> > >> > > >
> > >> > > > I was hoping to grab some sample MDL data to demonstrate
some
> > >> commands,
> > >> > > but
> > >> > > > all the links on this website aren't working:
> > >> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> > >> > > >
> > >> > > > So instead of using HRAP and NDFD data, I'll use sample
GFS and
> > NAM
> > >> > data.
> > >> > > >
> > >> > > > # Get sample GFS and NAM files
> > >> > > > wget
> > >> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > >> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > >> > > > wget
> > >> > > > http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > >> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> > >> > > >
> > >> > > > # Run regrid_data_plane to regrid 6-hour precip from GFS
to NAM
> > nest
> > >> > > domain
> > >> > > > met-6.0/bin/regrid_data_plane \
> > >> > > > gfs.t00z.pgrb2.0p25.f012 nam.t12z.conusnest.hiresf12.
> > tm00.grib2
> > >> > > > gfs_APCP_06_nam_nest_domain.nc \
> > >> > > > -field 'name="APCP"; level="A06";' -interp_mthd BUDGET
> > >> -interp_pnts
> > >> > 2
> > >> > > >
> > >> > > > # Run plot_data_plane to plot the 6-hour precip both
> before/after
> > >> the
> > >> > > > interpolation
> > >> > > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> > >> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
> level="A06";' \
> > >> > > > -color_table met-6.0/share/met/colortables/
> wrf_precip2.ctable
> > >> > > >
> > >> > > > met-6.0/bin/plot_data_plane
gfs_APCP_06_nam_nest_domain.nc
> > >> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
level="(*,*)";'
> \
> > >> > > > -color_table met-6.0/share/met/colortables/
> wrf_precip2.ctable
> > >> > > >
> > >> > > > The resulting plots (converted to png) are attached to
this
> email.
> > >> > > >
> > >> > > > In the regridding step, I defined the "to_grid" by
providing the
> > >> path
> > >> > to
> > >> > > a
> > >> > > > file that's already on the grid I'd like to use. MET
reads that
> > >> file,
> > >> > > > grabs the grid definition, and uses that in the
regridding step.
> > >> There
> > >> > > are
> > >> > > > multiple ways of defining "to_grid" (as described in
> > >> > > > met-6.0/data/config/README and the User's Guide), but I
find
> > >> providing
> > >> > a
> > >> > > > path to a sample file is pretty convenient.
> > >> > > >
> > >> > > > Or, even more convenient, if I'm running grid_stat for
example
> and
> > >> pass
> > >> > > in
> > >> > > > the GFS as the forecast and NAM as the observation, I can
just
> say
> > >> > > "to_grid
> > >> > > > = OBS;" and grid_stat will automatically regrid the
forecast
> data
> > >> (GFS)
> > >> > > to
> > >> > > > the observation grid (NAM).
> > >> > > >
> > >> > > > So, long story short, the following would work:
> > >> > > >
> > >> > > > regrid = {
> > >> > > > to_grid = "/path/to/sample/NDFD/data/file"; // Or
set
> this
> > to
> > >> > FCST
> > >> > > > or
> > >> > > > OBS if that's easier
> > >> > > > method = BUDGET;
> > >> > > > width = 2; // Budget should always have a width
of 2
> to
> > >> use
> > >> > > the 4
> > >> > > > closest points
> > >> > > > vld_thresh = 0.5;
> > >> > > > }
> > >> > > >
> > >> > > > Hope that helps get you going.
> > >> > > >
> > >> > > > Thanks,
> > >> > > > John
> > >> > > >
> > >> > > >
> > >> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA
Affiliate
> via
> > >> RT <
> > >> > > > met_help at ucar.edu> wrote:
> > >> > > >
> > >> > > > >
> > >> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted upon.
> > >> > > > > Transaction: Ticket created by daniel.plumb at noaa.gov
> > >> > > > > Queue: met_help
> > >> > > > > Subject: Forecast Regrid
> > >> > > > > Owner: Nobody
> > >> > > > > Requestors: daniel.plumb at noaa.gov
> > >> > > > > Status: new
> > >> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > >> > Ticket/Display.html?id=81369
> > >> > > >
> > >> > > > >
> > >> > > > >
> > >> > > > > Hello,
> > >> > > > >
> > >> > > > > I am a new user to the MET program, and attempting to
regrid a
> > >> > 1051x813
> > >> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54
km using
> > the
> > >> > > budget
> > >> > > > > regridding method. I am wondering if you have a simple,
> > >> > straightforward
> > >> > > > > example of how to do this in MET. So far, I have
managed to
> > figure
> > >> > out
> > >> > > > this
> > >> > > > > is needed in the config file:
> > >> > > > >
> > >> > > > > regrid = {
> > >> > > > > to_grid = NONE;
> > >> > > > > method = BUDGET;
> > >> > > > > width = 1;
> > >> > > > > vld_thresh = 0.5;
> > >> > > > > }
> > >> > > > >
> > >> > > > > Any other help would be much appreciated.
> > >> > > > >
> > >> > > > > Thanks in advance!
> > >> > > > >
> > >> > > > > Dan
> > >> > > > >
> > >> > > > > --
> > >> > > > > Daniel W. Plumb, Jr.
> > >> > > > > Scientific Programmer/Meteorologist
> > >> > > > > Ace Info Solutions, Inc. (contractor)
> > >> > > > > National Weather Service (NWS)- Meteorological
Development Lab
> > >> (MDL)
> > >> > > > >
> > >> > > > > Location:
> > >> > > > > SSMC2 – Room 10124
> > >> > > > > 1325 East West Highway
> > >> > > > > Silver Spring, MD 20910
> > >> > > > > Tel: (301) 427-9453
> > >> > > > >
> > >> > > > >
> > >> > > >
> > >> > > >
> > >> > >
> > >> > >
> > >> > > --
> > >> > > Daniel W. Plumb, Jr.
> > >> > > Scientific Programmer/Meteorologist
> > >> > > Ace Info Solutions, Inc. (contractor)
> > >> > > National Weather Service (NWS)- Meteorological Development
Lab
> (MDL)
> > >> > >
> > >> > > Location:
> > >> > > SSMC2 – Room 10124
> > >> > > 1325 East West Highway
> > >> > > Silver Spring, MD 20910
> > >> > > Tel: (301) 427-9453
> > >> > >
> > >> > >
> > >> >
> > >> >
> > >>
> > >>
> > >
> > >
> > > --
> > > Daniel W. Plumb, Jr.
> > > Scientific Programmer/Meteorologist
> > > Ace Info Solutions, Inc. (contractor)
> > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > >
> > > Location:
> > > SSMC2 – Room 10124
> > > 1325 East West Highway
> > > Silver Spring, MD 20910
> > > Tel: (301) 427-9453
> > >
> >
> >
> >
> > --
> > Daniel W. Plumb, Jr.
> > Scientific Programmer/Meteorologist
> > Ace Info Solutions, Inc. (contractor)
> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> >
> > Location:
> > SSMC2 – Room 10124
> > 1325 East West Highway
> > Silver Spring, MD 20910
> > Tel: (301) 427-9453
> >
> >
>
>
--
Daniel W. Plumb, Jr.
Scientific Programmer/Meteorologist
Ace Info Solutions, Inc. (contractor)
National Weather Service (NWS)- Meteorological Development Lab (MDL)
Location:
SSMC2 – Room 10124
1325 East West Highway
Silver Spring, MD 20910
Tel: (301) 427-9453
------------------------------------------------
Subject: Forecast Regrid
From: John Halley Gotway
Time: Fri Jul 28 16:37:42 2017
Dan,
I worked on some code in the development version of MET to parse polar
stereographic grids from CF-compliant NetCDF files, basically adding
in a
definition for this function:
NcCfFile::get_grid_mapping_polar_stereographic()
The good news is that it *seemed* to work. I ran the following
plot_data_plane command:
plot_data_plane qpf_06hr_2017072512f006Z.nc qpf.ps
'name="Forecast_total_precipitation"; level="(0,*,*)";'
And I've attached the resulting image. The image certainly looks
reasonable.
But I'd really like to step through this code with Randy Bullock, the
author of the grid/projection-related code in MET. Once we review the
changes next week, we can push them out to you.
Since this is a pretty isolated change... so one option would be a
bugfix
for met-6.0. Technically this is a "new feature" rather than a
"bugfix",
but that'd be the quickest way to get you this functionality.
Please take a look at this image and let me know what you think.
Thanks,
John
On Fri, Jul 28, 2017 at 3:40 PM, Daniel Plumb - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
> Hi Julie,
>
> That would be very helpful if I could assess MET v6.1. Is this
currently
> installed on WCOSS?
>
> Thanks,
> Dan
>
> On Fri, Jul 28, 2017 at 5:34 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Dan.
> >
> > Thank you for sending your second email. We never received your
first
> one,
> > possibly because the attachments were too large.
> >
> > Unfortunately, MET does not currently support polar stereographic
grids
> for
> > CF compliant files, but it does support polar stereographic grids
for
> > GRIB1, GRIB2, and MET-NetCDF format files.
> >
> > We are aware that this functionality is desired and are planning
on
> adding
> > support for it in METv6.1. When we have a development version of
6.1
> > ready, you could try it out if you'd like.
> >
> > I wish this were more helpful to you right now...
> >
> > Julie
> >
> >
> >
> > On Fri, Jul 28, 2017 at 2:05 PM, Daniel Plumb - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > >
> > > Please let me know if you were able to get all of my
attachments. I'm
> not
> > > sure if they all went through.
> > >
> > > Thanks,
> > > Dan
> > >
> > > On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA Affiliate <
> > > daniel.plumb at noaa.gov> wrote:
> > >
> > > > Hi Julie and John,
> > > >
> > > > Thank you for the help so far. I am able to run
regrid_data_plane
> > running
> > > > now. However, it is not giving me the desired output I would
like
> right
> > > > now. I have my command below and have attached the input
dataset
> > > > (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the
dataset
> with
> > > > the desired grid specs I want (ndfdcogb2p517072700qpfd13_006
> > > > <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2
format*):
> > > >
> > > > regrid_data_plane qpf_06hr_2017072512f006Z.nc
> > > > ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field
'name="APCP";
> > > > level="A06";' -method BUDGET -width 2
> > > >
> > > > DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> > > > ERROR :
> > > > ERROR : NcCfFile::get_grid_mapping_polar_stereographic() ->
Polar
> > > > stereographic grid not handled in MET.
> > > > ERROR :
> > > >
> > > > Is MET not able to handle polar stereographic grids? Or maybe
the
> > NetCDF
> > > > library is too old on WCOSS (dated: Dec 5 2012)?
> > > >
> > > > Thanks,
> > > > Dan
> > > >
> > > >
> > > >
> > > >
> > > > On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> I only have developer access to WCOSS, but I don't have
access to
> > WCOSS
> > > >> (gyre) at this time, as I believe it is still down for
developers.
> > > >>
> > > >> Please try running "module purge" and then run:
> > > >>
> > > >> - module use /usrx/local/dev/modulefiles
> > > >> - module load met
> > > >>
> > > >> and see if that resolves the problem. I believe there is a
conflict
> > > with
> > > >> the intel module you already have loaded. I know purging
isn't the
> > best
> > > >> method to solve this problem, however, I do believe that
should
> work.
> > > >>
> > > >> Please let me know if you have any questions.
> > > >>
> > > >> Thanks,
> > > >>
> > > >> Julie
> > > >>
> > > >>
> > > >>
> > > >>
> > > >>
> > > >> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT <
> > > >> met_help at ucar.edu> wrote:
> > > >>
> > > >> >
> > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > > >> >
> > > >> > Dan,
> > > >> >
> > > >> > I don't have access to WCOSS to replicate this problem, but
by
> > > colleague
> > > >> > Julie does.
> > > >> >
> > > >> > Julie, could you please take a look at this?
> > > >> >
> > > >> > Thanks,
> > > >> > John
> > > >> >
> > > >> >
> > > >> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA
Affiliate via
> > RT
> > > <
> > > >> > met_help at ucar.edu> wrote:
> > > >> >
> > > >> > >
> > > >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > > >> > >
> > > >> > > Hi John,
> > > >> > >
> > > >> > > Thank you very much for the reply and examples.
> > > >> > >
> > > >> > > I am trying to replicate you example for
regrid_data_plane
> (which
> > is
> > > >> > > exactly the example I needed), but there seems to be
something
> > wrong
> > > >> with
> > > >> > > MET on WCOSS tide. When I run "module load NetCDF
ics/15.0.3
> > udunits
> > > >> > ncview
> > > >> > > met/6.0 grib_util/v1.0.3", it is giving me an error
saying:
> > > >> > >
> > > >> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3' conflicts
with
> the
> > > >> > currently
> > > >> > > loaded module(s) 'ics/12.1'
> > > >> > > ics/15.0.3(16):ERROR:102: Tcl command execution failed:
conflict
> > > ics"
> > > >> > >
> > > >> > > I'm looking at the instructions under
> > > >> > >
http://www.dtcenter.org/met/users/support/online_tutorial/
> > > >> > ncwcp_201704.php
> > > >> > > .
> > > >> > >
> > > >> > > Thanks,
> > > >> > > Dan
> > > >> > >
> > > >> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via
RT <
> > > >> > > met_help at ucar.edu> wrote:
> > > >> > >
> > > >> > > > Hi Dan,
> > > >> > > >
> > > >> > > > Tara gave me a heads up that you've started using MET
at MDL.
> > > We're
> > > >> > > happy
> > > >> > > > to answer any questions you have or issues that come
up.
> > > >> > > >
> > > >> > > > I see that you'd like to use the regridding feature
within the
> > MET
> > > >> > config
> > > >> > > > files. I'd actually suggest starting by running the
> > > >> regrid_data_plane
> > > >> > > tool
> > > >> > > > first, and once you understand what it's doing, move
those
> > > settings
> > > >> > into
> > > >> > > > the config file. Both the regrid_data_plane tool and
the
> regrid
> > > >> > section
> > > >> > > of
> > > >> > > > the MET config files call the same underlying library
code
> > anyway.
> > > >> > > >
> > > >> > > > I was hoping to grab some sample MDL data to
demonstrate some
> > > >> commands,
> > > >> > > but
> > > >> > > > all the links on this website aren't working:
> > > >> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> > > >> > > >
> > > >> > > > So instead of using HRAP and NDFD data, I'll use sample
GFS
> and
> > > NAM
> > > >> > data.
> > > >> > > >
> > > >> > > > # Get sample GFS and NAM files
> > > >> > > > wget
> > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > > >> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > > >> > > > wget
> > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > > >> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> > > >> > > >
> > > >> > > > # Run regrid_data_plane to regrid 6-hour precip from
GFS to
> NAM
> > > nest
> > > >> > > domain
> > > >> > > > met-6.0/bin/regrid_data_plane \
> > > >> > > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.
> > > tm00.grib2
> > > >> > > > gfs_APCP_06_nam_nest_domain.nc \
> > > >> > > > -field 'name="APCP"; level="A06";' -interp_mthd
BUDGET
> > > >> -interp_pnts
> > > >> > 2
> > > >> > > >
> > > >> > > > # Run plot_data_plane to plot the 6-hour precip both
> > before/after
> > > >> the
> > > >> > > > interpolation
> > > >> > > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> > > >> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
> > level="A06";' \
> > > >> > > > -color_table met-6.0/share/met/colortables/
> > wrf_precip2.ctable
> > > >> > > >
> > > >> > > > met-6.0/bin/plot_data_plane
gfs_APCP_06_nam_nest_domain.nc
> > > >> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
> level="(*,*)";'
> > \
> > > >> > > > -color_table met-6.0/share/met/colortables/
> > wrf_precip2.ctable
> > > >> > > >
> > > >> > > > The resulting plots (converted to png) are attached to
this
> > email.
> > > >> > > >
> > > >> > > > In the regridding step, I defined the "to_grid" by
providing
> the
> > > >> path
> > > >> > to
> > > >> > > a
> > > >> > > > file that's already on the grid I'd like to use. MET
reads
> that
> > > >> file,
> > > >> > > > grabs the grid definition, and uses that in the
regridding
> step.
> > > >> There
> > > >> > > are
> > > >> > > > multiple ways of defining "to_grid" (as described in
> > > >> > > > met-6.0/data/config/README and the User's Guide), but I
find
> > > >> providing
> > > >> > a
> > > >> > > > path to a sample file is pretty convenient.
> > > >> > > >
> > > >> > > > Or, even more convenient, if I'm running grid_stat for
example
> > and
> > > >> pass
> > > >> > > in
> > > >> > > > the GFS as the forecast and NAM as the observation, I
can just
> > say
> > > >> > > "to_grid
> > > >> > > > = OBS;" and grid_stat will automatically regrid the
forecast
> > data
> > > >> (GFS)
> > > >> > > to
> > > >> > > > the observation grid (NAM).
> > > >> > > >
> > > >> > > > So, long story short, the following would work:
> > > >> > > >
> > > >> > > > regrid = {
> > > >> > > > to_grid = "/path/to/sample/NDFD/data/file"; // Or
set
> > this
> > > to
> > > >> > FCST
> > > >> > > > or
> > > >> > > > OBS if that's easier
> > > >> > > > method = BUDGET;
> > > >> > > > width = 2; // Budget should always have a
width of 2
> > to
> > > >> use
> > > >> > > the 4
> > > >> > > > closest points
> > > >> > > > vld_thresh = 0.5;
> > > >> > > > }
> > > >> > > >
> > > >> > > > Hope that helps get you going.
> > > >> > > >
> > > >> > > > Thanks,
> > > >> > > > John
> > > >> > > >
> > > >> > > >
> > > >> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA
Affiliate
> > via
> > > >> RT <
> > > >> > > > met_help at ucar.edu> wrote:
> > > >> > > >
> > > >> > > > >
> > > >> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted
upon.
> > > >> > > > > Transaction: Ticket created by daniel.plumb at noaa.gov
> > > >> > > > > Queue: met_help
> > > >> > > > > Subject: Forecast Regrid
> > > >> > > > > Owner: Nobody
> > > >> > > > > Requestors: daniel.plumb at noaa.gov
> > > >> > > > > Status: new
> > > >> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > >> > Ticket/Display.html?id=81369
> > > >> > > >
> > > >> > > > >
> > > >> > > > >
> > > >> > > > > Hello,
> > > >> > > > >
> > > >> > > > > I am a new user to the MET program, and attempting to
> regrid a
> > > >> > 1051x813
> > > >> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at 2.54
km
> using
> > > the
> > > >> > > budget
> > > >> > > > > regridding method. I am wondering if you have a
simple,
> > > >> > straightforward
> > > >> > > > > example of how to do this in MET. So far, I have
managed to
> > > figure
> > > >> > out
> > > >> > > > this
> > > >> > > > > is needed in the config file:
> > > >> > > > >
> > > >> > > > > regrid = {
> > > >> > > > > to_grid = NONE;
> > > >> > > > > method = BUDGET;
> > > >> > > > > width = 1;
> > > >> > > > > vld_thresh = 0.5;
> > > >> > > > > }
> > > >> > > > >
> > > >> > > > > Any other help would be much appreciated.
> > > >> > > > >
> > > >> > > > > Thanks in advance!
> > > >> > > > >
> > > >> > > > > Dan
> > > >> > > > >
> > > >> > > > > --
> > > >> > > > > Daniel W. Plumb, Jr.
> > > >> > > > > Scientific Programmer/Meteorologist
> > > >> > > > > Ace Info Solutions, Inc. (contractor)
> > > >> > > > > National Weather Service (NWS)- Meteorological
Development
> Lab
> > > >> (MDL)
> > > >> > > > >
> > > >> > > > > Location:
> > > >> > > > > SSMC2 – Room 10124
> > > >> > > > > 1325 East West Highway
> > > >> > > > > Silver Spring, MD 20910
> > > >> > > > > Tel: (301) 427-9453
> > > >> > > > >
> > > >> > > > >
> > > >> > > >
> > > >> > > >
> > > >> > >
> > > >> > >
> > > >> > > --
> > > >> > > Daniel W. Plumb, Jr.
> > > >> > > Scientific Programmer/Meteorologist
> > > >> > > Ace Info Solutions, Inc. (contractor)
> > > >> > > National Weather Service (NWS)- Meteorological
Development Lab
> > (MDL)
> > > >> > >
> > > >> > > Location:
> > > >> > > SSMC2 – Room 10124
> > > >> > > 1325 East West Highway
> > > >> > > Silver Spring, MD 20910
> > > >> > > Tel: (301) 427-9453
> > > >> > >
> > > >> > >
> > > >> >
> > > >> >
> > > >>
> > > >>
> > > >
> > > >
> > > > --
> > > > Daniel W. Plumb, Jr.
> > > > Scientific Programmer/Meteorologist
> > > > Ace Info Solutions, Inc. (contractor)
> > > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > > >
> > > > Location:
> > > > SSMC2 – Room 10124
> > > > 1325 East West Highway
> > > > Silver Spring, MD 20910
> > > > Tel: (301) 427-9453
> > > >
> > >
> > >
> > >
> > > --
> > > Daniel W. Plumb, Jr.
> > > Scientific Programmer/Meteorologist
> > > Ace Info Solutions, Inc. (contractor)
> > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > >
> > > Location:
> > > SSMC2 – Room 10124
> > > 1325 East West Highway
> > > Silver Spring, MD 20910
> > > Tel: (301) 427-9453
> > >
> > >
> >
> >
>
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
>
------------------------------------------------
Subject: Forecast Regrid
From: Daniel Plumb - NOAA Affiliate
Time: Mon Jul 31 08:37:43 2017
Hi John,
This image looks accurate! I look forward to the new code.
Thanks in advance,
Dan
On Fri, Jul 28, 2017 at 6:37 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Dan,
>
> I worked on some code in the development version of MET to parse
polar
> stereographic grids from CF-compliant NetCDF files, basically adding
in a
> definition for this function:
> NcCfFile::get_grid_mapping_polar_stereographic()
>
> The good news is that it *seemed* to work. I ran the following
> plot_data_plane command:
> plot_data_plane qpf_06hr_2017072512f006Z.nc qpf.ps
> 'name="Forecast_total_precipitation"; level="(0,*,*)";'
>
> And I've attached the resulting image. The image certainly looks
> reasonable.
>
> But I'd really like to step through this code with Randy Bullock,
the
> author of the grid/projection-related code in MET. Once we review
the
> changes next week, we can push them out to you.
>
> Since this is a pretty isolated change... so one option would be a
bugfix
> for met-6.0. Technically this is a "new feature" rather than a
"bugfix",
> but that'd be the quickest way to get you this functionality.
>
> Please take a look at this image and let me know what you think.
>
> Thanks,
> John
>
>
>
> On Fri, Jul 28, 2017 at 3:40 PM, Daniel Plumb - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> >
> > Hi Julie,
> >
> > That would be very helpful if I could assess MET v6.1. Is this
currently
> > installed on WCOSS?
> >
> > Thanks,
> > Dan
> >
> > On Fri, Jul 28, 2017 at 5:34 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Dan.
> > >
> > > Thank you for sending your second email. We never received your
first
> > one,
> > > possibly because the attachments were too large.
> > >
> > > Unfortunately, MET does not currently support polar
stereographic grids
> > for
> > > CF compliant files, but it does support polar stereographic
grids for
> > > GRIB1, GRIB2, and MET-NetCDF format files.
> > >
> > > We are aware that this functionality is desired and are planning
on
> > adding
> > > support for it in METv6.1. When we have a development version
of 6.1
> > > ready, you could try it out if you'd like.
> > >
> > > I wish this were more helpful to you right now...
> > >
> > > Julie
> > >
> > >
> > >
> > > On Fri, Jul 28, 2017 at 2:05 PM, Daniel Plumb - NOAA Affiliate
via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
>
> > > >
> > > > Please let me know if you were able to get all of my
attachments. I'm
> > not
> > > > sure if they all went through.
> > > >
> > > > Thanks,
> > > > Dan
> > > >
> > > > On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA Affiliate
<
> > > > daniel.plumb at noaa.gov> wrote:
> > > >
> > > > > Hi Julie and John,
> > > > >
> > > > > Thank you for the help so far. I am able to run
regrid_data_plane
> > > running
> > > > > now. However, it is not giving me the desired output I would
like
> > right
> > > > > now. I have my command below and have attached the input
dataset
> > > > > (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the
dataset
> > with
> > > > > the desired grid specs I want (ndfdcogb2p517072700qpfd13_006
> > > > > <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2
format*):
> > > > >
> > > > > regrid_data_plane qpf_06hr_2017072512f006Z.nc
> > > > > ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field
'name="APCP";
> > > > > level="A06";' -method BUDGET -width 2
> > > > >
> > > > > DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> > > > > ERROR :
> > > > > ERROR : NcCfFile::get_grid_mapping_polar_stereographic() ->
Polar
> > > > > stereographic grid not handled in MET.
> > > > > ERROR :
> > > > >
> > > > > Is MET not able to handle polar stereographic grids? Or
maybe the
> > > NetCDF
> > > > > library is too old on WCOSS (dated: Dec 5 2012)?
> > > > >
> > > > > Thanks,
> > > > > Dan
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > >> I only have developer access to WCOSS, but I don't have
access to
> > > WCOSS
> > > > >> (gyre) at this time, as I believe it is still down for
developers.
> > > > >>
> > > > >> Please try running "module purge" and then run:
> > > > >>
> > > > >> - module use /usrx/local/dev/modulefiles
> > > > >> - module load met
> > > > >>
> > > > >> and see if that resolves the problem. I believe there is a
> conflict
> > > > with
> > > > >> the intel module you already have loaded. I know purging
isn't
> the
> > > best
> > > > >> method to solve this problem, however, I do believe that
should
> > work.
> > > > >>
> > > > >> Please let me know if you have any questions.
> > > > >>
> > > > >> Thanks,
> > > > >>
> > > > >> Julie
> > > > >>
> > > > >>
> > > > >>
> > > > >>
> > > > >>
> > > > >> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via RT
<
> > > > >> met_help at ucar.edu> wrote:
> > > > >>
> > > > >> >
> > > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > > > >> >
> > > > >> > Dan,
> > > > >> >
> > > > >> > I don't have access to WCOSS to replicate this problem,
but by
> > > > colleague
> > > > >> > Julie does.
> > > > >> >
> > > > >> > Julie, could you please take a look at this?
> > > > >> >
> > > > >> > Thanks,
> > > > >> > John
> > > > >> >
> > > > >> >
> > > > >> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA
Affiliate
> via
> > > RT
> > > > <
> > > > >> > met_help at ucar.edu> wrote:
> > > > >> >
> > > > >> > >
> > > > >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
> >
> > > > >> > >
> > > > >> > > Hi John,
> > > > >> > >
> > > > >> > > Thank you very much for the reply and examples.
> > > > >> > >
> > > > >> > > I am trying to replicate you example for
regrid_data_plane
> > (which
> > > is
> > > > >> > > exactly the example I needed), but there seems to be
something
> > > wrong
> > > > >> with
> > > > >> > > MET on WCOSS tide. When I run "module load NetCDF
ics/15.0.3
> > > udunits
> > > > >> > ncview
> > > > >> > > met/6.0 grib_util/v1.0.3", it is giving me an error
saying:
> > > > >> > >
> > > > >> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3'
conflicts with
> > the
> > > > >> > currently
> > > > >> > > loaded module(s) 'ics/12.1'
> > > > >> > > ics/15.0.3(16):ERROR:102: Tcl command execution failed:
> conflict
> > > > ics"
> > > > >> > >
> > > > >> > > I'm looking at the instructions under
> > > > >> > >
http://www.dtcenter.org/met/users/support/online_tutorial/
> > > > >> > ncwcp_201704.php
> > > > >> > > .
> > > > >> > >
> > > > >> > > Thanks,
> > > > >> > > Dan
> > > > >> > >
> > > > >> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway via
RT <
> > > > >> > > met_help at ucar.edu> wrote:
> > > > >> > >
> > > > >> > > > Hi Dan,
> > > > >> > > >
> > > > >> > > > Tara gave me a heads up that you've started using MET
at
> MDL.
> > > > We're
> > > > >> > > happy
> > > > >> > > > to answer any questions you have or issues that come
up.
> > > > >> > > >
> > > > >> > > > I see that you'd like to use the regridding feature
within
> the
> > > MET
> > > > >> > config
> > > > >> > > > files. I'd actually suggest starting by running the
> > > > >> regrid_data_plane
> > > > >> > > tool
> > > > >> > > > first, and once you understand what it's doing, move
those
> > > > settings
> > > > >> > into
> > > > >> > > > the config file. Both the regrid_data_plane tool and
the
> > regrid
> > > > >> > section
> > > > >> > > of
> > > > >> > > > the MET config files call the same underlying library
code
> > > anyway.
> > > > >> > > >
> > > > >> > > > I was hoping to grab some sample MDL data to
demonstrate
> some
> > > > >> commands,
> > > > >> > > but
> > > > >> > > > all the links on this website aren't working:
> > > > >> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> > > > >> > > >
> > > > >> > > > So instead of using HRAP and NDFD data, I'll use
sample GFS
> > and
> > > > NAM
> > > > >> > data.
> > > > >> > > >
> > > > >> > > > # Get sample GFS and NAM files
> > > > >> > > > wget
> > > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > > > >> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > > > >> > > > wget
> > > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > > > >> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> > > > >> > > >
> > > > >> > > > # Run regrid_data_plane to regrid 6-hour precip from
GFS to
> > NAM
> > > > nest
> > > > >> > > domain
> > > > >> > > > met-6.0/bin/regrid_data_plane \
> > > > >> > > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.
> > > > tm00.grib2
> > > > >> > > > gfs_APCP_06_nam_nest_domain.nc \
> > > > >> > > > -field 'name="APCP"; level="A06";' -interp_mthd
BUDGET
> > > > >> -interp_pnts
> > > > >> > 2
> > > > >> > > >
> > > > >> > > > # Run plot_data_plane to plot the 6-hour precip both
> > > before/after
> > > > >> the
> > > > >> > > > interpolation
> > > > >> > > > met-6.0/bin/plot_data_plane gfs.t00z.pgrb2.0p25.f012
> > > > >> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
> > > level="A06";' \
> > > > >> > > > -color_table met-6.0/share/met/colortables/
> > > wrf_precip2.ctable
> > > > >> > > >
> > > > >> > > > met-6.0/bin/plot_data_plane
gfs_APCP_06_nam_nest_domain.nc
> > > > >> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
> > level="(*,*)";'
> > > \
> > > > >> > > > -color_table met-6.0/share/met/colortables/
> > > wrf_precip2.ctable
> > > > >> > > >
> > > > >> > > > The resulting plots (converted to png) are attached
to this
> > > email.
> > > > >> > > >
> > > > >> > > > In the regridding step, I defined the "to_grid" by
providing
> > the
> > > > >> path
> > > > >> > to
> > > > >> > > a
> > > > >> > > > file that's already on the grid I'd like to use. MET
reads
> > that
> > > > >> file,
> > > > >> > > > grabs the grid definition, and uses that in the
regridding
> > step.
> > > > >> There
> > > > >> > > are
> > > > >> > > > multiple ways of defining "to_grid" (as described in
> > > > >> > > > met-6.0/data/config/README and the User's Guide), but
I find
> > > > >> providing
> > > > >> > a
> > > > >> > > > path to a sample file is pretty convenient.
> > > > >> > > >
> > > > >> > > > Or, even more convenient, if I'm running grid_stat
for
> example
> > > and
> > > > >> pass
> > > > >> > > in
> > > > >> > > > the GFS as the forecast and NAM as the observation, I
can
> just
> > > say
> > > > >> > > "to_grid
> > > > >> > > > = OBS;" and grid_stat will automatically regrid the
forecast
> > > data
> > > > >> (GFS)
> > > > >> > > to
> > > > >> > > > the observation grid (NAM).
> > > > >> > > >
> > > > >> > > > So, long story short, the following would work:
> > > > >> > > >
> > > > >> > > > regrid = {
> > > > >> > > > to_grid = "/path/to/sample/NDFD/data/file"; //
Or set
> > > this
> > > > to
> > > > >> > FCST
> > > > >> > > > or
> > > > >> > > > OBS if that's easier
> > > > >> > > > method = BUDGET;
> > > > >> > > > width = 2; // Budget should always have a
width
> of 2
> > > to
> > > > >> use
> > > > >> > > the 4
> > > > >> > > > closest points
> > > > >> > > > vld_thresh = 0.5;
> > > > >> > > > }
> > > > >> > > >
> > > > >> > > > Hope that helps get you going.
> > > > >> > > >
> > > > >> > > > Thanks,
> > > > >> > > > John
> > > > >> > > >
> > > > >> > > >
> > > > >> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb - NOAA
> Affiliate
> > > via
> > > > >> RT <
> > > > >> > > > met_help at ucar.edu> wrote:
> > > > >> > > >
> > > > >> > > > >
> > > > >> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted
upon.
> > > > >> > > > > Transaction: Ticket created by
daniel.plumb at noaa.gov
> > > > >> > > > > Queue: met_help
> > > > >> > > > > Subject: Forecast Regrid
> > > > >> > > > > Owner: Nobody
> > > > >> > > > > Requestors: daniel.plumb at noaa.gov
> > > > >> > > > > Status: new
> > > > >> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > >> > Ticket/Display.html?id=81369
> > > > >> > > >
> > > > >> > > > >
> > > > >> > > > >
> > > > >> > > > > Hello,
> > > > >> > > > >
> > > > >> > > > > I am a new user to the MET program, and attempting
to
> > regrid a
> > > > >> > 1051x813
> > > > >> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at
2.54 km
> > using
> > > > the
> > > > >> > > budget
> > > > >> > > > > regridding method. I am wondering if you have a
simple,
> > > > >> > straightforward
> > > > >> > > > > example of how to do this in MET. So far, I have
managed
> to
> > > > figure
> > > > >> > out
> > > > >> > > > this
> > > > >> > > > > is needed in the config file:
> > > > >> > > > >
> > > > >> > > > > regrid = {
> > > > >> > > > > to_grid = NONE;
> > > > >> > > > > method = BUDGET;
> > > > >> > > > > width = 1;
> > > > >> > > > > vld_thresh = 0.5;
> > > > >> > > > > }
> > > > >> > > > >
> > > > >> > > > > Any other help would be much appreciated.
> > > > >> > > > >
> > > > >> > > > > Thanks in advance!
> > > > >> > > > >
> > > > >> > > > > Dan
> > > > >> > > > >
> > > > >> > > > > --
> > > > >> > > > > Daniel W. Plumb, Jr.
> > > > >> > > > > Scientific Programmer/Meteorologist
> > > > >> > > > > Ace Info Solutions, Inc. (contractor)
> > > > >> > > > > National Weather Service (NWS)- Meteorological
Development
> > Lab
> > > > >> (MDL)
> > > > >> > > > >
> > > > >> > > > > Location:
> > > > >> > > > > SSMC2 – Room 10124
> > > > >> > > > > 1325 East West Highway
> > > > >> > > > > Silver Spring, MD 20910
> > > > >> > > > > Tel: (301) 427-9453
> > > > >> > > > >
> > > > >> > > > >
> > > > >> > > >
> > > > >> > > >
> > > > >> > >
> > > > >> > >
> > > > >> > > --
> > > > >> > > Daniel W. Plumb, Jr.
> > > > >> > > Scientific Programmer/Meteorologist
> > > > >> > > Ace Info Solutions, Inc. (contractor)
> > > > >> > > National Weather Service (NWS)- Meteorological
Development Lab
> > > (MDL)
> > > > >> > >
> > > > >> > > Location:
> > > > >> > > SSMC2 – Room 10124
> > > > >> > > 1325 East West Highway
> > > > >> > > Silver Spring, MD 20910
> > > > >> > > Tel: (301) 427-9453
> > > > >> > >
> > > > >> > >
> > > > >> >
> > > > >> >
> > > > >>
> > > > >>
> > > > >
> > > > >
> > > > > --
> > > > > Daniel W. Plumb, Jr.
> > > > > Scientific Programmer/Meteorologist
> > > > > Ace Info Solutions, Inc. (contractor)
> > > > > National Weather Service (NWS)- Meteorological Development
Lab
> (MDL)
> > > > >
> > > > > Location:
> > > > > SSMC2 – Room 10124
> > > > > 1325 East West Highway
> > > > > Silver Spring, MD 20910
> > > > > Tel: (301) 427-9453
> > > > >
> > > >
> > > >
> > > >
> > > > --
> > > > Daniel W. Plumb, Jr.
> > > > Scientific Programmer/Meteorologist
> > > > Ace Info Solutions, Inc. (contractor)
> > > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > > >
> > > > Location:
> > > > SSMC2 – Room 10124
> > > > 1325 East West Highway
> > > > Silver Spring, MD 20910
> > > > Tel: (301) 427-9453
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Daniel W. Plumb, Jr.
> > Scientific Programmer/Meteorologist
> > Ace Info Solutions, Inc. (contractor)
> > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> >
> > Location:
> > SSMC2 – Room 10124
> > 1325 East West Highway
> > Silver Spring, MD 20910
> > Tel: (301) 427-9453
> >
> >
>
>
--
Daniel W. Plumb, Jr.
Scientific Programmer/Meteorologist
Ace Info Solutions, Inc. (contractor)
National Weather Service (NWS)- Meteorological Development Lab (MDL)
Location:
SSMC2 – Room 10124
1325 East West Highway
Silver Spring, MD 20910
Tel: (301) 427-9453
------------------------------------------------
Subject: Forecast Regrid
From: Julie Prestopnik
Time: Wed Aug 02 11:14:01 2017
Hi Dan.
I see that you were working on "tide", however, I only have developer
access and can only access "gyre" at this time. On "gyre", I was able
to
install the updated met-6.0 code in my area on WCOSS.
In a new shell, please run:
module use /awc/noscrub/Julie.Prestopnik/modulefiles
module load NetCDF ics/15.0.3 met/6.0
Then, you should be able to access the update MET tools.
I was able to successfully run the same plot_data_plane command that
John
ran. Please give the updated MET code a try and let us know how it
goes.
Thanks,
Julie
On Mon, Jul 31, 2017 at 8:37 AM, Daniel Plumb - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
>
> Hi John,
>
> This image looks accurate! I look forward to the new code.
>
> Thanks in advance,
> Dan
>
> On Fri, Jul 28, 2017 at 6:37 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Dan,
> >
> > I worked on some code in the development version of MET to parse
polar
> > stereographic grids from CF-compliant NetCDF files, basically
adding in a
> > definition for this function:
> > NcCfFile::get_grid_mapping_polar_stereographic()
> >
> > The good news is that it *seemed* to work. I ran the following
> > plot_data_plane command:
> > plot_data_plane qpf_06hr_2017072512f006Z.nc qpf.ps
> > 'name="Forecast_total_precipitation"; level="(0,*,*)";'
> >
> > And I've attached the resulting image. The image certainly looks
> > reasonable.
> >
> > But I'd really like to step through this code with Randy Bullock,
the
> > author of the grid/projection-related code in MET. Once we review
the
> > changes next week, we can push them out to you.
> >
> > Since this is a pretty isolated change... so one option would be a
bugfix
> > for met-6.0. Technically this is a "new feature" rather than a
"bugfix",
> > but that'd be the quickest way to get you this functionality.
> >
> > Please take a look at this image and let me know what you think.
> >
> > Thanks,
> > John
> >
> >
> >
> > On Fri, Jul 28, 2017 at 3:40 PM, Daniel Plumb - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > >
> > > Hi Julie,
> > >
> > > That would be very helpful if I could assess MET v6.1. Is this
> currently
> > > installed on WCOSS?
> > >
> > > Thanks,
> > > Dan
> > >
> > > On Fri, Jul 28, 2017 at 5:34 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Dan.
> > > >
> > > > Thank you for sending your second email. We never received
your
> first
> > > one,
> > > > possibly because the attachments were too large.
> > > >
> > > > Unfortunately, MET does not currently support polar
stereographic
> grids
> > > for
> > > > CF compliant files, but it does support polar stereographic
grids for
> > > > GRIB1, GRIB2, and MET-NetCDF format files.
> > > >
> > > > We are aware that this functionality is desired and are
planning on
> > > adding
> > > > support for it in METv6.1. When we have a development version
of 6.1
> > > > ready, you could try it out if you'd like.
> > > >
> > > > I wish this were more helpful to you right now...
> > > >
> > > > Julie
> > > >
> > > >
> > > >
> > > > On Fri, Jul 28, 2017 at 2:05 PM, Daniel Plumb - NOAA Affiliate
via
> RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369 >
> > > > >
> > > > > Please let me know if you were able to get all of my
attachments.
> I'm
> > > not
> > > > > sure if they all went through.
> > > > >
> > > > > Thanks,
> > > > > Dan
> > > > >
> > > > > On Fri, Jul 28, 2017 at 4:01 PM, Daniel Plumb - NOAA
Affiliate <
> > > > > daniel.plumb at noaa.gov> wrote:
> > > > >
> > > > > > Hi Julie and John,
> > > > > >
> > > > > > Thank you for the help so far. I am able to run
regrid_data_plane
> > > > running
> > > > > > now. However, it is not giving me the desired output I
would like
> > > right
> > > > > > now. I have my command below and have attached the input
dataset
> > > > > > (qpf_06hr_2017072512f006Z.nc - *in netcdf format*) and the
> dataset
> > > with
> > > > > > the desired grid specs I want
(ndfdcogb2p517072700qpfd13_006
> > > > > > <http://ndfdcogb2p517072700qpfd13_006.nc> -* in grib2
format*):
> > > > > >
> > > > > > regrid_data_plane qpf_06hr_2017072512f006Z.nc
> > > > > > ndfdcogb2p517072700qpfd13_006 NEW_2145x1377.nc -field
> 'name="APCP";
> > > > > > level="A06";' -method BUDGET -width 2
> > > > > >
> > > > > > DEBUG 1: Reading data file: qpf_06hr_2017072512f006Z.nc
> > > > > > ERROR :
> > > > > > ERROR : NcCfFile::get_grid_mapping_polar_stereographic()
->
> Polar
> > > > > > stereographic grid not handled in MET.
> > > > > > ERROR :
> > > > > >
> > > > > > Is MET not able to handle polar stereographic grids? Or
maybe the
> > > > NetCDF
> > > > > > library is too old on WCOSS (dated: Dec 5 2012)?
> > > > > >
> > > > > > Thanks,
> > > > > > Dan
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > On Thu, Jul 27, 2017 at 5:16 PM, Julie Prestopnik via RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > >> I only have developer access to WCOSS, but I don't have
access
> to
> > > > WCOSS
> > > > > >> (gyre) at this time, as I believe it is still down for
> developers.
> > > > > >>
> > > > > >> Please try running "module purge" and then run:
> > > > > >>
> > > > > >> - module use /usrx/local/dev/modulefiles
> > > > > >> - module load met
> > > > > >>
> > > > > >> and see if that resolves the problem. I believe there is
a
> > conflict
> > > > > with
> > > > > >> the intel module you already have loaded. I know purging
isn't
> > the
> > > > best
> > > > > >> method to solve this problem, however, I do believe that
should
> > > work.
> > > > > >>
> > > > > >> Please let me know if you have any questions.
> > > > > >>
> > > > > >> Thanks,
> > > > > >>
> > > > > >> Julie
> > > > > >>
> > > > > >>
> > > > > >>
> > > > > >>
> > > > > >>
> > > > > >> On Thu, Jul 27, 2017 at 3:02 PM, John Halley Gotway via
RT <
> > > > > >> met_help at ucar.edu> wrote:
> > > > > >>
> > > > > >> >
> > > > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81369
> >
> > > > > >> >
> > > > > >> > Dan,
> > > > > >> >
> > > > > >> > I don't have access to WCOSS to replicate this problem,
but by
> > > > > colleague
> > > > > >> > Julie does.
> > > > > >> >
> > > > > >> > Julie, could you please take a look at this?
> > > > > >> >
> > > > > >> > Thanks,
> > > > > >> > John
> > > > > >> >
> > > > > >> >
> > > > > >> > On Thu, Jul 27, 2017 at 2:53 PM, Daniel Plumb - NOAA
Affiliate
> > via
> > > > RT
> > > > > <
> > > > > >> > met_help at ucar.edu> wrote:
> > > > > >> >
> > > > > >> > >
> > > > > >> > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81369
> > >
> > > > > >> > >
> > > > > >> > > Hi John,
> > > > > >> > >
> > > > > >> > > Thank you very much for the reply and examples.
> > > > > >> > >
> > > > > >> > > I am trying to replicate you example for
regrid_data_plane
> > > (which
> > > > is
> > > > > >> > > exactly the example I needed), but there seems to be
> something
> > > > wrong
> > > > > >> with
> > > > > >> > > MET on WCOSS tide. When I run "module load NetCDF
ics/15.0.3
> > > > udunits
> > > > > >> > ncview
> > > > > >> > > met/6.0 grib_util/v1.0.3", it is giving me an error
saying:
> > > > > >> > >
> > > > > >> > > "ics/15.0.3(16):ERROR:150: Module 'ics/15.0.3'
conflicts
> with
> > > the
> > > > > >> > currently
> > > > > >> > > loaded module(s) 'ics/12.1'
> > > > > >> > > ics/15.0.3(16):ERROR:102: Tcl command execution
failed:
> > conflict
> > > > > ics"
> > > > > >> > >
> > > > > >> > > I'm looking at the instructions under
> > > > > >> > >
http://www.dtcenter.org/met/users/support/online_tutorial/
> > > > > >> > ncwcp_201704.php
> > > > > >> > > .
> > > > > >> > >
> > > > > >> > > Thanks,
> > > > > >> > > Dan
> > > > > >> > >
> > > > > >> > > On Thu, Jul 27, 2017 at 2:50 PM, John Halley Gotway
via RT <
> > > > > >> > > met_help at ucar.edu> wrote:
> > > > > >> > >
> > > > > >> > > > Hi Dan,
> > > > > >> > > >
> > > > > >> > > > Tara gave me a heads up that you've started using
MET at
> > MDL.
> > > > > We're
> > > > > >> > > happy
> > > > > >> > > > to answer any questions you have or issues that
come up.
> > > > > >> > > >
> > > > > >> > > > I see that you'd like to use the regridding feature
within
> > the
> > > > MET
> > > > > >> > config
> > > > > >> > > > files. I'd actually suggest starting by running
the
> > > > > >> regrid_data_plane
> > > > > >> > > tool
> > > > > >> > > > first, and once you understand what it's doing,
move those
> > > > > settings
> > > > > >> > into
> > > > > >> > > > the config file. Both the regrid_data_plane tool
and the
> > > regrid
> > > > > >> > section
> > > > > >> > > of
> > > > > >> > > > the MET config files call the same underlying
library code
> > > > anyway.
> > > > > >> > > >
> > > > > >> > > > I was hoping to grab some sample MDL data to
demonstrate
> > some
> > > > > >> commands,
> > > > > >> > > but
> > > > > >> > > > all the links on this website aren't working:
> > > > > >> > > > http://www.nws.noaa.gov/mdl/hrqpf/download.php
> > > > > >> > > >
> > > > > >> > > > So instead of using HRAP and NDFD data, I'll use
sample
> GFS
> > > and
> > > > > NAM
> > > > > >> > data.
> > > > > >> > > >
> > > > > >> > > > # Get sample GFS and NAM files
> > > > > >> > > > wget
> > > > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.
> > > > > >> > > > 2017072700/gfs.t00z.pgrb2.0p25.f012
> > > > > >> > > > wget
> > > > > >> > > >
http://www.ftp.ncep.noaa.gov/data/nccf/com/nam/prod/nam.
> > > > > >> > > > 20170727/nam.t12z.conusnest.hiresf12.tm00.grib2
> > > > > >> > > >
> > > > > >> > > > # Run regrid_data_plane to regrid 6-hour precip
from GFS
> to
> > > NAM
> > > > > nest
> > > > > >> > > domain
> > > > > >> > > > met-6.0/bin/regrid_data_plane \
> > > > > >> > > > gfs.t00z.pgrb2.0p25.f012
nam.t12z.conusnest.hiresf12.
> > > > > tm00.grib2
> > > > > >> > > > gfs_APCP_06_nam_nest_domain.nc \
> > > > > >> > > > -field 'name="APCP"; level="A06";' -interp_mthd
BUDGET
> > > > > >> -interp_pnts
> > > > > >> > 2
> > > > > >> > > >
> > > > > >> > > > # Run plot_data_plane to plot the 6-hour precip
both
> > > > before/after
> > > > > >> the
> > > > > >> > > > interpolation
> > > > > >> > > > met-6.0/bin/plot_data_plane
gfs.t00z.pgrb2.0p25.f012
> > > > > >> > > > gfs.t00z.pgrb2.0p25.f012_APCP_06.ps 'name="APCP";
> > > > level="A06";' \
> > > > > >> > > > -color_table met-6.0/share/met/colortables/
> > > > wrf_precip2.ctable
> > > > > >> > > >
> > > > > >> > > > met-6.0/bin/plot_data_plane
gfs_APCP_06_nam_nest_domain.nc
> > > > > >> > > > gfs_APCP_06_nam_nest_domain.ps 'name="APCP_06";
> > > level="(*,*)";'
> > > > \
> > > > > >> > > > -color_table met-6.0/share/met/colortables/
> > > > wrf_precip2.ctable
> > > > > >> > > >
> > > > > >> > > > The resulting plots (converted to png) are attached
to
> this
> > > > email.
> > > > > >> > > >
> > > > > >> > > > In the regridding step, I defined the "to_grid" by
> providing
> > > the
> > > > > >> path
> > > > > >> > to
> > > > > >> > > a
> > > > > >> > > > file that's already on the grid I'd like to use.
MET
> reads
> > > that
> > > > > >> file,
> > > > > >> > > > grabs the grid definition, and uses that in the
regridding
> > > step.
> > > > > >> There
> > > > > >> > > are
> > > > > >> > > > multiple ways of defining "to_grid" (as described
in
> > > > > >> > > > met-6.0/data/config/README and the User's Guide),
but I
> find
> > > > > >> providing
> > > > > >> > a
> > > > > >> > > > path to a sample file is pretty convenient.
> > > > > >> > > >
> > > > > >> > > > Or, even more convenient, if I'm running grid_stat
for
> > example
> > > > and
> > > > > >> pass
> > > > > >> > > in
> > > > > >> > > > the GFS as the forecast and NAM as the observation,
I can
> > just
> > > > say
> > > > > >> > > "to_grid
> > > > > >> > > > = OBS;" and grid_stat will automatically regrid the
> forecast
> > > > data
> > > > > >> (GFS)
> > > > > >> > > to
> > > > > >> > > > the observation grid (NAM).
> > > > > >> > > >
> > > > > >> > > > So, long story short, the following would work:
> > > > > >> > > >
> > > > > >> > > > regrid = {
> > > > > >> > > > to_grid = "/path/to/sample/NDFD/data/file";
// Or
> set
> > > > this
> > > > > to
> > > > > >> > FCST
> > > > > >> > > > or
> > > > > >> > > > OBS if that's easier
> > > > > >> > > > method = BUDGET;
> > > > > >> > > > width = 2; // Budget should always have a
width
> > of 2
> > > > to
> > > > > >> use
> > > > > >> > > the 4
> > > > > >> > > > closest points
> > > > > >> > > > vld_thresh = 0.5;
> > > > > >> > > > }
> > > > > >> > > >
> > > > > >> > > > Hope that helps get you going.
> > > > > >> > > >
> > > > > >> > > > Thanks,
> > > > > >> > > > John
> > > > > >> > > >
> > > > > >> > > >
> > > > > >> > > > On Thu, Jul 27, 2017 at 9:03 AM, Daniel Plumb -
NOAA
> > Affiliate
> > > > via
> > > > > >> RT <
> > > > > >> > > > met_help at ucar.edu> wrote:
> > > > > >> > > >
> > > > > >> > > > >
> > > > > >> > > > > Thu Jul 27 09:03:13 2017: Request 81369 was acted
upon.
> > > > > >> > > > > Transaction: Ticket created by
daniel.plumb at noaa.gov
> > > > > >> > > > > Queue: met_help
> > > > > >> > > > > Subject: Forecast Regrid
> > > > > >> > > > > Owner: Nobody
> > > > > >> > > > > Requestors: daniel.plumb at noaa.gov
> > > > > >> > > > > Status: new
> > > > > >> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > >> > Ticket/Display.html?id=81369
> > > > > >> > > >
> > > > > >> > > > >
> > > > > >> > > > >
> > > > > >> > > > > Hello,
> > > > > >> > > > >
> > > > > >> > > > > I am a new user to the MET program, and
attempting to
> > > regrid a
> > > > > >> > 1051x813
> > > > > >> > > > > HRAP grid at 4.762km to a 2145x1377 NDFD grid at
2.54 km
> > > using
> > > > > the
> > > > > >> > > budget
> > > > > >> > > > > regridding method. I am wondering if you have a
simple,
> > > > > >> > straightforward
> > > > > >> > > > > example of how to do this in MET. So far, I have
managed
> > to
> > > > > figure
> > > > > >> > out
> > > > > >> > > > this
> > > > > >> > > > > is needed in the config file:
> > > > > >> > > > >
> > > > > >> > > > > regrid = {
> > > > > >> > > > > to_grid = NONE;
> > > > > >> > > > > method = BUDGET;
> > > > > >> > > > > width = 1;
> > > > > >> > > > > vld_thresh = 0.5;
> > > > > >> > > > > }
> > > > > >> > > > >
> > > > > >> > > > > Any other help would be much appreciated.
> > > > > >> > > > >
> > > > > >> > > > > Thanks in advance!
> > > > > >> > > > >
> > > > > >> > > > > Dan
> > > > > >> > > > >
> > > > > >> > > > > --
> > > > > >> > > > > Daniel W. Plumb, Jr.
> > > > > >> > > > > Scientific Programmer/Meteorologist
> > > > > >> > > > > Ace Info Solutions, Inc. (contractor)
> > > > > >> > > > > National Weather Service (NWS)- Meteorological
> Development
> > > Lab
> > > > > >> (MDL)
> > > > > >> > > > >
> > > > > >> > > > > Location:
> > > > > >> > > > > SSMC2 – Room 10124
> > > > > >> > > > > 1325 East West Highway
> > > > > >> > > > > Silver Spring, MD 20910
> > > > > >> > > > > Tel: (301) 427-9453
> > > > > >> > > > >
> > > > > >> > > > >
> > > > > >> > > >
> > > > > >> > > >
> > > > > >> > >
> > > > > >> > >
> > > > > >> > > --
> > > > > >> > > Daniel W. Plumb, Jr.
> > > > > >> > > Scientific Programmer/Meteorologist
> > > > > >> > > Ace Info Solutions, Inc. (contractor)
> > > > > >> > > National Weather Service (NWS)- Meteorological
Development
> Lab
> > > > (MDL)
> > > > > >> > >
> > > > > >> > > Location:
> > > > > >> > > SSMC2 – Room 10124
> > > > > >> > > 1325 East West Highway
> > > > > >> > > Silver Spring, MD 20910
> > > > > >> > > Tel: (301) 427-9453
> > > > > >> > >
> > > > > >> > >
> > > > > >> >
> > > > > >> >
> > > > > >>
> > > > > >>
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Daniel W. Plumb, Jr.
> > > > > > Scientific Programmer/Meteorologist
> > > > > > Ace Info Solutions, Inc. (contractor)
> > > > > > National Weather Service (NWS)- Meteorological Development
Lab
> > (MDL)
> > > > > >
> > > > > > Location:
> > > > > > SSMC2 – Room 10124
> > > > > > 1325 East West Highway
> > > > > > Silver Spring, MD 20910
> > > > > > Tel: (301) 427-9453
> > > > > >
> > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Daniel W. Plumb, Jr.
> > > > > Scientific Programmer/Meteorologist
> > > > > Ace Info Solutions, Inc. (contractor)
> > > > > National Weather Service (NWS)- Meteorological Development
Lab
> (MDL)
> > > > >
> > > > > Location:
> > > > > SSMC2 – Room 10124
> > > > > 1325 East West Highway
> > > > > Silver Spring, MD 20910
> > > > > Tel: (301) 427-9453
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > > --
> > > Daniel W. Plumb, Jr.
> > > Scientific Programmer/Meteorologist
> > > Ace Info Solutions, Inc. (contractor)
> > > National Weather Service (NWS)- Meteorological Development Lab
(MDL)
> > >
> > > Location:
> > > SSMC2 – Room 10124
> > > 1325 East West Highway
> > > Silver Spring, MD 20910
> > > Tel: (301) 427-9453
> > >
> > >
> >
> >
>
>
> --
> Daniel W. Plumb, Jr.
> Scientific Programmer/Meteorologist
> Ace Info Solutions, Inc. (contractor)
> National Weather Service (NWS)- Meteorological Development Lab (MDL)
>
> Location:
> SSMC2 – Room 10124
> 1325 East West Highway
> Silver Spring, MD 20910
> Tel: (301) 427-9453
>
>
------------------------------------------------
More information about the Met_help
mailing list