[Met_help] [rt.rap.ucar.edu #73975] History for Regridding problem
John Halley Gotway via RT
met_help at ucar.edu
Fri Nov 6 17:24:15 MST 2015
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I'm a new MET user, and I'm attempting to regrid RTMA data to an NMMB output grid. I used the following tool:
./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2 nmbprs_d01.12 regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN -width 2 -v 1
which worked. But, when I subtracted the fields from one dataset to the other, there was a clear offset on the values. So, the regridding was done indeed, but the gridpoint values did not seem to be assigned correctly with respect to their counterparts.
I attached the plot, so you can see the offset problem. Any help will be greatly appreciated.
Best Regards,
Jose
Jose M. Garcia-Rivera, PhD
NWP Scientist
I.M. Systems Group
This e-mail, including attachments, may include confidential and/or proprietary information, and may be used only by the person or entity to which it is addressed. If the reader of this e-mail is not the intended recipient or his or her authorized agent, the reader is hereby notified that any dissemination, distribution or copying of this e-mail is prohibited. If you have received this e-mail in error, please notify the sender by replying to this message and delete this e-mail immediately.
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Tue Nov 03 14:03:55 2015
Jose,
Thanks for letting us know about the issue you're seeing. I'd like to
replicate the behavior your describe using the exact data you used.
Once I
can replicate it, I'll be able to investigate further.
Are you able to send me the following 2 files you used:
rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
You can post them to our anonymous ftp site, following these
instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Please write back to let me know when you've posted the data and I'll
go
grab it.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> Transaction: Ticket created by GarciaJ at imsg.com
> Queue: met_help
> Subject: Regridding problem
> Owner: Nobody
> Requestors: GarciaJ at imsg.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
>
> Hello,
>
> I'm a new MET user, and I'm attempting to regrid RTMA data to an
NMMB
> output grid. I used the following tool:
>
>
>
> ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2 nmbprs_d01.12
> regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN -width 2
-v 1
>
> which worked. But, when I subtracted the fields from one dataset to
the
> other, there was a clear offset on the values. So, the regridding
was done
> indeed, but the gridpoint values did not seem to be assigned
correctly with
> respect to their counterparts.
>
> I attached the plot, so you can see the offset problem. Any help
will be
> greatly appreciated.
>
> Best Regards,
> Jose
>
> Jose M. Garcia-Rivera, PhD
> NWP Scientist
> I.M. Systems Group
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Tue Nov 03 14:16:05 2015
Jose,
I wonder about the files you're using. I grabbed this RTMA file which
matches the naming convention you used:
ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/rtma/prod/rtma2p5.20151102/rtma2p5.t23z.2dvaranl_nwrfc.grb2
But the "nwrfc" means that it's defined over the northwest region of
the
United States, as shown in the attached plot.
Perhaps you actually used the full CONUS "ndfd" file (image also
attached)?
ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/rtma/prod/rtma2p5.20151102/rtma2p5.t23z.2dvaranl_ndfd.grb2
Thanks,
John
On Tue, Nov 3, 2015 at 2:03 PM, The RT System itself via RT <
met_help at ucar.edu> wrote:
>
> Tue Nov 03 14:03:55 2015: Request 73975 was acted upon.
> Transaction: Given to johnhg (John Halley Gotway) by RT_System
> Queue: met_help
> Subject: Regridding problem
> Owner: johnhg
> Requestors: GarciaJ at imsg.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
>
> This transaction appears to have no content
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Tue Nov 03 15:07:47 2015
Dear John,
Thank you so much for the quick reply.
I followed the instructions, and ftp'ed the files to "garcia-
rivera_data".
The nmbprs file is the .grb NMMB output file.
The rtma file is the observational.
The nmm_T2_bias.ncl subtracts the fields from both outputs in the same
grid, and produces the plot.
Any help will be greatly appreciated.
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Tuesday, November 03, 2015 4:04 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
Thanks for letting us know about the issue you're seeing. I'd like to
replicate the behavior your describe using the exact data you used.
Once I can replicate it, I'll be able to investigate further.
Are you able to send me the following 2 files you used:
rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
You can post them to our anonymous ftp site, following these
instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Please write back to let me know when you've posted the data and I'll
go grab it.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> Transaction: Ticket created by GarciaJ at imsg.com
> Queue: met_help
> Subject: Regridding problem
> Owner: Nobody
> Requestors: GarciaJ at imsg.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >
>
>
> Hello,
>
> I'm a new MET user, and I'm attempting to regrid RTMA data to an
NMMB
> output grid. I used the following tool:
>
>
>
> ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2 nmbprs_d01.12
> regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN -width 2
> -v 1
>
> which worked. But, when I subtracted the fields from one dataset to
> the other, there was a clear offset on the values. So, the
regridding
> was done indeed, but the gridpoint values did not seem to be
assigned
> correctly with respect to their counterparts.
>
> I attached the plot, so you can see the offset problem. Any help
will
> be greatly appreciated.
>
> Best Regards,
> Jose
>
> Jose M. Garcia-Rivera, PhD
> NWP Scientist
> I.M. Systems Group
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Tue Nov 03 15:29:07 2015
Jose,
I believe the problem is in our handling the domain on which your
model
data resides. I ran the following commands to plot temperature from
your
model data file:
met-5.1/bin/plot_data_plane \
nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
'name="TMP"; level="Z0";'
convert -rotate 90 -background white -flatten \
nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
The resulting image is attached and shows the shift you described.
So it isn't the regridding step that's the problem... its our handling
of
that domain. I plotted the same data using IDV and it puts things in
the
correct place (also attached).
One issue is that the 2-m temperature data in GRIB record 296 contains
a
constant value of 273. The surface temperature data in GRIB record
276
contains more reasonable values. So that's the field I'm
investigating.
But you might want to check why GRIB record 296 looks so odd.
I need to head out of the office, but will investigate further
tomorrow.
Thanks,
John
On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thank you so much for the quick reply.
>
> I followed the instructions, and ftp'ed the files to "garcia-
rivera_data".
>
> The nmbprs file is the .grb NMMB output file.
> The rtma file is the observational.
> The nmm_T2_bias.ncl subtracts the fields from both outputs in the
same
> grid, and produces the plot.
>
> Any help will be greatly appreciated.
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Tuesday, November 03, 2015 4:04 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Thanks for letting us know about the issue you're seeing. I'd like
to
> replicate the behavior your describe using the exact data you used.
Once I
> can replicate it, I'll be able to investigate further.
>
> Are you able to send me the following 2 files you used:
> rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>
> You can post them to our anonymous ftp site, following these
instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please write back to let me know when you've posted the data and
I'll go
> grab it.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
>
> On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > Transaction: Ticket created by GarciaJ at imsg.com
> > Queue: met_help
> > Subject: Regridding problem
> > Owner: Nobody
> > Requestors: GarciaJ at imsg.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >
> >
> >
> > Hello,
> >
> > I'm a new MET user, and I'm attempting to regrid RTMA data to an
NMMB
> > output grid. I used the following tool:
> >
> >
> >
> > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2 nmbprs_d01.12
> > regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN -width
2
> > -v 1
> >
> > which worked. But, when I subtracted the fields from one dataset
to
> > the other, there was a clear offset on the values. So, the
regridding
> > was done indeed, but the gridpoint values did not seem to be
assigned
> > correctly with respect to their counterparts.
> >
> > I attached the plot, so you can see the offset problem. Any help
will
> > be greatly appreciated.
> >
> > Best Regards,
> > Jose
> >
> > Jose M. Garcia-Rivera, PhD
> > NWP Scientist
> > I.M. Systems Group
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: Regridding problem
From: Jose Garcia-Rivera
Time: Wed Nov 04 08:48:05 2015
Hi John,
Thanks for all the help. I want to add this pdf, showing something odd
with the regridding. The domains did not match perfectly on the north
side, it's like a skewing effect. The south sides match perfectly, but
note the differences in the regridded RTMA file and the model file in
Michigan, for example. Maybe that has something to do with it.
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Tuesday, November 03, 2015 5:29 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
I believe the problem is in our handling the domain on which your
model data resides. I ran the following commands to plot temperature
from your model data file:
met-5.1/bin/plot_data_plane \
nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
'name="TMP"; level="Z0";'
convert -rotate 90 -background white -flatten \
nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
The resulting image is attached and shows the shift you described.
So it isn't the regridding step that's the problem... its our handling
of that domain. I plotted the same data using IDV and it puts things
in the correct place (also attached).
One issue is that the 2-m temperature data in GRIB record 296 contains
a constant value of 273. The surface temperature data in GRIB record
276 contains more reasonable values. So that's the field I'm
investigating.
But you might want to check why GRIB record 296 looks so odd.
I need to head out of the office, but will investigate further
tomorrow.
Thanks,
John
On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thank you so much for the quick reply.
>
> I followed the instructions, and ftp'ed the files to "garcia-
rivera_data".
>
> The nmbprs file is the .grb NMMB output file.
> The rtma file is the observational.
> The nmm_T2_bias.ncl subtracts the fields from both outputs in the
same
> grid, and produces the plot.
>
> Any help will be greatly appreciated.
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Tuesday, November 03, 2015 4:04 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Thanks for letting us know about the issue you're seeing. I'd like
to
> replicate the behavior your describe using the exact data you used.
> Once I can replicate it, I'll be able to investigate further.
>
> Are you able to send me the following 2 files you used:
> rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>
> You can post them to our anonymous ftp site, following these
instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please write back to let me know when you've posted the data and
I'll
> go grab it.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
>
> On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > Transaction: Ticket created by GarciaJ at imsg.com
> > Queue: met_help
> > Subject: Regridding problem
> > Owner: Nobody
> > Requestors: GarciaJ at imsg.com
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >
> >
> >
> > Hello,
> >
> > I'm a new MET user, and I'm attempting to regrid RTMA data to an
> > NMMB output grid. I used the following tool:
> >
> >
> >
> > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2 nmbprs_d01.12
> > regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN -width
2
> > -v 1
> >
> > which worked. But, when I subtracted the fields from one dataset
to
> > the other, there was a clear offset on the values. So, the
> > regridding was done indeed, but the gridpoint values did not seem
to
> > be assigned correctly with respect to their counterparts.
> >
> > I attached the plot, so you can see the offset problem. Any help
> > will be greatly appreciated.
> >
> > Best Regards,
> > Jose
> >
> > Jose M. Garcia-Rivera, PhD
> > NWP Scientist
> > I.M. Systems Group
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Wed Nov 04 10:28:43 2015
Jose,
I believe I've figured out the problem. The grid code in MET
currently
only supports Lambert Conformal grids where the grid spacing is the
same in
the x and y directions.
Your grid is spaced 3.317km in the x-direction and 2.97km in the
y-direction. The MET code is reading the Dx value and using it for
both x
and y. This causes the skewing of the data to the north that you're
seeing.
Really, there are 2 problems here.
First, if we're requiring that Dx and Dy be the same for Lambert
Conformal
grids, MET should error out with a descriptive error message when they
are
not.
Second, obviously it would be preferable for MET to support Lambert
Conformal grids where Dx and Dy are not the same! I'll ask around
here for
more information on that.
Thanks for paying close attention to your data and uncovering this
shortcoming in MET.
John
On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Thanks for all the help. I want to add this pdf, showing something
odd
> with the regridding. The domains did not match perfectly on the
north side,
> it's like a skewing effect. The south sides match perfectly, but
note the
> differences in the regridded RTMA file and the model file in
Michigan, for
> example. Maybe that has something to do with it.
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Tuesday, November 03, 2015 5:29 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> I believe the problem is in our handling the domain on which your
model
> data resides. I ran the following commands to plot temperature from
your
> model data file:
>
> met-5.1/bin/plot_data_plane \
> nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> 'name="TMP"; level="Z0";'
>
> convert -rotate 90 -background white -flatten \
> nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>
> The resulting image is attached and shows the shift you described.
>
> So it isn't the regridding step that's the problem... its our
handling of
> that domain. I plotted the same data using IDV and it puts things
in the
> correct place (also attached).
>
> One issue is that the 2-m temperature data in GRIB record 296
contains a
> constant value of 273. The surface temperature data in GRIB record
276
> contains more reasonable values. So that's the field I'm
investigating.
> But you might want to check why GRIB record 296 looks so odd.
>
> I need to head out of the office, but will investigate further
tomorrow.
>
> Thanks,
> John
>
>
>
> On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > Thank you so much for the quick reply.
> >
> > I followed the instructions, and ftp'ed the files to
> "garcia-rivera_data".
> >
> > The nmbprs file is the .grb NMMB output file.
> > The rtma file is the observational.
> > The nmm_T2_bias.ncl subtracts the fields from both outputs in the
same
> > grid, and produces the plot.
> >
> > Any help will be greatly appreciated.
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Tuesday, November 03, 2015 4:04 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Thanks for letting us know about the issue you're seeing. I'd
like to
> > replicate the behavior your describe using the exact data you
used.
> > Once I can replicate it, I'll be able to investigate further.
> >
> > Are you able to send me the following 2 files you used:
> > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> >
> > You can post them to our anonymous ftp site, following these
> instructions:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Please write back to let me know when you've posted the data and
I'll
> > go grab it.
> >
> > Thanks,
> > John Halley Gotway
> > met_help at ucar.edu
> >
> >
> >
> > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > > Transaction: Ticket created by GarciaJ at imsg.com
> > > Queue: met_help
> > > Subject: Regridding problem
> > > Owner: Nobody
> > > Requestors: GarciaJ at imsg.com
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >
> > >
> > >
> > > Hello,
> > >
> > > I'm a new MET user, and I'm attempting to regrid RTMA data to an
> > > NMMB output grid. I used the following tool:
> > >
> > >
> > >
> > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
nmbprs_d01.12
> > > regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN
-width 2
> > > -v 1
> > >
> > > which worked. But, when I subtracted the fields from one dataset
to
> > > the other, there was a clear offset on the values. So, the
> > > regridding was done indeed, but the gridpoint values did not
seem to
> > > be assigned correctly with respect to their counterparts.
> > >
> > > I attached the plot, so you can see the offset problem. Any help
> > > will be greatly appreciated.
> > >
> > > Best Regards,
> > > Jose
> > >
> > > Jose M. Garcia-Rivera, PhD
> > > NWP Scientist
> > > I.M. Systems Group
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Wed Nov 04 10:34:26 2015
Hi John,
Thanks so much for the information. So, that means that I'll have the
same problem when running any other analyses with this particular grid
using MET?
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, November 04, 2015 12:29 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
I believe I've figured out the problem. The grid code in MET
currently only supports Lambert Conformal grids where the grid spacing
is the same in the x and y directions.
Your grid is spaced 3.317km in the x-direction and 2.97km in the y-
direction. The MET code is reading the Dx value and using it for both
x and y. This causes the skewing of the data to the north that you're
seeing.
Really, there are 2 problems here.
First, if we're requiring that Dx and Dy be the same for Lambert
Conformal grids, MET should error out with a descriptive error message
when they are not.
Second, obviously it would be preferable for MET to support Lambert
Conformal grids where Dx and Dy are not the same! I'll ask around
here for more information on that.
Thanks for paying close attention to your data and uncovering this
shortcoming in MET.
John
On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Thanks for all the help. I want to add this pdf, showing something
odd
> with the regridding. The domains did not match perfectly on the
north
> side, it's like a skewing effect. The south sides match perfectly,
but
> note the differences in the regridded RTMA file and the model file
in
> Michigan, for example. Maybe that has something to do with it.
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Tuesday, November 03, 2015 5:29 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> I believe the problem is in our handling the domain on which your
> model data resides. I ran the following commands to plot
temperature
> from your model data file:
>
> met-5.1/bin/plot_data_plane \
> nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> 'name="TMP"; level="Z0";'
>
> convert -rotate 90 -background white -flatten \
> nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>
> The resulting image is attached and shows the shift you described.
>
> So it isn't the regridding step that's the problem... its our
handling
> of that domain. I plotted the same data using IDV and it puts
things
> in the correct place (also attached).
>
> One issue is that the 2-m temperature data in GRIB record 296
contains
> a constant value of 273. The surface temperature data in GRIB
record
> 276 contains more reasonable values. So that's the field I'm
investigating.
> But you might want to check why GRIB record 296 looks so odd.
>
> I need to head out of the office, but will investigate further
tomorrow.
>
> Thanks,
> John
>
>
>
> On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > Thank you so much for the quick reply.
> >
> > I followed the instructions, and ftp'ed the files to
> "garcia-rivera_data".
> >
> > The nmbprs file is the .grb NMMB output file.
> > The rtma file is the observational.
> > The nmm_T2_bias.ncl subtracts the fields from both outputs in the
> > same grid, and produces the plot.
> >
> > Any help will be greatly appreciated.
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Tuesday, November 03, 2015 4:04 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Thanks for letting us know about the issue you're seeing. I'd
like
> > to replicate the behavior your describe using the exact data you
used.
> > Once I can replicate it, I'll be able to investigate further.
> >
> > Are you able to send me the following 2 files you used:
> > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> >
> > You can post them to our anonymous ftp site, following these
> instructions:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Please write back to let me know when you've posted the data and
> > I'll go grab it.
> >
> > Thanks,
> > John Halley Gotway
> > met_help at ucar.edu
> >
> >
> >
> > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > > Transaction: Ticket created by GarciaJ at imsg.com
> > > Queue: met_help
> > > Subject: Regridding problem
> > > Owner: Nobody
> > > Requestors: GarciaJ at imsg.com
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >
> > >
> > >
> > > Hello,
> > >
> > > I'm a new MET user, and I'm attempting to regrid RTMA data to an
> > > NMMB output grid. I used the following tool:
> > >
> > >
> > >
> > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
nmbprs_d01.12
> > > regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN
-width
> > > 2 -v 1
> > >
> > > which worked. But, when I subtracted the fields from one dataset
> > > to the other, there was a clear offset on the values. So, the
> > > regridding was done indeed, but the gridpoint values did not
seem
> > > to be assigned correctly with respect to their counterparts.
> > >
> > > I attached the plot, so you can see the offset problem. Any help
> > > will be greatly appreciated.
> > >
> > > Best Regards,
> > > Jose
> > >
> > > Jose M. Garcia-Rivera, PhD
> > > NWP Scientist
> > > I.M. Systems Group
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Wed Nov 04 11:01:35 2015
Jose,
Yes, that's correct. Not supporting differing Lambert Conformal Dx
and Dy
values is a short-coming in MET. However, we're investigating the
impact/cost of adding support for it. We may be able to add support
by
patching the build.
I'll write you back once we figure out how to proceed.
Thanks,
John
On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Thanks so much for the information. So, that means that I'll have
the same
> problem when running any other analyses with this particular grid
using MET?
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, November 04, 2015 12:29 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> I believe I've figured out the problem. The grid code in MET
currently
> only supports Lambert Conformal grids where the grid spacing is the
same in
> the x and y directions.
>
> Your grid is spaced 3.317km in the x-direction and 2.97km in the
> y-direction. The MET code is reading the Dx value and using it for
both x
> and y. This causes the skewing of the data to the north that you're
seeing.
>
> Really, there are 2 problems here.
>
> First, if we're requiring that Dx and Dy be the same for Lambert
Conformal
> grids, MET should error out with a descriptive error message when
they are
> not.
>
> Second, obviously it would be preferable for MET to support Lambert
> Conformal grids where Dx and Dy are not the same! I'll ask around
here for
> more information on that.
>
> Thanks for paying close attention to your data and uncovering this
> shortcoming in MET.
>
> John
>
> On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Hi John,
> >
> > Thanks for all the help. I want to add this pdf, showing something
odd
> > with the regridding. The domains did not match perfectly on the
north
> > side, it's like a skewing effect. The south sides match perfectly,
but
> > note the differences in the regridded RTMA file and the model file
in
> > Michigan, for example. Maybe that has something to do with it.
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Tuesday, November 03, 2015 5:29 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > I believe the problem is in our handling the domain on which your
> > model data resides. I ran the following commands to plot
temperature
> > from your model data file:
> >
> > met-5.1/bin/plot_data_plane \
> > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > 'name="TMP"; level="Z0";'
> >
> > convert -rotate 90 -background white -flatten \
> > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> >
> > The resulting image is attached and shows the shift you described.
> >
> > So it isn't the regridding step that's the problem... its our
handling
> > of that domain. I plotted the same data using IDV and it puts
things
> > in the correct place (also attached).
> >
> > One issue is that the 2-m temperature data in GRIB record 296
contains
> > a constant value of 273. The surface temperature data in GRIB
record
> > 276 contains more reasonable values. So that's the field I'm
> investigating.
> > But you might want to check why GRIB record 296 looks so odd.
> >
> > I need to head out of the office, but will investigate further
tomorrow.
> >
> > Thanks,
> > John
> >
> >
> >
> > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Dear John,
> > >
> > > Thank you so much for the quick reply.
> > >
> > > I followed the instructions, and ftp'ed the files to
> > "garcia-rivera_data".
> > >
> > > The nmbprs file is the .grb NMMB output file.
> > > The rtma file is the observational.
> > > The nmm_T2_bias.ncl subtracts the fields from both outputs in
the
> > > same grid, and produces the plot.
> > >
> > > Any help will be greatly appreciated.
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > Thanks for letting us know about the issue you're seeing. I'd
like
> > > to replicate the behavior your describe using the exact data you
used.
> > > Once I can replicate it, I'll be able to investigate further.
> > >
> > > Are you able to send me the following 2 files you used:
> > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > >
> > > You can post them to our anonymous ftp site, following these
> > instructions:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Please write back to let me know when you've posted the data and
> > > I'll go grab it.
> > >
> > > Thanks,
> > > John Halley Gotway
> > > met_help at ucar.edu
> > >
> > >
> > >
> > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > Queue: met_help
> > > > Subject: Regridding problem
> > > > Owner: Nobody
> > > > Requestors: GarciaJ at imsg.com
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > >
> > > >
> > > >
> > > > Hello,
> > > >
> > > > I'm a new MET user, and I'm attempting to regrid RTMA data to
an
> > > > NMMB output grid. I used the following tool:
> > > >
> > > >
> > > >
> > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
nmbprs_d01.12
> > > > regridded.nc -field 'name="TMP"; level="Z2";' -method BILIN
-width
> > > > 2 -v 1
> > > >
> > > > which worked. But, when I subtracted the fields from one
dataset
> > > > to the other, there was a clear offset on the values. So, the
> > > > regridding was done indeed, but the gridpoint values did not
seem
> > > > to be assigned correctly with respect to their counterparts.
> > > >
> > > > I attached the plot, so you can see the offset problem. Any
help
> > > > will be greatly appreciated.
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > > Jose M. Garcia-Rivera, PhD
> > > > NWP Scientist
> > > > I.M. Systems Group
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Wed Nov 04 11:44:09 2015
Dear John,
Thank you for your insights. One last question: If we run the model
with dx = dy, should that solve the problem?
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, November 04, 2015 1:02 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
Yes, that's correct. Not supporting differing Lambert Conformal Dx
and Dy values is a short-coming in MET. However, we're investigating
the impact/cost of adding support for it. We may be able to add
support by patching the build.
I'll write you back once we figure out how to proceed.
Thanks,
John
On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Thanks so much for the information. So, that means that I'll have
the
> same problem when running any other analyses with this particular
grid using MET?
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, November 04, 2015 12:29 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> I believe I've figured out the problem. The grid code in MET
> currently only supports Lambert Conformal grids where the grid
spacing
> is the same in the x and y directions.
>
> Your grid is spaced 3.317km in the x-direction and 2.97km in the
> y-direction. The MET code is reading the Dx value and using it for
> both x and y. This causes the skewing of the data to the north that
you're seeing.
>
> Really, there are 2 problems here.
>
> First, if we're requiring that Dx and Dy be the same for Lambert
> Conformal grids, MET should error out with a descriptive error
message
> when they are not.
>
> Second, obviously it would be preferable for MET to support Lambert
> Conformal grids where Dx and Dy are not the same! I'll ask around
> here for more information on that.
>
> Thanks for paying close attention to your data and uncovering this
> shortcoming in MET.
>
> John
>
> On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Hi John,
> >
> > Thanks for all the help. I want to add this pdf, showing something
> > odd with the regridding. The domains did not match perfectly on
the
> > north side, it's like a skewing effect. The south sides match
> > perfectly, but note the differences in the regridded RTMA file and
> > the model file in Michigan, for example. Maybe that has something
to do with it.
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Tuesday, November 03, 2015 5:29 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > I believe the problem is in our handling the domain on which your
> > model data resides. I ran the following commands to plot
> > temperature from your model data file:
> >
> > met-5.1/bin/plot_data_plane \
> > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > 'name="TMP"; level="Z0";'
> >
> > convert -rotate 90 -background white -flatten \
> > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> >
> > The resulting image is attached and shows the shift you described.
> >
> > So it isn't the regridding step that's the problem... its our
> > handling of that domain. I plotted the same data using IDV and it
> > puts things in the correct place (also attached).
> >
> > One issue is that the 2-m temperature data in GRIB record 296
> > contains a constant value of 273. The surface temperature data in
> > GRIB record
> > 276 contains more reasonable values. So that's the field I'm
> investigating.
> > But you might want to check why GRIB record 296 looks so odd.
> >
> > I need to head out of the office, but will investigate further
tomorrow.
> >
> > Thanks,
> > John
> >
> >
> >
> > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Dear John,
> > >
> > > Thank you so much for the quick reply.
> > >
> > > I followed the instructions, and ftp'ed the files to
> > "garcia-rivera_data".
> > >
> > > The nmbprs file is the .grb NMMB output file.
> > > The rtma file is the observational.
> > > The nmm_T2_bias.ncl subtracts the fields from both outputs in
the
> > > same grid, and produces the plot.
> > >
> > > Any help will be greatly appreciated.
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > Thanks for letting us know about the issue you're seeing. I'd
> > > like to replicate the behavior your describe using the exact
data you used.
> > > Once I can replicate it, I'll be able to investigate further.
> > >
> > > Are you able to send me the following 2 files you used:
> > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > >
> > > You can post them to our anonymous ftp site, following these
> > instructions:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Please write back to let me know when you've posted the data and
> > > I'll go grab it.
> > >
> > > Thanks,
> > > John Halley Gotway
> > > met_help at ucar.edu
> > >
> > >
> > >
> > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > Queue: met_help
> > > > Subject: Regridding problem
> > > > Owner: Nobody
> > > > Requestors: GarciaJ at imsg.com
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > >
> > > >
> > > >
> > > > Hello,
> > > >
> > > > I'm a new MET user, and I'm attempting to regrid RTMA data to
an
> > > > NMMB output grid. I used the following tool:
> > > >
> > > >
> > > >
> > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > > > nmbprs_d01.12 regridded.nc -field 'name="TMP"; level="Z2";'
> > > > -method BILIN -width
> > > > 2 -v 1
> > > >
> > > > which worked. But, when I subtracted the fields from one
dataset
> > > > to the other, there was a clear offset on the values. So, the
> > > > regridding was done indeed, but the gridpoint values did not
> > > > seem to be assigned correctly with respect to their
counterparts.
> > > >
> > > > I attached the plot, so you can see the offset problem. Any
help
> > > > will be greatly appreciated.
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > > Jose M. Garcia-Rivera, PhD
> > > > NWP Scientist
> > > > I.M. Systems Group
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
> > > > person or entity to which it is addressed. If the reader of
this
> > > > e-mail is not the intended recipient or his or her authorized
> > > > agent, the reader is hereby notified that any dissemination,
> > > > distribution or copying of this e-mail is prohibited. If you
> > > > have received this e-mail in error, please notify the sender
by
> > > > replying to this message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Wed Nov 04 11:51:11 2015
Jose, yes, running with dx = dy would solve the problem.
John
On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thank you for your insights. One last question: If we run the model
with
> dx = dy, should that solve the problem?
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, November 04, 2015 1:02 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Yes, that's correct. Not supporting differing Lambert Conformal Dx
and Dy
> values is a short-coming in MET. However, we're investigating the
> impact/cost of adding support for it. We may be able to add support
by
> patching the build.
>
> I'll write you back once we figure out how to proceed.
>
> Thanks,
> John
>
> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Hi John,
> >
> > Thanks so much for the information. So, that means that I'll have
the
> > same problem when running any other analyses with this particular
grid
> using MET?
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Wednesday, November 04, 2015 12:29 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > I believe I've figured out the problem. The grid code in MET
> > currently only supports Lambert Conformal grids where the grid
spacing
> > is the same in the x and y directions.
> >
> > Your grid is spaced 3.317km in the x-direction and 2.97km in the
> > y-direction. The MET code is reading the Dx value and using it
for
> > both x and y. This causes the skewing of the data to the north
that
> you're seeing.
> >
> > Really, there are 2 problems here.
> >
> > First, if we're requiring that Dx and Dy be the same for Lambert
> > Conformal grids, MET should error out with a descriptive error
message
> > when they are not.
> >
> > Second, obviously it would be preferable for MET to support
Lambert
> > Conformal grids where Dx and Dy are not the same! I'll ask around
> > here for more information on that.
> >
> > Thanks for paying close attention to your data and uncovering this
> > shortcoming in MET.
> >
> > John
> >
> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Hi John,
> > >
> > > Thanks for all the help. I want to add this pdf, showing
something
> > > odd with the regridding. The domains did not match perfectly on
the
> > > north side, it's like a skewing effect. The south sides match
> > > perfectly, but note the differences in the regridded RTMA file
and
> > > the model file in Michigan, for example. Maybe that has
something to
> do with it.
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > I believe the problem is in our handling the domain on which
your
> > > model data resides. I ran the following commands to plot
> > > temperature from your model data file:
> > >
> > > met-5.1/bin/plot_data_plane \
> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > > 'name="TMP"; level="Z0";'
> > >
> > > convert -rotate 90 -background white -flatten \
> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > >
> > > The resulting image is attached and shows the shift you
described.
> > >
> > > So it isn't the regridding step that's the problem... its our
> > > handling of that domain. I plotted the same data using IDV and
it
> > > puts things in the correct place (also attached).
> > >
> > > One issue is that the 2-m temperature data in GRIB record 296
> > > contains a constant value of 273. The surface temperature data
in
> > > GRIB record
> > > 276 contains more reasonable values. So that's the field I'm
> > investigating.
> > > But you might want to check why GRIB record 296 looks so odd.
> > >
> > > I need to head out of the office, but will investigate further
> tomorrow.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > >
> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> > > >
> > > > Dear John,
> > > >
> > > > Thank you so much for the quick reply.
> > > >
> > > > I followed the instructions, and ftp'ed the files to
> > > "garcia-rivera_data".
> > > >
> > > > The nmbprs file is the .grb NMMB output file.
> > > > The rtma file is the observational.
> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs in
the
> > > > same grid, and produces the plot.
> > > >
> > > > Any help will be greatly appreciated.
> > > >
> > > > Best Regards,
> > > > Jose
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >
> > > > Jose,
> > > >
> > > > Thanks for letting us know about the issue you're seeing. I'd
> > > > like to replicate the behavior your describe using the exact
data
> you used.
> > > > Once I can replicate it, I'll be able to investigate further.
> > > >
> > > > Are you able to send me the following 2 files you used:
> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > > >
> > > > You can post them to our anonymous ftp site, following these
> > > instructions:
> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Please write back to let me know when you've posted the data
and
> > > > I'll go grab it.
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > > met_help at ucar.edu
> > > >
> > > >
> > > >
> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > > Queue: met_help
> > > > > Subject: Regridding problem
> > > > > Owner: Nobody
> > > > > Requestors: GarciaJ at imsg.com
> > > > > Status: new
> > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > > >
> > > > >
> > > > >
> > > > > Hello,
> > > > >
> > > > > I'm a new MET user, and I'm attempting to regrid RTMA data
to an
> > > > > NMMB output grid. I used the following tool:
> > > > >
> > > > >
> > > > >
> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP"; level="Z2";'
> > > > > -method BILIN -width
> > > > > 2 -v 1
> > > > >
> > > > > which worked. But, when I subtracted the fields from one
dataset
> > > > > to the other, there was a clear offset on the values. So,
the
> > > > > regridding was done indeed, but the gridpoint values did not
> > > > > seem to be assigned correctly with respect to their
counterparts.
> > > > >
> > > > > I attached the plot, so you can see the offset problem. Any
help
> > > > > will be greatly appreciated.
> > > > >
> > > > > Best Regards,
> > > > > Jose
> > > > >
> > > > > Jose M. Garcia-Rivera, PhD
> > > > > NWP Scientist
> > > > > I.M. Systems Group
> > > > > This e-mail, including attachments, may include confidential
> > > > > and/or proprietary information, and may be used only by the
> > > > > person or entity to which it is addressed. If the reader of
this
> > > > > e-mail is not the intended recipient or his or her
authorized
> > > > > agent, the reader is hereby notified that any dissemination,
> > > > > distribution or copying of this e-mail is prohibited. If you
> > > > > have received this e-mail in error, please notify the sender
by
> > > > > replying to this message and delete this
> > > > e-mail immediately.
> > > > >
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Wed Nov 04 14:28:09 2015
Jose,
We're still trying to decide how to address this issue. We definitely
plan
on adding support for grids where dx != dy, but the question is
whether we
should add it as a patch for met-5.1 or just include it in the next
release
of the code.
The enhancement will touch several files, and we worry about the
potential
for unintended consequences of the changes.
Is supporting lambert conformal and polar stereographic grids where dx
!=
dy a requirement for your work? Is this shortcoming standing in your
way
of using MET?
Thanks,
John
On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Jose, yes, running with dx = dy would solve the problem.
>
> John
>
> On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>>
>> Dear John,
>>
>> Thank you for your insights. One last question: If we run the model
with
>> dx = dy, should that solve the problem?
>>
>> Best Regards,
>> Jose
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Wednesday, November 04, 2015 1:02 PM
>> To: Jose Garcia-Rivera
>> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>>
>> Jose,
>>
>> Yes, that's correct. Not supporting differing Lambert Conformal Dx
and
>> Dy values is a short-coming in MET. However, we're investigating
the
>> impact/cost of adding support for it. We may be able to add
support by
>> patching the build.
>>
>> I'll write you back once we figure out how to proceed.
>>
>> Thanks,
>> John
>>
>> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >
>> > Hi John,
>> >
>> > Thanks so much for the information. So, that means that I'll have
the
>> > same problem when running any other analyses with this particular
grid
>> using MET?
>> >
>> > Best Regards,
>> > Jose
>> >
>> > -----Original Message-----
>> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > Sent: Wednesday, November 04, 2015 12:29 PM
>> > To: Jose Garcia-Rivera
>> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >
>> > Jose,
>> >
>> > I believe I've figured out the problem. The grid code in MET
>> > currently only supports Lambert Conformal grids where the grid
spacing
>> > is the same in the x and y directions.
>> >
>> > Your grid is spaced 3.317km in the x-direction and 2.97km in the
>> > y-direction. The MET code is reading the Dx value and using it
for
>> > both x and y. This causes the skewing of the data to the north
that
>> you're seeing.
>> >
>> > Really, there are 2 problems here.
>> >
>> > First, if we're requiring that Dx and Dy be the same for Lambert
>> > Conformal grids, MET should error out with a descriptive error
message
>> > when they are not.
>> >
>> > Second, obviously it would be preferable for MET to support
Lambert
>> > Conformal grids where Dx and Dy are not the same! I'll ask
around
>> > here for more information on that.
>> >
>> > Thanks for paying close attention to your data and uncovering
this
>> > shortcoming in MET.
>> >
>> > John
>> >
>> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
>> > met_help at ucar.edu
>> > > wrote:
>> >
>> > >
>> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> > >
>> > > Hi John,
>> > >
>> > > Thanks for all the help. I want to add this pdf, showing
something
>> > > odd with the regridding. The domains did not match perfectly on
the
>> > > north side, it's like a skewing effect. The south sides match
>> > > perfectly, but note the differences in the regridded RTMA file
and
>> > > the model file in Michigan, for example. Maybe that has
something to
>> do with it.
>> > >
>> > > Best Regards,
>> > > Jose
>> > >
>> > > -----Original Message-----
>> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > > Sent: Tuesday, November 03, 2015 5:29 PM
>> > > To: Jose Garcia-Rivera
>> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> > >
>> > > Jose,
>> > >
>> > > I believe the problem is in our handling the domain on which
your
>> > > model data resides. I ran the following commands to plot
>> > > temperature from your model data file:
>> > >
>> > > met-5.1/bin/plot_data_plane \
>> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
>> > > 'name="TMP"; level="Z0";'
>> > >
>> > > convert -rotate 90 -background white -flatten \
>> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>> > >
>> > > The resulting image is attached and shows the shift you
described.
>> > >
>> > > So it isn't the regridding step that's the problem... its our
>> > > handling of that domain. I plotted the same data using IDV and
it
>> > > puts things in the correct place (also attached).
>> > >
>> > > One issue is that the 2-m temperature data in GRIB record 296
>> > > contains a constant value of 273. The surface temperature data
in
>> > > GRIB record
>> > > 276 contains more reasonable values. So that's the field I'm
>> > investigating.
>> > > But you might want to check why GRIB record 296 looks so odd.
>> > >
>> > > I need to head out of the office, but will investigate further
>> tomorrow.
>> > >
>> > > Thanks,
>> > > John
>> > >
>> > >
>> > >
>> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
>> > > met_help at ucar.edu
>> > > > wrote:
>> > >
>> > > >
>> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
>> > > >
>> > > > Dear John,
>> > > >
>> > > > Thank you so much for the quick reply.
>> > > >
>> > > > I followed the instructions, and ftp'ed the files to
>> > > "garcia-rivera_data".
>> > > >
>> > > > The nmbprs file is the .grb NMMB output file.
>> > > > The rtma file is the observational.
>> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs in
the
>> > > > same grid, and produces the plot.
>> > > >
>> > > > Any help will be greatly appreciated.
>> > > >
>> > > > Best Regards,
>> > > > Jose
>> > > >
>> > > > -----Original Message-----
>> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > > > Sent: Tuesday, November 03, 2015 4:04 PM
>> > > > To: Jose Garcia-Rivera
>> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> > > >
>> > > > Jose,
>> > > >
>> > > > Thanks for letting us know about the issue you're seeing.
I'd
>> > > > like to replicate the behavior your describe using the exact
data
>> you used.
>> > > > Once I can replicate it, I'll be able to investigate further.
>> > > >
>> > > > Are you able to send me the following 2 files you used:
>> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>> > > >
>> > > > You can post them to our anonymous ftp site, following these
>> > > instructions:
>> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> > > >
>> > > > Please write back to let me know when you've posted the data
and
>> > > > I'll go grab it.
>> > > >
>> > > > Thanks,
>> > > > John Halley Gotway
>> > > > met_help at ucar.edu
>> > > >
>> > > >
>> > > >
>> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
>> > > > met_help at ucar.edu
>> > > > > wrote:
>> > > >
>> > > > >
>> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
>> > > > > Transaction: Ticket created by GarciaJ at imsg.com
>> > > > > Queue: met_help
>> > > > > Subject: Regridding problem
>> > > > > Owner: Nobody
>> > > > > Requestors: GarciaJ at imsg.com
>> > > > > Status: new
>> > > > > Ticket <URL:
>> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>> > > > > >
>> > > > >
>> > > > >
>> > > > > Hello,
>> > > > >
>> > > > > I'm a new MET user, and I'm attempting to regrid RTMA data
to an
>> > > > > NMMB output grid. I used the following tool:
>> > > > >
>> > > > >
>> > > > >
>> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
>> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP"; level="Z2";'
>> > > > > -method BILIN -width
>> > > > > 2 -v 1
>> > > > >
>> > > > > which worked. But, when I subtracted the fields from one
dataset
>> > > > > to the other, there was a clear offset on the values. So,
the
>> > > > > regridding was done indeed, but the gridpoint values did
not
>> > > > > seem to be assigned correctly with respect to their
counterparts.
>> > > > >
>> > > > > I attached the plot, so you can see the offset problem. Any
help
>> > > > > will be greatly appreciated.
>> > > > >
>> > > > > Best Regards,
>> > > > > Jose
>> > > > >
>> > > > > Jose M. Garcia-Rivera, PhD
>> > > > > NWP Scientist
>> > > > > I.M. Systems Group
>> > > > > This e-mail, including attachments, may include
confidential
>> > > > > and/or proprietary information, and may be used only by the
>> > > > > person or entity to which it is addressed. If the reader of
this
>> > > > > e-mail is not the intended recipient or his or her
authorized
>> > > > > agent, the reader is hereby notified that any
dissemination,
>> > > > > distribution or copying of this e-mail is prohibited. If
you
>> > > > > have received this e-mail in error, please notify the
sender by
>> > > > > replying to this message and delete this
>> > > > e-mail immediately.
>> > > > >
>> > > > >
>> > > >
>> > > > This e-mail, including attachments, may include confidential
>> > > > and/or proprietary information, and may be used only by the
person
>> > > > or entity to which it is addressed. If the reader of this e-
mail
>> > > > is not the intended recipient or his or her authorized agent,
the
>> > > > reader is hereby notified that any dissemination,
distribution or
>> > > > copying of this e-mail is prohibited. If you have received
this
>> > > > e-mail in error, please notify the sender by replying to this
>> > > > message and delete this
>> > > e-mail immediately.
>> > > >
>> > > >
>> > > >
>> > >
>> > > This e-mail, including attachments, may include confidential
and/or
>> > > proprietary information, and may be used only by the person or
>> > > entity to which it is addressed. If the reader of this e-mail
is not
>> > > the intended recipient or his or her authorized agent, the
reader is
>> > > hereby notified that any dissemination, distribution or copying
of
>> > > this e-mail is prohibited. If you have received this e-mail in
>> > > error, please notify the sender by replying to this message and
>> > > delete this
>> > e-mail immediately.
>> > >
>> > >
>> >
>> > This e-mail, including attachments, may include confidential
and/or
>> > proprietary information, and may be used only by the person or
entity
>> > to which it is addressed. If the reader of this e-mail is not the
>> > intended recipient or his or her authorized agent, the reader is
>> > hereby notified that any dissemination, distribution or copying
of
>> > this e-mail is prohibited. If you have received this e-mail in
error,
>> > please notify the sender by replying to this message and delete
this
>> e-mail immediately.
>> >
>> >
>> >
>>
>> This e-mail, including attachments, may include confidential and/or
>> proprietary information, and may be used only by the person or
entity to
>> which it is addressed. If the reader of this e-mail is not the
intended
>> recipient or his or her authorized agent, the reader is hereby
notified
>> that any dissemination, distribution or copying of this e-mail is
>> prohibited. If you have received this e-mail in error, please
notify the
>> sender by replying to this message and delete this e-mail
immediately.
>>
>>
>>
>
------------------------------------------------
Subject: Regridding problem
From: Jose Garcia-Rivera
Time: Wed Nov 04 14:35:48 2015
Dear John,
I just ran a simulation in which DX=DY, and the skewing was much
smaller, but it's still observed.
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Wednesday, November 04, 2015 4:28 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
We're still trying to decide how to address this issue. We definitely
plan on adding support for grids where dx != dy, but the question is
whether we should add it as a patch for met-5.1 or just include it in
the next release of the code.
The enhancement will touch several files, and we worry about the
potential for unintended consequences of the changes.
Is supporting lambert conformal and polar stereographic grids where dx
!= dy a requirement for your work? Is this shortcoming standing in
your way of using MET?
Thanks,
John
On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Jose, yes, running with dx = dy would solve the problem.
>
> John
>
> On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>>
>> Dear John,
>>
>> Thank you for your insights. One last question: If we run the model
>> with dx = dy, should that solve the problem?
>>
>> Best Regards,
>> Jose
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Wednesday, November 04, 2015 1:02 PM
>> To: Jose Garcia-Rivera
>> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>>
>> Jose,
>>
>> Yes, that's correct. Not supporting differing Lambert Conformal Dx
>> and Dy values is a short-coming in MET. However, we're
investigating
>> the impact/cost of adding support for it. We may be able to add
>> support by patching the build.
>>
>> I'll write you back once we figure out how to proceed.
>>
>> Thanks,
>> John
>>
>> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >
>> > Hi John,
>> >
>> > Thanks so much for the information. So, that means that I'll have
>> > the same problem when running any other analyses with this
>> > particular grid
>> using MET?
>> >
>> > Best Regards,
>> > Jose
>> >
>> > -----Original Message-----
>> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > Sent: Wednesday, November 04, 2015 12:29 PM
>> > To: Jose Garcia-Rivera
>> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >
>> > Jose,
>> >
>> > I believe I've figured out the problem. The grid code in MET
>> > currently only supports Lambert Conformal grids where the grid
>> > spacing is the same in the x and y directions.
>> >
>> > Your grid is spaced 3.317km in the x-direction and 2.97km in the
>> > y-direction. The MET code is reading the Dx value and using it
for
>> > both x and y. This causes the skewing of the data to the north
>> > that
>> you're seeing.
>> >
>> > Really, there are 2 problems here.
>> >
>> > First, if we're requiring that Dx and Dy be the same for Lambert
>> > Conformal grids, MET should error out with a descriptive error
>> > message when they are not.
>> >
>> > Second, obviously it would be preferable for MET to support
Lambert
>> > Conformal grids where Dx and Dy are not the same! I'll ask
around
>> > here for more information on that.
>> >
>> > Thanks for paying close attention to your data and uncovering
this
>> > shortcoming in MET.
>> >
>> > John
>> >
>> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
>> > met_help at ucar.edu
>> > > wrote:
>> >
>> > >
>> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> > >
>> > > Hi John,
>> > >
>> > > Thanks for all the help. I want to add this pdf, showing
>> > > something odd with the regridding. The domains did not match
>> > > perfectly on the north side, it's like a skewing effect. The
>> > > south sides match perfectly, but note the differences in the
>> > > regridded RTMA file and the model file in Michigan, for
example.
>> > > Maybe that has something to
>> do with it.
>> > >
>> > > Best Regards,
>> > > Jose
>> > >
>> > > -----Original Message-----
>> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > > Sent: Tuesday, November 03, 2015 5:29 PM
>> > > To: Jose Garcia-Rivera
>> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> > >
>> > > Jose,
>> > >
>> > > I believe the problem is in our handling the domain on which
your
>> > > model data resides. I ran the following commands to plot
>> > > temperature from your model data file:
>> > >
>> > > met-5.1/bin/plot_data_plane \
>> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
>> > > 'name="TMP"; level="Z0";'
>> > >
>> > > convert -rotate 90 -background white -flatten \
>> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>> > >
>> > > The resulting image is attached and shows the shift you
described.
>> > >
>> > > So it isn't the regridding step that's the problem... its our
>> > > handling of that domain. I plotted the same data using IDV and
>> > > it puts things in the correct place (also attached).
>> > >
>> > > One issue is that the 2-m temperature data in GRIB record 296
>> > > contains a constant value of 273. The surface temperature data
>> > > in GRIB record
>> > > 276 contains more reasonable values. So that's the field I'm
>> > investigating.
>> > > But you might want to check why GRIB record 296 looks so odd.
>> > >
>> > > I need to head out of the office, but will investigate further
>> tomorrow.
>> > >
>> > > Thanks,
>> > > John
>> > >
>> > >
>> > >
>> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
>> > > met_help at ucar.edu
>> > > > wrote:
>> > >
>> > > >
>> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
>> > > >
>> > > > Dear John,
>> > > >
>> > > > Thank you so much for the quick reply.
>> > > >
>> > > > I followed the instructions, and ftp'ed the files to
>> > > "garcia-rivera_data".
>> > > >
>> > > > The nmbprs file is the .grb NMMB output file.
>> > > > The rtma file is the observational.
>> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs in
>> > > > the same grid, and produces the plot.
>> > > >
>> > > > Any help will be greatly appreciated.
>> > > >
>> > > > Best Regards,
>> > > > Jose
>> > > >
>> > > > -----Original Message-----
>> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> > > > Sent: Tuesday, November 03, 2015 4:04 PM
>> > > > To: Jose Garcia-Rivera
>> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> > > >
>> > > > Jose,
>> > > >
>> > > > Thanks for letting us know about the issue you're seeing.
I'd
>> > > > like to replicate the behavior your describe using the exact
>> > > > data
>> you used.
>> > > > Once I can replicate it, I'll be able to investigate further.
>> > > >
>> > > > Are you able to send me the following 2 files you used:
>> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>> > > >
>> > > > You can post them to our anonymous ftp site, following these
>> > > instructions:
>> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> > > >
>> > > > Please write back to let me know when you've posted the data
>> > > > and I'll go grab it.
>> > > >
>> > > > Thanks,
>> > > > John Halley Gotway
>> > > > met_help at ucar.edu
>> > > >
>> > > >
>> > > >
>> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
>> > > > met_help at ucar.edu
>> > > > > wrote:
>> > > >
>> > > > >
>> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
>> > > > > Transaction: Ticket created by GarciaJ at imsg.com
>> > > > > Queue: met_help
>> > > > > Subject: Regridding problem
>> > > > > Owner: Nobody
>> > > > > Requestors: GarciaJ at imsg.com
>> > > > > Status: new
>> > > > > Ticket <URL:
>> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>> > > > > >
>> > > > >
>> > > > >
>> > > > > Hello,
>> > > > >
>> > > > > I'm a new MET user, and I'm attempting to regrid RTMA data
to
>> > > > > an NMMB output grid. I used the following tool:
>> > > > >
>> > > > >
>> > > > >
>> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
>> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP"; level="Z2";'
>> > > > > -method BILIN -width
>> > > > > 2 -v 1
>> > > > >
>> > > > > which worked. But, when I subtracted the fields from one
>> > > > > dataset to the other, there was a clear offset on the
values.
>> > > > > So, the regridding was done indeed, but the gridpoint
values
>> > > > > did not seem to be assigned correctly with respect to their
counterparts.
>> > > > >
>> > > > > I attached the plot, so you can see the offset problem. Any
>> > > > > help will be greatly appreciated.
>> > > > >
>> > > > > Best Regards,
>> > > > > Jose
>> > > > >
>> > > > > Jose M. Garcia-Rivera, PhD
>> > > > > NWP Scientist
>> > > > > I.M. Systems Group
>> > > > > This e-mail, including attachments, may include
confidential
>> > > > > and/or proprietary information, and may be used only by the
>> > > > > person or entity to which it is addressed. If the reader of
>> > > > > this e-mail is not the intended recipient or his or her
>> > > > > authorized agent, the reader is hereby notified that any
>> > > > > dissemination, distribution or copying of this e-mail is
>> > > > > prohibited. If you have received this e-mail in error,
please
>> > > > > notify the sender by replying to this message and delete
this
>> > > > e-mail immediately.
>> > > > >
>> > > > >
>> > > >
>> > > > This e-mail, including attachments, may include confidential
>> > > > and/or proprietary information, and may be used only by the
>> > > > person or entity to which it is addressed. If the reader of
>> > > > this e-mail is not the intended recipient or his or her
>> > > > authorized agent, the reader is hereby notified that any
>> > > > dissemination, distribution or copying of this e-mail is
>> > > > prohibited. If you have received this e-mail in error, please
>> > > > notify the sender by replying to this message and delete this
>> > > e-mail immediately.
>> > > >
>> > > >
>> > > >
>> > >
>> > > This e-mail, including attachments, may include confidential
>> > > and/or proprietary information, and may be used only by the
>> > > person or entity to which it is addressed. If the reader of
this
>> > > e-mail is not the intended recipient or his or her authorized
>> > > agent, the reader is hereby notified that any dissemination,
>> > > distribution or copying of this e-mail is prohibited. If you
have
>> > > received this e-mail in error, please notify the sender by
>> > > replying to this message and delete this
>> > e-mail immediately.
>> > >
>> > >
>> >
>> > This e-mail, including attachments, may include confidential
and/or
>> > proprietary information, and may be used only by the person or
>> > entity to which it is addressed. If the reader of this e-mail is
>> > not the intended recipient or his or her authorized agent, the
>> > reader is hereby notified that any dissemination, distribution or
>> > copying of this e-mail is prohibited. If you have received this
>> > e-mail in error, please notify the sender by replying to this
>> > message and delete this
>> e-mail immediately.
>> >
>> >
>> >
>>
>> This e-mail, including attachments, may include confidential and/or
>> proprietary information, and may be used only by the person or
entity
>> to which it is addressed. If the reader of this e-mail is not the
>> intended recipient or his or her authorized agent, the reader is
>> hereby notified that any dissemination, distribution or copying of
>> this e-mail is prohibited. If you have received this e-mail in
error,
>> please notify the sender by replying to this message and delete
this e-mail immediately.
>>
>>
>>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Wed Nov 04 15:03:32 2015
Jose,
Are you able to send me the GRIB files as well by posting it to our
anonymous ftp site?
I took a look at the images you sent, but the skew you mention isn't
jumping out at me. Before I saw that the values that should be over
the
great lakes were shifted north by a large extent. What should I be
looking
for now?
Rather than diffing the NMM and RTMA values, wouldn't it make more
sense
to...
(1) Run the RTMA through regrid_data_plane to regrid it
(2) Run plot_data_plane to plot the input to regrid_data_plane
(3) Run plot_data_plane to plot the output from regrid_data_plane
Compare images from the input/output to make sure the values are all
in the
right places.
Thanks,
John
On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> I just ran a simulation in which DX=DY, and the skewing was much
smaller,
> but it's still observed.
>
> Best Regards,
> Jose
>
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Wednesday, November 04, 2015 4:28 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> We're still trying to decide how to address this issue. We
definitely
> plan on adding support for grids where dx != dy, but the question is
> whether we should add it as a patch for met-5.1 or just include it
in the
> next release of the code.
>
> The enhancement will touch several files, and we worry about the
potential
> for unintended consequences of the changes.
>
> Is supporting lambert conformal and polar stereographic grids where
dx !=
> dy a requirement for your work? Is this shortcoming standing in
your way
> of using MET?
>
> Thanks,
> John
>
> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
<johnhg at ucar.edu>
> wrote:
>
> > Jose, yes, running with dx = dy would solve the problem.
> >
> > John
> >
> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >>
> >> Dear John,
> >>
> >> Thank you for your insights. One last question: If we run the
model
> >> with dx = dy, should that solve the problem?
> >>
> >> Best Regards,
> >> Jose
> >>
> >> -----Original Message-----
> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> Sent: Wednesday, November 04, 2015 1:02 PM
> >> To: Jose Garcia-Rivera
> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >>
> >> Jose,
> >>
> >> Yes, that's correct. Not supporting differing Lambert Conformal
Dx
> >> and Dy values is a short-coming in MET. However, we're
investigating
> >> the impact/cost of adding support for it. We may be able to add
> >> support by patching the build.
> >>
> >> I'll write you back once we figure out how to proceed.
> >>
> >> Thanks,
> >> John
> >>
> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >> >
> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >> >
> >> > Hi John,
> >> >
> >> > Thanks so much for the information. So, that means that I'll
have
> >> > the same problem when running any other analyses with this
> >> > particular grid
> >> using MET?
> >> >
> >> > Best Regards,
> >> > Jose
> >> >
> >> > -----Original Message-----
> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> >> > To: Jose Garcia-Rivera
> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >
> >> > Jose,
> >> >
> >> > I believe I've figured out the problem. The grid code in MET
> >> > currently only supports Lambert Conformal grids where the grid
> >> > spacing is the same in the x and y directions.
> >> >
> >> > Your grid is spaced 3.317km in the x-direction and 2.97km in
the
> >> > y-direction. The MET code is reading the Dx value and using it
for
> >> > both x and y. This causes the skewing of the data to the north
> >> > that
> >> you're seeing.
> >> >
> >> > Really, there are 2 problems here.
> >> >
> >> > First, if we're requiring that Dx and Dy be the same for
Lambert
> >> > Conformal grids, MET should error out with a descriptive error
> >> > message when they are not.
> >> >
> >> > Second, obviously it would be preferable for MET to support
Lambert
> >> > Conformal grids where Dx and Dy are not the same! I'll ask
around
> >> > here for more information on that.
> >> >
> >> > Thanks for paying close attention to your data and uncovering
this
> >> > shortcoming in MET.
> >> >
> >> > John
> >> >
> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> >> > met_help at ucar.edu
> >> > > wrote:
> >> >
> >> > >
> >> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> >> > >
> >> > > Hi John,
> >> > >
> >> > > Thanks for all the help. I want to add this pdf, showing
> >> > > something odd with the regridding. The domains did not match
> >> > > perfectly on the north side, it's like a skewing effect. The
> >> > > south sides match perfectly, but note the differences in the
> >> > > regridded RTMA file and the model file in Michigan, for
example.
> >> > > Maybe that has something to
> >> do with it.
> >> > >
> >> > > Best Regards,
> >> > > Jose
> >> > >
> >> > > -----Original Message-----
> >> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> >> > > To: Jose Garcia-Rivera
> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> > >
> >> > > Jose,
> >> > >
> >> > > I believe the problem is in our handling the domain on which
your
> >> > > model data resides. I ran the following commands to plot
> >> > > temperature from your model data file:
> >> > >
> >> > > met-5.1/bin/plot_data_plane \
> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> >> > > 'name="TMP"; level="Z0";'
> >> > >
> >> > > convert -rotate 90 -background white -flatten \
> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> >> > >
> >> > > The resulting image is attached and shows the shift you
described.
> >> > >
> >> > > So it isn't the regridding step that's the problem... its our
> >> > > handling of that domain. I plotted the same data using IDV
and
> >> > > it puts things in the correct place (also attached).
> >> > >
> >> > > One issue is that the 2-m temperature data in GRIB record 296
> >> > > contains a constant value of 273. The surface temperature
data
> >> > > in GRIB record
> >> > > 276 contains more reasonable values. So that's the field I'm
> >> > investigating.
> >> > > But you might want to check why GRIB record 296 looks so odd.
> >> > >
> >> > > I need to head out of the office, but will investigate
further
> >> tomorrow.
> >> > >
> >> > > Thanks,
> >> > > John
> >> > >
> >> > >
> >> > >
> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
> >> > > met_help at ucar.edu
> >> > > > wrote:
> >> > >
> >> > > >
> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >> > > >
> >> > > > Dear John,
> >> > > >
> >> > > > Thank you so much for the quick reply.
> >> > > >
> >> > > > I followed the instructions, and ftp'ed the files to
> >> > > "garcia-rivera_data".
> >> > > >
> >> > > > The nmbprs file is the .grb NMMB output file.
> >> > > > The rtma file is the observational.
> >> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs
in
> >> > > > the same grid, and produces the plot.
> >> > > >
> >> > > > Any help will be greatly appreciated.
> >> > > >
> >> > > > Best Regards,
> >> > > > Jose
> >> > > >
> >> > > > -----Original Message-----
> >> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> >> > > > To: Jose Garcia-Rivera
> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> > > >
> >> > > > Jose,
> >> > > >
> >> > > > Thanks for letting us know about the issue you're seeing.
I'd
> >> > > > like to replicate the behavior your describe using the
exact
> >> > > > data
> >> you used.
> >> > > > Once I can replicate it, I'll be able to investigate
further.
> >> > > >
> >> > > > Are you able to send me the following 2 files you used:
> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> >> > > >
> >> > > > You can post them to our anonymous ftp site, following
these
> >> > > instructions:
> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >> > > >
> >> > > > Please write back to let me know when you've posted the
data
> >> > > > and I'll go grab it.
> >> > > >
> >> > > > Thanks,
> >> > > > John Halley Gotway
> >> > > > met_help at ucar.edu
> >> > > >
> >> > > >
> >> > > >
> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT <
> >> > > > met_help at ucar.edu
> >> > > > > wrote:
> >> > > >
> >> > > > >
> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> >> > > > > Queue: met_help
> >> > > > > Subject: Regridding problem
> >> > > > > Owner: Nobody
> >> > > > > Requestors: GarciaJ at imsg.com
> >> > > > > Status: new
> >> > > > > Ticket <URL:
> >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> > > > > >
> >> > > > >
> >> > > > >
> >> > > > > Hello,
> >> > > > >
> >> > > > > I'm a new MET user, and I'm attempting to regrid RTMA
data to
> >> > > > > an NMMB output grid. I used the following tool:
> >> > > > >
> >> > > > >
> >> > > > >
> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
level="Z2";'
> >> > > > > -method BILIN -width
> >> > > > > 2 -v 1
> >> > > > >
> >> > > > > which worked. But, when I subtracted the fields from one
> >> > > > > dataset to the other, there was a clear offset on the
values.
> >> > > > > So, the regridding was done indeed, but the gridpoint
values
> >> > > > > did not seem to be assigned correctly with respect to
their
> counterparts.
> >> > > > >
> >> > > > > I attached the plot, so you can see the offset problem.
Any
> >> > > > > help will be greatly appreciated.
> >> > > > >
> >> > > > > Best Regards,
> >> > > > > Jose
> >> > > > >
> >> > > > > Jose M. Garcia-Rivera, PhD
> >> > > > > NWP Scientist
> >> > > > > I.M. Systems Group
> >> > > > > This e-mail, including attachments, may include
confidential
> >> > > > > and/or proprietary information, and may be used only by
the
> >> > > > > person or entity to which it is addressed. If the reader
of
> >> > > > > this e-mail is not the intended recipient or his or her
> >> > > > > authorized agent, the reader is hereby notified that any
> >> > > > > dissemination, distribution or copying of this e-mail is
> >> > > > > prohibited. If you have received this e-mail in error,
please
> >> > > > > notify the sender by replying to this message and delete
this
> >> > > > e-mail immediately.
> >> > > > >
> >> > > > >
> >> > > >
> >> > > > This e-mail, including attachments, may include
confidential
> >> > > > and/or proprietary information, and may be used only by the
> >> > > > person or entity to which it is addressed. If the reader of
> >> > > > this e-mail is not the intended recipient or his or her
> >> > > > authorized agent, the reader is hereby notified that any
> >> > > > dissemination, distribution or copying of this e-mail is
> >> > > > prohibited. If you have received this e-mail in error,
please
> >> > > > notify the sender by replying to this message and delete
this
> >> > > e-mail immediately.
> >> > > >
> >> > > >
> >> > > >
> >> > >
> >> > > This e-mail, including attachments, may include confidential
> >> > > and/or proprietary information, and may be used only by the
> >> > > person or entity to which it is addressed. If the reader of
this
> >> > > e-mail is not the intended recipient or his or her authorized
> >> > > agent, the reader is hereby notified that any dissemination,
> >> > > distribution or copying of this e-mail is prohibited. If you
have
> >> > > received this e-mail in error, please notify the sender by
> >> > > replying to this message and delete this
> >> > e-mail immediately.
> >> > >
> >> > >
> >> >
> >> > This e-mail, including attachments, may include confidential
and/or
> >> > proprietary information, and may be used only by the person or
> >> > entity to which it is addressed. If the reader of this e-mail
is
> >> > not the intended recipient or his or her authorized agent, the
> >> > reader is hereby notified that any dissemination, distribution
or
> >> > copying of this e-mail is prohibited. If you have received this
> >> > e-mail in error, please notify the sender by replying to this
> >> > message and delete this
> >> e-mail immediately.
> >> >
> >> >
> >> >
> >>
> >> This e-mail, including attachments, may include confidential
and/or
> >> proprietary information, and may be used only by the person or
entity
> >> to which it is addressed. If the reader of this e-mail is not the
> >> intended recipient or his or her authorized agent, the reader is
> >> hereby notified that any dissemination, distribution or copying
of
> >> this e-mail is prohibited. If you have received this e-mail in
error,
> >> please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >>
> >>
> >>
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Thu Nov 05 16:31:16 2015
Jose,
Just wanted to follow up. Do you still see a shifting problem in the
regridding? And is this a showstopper for you guys? Is adding
support for
differing dx/dy as a patch for 5.1 a high priority or can we delay
until
the next release?
Thanks,
John
On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Jose,
>
> Are you able to send me the GRIB files as well by posting it to our
> anonymous ftp site?
>
> I took a look at the images you sent, but the skew you mention isn't
> jumping out at me. Before I saw that the values that should be over
the
> great lakes were shifted north by a large extent. What should I be
looking
> for now?
>
> Rather than diffing the NMM and RTMA values, wouldn't it make more
sense
> to...
>
> (1) Run the RTMA through regrid_data_plane to regrid it
> (2) Run plot_data_plane to plot the input to regrid_data_plane
> (3) Run plot_data_plane to plot the output from regrid_data_plane
>
> Compare images from the input/output to make sure the values are all
in
> the right places.
>
> Thanks,
> John
>
>
> On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>>
>> Dear John,
>>
>> I just ran a simulation in which DX=DY, and the skewing was much
smaller,
>> but it's still observed.
>>
>> Best Regards,
>> Jose
>>
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Wednesday, November 04, 2015 4:28 PM
>> To: Jose Garcia-Rivera
>> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>>
>> Jose,
>>
>> We're still trying to decide how to address this issue. We
definitely
>> plan on adding support for grids where dx != dy, but the question
is
>> whether we should add it as a patch for met-5.1 or just include it
in the
>> next release of the code.
>>
>> The enhancement will touch several files, and we worry about the
>> potential for unintended consequences of the changes.
>>
>> Is supporting lambert conformal and polar stereographic grids where
dx !=
>> dy a requirement for your work? Is this shortcoming standing in
your way
>> of using MET?
>>
>> Thanks,
>> John
>>
>> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
<johnhg at ucar.edu>
>> wrote:
>>
>> > Jose, yes, running with dx = dy would solve the problem.
>> >
>> > John
>> >
>> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
>> > met_help at ucar.edu> wrote:
>> >
>> >>
>> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >>
>> >> Dear John,
>> >>
>> >> Thank you for your insights. One last question: If we run the
model
>> >> with dx = dy, should that solve the problem?
>> >>
>> >> Best Regards,
>> >> Jose
>> >>
>> >> -----Original Message-----
>> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> Sent: Wednesday, November 04, 2015 1:02 PM
>> >> To: Jose Garcia-Rivera
>> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >>
>> >> Jose,
>> >>
>> >> Yes, that's correct. Not supporting differing Lambert Conformal
Dx
>> >> and Dy values is a short-coming in MET. However, we're
investigating
>> >> the impact/cost of adding support for it. We may be able to add
>> >> support by patching the build.
>> >>
>> >> I'll write you back once we figure out how to proceed.
>> >>
>> >> Thanks,
>> >> John
>> >>
>> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
>> >> met_help at ucar.edu> wrote:
>> >>
>> >> >
>> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
>> >> >
>> >> > Hi John,
>> >> >
>> >> > Thanks so much for the information. So, that means that I'll
have
>> >> > the same problem when running any other analyses with this
>> >> > particular grid
>> >> using MET?
>> >> >
>> >> > Best Regards,
>> >> > Jose
>> >> >
>> >> > -----Original Message-----
>> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > Sent: Wednesday, November 04, 2015 12:29 PM
>> >> > To: Jose Garcia-Rivera
>> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> >
>> >> > Jose,
>> >> >
>> >> > I believe I've figured out the problem. The grid code in MET
>> >> > currently only supports Lambert Conformal grids where the grid
>> >> > spacing is the same in the x and y directions.
>> >> >
>> >> > Your grid is spaced 3.317km in the x-direction and 2.97km in
the
>> >> > y-direction. The MET code is reading the Dx value and using
it for
>> >> > both x and y. This causes the skewing of the data to the
north
>> >> > that
>> >> you're seeing.
>> >> >
>> >> > Really, there are 2 problems here.
>> >> >
>> >> > First, if we're requiring that Dx and Dy be the same for
Lambert
>> >> > Conformal grids, MET should error out with a descriptive error
>> >> > message when they are not.
>> >> >
>> >> > Second, obviously it would be preferable for MET to support
Lambert
>> >> > Conformal grids where Dx and Dy are not the same! I'll ask
around
>> >> > here for more information on that.
>> >> >
>> >> > Thanks for paying close attention to your data and uncovering
this
>> >> > shortcoming in MET.
>> >> >
>> >> > John
>> >> >
>> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
>> >> > met_help at ucar.edu
>> >> > > wrote:
>> >> >
>> >> > >
>> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >> > >
>> >> > > Hi John,
>> >> > >
>> >> > > Thanks for all the help. I want to add this pdf, showing
>> >> > > something odd with the regridding. The domains did not match
>> >> > > perfectly on the north side, it's like a skewing effect. The
>> >> > > south sides match perfectly, but note the differences in the
>> >> > > regridded RTMA file and the model file in Michigan, for
example.
>> >> > > Maybe that has something to
>> >> do with it.
>> >> > >
>> >> > > Best Regards,
>> >> > > Jose
>> >> > >
>> >> > > -----Original Message-----
>> >> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
>> >> > > To: Jose Garcia-Rivera
>> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> > >
>> >> > > Jose,
>> >> > >
>> >> > > I believe the problem is in our handling the domain on which
your
>> >> > > model data resides. I ran the following commands to plot
>> >> > > temperature from your model data file:
>> >> > >
>> >> > > met-5.1/bin/plot_data_plane \
>> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
>> >> > > 'name="TMP"; level="Z0";'
>> >> > >
>> >> > > convert -rotate 90 -background white -flatten \
>> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>> >> > >
>> >> > > The resulting image is attached and shows the shift you
described.
>> >> > >
>> >> > > So it isn't the regridding step that's the problem... its
our
>> >> > > handling of that domain. I plotted the same data using IDV
and
>> >> > > it puts things in the correct place (also attached).
>> >> > >
>> >> > > One issue is that the 2-m temperature data in GRIB record
296
>> >> > > contains a constant value of 273. The surface temperature
data
>> >> > > in GRIB record
>> >> > > 276 contains more reasonable values. So that's the field
I'm
>> >> > investigating.
>> >> > > But you might want to check why GRIB record 296 looks so
odd.
>> >> > >
>> >> > > I need to head out of the office, but will investigate
further
>> >> tomorrow.
>> >> > >
>> >> > > Thanks,
>> >> > > John
>> >> > >
>> >> > >
>> >> > >
>> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
>> >> > > met_help at ucar.edu
>> >> > > > wrote:
>> >> > >
>> >> > > >
>> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >> > > >
>> >> > > > Dear John,
>> >> > > >
>> >> > > > Thank you so much for the quick reply.
>> >> > > >
>> >> > > > I followed the instructions, and ftp'ed the files to
>> >> > > "garcia-rivera_data".
>> >> > > >
>> >> > > > The nmbprs file is the .grb NMMB output file.
>> >> > > > The rtma file is the observational.
>> >> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs
in
>> >> > > > the same grid, and produces the plot.
>> >> > > >
>> >> > > > Any help will be greatly appreciated.
>> >> > > >
>> >> > > > Best Regards,
>> >> > > > Jose
>> >> > > >
>> >> > > > -----Original Message-----
>> >> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
>> >> > > > To: Jose Garcia-Rivera
>> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> > > >
>> >> > > > Jose,
>> >> > > >
>> >> > > > Thanks for letting us know about the issue you're seeing.
I'd
>> >> > > > like to replicate the behavior your describe using the
exact
>> >> > > > data
>> >> you used.
>> >> > > > Once I can replicate it, I'll be able to investigate
further.
>> >> > > >
>> >> > > > Are you able to send me the following 2 files you used:
>> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>> >> > > >
>> >> > > > You can post them to our anonymous ftp site, following
these
>> >> > > instructions:
>> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> >> > > >
>> >> > > > Please write back to let me know when you've posted the
data
>> >> > > > and I'll go grab it.
>> >> > > >
>> >> > > > Thanks,
>> >> > > > John Halley Gotway
>> >> > > > met_help at ucar.edu
>> >> > > >
>> >> > > >
>> >> > > >
>> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT
<
>> >> > > > met_help at ucar.edu
>> >> > > > > wrote:
>> >> > > >
>> >> > > > >
>> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
>> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
>> >> > > > > Queue: met_help
>> >> > > > > Subject: Regridding problem
>> >> > > > > Owner: Nobody
>> >> > > > > Requestors: GarciaJ at imsg.com
>> >> > > > > Status: new
>> >> > > > > Ticket <URL:
>> >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>> >> > > > > >
>> >> > > > >
>> >> > > > >
>> >> > > > > Hello,
>> >> > > > >
>> >> > > > > I'm a new MET user, and I'm attempting to regrid RTMA
data to
>> >> > > > > an NMMB output grid. I used the following tool:
>> >> > > > >
>> >> > > > >
>> >> > > > >
>> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
>> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
level="Z2";'
>> >> > > > > -method BILIN -width
>> >> > > > > 2 -v 1
>> >> > > > >
>> >> > > > > which worked. But, when I subtracted the fields from one
>> >> > > > > dataset to the other, there was a clear offset on the
values.
>> >> > > > > So, the regridding was done indeed, but the gridpoint
values
>> >> > > > > did not seem to be assigned correctly with respect to
their
>> counterparts.
>> >> > > > >
>> >> > > > > I attached the plot, so you can see the offset problem.
Any
>> >> > > > > help will be greatly appreciated.
>> >> > > > >
>> >> > > > > Best Regards,
>> >> > > > > Jose
>> >> > > > >
>> >> > > > > Jose M. Garcia-Rivera, PhD
>> >> > > > > NWP Scientist
>> >> > > > > I.M. Systems Group
>> >> > > > > This e-mail, including attachments, may include
confidential
>> >> > > > > and/or proprietary information, and may be used only by
the
>> >> > > > > person or entity to which it is addressed. If the reader
of
>> >> > > > > this e-mail is not the intended recipient or his or her
>> >> > > > > authorized agent, the reader is hereby notified that any
>> >> > > > > dissemination, distribution or copying of this e-mail is
>> >> > > > > prohibited. If you have received this e-mail in error,
please
>> >> > > > > notify the sender by replying to this message and delete
this
>> >> > > > e-mail immediately.
>> >> > > > >
>> >> > > > >
>> >> > > >
>> >> > > > This e-mail, including attachments, may include
confidential
>> >> > > > and/or proprietary information, and may be used only by
the
>> >> > > > person or entity to which it is addressed. If the reader
of
>> >> > > > this e-mail is not the intended recipient or his or her
>> >> > > > authorized agent, the reader is hereby notified that any
>> >> > > > dissemination, distribution or copying of this e-mail is
>> >> > > > prohibited. If you have received this e-mail in error,
please
>> >> > > > notify the sender by replying to this message and delete
this
>> >> > > e-mail immediately.
>> >> > > >
>> >> > > >
>> >> > > >
>> >> > >
>> >> > > This e-mail, including attachments, may include confidential
>> >> > > and/or proprietary information, and may be used only by the
>> >> > > person or entity to which it is addressed. If the reader of
this
>> >> > > e-mail is not the intended recipient or his or her
authorized
>> >> > > agent, the reader is hereby notified that any dissemination,
>> >> > > distribution or copying of this e-mail is prohibited. If you
have
>> >> > > received this e-mail in error, please notify the sender by
>> >> > > replying to this message and delete this
>> >> > e-mail immediately.
>> >> > >
>> >> > >
>> >> >
>> >> > This e-mail, including attachments, may include confidential
and/or
>> >> > proprietary information, and may be used only by the person or
>> >> > entity to which it is addressed. If the reader of this e-mail
is
>> >> > not the intended recipient or his or her authorized agent, the
>> >> > reader is hereby notified that any dissemination, distribution
or
>> >> > copying of this e-mail is prohibited. If you have received
this
>> >> > e-mail in error, please notify the sender by replying to this
>> >> > message and delete this
>> >> e-mail immediately.
>> >> >
>> >> >
>> >> >
>> >>
>> >> This e-mail, including attachments, may include confidential
and/or
>> >> proprietary information, and may be used only by the person or
entity
>> >> to which it is addressed. If the reader of this e-mail is not
the
>> >> intended recipient or his or her authorized agent, the reader is
>> >> hereby notified that any dissemination, distribution or copying
of
>> >> this e-mail is prohibited. If you have received this e-mail in
error,
>> >> please notify the sender by replying to this message and delete
this
>> e-mail immediately.
>> >>
>> >>
>> >>
>> >
>>
>> This e-mail, including attachments, may include confidential and/or
>> proprietary information, and may be used only by the person or
entity to
>> which it is addressed. If the reader of this e-mail is not the
intended
>> recipient or his or her authorized agent, the reader is hereby
notified
>> that any dissemination, distribution or copying of this e-mail is
>> prohibited. If you have received this e-mail in error, please
notify the
>> sender by replying to this message and delete this e-mail
immediately.
>>
>>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Fri Nov 06 08:54:40 2015
Dear John,
We want to use MET, but we need to get the regridding right. I was
wondering if you know how to regrid the RTMA file to match an NMMB
grid by using UPP (as in previous MET versions)?
Any help will be greatly appreciated,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, November 05, 2015 6:31 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
Just wanted to follow up. Do you still see a shifting problem in the
regridding? And is this a showstopper for you guys? Is adding
support for differing dx/dy as a patch for 5.1 a high priority or can
we delay until the next release?
Thanks,
John
On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Jose,
>
> Are you able to send me the GRIB files as well by posting it to our
> anonymous ftp site?
>
> I took a look at the images you sent, but the skew you mention isn't
> jumping out at me. Before I saw that the values that should be over
> the great lakes were shifted north by a large extent. What should I
> be looking for now?
>
> Rather than diffing the NMM and RTMA values, wouldn't it make more
> sense to...
>
> (1) Run the RTMA through regrid_data_plane to regrid it
> (2) Run plot_data_plane to plot the input to regrid_data_plane
> (3) Run plot_data_plane to plot the output from regrid_data_plane
>
> Compare images from the input/output to make sure the values are all
> in the right places.
>
> Thanks,
> John
>
>
> On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>>
>> Dear John,
>>
>> I just ran a simulation in which DX=DY, and the skewing was much
>> smaller, but it's still observed.
>>
>> Best Regards,
>> Jose
>>
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Wednesday, November 04, 2015 4:28 PM
>> To: Jose Garcia-Rivera
>> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>>
>> Jose,
>>
>> We're still trying to decide how to address this issue. We
>> definitely plan on adding support for grids where dx != dy, but the
>> question is whether we should add it as a patch for met-5.1 or just
>> include it in the next release of the code.
>>
>> The enhancement will touch several files, and we worry about the
>> potential for unintended consequences of the changes.
>>
>> Is supporting lambert conformal and polar stereographic grids where
>> dx != dy a requirement for your work? Is this shortcoming standing
>> in your way of using MET?
>>
>> Thanks,
>> John
>>
>> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
<johnhg at ucar.edu>
>> wrote:
>>
>> > Jose, yes, running with dx = dy would solve the problem.
>> >
>> > John
>> >
>> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
>> > met_help at ucar.edu> wrote:
>> >
>> >>
>> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >>
>> >> Dear John,
>> >>
>> >> Thank you for your insights. One last question: If we run the
>> >> model with dx = dy, should that solve the problem?
>> >>
>> >> Best Regards,
>> >> Jose
>> >>
>> >> -----Original Message-----
>> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> Sent: Wednesday, November 04, 2015 1:02 PM
>> >> To: Jose Garcia-Rivera
>> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >>
>> >> Jose,
>> >>
>> >> Yes, that's correct. Not supporting differing Lambert Conformal
>> >> Dx and Dy values is a short-coming in MET. However, we're
>> >> investigating the impact/cost of adding support for it. We may
be
>> >> able to add support by patching the build.
>> >>
>> >> I'll write you back once we figure out how to proceed.
>> >>
>> >> Thanks,
>> >> John
>> >>
>> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
>> >> met_help at ucar.edu> wrote:
>> >>
>> >> >
>> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
>> >> >
>> >> > Hi John,
>> >> >
>> >> > Thanks so much for the information. So, that means that I'll
>> >> > have the same problem when running any other analyses with
this
>> >> > particular grid
>> >> using MET?
>> >> >
>> >> > Best Regards,
>> >> > Jose
>> >> >
>> >> > -----Original Message-----
>> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > Sent: Wednesday, November 04, 2015 12:29 PM
>> >> > To: Jose Garcia-Rivera
>> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> >
>> >> > Jose,
>> >> >
>> >> > I believe I've figured out the problem. The grid code in MET
>> >> > currently only supports Lambert Conformal grids where the grid
>> >> > spacing is the same in the x and y directions.
>> >> >
>> >> > Your grid is spaced 3.317km in the x-direction and 2.97km in
the
>> >> > y-direction. The MET code is reading the Dx value and using
it
>> >> > for both x and y. This causes the skewing of the data to the
>> >> > north that
>> >> you're seeing.
>> >> >
>> >> > Really, there are 2 problems here.
>> >> >
>> >> > First, if we're requiring that Dx and Dy be the same for
Lambert
>> >> > Conformal grids, MET should error out with a descriptive error
>> >> > message when they are not.
>> >> >
>> >> > Second, obviously it would be preferable for MET to support
>> >> > Lambert Conformal grids where Dx and Dy are not the same!
I'll
>> >> > ask around here for more information on that.
>> >> >
>> >> > Thanks for paying close attention to your data and uncovering
>> >> > this shortcoming in MET.
>> >> >
>> >> > John
>> >> >
>> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
>> >> > met_help at ucar.edu
>> >> > > wrote:
>> >> >
>> >> > >
>> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>> >> > > >
>> >> > >
>> >> > > Hi John,
>> >> > >
>> >> > > Thanks for all the help. I want to add this pdf, showing
>> >> > > something odd with the regridding. The domains did not match
>> >> > > perfectly on the north side, it's like a skewing effect. The
>> >> > > south sides match perfectly, but note the differences in the
>> >> > > regridded RTMA file and the model file in Michigan, for
example.
>> >> > > Maybe that has something to
>> >> do with it.
>> >> > >
>> >> > > Best Regards,
>> >> > > Jose
>> >> > >
>> >> > > -----Original Message-----
>> >> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
>> >> > > To: Jose Garcia-Rivera
>> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> > >
>> >> > > Jose,
>> >> > >
>> >> > > I believe the problem is in our handling the domain on which
>> >> > > your model data resides. I ran the following commands to
plot
>> >> > > temperature from your model data file:
>> >> > >
>> >> > > met-5.1/bin/plot_data_plane \
>> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
>> >> > > 'name="TMP"; level="Z0";'
>> >> > >
>> >> > > convert -rotate 90 -background white -flatten \
>> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
>> >> > >
>> >> > > The resulting image is attached and shows the shift you
described.
>> >> > >
>> >> > > So it isn't the regridding step that's the problem... its
our
>> >> > > handling of that domain. I plotted the same data using IDV
>> >> > > and it puts things in the correct place (also attached).
>> >> > >
>> >> > > One issue is that the 2-m temperature data in GRIB record
296
>> >> > > contains a constant value of 273. The surface temperature
>> >> > > data in GRIB record
>> >> > > 276 contains more reasonable values. So that's the field
I'm
>> >> > investigating.
>> >> > > But you might want to check why GRIB record 296 looks so
odd.
>> >> > >
>> >> > > I need to head out of the office, but will investigate
further
>> >> tomorrow.
>> >> > >
>> >> > > Thanks,
>> >> > > John
>> >> > >
>> >> > >
>> >> > >
>> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT <
>> >> > > met_help at ucar.edu
>> >> > > > wrote:
>> >> > >
>> >> > > >
>> >> > > > <URL:
>> >> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>> >> > > >
>> >> > > > Dear John,
>> >> > > >
>> >> > > > Thank you so much for the quick reply.
>> >> > > >
>> >> > > > I followed the instructions, and ftp'ed the files to
>> >> > > "garcia-rivera_data".
>> >> > > >
>> >> > > > The nmbprs file is the .grb NMMB output file.
>> >> > > > The rtma file is the observational.
>> >> > > > The nmm_T2_bias.ncl subtracts the fields from both outputs
>> >> > > > in the same grid, and produces the plot.
>> >> > > >
>> >> > > > Any help will be greatly appreciated.
>> >> > > >
>> >> > > > Best Regards,
>> >> > > > Jose
>> >> > > >
>> >> > > > -----Original Message-----
>> >> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
>> >> > > > To: Jose Garcia-Rivera
>> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>> >> > > >
>> >> > > > Jose,
>> >> > > >
>> >> > > > Thanks for letting us know about the issue you're seeing.
>> >> > > > I'd like to replicate the behavior your describe using the
>> >> > > > exact data
>> >> you used.
>> >> > > > Once I can replicate it, I'll be able to investigate
further.
>> >> > > >
>> >> > > > Are you able to send me the following 2 files you used:
>> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
>> >> > > >
>> >> > > > You can post them to our anonymous ftp site, following
these
>> >> > > instructions:
>> >> > > >
>> >> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> >> > > >
>> >> > > > Please write back to let me know when you've posted the
data
>> >> > > > and I'll go grab it.
>> >> > > >
>> >> > > > Thanks,
>> >> > > > John Halley Gotway
>> >> > > > met_help at ucar.edu
>> >> > > >
>> >> > > >
>> >> > > >
>> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via RT
<
>> >> > > > met_help at ucar.edu
>> >> > > > > wrote:
>> >> > > >
>> >> > > > >
>> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted upon.
>> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
>> >> > > > > Queue: met_help
>> >> > > > > Subject: Regridding problem
>> >> > > > > Owner: Nobody
>> >> > > > > Requestors: GarciaJ at imsg.com
>> >> > > > > Status: new
>> >> > > > > Ticket <URL:
>> >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>> >> > > > > >
>> >> > > > >
>> >> > > > >
>> >> > > > > Hello,
>> >> > > > >
>> >> > > > > I'm a new MET user, and I'm attempting to regrid RTMA
data
>> >> > > > > to an NMMB output grid. I used the following tool:
>> >> > > > >
>> >> > > > >
>> >> > > > >
>> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
>> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
level="Z2";'
>> >> > > > > -method BILIN -width
>> >> > > > > 2 -v 1
>> >> > > > >
>> >> > > > > which worked. But, when I subtracted the fields from one
>> >> > > > > dataset to the other, there was a clear offset on the
values.
>> >> > > > > So, the regridding was done indeed, but the gridpoint
>> >> > > > > values did not seem to be assigned correctly with
respect
>> >> > > > > to their
>> counterparts.
>> >> > > > >
>> >> > > > > I attached the plot, so you can see the offset problem.
>> >> > > > > Any help will be greatly appreciated.
>> >> > > > >
>> >> > > > > Best Regards,
>> >> > > > > Jose
>> >> > > > >
>> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M. Systems
>> >> > > > > Group This e-mail, including attachments, may include
>> >> > > > > confidential and/or proprietary information, and may be
>> >> > > > > used only by the person or entity to which it is
>> >> > > > > addressed. If the reader of this e-mail is not the
>> >> > > > > intended recipient or his or her authorized agent, the
>> >> > > > > reader is hereby notified that any dissemination,
>> >> > > > > distribution or copying of this e-mail is prohibited. If
>> >> > > > > you have received this e-mail in error, please notify
the
>> >> > > > > sender by replying to this message and delete this
>> >> > > > e-mail immediately.
>> >> > > > >
>> >> > > > >
>> >> > > >
>> >> > > > This e-mail, including attachments, may include
confidential
>> >> > > > and/or proprietary information, and may be used only by
the
>> >> > > > person or entity to which it is addressed. If the reader
of
>> >> > > > this e-mail is not the intended recipient or his or her
>> >> > > > authorized agent, the reader is hereby notified that any
>> >> > > > dissemination, distribution or copying of this e-mail is
>> >> > > > prohibited. If you have received this e-mail in error,
>> >> > > > please notify the sender by replying to this message and
>> >> > > > delete this
>> >> > > e-mail immediately.
>> >> > > >
>> >> > > >
>> >> > > >
>> >> > >
>> >> > > This e-mail, including attachments, may include confidential
>> >> > > and/or proprietary information, and may be used only by the
>> >> > > person or entity to which it is addressed. If the reader of
>> >> > > this e-mail is not the intended recipient or his or her
>> >> > > authorized agent, the reader is hereby notified that any
>> >> > > dissemination, distribution or copying of this e-mail is
>> >> > > prohibited. If you have received this e-mail in error,
please
>> >> > > notify the sender by replying to this message and delete
this
>> >> > e-mail immediately.
>> >> > >
>> >> > >
>> >> >
>> >> > This e-mail, including attachments, may include confidential
>> >> > and/or proprietary information, and may be used only by the
>> >> > person or entity to which it is addressed. If the reader of
this
>> >> > e-mail is not the intended recipient or his or her authorized
>> >> > agent, the reader is hereby notified that any dissemination,
>> >> > distribution or copying of this e-mail is prohibited. If you
>> >> > have received this e-mail in error, please notify the sender
by
>> >> > replying to this message and delete this
>> >> e-mail immediately.
>> >> >
>> >> >
>> >> >
>> >>
>> >> This e-mail, including attachments, may include confidential
>> >> and/or proprietary information, and may be used only by the
person
>> >> or entity to which it is addressed. If the reader of this e-mail
>> >> is not the intended recipient or his or her authorized agent,
the
>> >> reader is hereby notified that any dissemination, distribution
or
>> >> copying of this e-mail is prohibited. If you have received this
>> >> e-mail in error, please notify the sender by replying to this
>> >> message and delete this
>> e-mail immediately.
>> >>
>> >>
>> >>
>> >
>>
>> This e-mail, including attachments, may include confidential and/or
>> proprietary information, and may be used only by the person or
entity
>> to which it is addressed. If the reader of this e-mail is not the
>> intended recipient or his or her authorized agent, the reader is
>> hereby notified that any dissemination, distribution or copying of
>> this e-mail is prohibited. If you have received this e-mail in
error,
>> please notify the sender by replying to this message and delete
this e-mail immediately.
>>
>>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Fri Nov 06 09:49:16 2015
Jose,
Sure, I can help with you that, but I'll need to know the "target"
grid
definition. You said you reran with an NMMB grid with dx = dy. And
you
sent images of the output, but I don't think you sent me a data file.
If you send me that GRIB file, I can see the desired grid info and
send you
a copygb command that'll put the RTMA data onto that grid.
Thanks,
John
On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> We want to use MET, but we need to get the regridding right. I was
> wondering if you know how to regrid the RTMA file to match an NMMB
grid by
> using UPP (as in previous MET versions)?
>
> Any help will be greatly appreciated,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, November 05, 2015 6:31 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Just wanted to follow up. Do you still see a shifting problem in
the
> regridding? And is this a showstopper for you guys? Is adding
support for
> differing dx/dy as a patch for 5.1 a high priority or can we delay
until
> the next release?
>
> Thanks,
> John
>
> On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway <johnhg at ucar.edu>
> wrote:
>
> > Jose,
> >
> > Are you able to send me the GRIB files as well by posting it to
our
> > anonymous ftp site?
> >
> > I took a look at the images you sent, but the skew you mention
isn't
> > jumping out at me. Before I saw that the values that should be
over
> > the great lakes were shifted north by a large extent. What should
I
> > be looking for now?
> >
> > Rather than diffing the NMM and RTMA values, wouldn't it make more
> > sense to...
> >
> > (1) Run the RTMA through regrid_data_plane to regrid it
> > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > (3) Run plot_data_plane to plot the output from regrid_data_plane
> >
> > Compare images from the input/output to make sure the values are
all
> > in the right places.
> >
> > Thanks,
> > John
> >
> >
> > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >>
> >> Dear John,
> >>
> >> I just ran a simulation in which DX=DY, and the skewing was much
> >> smaller, but it's still observed.
> >>
> >> Best Regards,
> >> Jose
> >>
> >>
> >> -----Original Message-----
> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> Sent: Wednesday, November 04, 2015 4:28 PM
> >> To: Jose Garcia-Rivera
> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >>
> >> Jose,
> >>
> >> We're still trying to decide how to address this issue. We
> >> definitely plan on adding support for grids where dx != dy, but
the
> >> question is whether we should add it as a patch for met-5.1 or
just
> >> include it in the next release of the code.
> >>
> >> The enhancement will touch several files, and we worry about the
> >> potential for unintended consequences of the changes.
> >>
> >> Is supporting lambert conformal and polar stereographic grids
where
> >> dx != dy a requirement for your work? Is this shortcoming
standing
> >> in your way of using MET?
> >>
> >> Thanks,
> >> John
> >>
> >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
<johnhg at ucar.edu>
> >> wrote:
> >>
> >> > Jose, yes, running with dx = dy would solve the problem.
> >> >
> >> > John
> >> >
> >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> >> > met_help at ucar.edu> wrote:
> >> >
> >> >>
> >> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> >> >>
> >> >> Dear John,
> >> >>
> >> >> Thank you for your insights. One last question: If we run the
> >> >> model with dx = dy, should that solve the problem?
> >> >>
> >> >> Best Regards,
> >> >> Jose
> >> >>
> >> >> -----Original Message-----
> >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> >> >> To: Jose Garcia-Rivera
> >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >>
> >> >> Jose,
> >> >>
> >> >> Yes, that's correct. Not supporting differing Lambert
Conformal
> >> >> Dx and Dy values is a short-coming in MET. However, we're
> >> >> investigating the impact/cost of adding support for it. We
may be
> >> >> able to add support by patching the build.
> >> >>
> >> >> I'll write you back once we figure out how to proceed.
> >> >>
> >> >> Thanks,
> >> >> John
> >> >>
> >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> >> >> met_help at ucar.edu> wrote:
> >> >>
> >> >> >
> >> >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >> >> >
> >> >> > Hi John,
> >> >> >
> >> >> > Thanks so much for the information. So, that means that I'll
> >> >> > have the same problem when running any other analyses with
this
> >> >> > particular grid
> >> >> using MET?
> >> >> >
> >> >> > Best Regards,
> >> >> > Jose
> >> >> >
> >> >> > -----Original Message-----
> >> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> >> >> > To: Jose Garcia-Rivera
> >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> >
> >> >> > Jose,
> >> >> >
> >> >> > I believe I've figured out the problem. The grid code in
MET
> >> >> > currently only supports Lambert Conformal grids where the
grid
> >> >> > spacing is the same in the x and y directions.
> >> >> >
> >> >> > Your grid is spaced 3.317km in the x-direction and 2.97km in
the
> >> >> > y-direction. The MET code is reading the Dx value and using
it
> >> >> > for both x and y. This causes the skewing of the data to
the
> >> >> > north that
> >> >> you're seeing.
> >> >> >
> >> >> > Really, there are 2 problems here.
> >> >> >
> >> >> > First, if we're requiring that Dx and Dy be the same for
Lambert
> >> >> > Conformal grids, MET should error out with a descriptive
error
> >> >> > message when they are not.
> >> >> >
> >> >> > Second, obviously it would be preferable for MET to support
> >> >> > Lambert Conformal grids where Dx and Dy are not the same!
I'll
> >> >> > ask around here for more information on that.
> >> >> >
> >> >> > Thanks for paying close attention to your data and
uncovering
> >> >> > this shortcoming in MET.
> >> >> >
> >> >> > John
> >> >> >
> >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> >> >> > met_help at ucar.edu
> >> >> > > wrote:
> >> >> >
> >> >> > >
> >> >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> >> > > >
> >> >> > >
> >> >> > > Hi John,
> >> >> > >
> >> >> > > Thanks for all the help. I want to add this pdf, showing
> >> >> > > something odd with the regridding. The domains did not
match
> >> >> > > perfectly on the north side, it's like a skewing effect.
The
> >> >> > > south sides match perfectly, but note the differences in
the
> >> >> > > regridded RTMA file and the model file in Michigan, for
example.
> >> >> > > Maybe that has something to
> >> >> do with it.
> >> >> > >
> >> >> > > Best Regards,
> >> >> > > Jose
> >> >> > >
> >> >> > > -----Original Message-----
> >> >> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> >> >> > > To: Jose Garcia-Rivera
> >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> > >
> >> >> > > Jose,
> >> >> > >
> >> >> > > I believe the problem is in our handling the domain on
which
> >> >> > > your model data resides. I ran the following commands to
plot
> >> >> > > temperature from your model data file:
> >> >> > >
> >> >> > > met-5.1/bin/plot_data_plane \
> >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> >> >> > > 'name="TMP"; level="Z0";'
> >> >> > >
> >> >> > > convert -rotate 90 -background white -flatten \
> >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> >> >> > >
> >> >> > > The resulting image is attached and shows the shift you
> described.
> >> >> > >
> >> >> > > So it isn't the regridding step that's the problem... its
our
> >> >> > > handling of that domain. I plotted the same data using
IDV
> >> >> > > and it puts things in the correct place (also attached).
> >> >> > >
> >> >> > > One issue is that the 2-m temperature data in GRIB record
296
> >> >> > > contains a constant value of 273. The surface temperature
> >> >> > > data in GRIB record
> >> >> > > 276 contains more reasonable values. So that's the field
I'm
> >> >> > investigating.
> >> >> > > But you might want to check why GRIB record 296 looks so
odd.
> >> >> > >
> >> >> > > I need to head out of the office, but will investigate
further
> >> >> tomorrow.
> >> >> > >
> >> >> > > Thanks,
> >> >> > > John
> >> >> > >
> >> >> > >
> >> >> > >
> >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT
<
> >> >> > > met_help at ucar.edu
> >> >> > > > wrote:
> >> >> > >
> >> >> > > >
> >> >> > > > <URL:
> >> >> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> >> >> > > >
> >> >> > > > Dear John,
> >> >> > > >
> >> >> > > > Thank you so much for the quick reply.
> >> >> > > >
> >> >> > > > I followed the instructions, and ftp'ed the files to
> >> >> > > "garcia-rivera_data".
> >> >> > > >
> >> >> > > > The nmbprs file is the .grb NMMB output file.
> >> >> > > > The rtma file is the observational.
> >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
outputs
> >> >> > > > in the same grid, and produces the plot.
> >> >> > > >
> >> >> > > > Any help will be greatly appreciated.
> >> >> > > >
> >> >> > > > Best Regards,
> >> >> > > > Jose
> >> >> > > >
> >> >> > > > -----Original Message-----
> >> >> > > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> >> >> > > > To: Jose Garcia-Rivera
> >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> > > >
> >> >> > > > Jose,
> >> >> > > >
> >> >> > > > Thanks for letting us know about the issue you're
seeing.
> >> >> > > > I'd like to replicate the behavior your describe using
the
> >> >> > > > exact data
> >> >> you used.
> >> >> > > > Once I can replicate it, I'll be able to investigate
further.
> >> >> > > >
> >> >> > > > Are you able to send me the following 2 files you used:
> >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> >> >> > > >
> >> >> > > > You can post them to our anonymous ftp site, following
these
> >> >> > > instructions:
> >> >> > > >
> >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >> >> > > >
> >> >> > > > Please write back to let me know when you've posted the
data
> >> >> > > > and I'll go grab it.
> >> >> > > >
> >> >> > > > Thanks,
> >> >> > > > John Halley Gotway
> >> >> > > > met_help at ucar.edu
> >> >> > > >
> >> >> > > >
> >> >> > > >
> >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via
RT <
> >> >> > > > met_help at ucar.edu
> >> >> > > > > wrote:
> >> >> > > >
> >> >> > > > >
> >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> >> >> > > > > Queue: met_help
> >> >> > > > > Subject: Regridding problem
> >> >> > > > > Owner: Nobody
> >> >> > > > > Requestors: GarciaJ at imsg.com
> >> >> > > > > Status: new
> >> >> > > > > Ticket <URL:
> >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> >> > > > > >
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > Hello,
> >> >> > > > >
> >> >> > > > > I'm a new MET user, and I'm attempting to regrid RTMA
data
> >> >> > > > > to an NMMB output grid. I used the following tool:
> >> >> > > > >
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
level="Z2";'
> >> >> > > > > -method BILIN -width
> >> >> > > > > 2 -v 1
> >> >> > > > >
> >> >> > > > > which worked. But, when I subtracted the fields from
one
> >> >> > > > > dataset to the other, there was a clear offset on the
values.
> >> >> > > > > So, the regridding was done indeed, but the gridpoint
> >> >> > > > > values did not seem to be assigned correctly with
respect
> >> >> > > > > to their
> >> counterparts.
> >> >> > > > >
> >> >> > > > > I attached the plot, so you can see the offset
problem.
> >> >> > > > > Any help will be greatly appreciated.
> >> >> > > > >
> >> >> > > > > Best Regards,
> >> >> > > > > Jose
> >> >> > > > >
> >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M. Systems
> >> >> > > > > Group This e-mail, including attachments, may include
> >> >> > > > > confidential and/or proprietary information, and may
be
> >> >> > > > > used only by the person or entity to which it is
> >> >> > > > > addressed. If the reader of this e-mail is not the
> >> >> > > > > intended recipient or his or her authorized agent, the
> >> >> > > > > reader is hereby notified that any dissemination,
> >> >> > > > > distribution or copying of this e-mail is prohibited.
If
> >> >> > > > > you have received this e-mail in error, please notify
the
> >> >> > > > > sender by replying to this message and delete this
> >> >> > > > e-mail immediately.
> >> >> > > > >
> >> >> > > > >
> >> >> > > >
> >> >> > > > This e-mail, including attachments, may include
confidential
> >> >> > > > and/or proprietary information, and may be used only by
the
> >> >> > > > person or entity to which it is addressed. If the reader
of
> >> >> > > > this e-mail is not the intended recipient or his or her
> >> >> > > > authorized agent, the reader is hereby notified that any
> >> >> > > > dissemination, distribution or copying of this e-mail is
> >> >> > > > prohibited. If you have received this e-mail in error,
> >> >> > > > please notify the sender by replying to this message and
> >> >> > > > delete this
> >> >> > > e-mail immediately.
> >> >> > > >
> >> >> > > >
> >> >> > > >
> >> >> > >
> >> >> > > This e-mail, including attachments, may include
confidential
> >> >> > > and/or proprietary information, and may be used only by
the
> >> >> > > person or entity to which it is addressed. If the reader
of
> >> >> > > this e-mail is not the intended recipient or his or her
> >> >> > > authorized agent, the reader is hereby notified that any
> >> >> > > dissemination, distribution or copying of this e-mail is
> >> >> > > prohibited. If you have received this e-mail in error,
please
> >> >> > > notify the sender by replying to this message and delete
this
> >> >> > e-mail immediately.
> >> >> > >
> >> >> > >
> >> >> >
> >> >> > This e-mail, including attachments, may include confidential
> >> >> > and/or proprietary information, and may be used only by the
> >> >> > person or entity to which it is addressed. If the reader of
this
> >> >> > e-mail is not the intended recipient or his or her
authorized
> >> >> > agent, the reader is hereby notified that any dissemination,
> >> >> > distribution or copying of this e-mail is prohibited. If you
> >> >> > have received this e-mail in error, please notify the sender
by
> >> >> > replying to this message and delete this
> >> >> e-mail immediately.
> >> >> >
> >> >> >
> >> >> >
> >> >>
> >> >> This e-mail, including attachments, may include confidential
> >> >> and/or proprietary information, and may be used only by the
person
> >> >> or entity to which it is addressed. If the reader of this e-
mail
> >> >> is not the intended recipient or his or her authorized agent,
the
> >> >> reader is hereby notified that any dissemination, distribution
or
> >> >> copying of this e-mail is prohibited. If you have received
this
> >> >> e-mail in error, please notify the sender by replying to this
> >> >> message and delete this
> >> e-mail immediately.
> >> >>
> >> >>
> >> >>
> >> >
> >>
> >> This e-mail, including attachments, may include confidential
and/or
> >> proprietary information, and may be used only by the person or
entity
> >> to which it is addressed. If the reader of this e-mail is not the
> >> intended recipient or his or her authorized agent, the reader is
> >> hereby notified that any dissemination, distribution or copying
of
> >> this e-mail is prohibited. If you have received this e-mail in
error,
> >> please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >>
> >>
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Fri Nov 06 09:51:44 2015
Dear John,
Thanks for all the help. The target grid would be the original grid I
sent you.
This is the same grid in copy_gridnav.txt:
255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64 38947
38947 44113 -69980
Best Regards,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Friday, November 06, 2015 11:49 AM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
Sure, I can help with you that, but I'll need to know the "target"
grid definition. You said you reran with an NMMB grid with dx = dy.
And you sent images of the output, but I don't think you sent me a
data file.
If you send me that GRIB file, I can see the desired grid info and
send you a copygb command that'll put the RTMA data onto that grid.
Thanks,
John
On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> We want to use MET, but we need to get the regridding right. I was
> wondering if you know how to regrid the RTMA file to match an NMMB
> grid by using UPP (as in previous MET versions)?
>
> Any help will be greatly appreciated,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, November 05, 2015 6:31 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Just wanted to follow up. Do you still see a shifting problem in
the
> regridding? And is this a showstopper for you guys? Is adding
> support for differing dx/dy as a patch for 5.1 a high priority or
can
> we delay until the next release?
>
> Thanks,
> John
>
> On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway <johnhg at ucar.edu>
> wrote:
>
> > Jose,
> >
> > Are you able to send me the GRIB files as well by posting it to
our
> > anonymous ftp site?
> >
> > I took a look at the images you sent, but the skew you mention
isn't
> > jumping out at me. Before I saw that the values that should be
over
> > the great lakes were shifted north by a large extent. What should
I
> > be looking for now?
> >
> > Rather than diffing the NMM and RTMA values, wouldn't it make more
> > sense to...
> >
> > (1) Run the RTMA through regrid_data_plane to regrid it
> > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > (3) Run plot_data_plane to plot the output from regrid_data_plane
> >
> > Compare images from the input/output to make sure the values are
all
> > in the right places.
> >
> > Thanks,
> > John
> >
> >
> > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >>
> >> Dear John,
> >>
> >> I just ran a simulation in which DX=DY, and the skewing was much
> >> smaller, but it's still observed.
> >>
> >> Best Regards,
> >> Jose
> >>
> >>
> >> -----Original Message-----
> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> Sent: Wednesday, November 04, 2015 4:28 PM
> >> To: Jose Garcia-Rivera
> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >>
> >> Jose,
> >>
> >> We're still trying to decide how to address this issue. We
> >> definitely plan on adding support for grids where dx != dy, but
the
> >> question is whether we should add it as a patch for met-5.1 or
just
> >> include it in the next release of the code.
> >>
> >> The enhancement will touch several files, and we worry about the
> >> potential for unintended consequences of the changes.
> >>
> >> Is supporting lambert conformal and polar stereographic grids
where
> >> dx != dy a requirement for your work? Is this shortcoming
standing
> >> in your way of using MET?
> >>
> >> Thanks,
> >> John
> >>
> >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> >> <johnhg at ucar.edu>
> >> wrote:
> >>
> >> > Jose, yes, running with dx = dy would solve the problem.
> >> >
> >> > John
> >> >
> >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> >> > met_help at ucar.edu> wrote:
> >> >
> >> >>
> >> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> >> >>
> >> >> Dear John,
> >> >>
> >> >> Thank you for your insights. One last question: If we run the
> >> >> model with dx = dy, should that solve the problem?
> >> >>
> >> >> Best Regards,
> >> >> Jose
> >> >>
> >> >> -----Original Message-----
> >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> >> >> To: Jose Garcia-Rivera
> >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >>
> >> >> Jose,
> >> >>
> >> >> Yes, that's correct. Not supporting differing Lambert
Conformal
> >> >> Dx and Dy values is a short-coming in MET. However, we're
> >> >> investigating the impact/cost of adding support for it. We
may
> >> >> be able to add support by patching the build.
> >> >>
> >> >> I'll write you back once we figure out how to proceed.
> >> >>
> >> >> Thanks,
> >> >> John
> >> >>
> >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> >> >> met_help at ucar.edu> wrote:
> >> >>
> >> >> >
> >> >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> >> > >
> >> >> >
> >> >> > Hi John,
> >> >> >
> >> >> > Thanks so much for the information. So, that means that I'll
> >> >> > have the same problem when running any other analyses with
> >> >> > this particular grid
> >> >> using MET?
> >> >> >
> >> >> > Best Regards,
> >> >> > Jose
> >> >> >
> >> >> > -----Original Message-----
> >> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> >> >> > To: Jose Garcia-Rivera
> >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> >
> >> >> > Jose,
> >> >> >
> >> >> > I believe I've figured out the problem. The grid code in
MET
> >> >> > currently only supports Lambert Conformal grids where the
grid
> >> >> > spacing is the same in the x and y directions.
> >> >> >
> >> >> > Your grid is spaced 3.317km in the x-direction and 2.97km in
> >> >> > the y-direction. The MET code is reading the Dx value and
> >> >> > using it for both x and y. This causes the skewing of the
> >> >> > data to the north that
> >> >> you're seeing.
> >> >> >
> >> >> > Really, there are 2 problems here.
> >> >> >
> >> >> > First, if we're requiring that Dx and Dy be the same for
> >> >> > Lambert Conformal grids, MET should error out with a
> >> >> > descriptive error message when they are not.
> >> >> >
> >> >> > Second, obviously it would be preferable for MET to support
> >> >> > Lambert Conformal grids where Dx and Dy are not the same!
> >> >> > I'll ask around here for more information on that.
> >> >> >
> >> >> > Thanks for paying close attention to your data and
uncovering
> >> >> > this shortcoming in MET.
> >> >> >
> >> >> > John
> >> >> >
> >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT <
> >> >> > met_help at ucar.edu
> >> >> > > wrote:
> >> >> >
> >> >> > >
> >> >> > > <URL:
> >> >> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> >> > > >
> >> >> > >
> >> >> > > Hi John,
> >> >> > >
> >> >> > > Thanks for all the help. I want to add this pdf, showing
> >> >> > > something odd with the regridding. The domains did not
match
> >> >> > > perfectly on the north side, it's like a skewing effect.
The
> >> >> > > south sides match perfectly, but note the differences in
the
> >> >> > > regridded RTMA file and the model file in Michigan, for
example.
> >> >> > > Maybe that has something to
> >> >> do with it.
> >> >> > >
> >> >> > > Best Regards,
> >> >> > > Jose
> >> >> > >
> >> >> > > -----Original Message-----
> >> >> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> >> >> > > To: Jose Garcia-Rivera
> >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> > >
> >> >> > > Jose,
> >> >> > >
> >> >> > > I believe the problem is in our handling the domain on
which
> >> >> > > your model data resides. I ran the following commands to
> >> >> > > plot temperature from your model data file:
> >> >> > >
> >> >> > > met-5.1/bin/plot_data_plane \
> >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> >> >> > > 'name="TMP"; level="Z0";'
> >> >> > >
> >> >> > > convert -rotate 90 -background white -flatten \
> >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> >> >> > >
> >> >> > > The resulting image is attached and shows the shift you
> described.
> >> >> > >
> >> >> > > So it isn't the regridding step that's the problem... its
> >> >> > > our handling of that domain. I plotted the same data
using
> >> >> > > IDV and it puts things in the correct place (also
attached).
> >> >> > >
> >> >> > > One issue is that the 2-m temperature data in GRIB record
> >> >> > > 296 contains a constant value of 273. The surface
> >> >> > > temperature data in GRIB record
> >> >> > > 276 contains more reasonable values. So that's the field
> >> >> > > I'm
> >> >> > investigating.
> >> >> > > But you might want to check why GRIB record 296 looks so
odd.
> >> >> > >
> >> >> > > I need to head out of the office, but will investigate
> >> >> > > further
> >> >> tomorrow.
> >> >> > >
> >> >> > > Thanks,
> >> >> > > John
> >> >> > >
> >> >> > >
> >> >> > >
> >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via RT
<
> >> >> > > met_help at ucar.edu
> >> >> > > > wrote:
> >> >> > >
> >> >> > > >
> >> >> > > > <URL:
> >> >> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> >> >> > > >
> >> >> > > > Dear John,
> >> >> > > >
> >> >> > > > Thank you so much for the quick reply.
> >> >> > > >
> >> >> > > > I followed the instructions, and ftp'ed the files to
> >> >> > > "garcia-rivera_data".
> >> >> > > >
> >> >> > > > The nmbprs file is the .grb NMMB output file.
> >> >> > > > The rtma file is the observational.
> >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
outputs
> >> >> > > > in the same grid, and produces the plot.
> >> >> > > >
> >> >> > > > Any help will be greatly appreciated.
> >> >> > > >
> >> >> > > > Best Regards,
> >> >> > > > Jose
> >> >> > > >
> >> >> > > > -----Original Message-----
> >> >> > > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> >> >> > > > To: Jose Garcia-Rivera
> >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >> >> > > >
> >> >> > > > Jose,
> >> >> > > >
> >> >> > > > Thanks for letting us know about the issue you're
seeing.
> >> >> > > > I'd like to replicate the behavior your describe using
the
> >> >> > > > exact data
> >> >> you used.
> >> >> > > > Once I can replicate it, I'll be able to investigate
further.
> >> >> > > >
> >> >> > > > Are you able to send me the following 2 files you used:
> >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> >> >> > > >
> >> >> > > > You can post them to our anonymous ftp site, following
> >> >> > > > these
> >> >> > > instructions:
> >> >> > > >
> >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >> >> > > >
> >> >> > > > Please write back to let me know when you've posted the
> >> >> > > > data and I'll go grab it.
> >> >> > > >
> >> >> > > > Thanks,
> >> >> > > > John Halley Gotway
> >> >> > > > met_help at ucar.edu
> >> >> > > >
> >> >> > > >
> >> >> > > >
> >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via
RT
> >> >> > > > < met_help at ucar.edu
> >> >> > > > > wrote:
> >> >> > > >
> >> >> > > > >
> >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> >> >> > > > > Queue: met_help
> >> >> > > > > Subject: Regridding problem
> >> >> > > > > Owner: Nobody
> >> >> > > > > Requestors: GarciaJ at imsg.com
> >> >> > > > > Status: new
> >> >> > > > > Ticket <URL:
> >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> >> >> > > > > >
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > Hello,
> >> >> > > > >
> >> >> > > > > I'm a new MET user, and I'm attempting to regrid RTMA
> >> >> > > > > data to an NMMB output grid. I used the following
tool:
> >> >> > > > >
> >> >> > > > >
> >> >> > > > >
> >> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
level="Z2";'
> >> >> > > > > -method BILIN -width
> >> >> > > > > 2 -v 1
> >> >> > > > >
> >> >> > > > > which worked. But, when I subtracted the fields from
one
> >> >> > > > > dataset to the other, there was a clear offset on the
values.
> >> >> > > > > So, the regridding was done indeed, but the gridpoint
> >> >> > > > > values did not seem to be assigned correctly with
> >> >> > > > > respect to their
> >> counterparts.
> >> >> > > > >
> >> >> > > > > I attached the plot, so you can see the offset
problem.
> >> >> > > > > Any help will be greatly appreciated.
> >> >> > > > >
> >> >> > > > > Best Regards,
> >> >> > > > > Jose
> >> >> > > > >
> >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M. Systems
> >> >> > > > > Group This e-mail, including attachments, may include
> >> >> > > > > confidential and/or proprietary information, and may
be
> >> >> > > > > used only by the person or entity to which it is
> >> >> > > > > addressed. If the reader of this e-mail is not the
> >> >> > > > > intended recipient or his or her authorized agent, the
> >> >> > > > > reader is hereby notified that any dissemination,
> >> >> > > > > distribution or copying of this e-mail is prohibited.
If
> >> >> > > > > you have received this e-mail in error, please notify
> >> >> > > > > the sender by replying to this message and delete this
> >> >> > > > e-mail immediately.
> >> >> > > > >
> >> >> > > > >
> >> >> > > >
> >> >> > > > This e-mail, including attachments, may include
> >> >> > > > confidential and/or proprietary information, and may be
> >> >> > > > used only by the person or entity to which it is
> >> >> > > > addressed. If the reader of this e-mail is not the
> >> >> > > > intended recipient or his or her authorized agent, the
> >> >> > > > reader is hereby notified that any dissemination,
> >> >> > > > distribution or copying of this e-mail is prohibited. If
> >> >> > > > you have received this e-mail in error, please notify
the
> >> >> > > > sender by replying to this message and delete this
> >> >> > > e-mail immediately.
> >> >> > > >
> >> >> > > >
> >> >> > > >
> >> >> > >
> >> >> > > This e-mail, including attachments, may include
confidential
> >> >> > > and/or proprietary information, and may be used only by
the
> >> >> > > person or entity to which it is addressed. If the reader
of
> >> >> > > this e-mail is not the intended recipient or his or her
> >> >> > > authorized agent, the reader is hereby notified that any
> >> >> > > dissemination, distribution or copying of this e-mail is
> >> >> > > prohibited. If you have received this e-mail in error,
> >> >> > > please notify the sender by replying to this message and
> >> >> > > delete this
> >> >> > e-mail immediately.
> >> >> > >
> >> >> > >
> >> >> >
> >> >> > This e-mail, including attachments, may include confidential
> >> >> > and/or proprietary information, and may be used only by the
> >> >> > person or entity to which it is addressed. If the reader of
> >> >> > this e-mail is not the intended recipient or his or her
> >> >> > authorized agent, the reader is hereby notified that any
> >> >> > dissemination, distribution or copying of this e-mail is
> >> >> > prohibited. If you have received this e-mail in error,
please
> >> >> > notify the sender by replying to this message and delete
this
> >> >> e-mail immediately.
> >> >> >
> >> >> >
> >> >> >
> >> >>
> >> >> This e-mail, including attachments, may include confidential
> >> >> and/or proprietary information, and may be used only by the
> >> >> person or entity to which it is addressed. If the reader of
this
> >> >> e-mail is not the intended recipient or his or her authorized
> >> >> agent, the reader is hereby notified that any dissemination,
> >> >> distribution or copying of this e-mail is prohibited. If you
> >> >> have received this e-mail in error, please notify the sender
by
> >> >> replying to this message and delete this
> >> e-mail immediately.
> >> >>
> >> >>
> >> >>
> >> >
> >>
> >> This e-mail, including attachments, may include confidential
and/or
> >> proprietary information, and may be used only by the person or
> >> entity to which it is addressed. If the reader of this e-mail is
> >> not the intended recipient or his or her authorized agent, the
> >> reader is hereby notified that any dissemination, distribution or
> >> copying of this e-mail is prohibited. If you have received this
> >> e-mail in error, please notify the sender by replying to this
> >> message and delete this
> e-mail immediately.
> >>
> >>
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Fri Nov 06 10:14:53 2015
Jose,
The regridding tool is named copygb, but it only works on GRIB1 files.
The
RTMA data is in GRIB2. Here's one way of doing it:
# Convert RTMA from GRIB2 to GRIB1
cnvgrib -g21
rtma2p5.t00z.2dvaranl_ndfd.grb2rtma2p5.t00z.2dvaranl_ndfd.grb
# Run copygb to regrid GRIB1 file
copygb \
-xg'255 3 400 400 33302 -84951 8 -77459 3317 2970 0 64 38947 38947'
\
rtma2p5.t00z.2dvaranl_ndfd.grb
rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
FYI, to do this I just followed an example from the MET online
tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run3.php
Run wgrib on the output to verify that it's on the intended grid:
wgrib -d 1 -V rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
...
Lambert Conf: Lat1 33.302000 Lon1 -84.951000 Lov -77.459000
Latin1 38.947000 Latin2 38.947000 LatSP 0.000000 LonSP 0.000000
North Pole (400 x 400) Dx 3.317000 Dy 2.970000 scan 64 mode 8
...
This data can now be run through MET. For example, running it through
grid-stat will give you statistics computed point-by-point and those
should
be correct. But the underlying projection is still wrong because dx
!=
dy. So any plots created by MET will show that skew to the north.
It may also be possible to do the regridding using the "wgrib2" tool
directly on the GRIB2 file. But I haven't done that as much. You
could
look into that if you want to skip the GRIB2 to GRIB1 conversion step.
Thanks,
John
On Fri, Nov 6, 2015 at 9:51 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thanks for all the help. The target grid would be the original grid
I sent
> you.
>
> This is the same grid in copy_gridnav.txt:
>
> 255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64 38947
38947
> 44113 -69980
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, November 06, 2015 11:49 AM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Sure, I can help with you that, but I'll need to know the "target"
grid
> definition. You said you reran with an NMMB grid with dx = dy. And
you
> sent images of the output, but I don't think you sent me a data
file.
>
> If you send me that GRIB file, I can see the desired grid info and
send
> you a copygb command that'll put the RTMA data onto that grid.
>
> Thanks,
> John
>
> On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > We want to use MET, but we need to get the regridding right. I was
> > wondering if you know how to regrid the RTMA file to match an NMMB
> > grid by using UPP (as in previous MET versions)?
> >
> > Any help will be greatly appreciated,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Thursday, November 05, 2015 6:31 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Just wanted to follow up. Do you still see a shifting problem in
the
> > regridding? And is this a showstopper for you guys? Is adding
> > support for differing dx/dy as a patch for 5.1 a high priority or
can
> > we delay until the next release?
> >
> > Thanks,
> > John
> >
> > On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway
<johnhg at ucar.edu>
> > wrote:
> >
> > > Jose,
> > >
> > > Are you able to send me the GRIB files as well by posting it to
our
> > > anonymous ftp site?
> > >
> > > I took a look at the images you sent, but the skew you mention
isn't
> > > jumping out at me. Before I saw that the values that should be
over
> > > the great lakes were shifted north by a large extent. What
should I
> > > be looking for now?
> > >
> > > Rather than diffing the NMM and RTMA values, wouldn't it make
more
> > > sense to...
> > >
> > > (1) Run the RTMA through regrid_data_plane to regrid it
> > > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > > (3) Run plot_data_plane to plot the output from
regrid_data_plane
> > >
> > > Compare images from the input/output to make sure the values are
all
> > > in the right places.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >>
> > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >>
> > >> Dear John,
> > >>
> > >> I just ran a simulation in which DX=DY, and the skewing was
much
> > >> smaller, but it's still observed.
> > >>
> > >> Best Regards,
> > >> Jose
> > >>
> > >>
> > >> -----Original Message-----
> > >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> Sent: Wednesday, November 04, 2015 4:28 PM
> > >> To: Jose Garcia-Rivera
> > >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >>
> > >> Jose,
> > >>
> > >> We're still trying to decide how to address this issue. We
> > >> definitely plan on adding support for grids where dx != dy, but
the
> > >> question is whether we should add it as a patch for met-5.1 or
just
> > >> include it in the next release of the code.
> > >>
> > >> The enhancement will touch several files, and we worry about
the
> > >> potential for unintended consequences of the changes.
> > >>
> > >> Is supporting lambert conformal and polar stereographic grids
where
> > >> dx != dy a requirement for your work? Is this shortcoming
standing
> > >> in your way of using MET?
> > >>
> > >> Thanks,
> > >> John
> > >>
> > >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> > >> <johnhg at ucar.edu>
> > >> wrote:
> > >>
> > >> > Jose, yes, running with dx = dy would solve the problem.
> > >> >
> > >> > John
> > >> >
> > >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> > >> > met_help at ucar.edu> wrote:
> > >> >
> > >> >>
> > >> >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >> >>
> > >> >> Dear John,
> > >> >>
> > >> >> Thank you for your insights. One last question: If we run
the
> > >> >> model with dx = dy, should that solve the problem?
> > >> >>
> > >> >> Best Regards,
> > >> >> Jose
> > >> >>
> > >> >> -----Original Message-----
> > >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> > >> >> To: Jose Garcia-Rivera
> > >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >>
> > >> >> Jose,
> > >> >>
> > >> >> Yes, that's correct. Not supporting differing Lambert
Conformal
> > >> >> Dx and Dy values is a short-coming in MET. However, we're
> > >> >> investigating the impact/cost of adding support for it. We
may
> > >> >> be able to add support by patching the build.
> > >> >>
> > >> >> I'll write you back once we figure out how to proceed.
> > >> >>
> > >> >> Thanks,
> > >> >> John
> > >> >>
> > >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> > >> >> met_help at ucar.edu> wrote:
> > >> >>
> > >> >> >
> > >> >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > >
> > >> >> >
> > >> >> > Hi John,
> > >> >> >
> > >> >> > Thanks so much for the information. So, that means that
I'll
> > >> >> > have the same problem when running any other analyses with
> > >> >> > this particular grid
> > >> >> using MET?
> > >> >> >
> > >> >> > Best Regards,
> > >> >> > Jose
> > >> >> >
> > >> >> > -----Original Message-----
> > >> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> > >> >> > To: Jose Garcia-Rivera
> > >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >> >
> > >> >> > Jose,
> > >> >> >
> > >> >> > I believe I've figured out the problem. The grid code in
MET
> > >> >> > currently only supports Lambert Conformal grids where the
grid
> > >> >> > spacing is the same in the x and y directions.
> > >> >> >
> > >> >> > Your grid is spaced 3.317km in the x-direction and 2.97km
in
> > >> >> > the y-direction. The MET code is reading the Dx value and
> > >> >> > using it for both x and y. This causes the skewing of the
> > >> >> > data to the north that
> > >> >> you're seeing.
> > >> >> >
> > >> >> > Really, there are 2 problems here.
> > >> >> >
> > >> >> > First, if we're requiring that Dx and Dy be the same for
> > >> >> > Lambert Conformal grids, MET should error out with a
> > >> >> > descriptive error message when they are not.
> > >> >> >
> > >> >> > Second, obviously it would be preferable for MET to
support
> > >> >> > Lambert Conformal grids where Dx and Dy are not the same!
> > >> >> > I'll ask around here for more information on that.
> > >> >> >
> > >> >> > Thanks for paying close attention to your data and
uncovering
> > >> >> > this shortcoming in MET.
> > >> >> >
> > >> >> > John
> > >> >> >
> > >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT
<
> > >> >> > met_help at ucar.edu
> > >> >> > > wrote:
> > >> >> >
> > >> >> > >
> > >> >> > > <URL:
> > >> >> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > > >
> > >> >> > >
> > >> >> > > Hi John,
> > >> >> > >
> > >> >> > > Thanks for all the help. I want to add this pdf, showing
> > >> >> > > something odd with the regridding. The domains did not
match
> > >> >> > > perfectly on the north side, it's like a skewing effect.
The
> > >> >> > > south sides match perfectly, but note the differences in
the
> > >> >> > > regridded RTMA file and the model file in Michigan, for
> example.
> > >> >> > > Maybe that has something to
> > >> >> do with it.
> > >> >> > >
> > >> >> > > Best Regards,
> > >> >> > > Jose
> > >> >> > >
> > >> >> > > -----Original Message-----
> > >> >> > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > >> >> > > To: Jose Garcia-Rivera
> > >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >> > >
> > >> >> > > Jose,
> > >> >> > >
> > >> >> > > I believe the problem is in our handling the domain on
which
> > >> >> > > your model data resides. I ran the following commands
to
> > >> >> > > plot temperature from your model data file:
> > >> >> > >
> > >> >> > > met-5.1/bin/plot_data_plane \
> > >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > >> >> > > 'name="TMP"; level="Z0";'
> > >> >> > >
> > >> >> > > convert -rotate 90 -background white -flatten \
> > >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > >> >> > >
> > >> >> > > The resulting image is attached and shows the shift you
> > described.
> > >> >> > >
> > >> >> > > So it isn't the regridding step that's the problem...
its
> > >> >> > > our handling of that domain. I plotted the same data
using
> > >> >> > > IDV and it puts things in the correct place (also
attached).
> > >> >> > >
> > >> >> > > One issue is that the 2-m temperature data in GRIB
record
> > >> >> > > 296 contains a constant value of 273. The surface
> > >> >> > > temperature data in GRIB record
> > >> >> > > 276 contains more reasonable values. So that's the
field
> > >> >> > > I'm
> > >> >> > investigating.
> > >> >> > > But you might want to check why GRIB record 296 looks so
odd.
> > >> >> > >
> > >> >> > > I need to head out of the office, but will investigate
> > >> >> > > further
> > >> >> tomorrow.
> > >> >> > >
> > >> >> > > Thanks,
> > >> >> > > John
> > >> >> > >
> > >> >> > >
> > >> >> > >
> > >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via
RT <
> > >> >> > > met_help at ucar.edu
> > >> >> > > > wrote:
> > >> >> > >
> > >> >> > > >
> > >> >> > > > <URL:
> > >> >> > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >> >> > > >
> > >> >> > > > Dear John,
> > >> >> > > >
> > >> >> > > > Thank you so much for the quick reply.
> > >> >> > > >
> > >> >> > > > I followed the instructions, and ftp'ed the files to
> > >> >> > > "garcia-rivera_data".
> > >> >> > > >
> > >> >> > > > The nmbprs file is the .grb NMMB output file.
> > >> >> > > > The rtma file is the observational.
> > >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
outputs
> > >> >> > > > in the same grid, and produces the plot.
> > >> >> > > >
> > >> >> > > > Any help will be greatly appreciated.
> > >> >> > > >
> > >> >> > > > Best Regards,
> > >> >> > > > Jose
> > >> >> > > >
> > >> >> > > > -----Original Message-----
> > >> >> > > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > >> >> > > > To: Jose Garcia-Rivera
> > >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > >> >> > > >
> > >> >> > > > Jose,
> > >> >> > > >
> > >> >> > > > Thanks for letting us know about the issue you're
seeing.
> > >> >> > > > I'd like to replicate the behavior your describe using
the
> > >> >> > > > exact data
> > >> >> you used.
> > >> >> > > > Once I can replicate it, I'll be able to investigate
further.
> > >> >> > > >
> > >> >> > > > Are you able to send me the following 2 files you
used:
> > >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > >> >> > > >
> > >> >> > > > You can post them to our anonymous ftp site, following
> > >> >> > > > these
> > >> >> > > instructions:
> > >> >> > > >
> > >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >> >> > > >
> > >> >> > > > Please write back to let me know when you've posted
the
> > >> >> > > > data and I'll go grab it.
> > >> >> > > >
> > >> >> > > > Thanks,
> > >> >> > > > John Halley Gotway
> > >> >> > > > met_help at ucar.edu
> > >> >> > > >
> > >> >> > > >
> > >> >> > > >
> > >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via
RT
> > >> >> > > > < met_help at ucar.edu
> > >> >> > > > > wrote:
> > >> >> > > >
> > >> >> > > > >
> > >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> > >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > >> >> > > > > Queue: met_help
> > >> >> > > > > Subject: Regridding problem
> > >> >> > > > > Owner: Nobody
> > >> >> > > > > Requestors: GarciaJ at imsg.com
> > >> >> > > > > Status: new
> > >> >> > > > > Ticket <URL:
> > >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > > > > >
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > Hello,
> > >> >> > > > >
> > >> >> > > > > I'm a new MET user, and I'm attempting to regrid
RTMA
> > >> >> > > > > data to an NMMB output grid. I used the following
tool:
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
> level="Z2";'
> > >> >> > > > > -method BILIN -width
> > >> >> > > > > 2 -v 1
> > >> >> > > > >
> > >> >> > > > > which worked. But, when I subtracted the fields from
one
> > >> >> > > > > dataset to the other, there was a clear offset on
the
> values.
> > >> >> > > > > So, the regridding was done indeed, but the
gridpoint
> > >> >> > > > > values did not seem to be assigned correctly with
> > >> >> > > > > respect to their
> > >> counterparts.
> > >> >> > > > >
> > >> >> > > > > I attached the plot, so you can see the offset
problem.
> > >> >> > > > > Any help will be greatly appreciated.
> > >> >> > > > >
> > >> >> > > > > Best Regards,
> > >> >> > > > > Jose
> > >> >> > > > >
> > >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M.
Systems
> > >> >> > > > > Group This e-mail, including attachments, may
include
> > >> >> > > > > confidential and/or proprietary information, and may
be
> > >> >> > > > > used only by the person or entity to which it is
> > >> >> > > > > addressed. If the reader of this e-mail is not the
> > >> >> > > > > intended recipient or his or her authorized agent,
the
> > >> >> > > > > reader is hereby notified that any dissemination,
> > >> >> > > > > distribution or copying of this e-mail is
prohibited. If
> > >> >> > > > > you have received this e-mail in error, please
notify
> > >> >> > > > > the sender by replying to this message and delete
this
> > >> >> > > > e-mail immediately.
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > >
> > >> >> > > > This e-mail, including attachments, may include
> > >> >> > > > confidential and/or proprietary information, and may
be
> > >> >> > > > used only by the person or entity to which it is
> > >> >> > > > addressed. If the reader of this e-mail is not the
> > >> >> > > > intended recipient or his or her authorized agent, the
> > >> >> > > > reader is hereby notified that any dissemination,
> > >> >> > > > distribution or copying of this e-mail is prohibited.
If
> > >> >> > > > you have received this e-mail in error, please notify
the
> > >> >> > > > sender by replying to this message and delete this
> > >> >> > > e-mail immediately.
> > >> >> > > >
> > >> >> > > >
> > >> >> > > >
> > >> >> > >
> > >> >> > > This e-mail, including attachments, may include
confidential
> > >> >> > > and/or proprietary information, and may be used only by
the
> > >> >> > > person or entity to which it is addressed. If the reader
of
> > >> >> > > this e-mail is not the intended recipient or his or her
> > >> >> > > authorized agent, the reader is hereby notified that any
> > >> >> > > dissemination, distribution or copying of this e-mail is
> > >> >> > > prohibited. If you have received this e-mail in error,
> > >> >> > > please notify the sender by replying to this message and
> > >> >> > > delete this
> > >> >> > e-mail immediately.
> > >> >> > >
> > >> >> > >
> > >> >> >
> > >> >> > This e-mail, including attachments, may include
confidential
> > >> >> > and/or proprietary information, and may be used only by
the
> > >> >> > person or entity to which it is addressed. If the reader
of
> > >> >> > this e-mail is not the intended recipient or his or her
> > >> >> > authorized agent, the reader is hereby notified that any
> > >> >> > dissemination, distribution or copying of this e-mail is
> > >> >> > prohibited. If you have received this e-mail in error,
please
> > >> >> > notify the sender by replying to this message and delete
this
> > >> >> e-mail immediately.
> > >> >> >
> > >> >> >
> > >> >> >
> > >> >>
> > >> >> This e-mail, including attachments, may include confidential
> > >> >> and/or proprietary information, and may be used only by the
> > >> >> person or entity to which it is addressed. If the reader of
this
> > >> >> e-mail is not the intended recipient or his or her
authorized
> > >> >> agent, the reader is hereby notified that any dissemination,
> > >> >> distribution or copying of this e-mail is prohibited. If you
> > >> >> have received this e-mail in error, please notify the sender
by
> > >> >> replying to this message and delete this
> > >> e-mail immediately.
> > >> >>
> > >> >>
> > >> >>
> > >> >
> > >>
> > >> This e-mail, including attachments, may include confidential
and/or
> > >> proprietary information, and may be used only by the person or
> > >> entity to which it is addressed. If the reader of this e-mail
is
> > >> not the intended recipient or his or her authorized agent, the
> > >> reader is hereby notified that any dissemination, distribution
or
> > >> copying of this e-mail is prohibited. If you have received this
> > >> e-mail in error, please notify the sender by replying to this
> > >> message and delete this
> > e-mail immediately.
> > >>
> > >>
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Fri Nov 06 11:05:26 2015
Hi John,
Where is the cnvgrib utility?
Thanks,
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Friday, November 06, 2015 12:15 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
The regridding tool is named copygb, but it only works on GRIB1 files.
The RTMA data is in GRIB2. Here's one way of doing it:
# Convert RTMA from GRIB2 to GRIB1
cnvgrib -g21
rtma2p5.t00z.2dvaranl_ndfd.grb2rtma2p5.t00z.2dvaranl_ndfd.grb
# Run copygb to regrid GRIB1 file
copygb \
-xg'255 3 400 400 33302 -84951 8 -77459 3317 2970 0 64 38947 38947'
\
rtma2p5.t00z.2dvaranl_ndfd.grb
rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
FYI, to do this I just followed an example from the MET online
tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run3.php
Run wgrib on the output to verify that it's on the intended grid:
wgrib -d 1 -V rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
...
Lambert Conf: Lat1 33.302000 Lon1 -84.951000 Lov -77.459000
Latin1 38.947000 Latin2 38.947000 LatSP 0.000000 LonSP 0.000000
North Pole (400 x 400) Dx 3.317000 Dy 2.970000 scan 64 mode 8
...
This data can now be run through MET. For example, running it through
grid-stat will give you statistics computed point-by-point and those
should be correct. But the underlying projection is still wrong
because dx != dy. So any plots created by MET will show that skew to
the north.
It may also be possible to do the regridding using the "wgrib2" tool
directly on the GRIB2 file. But I haven't done that as much. You
could look into that if you want to skip the GRIB2 to GRIB1 conversion
step.
Thanks,
John
On Fri, Nov 6, 2015 at 9:51 AM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thanks for all the help. The target grid would be the original grid
I
> sent you.
>
> This is the same grid in copy_gridnav.txt:
>
> 255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64 38947
38947
> 44113 -69980
>
> Best Regards,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, November 06, 2015 11:49 AM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> Sure, I can help with you that, but I'll need to know the "target"
> grid definition. You said you reran with an NMMB grid with dx = dy.
> And you sent images of the output, but I don't think you sent me a
data file.
>
> If you send me that GRIB file, I can see the desired grid info and
> send you a copygb command that'll put the RTMA data onto that grid.
>
> Thanks,
> John
>
> On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > We want to use MET, but we need to get the regridding right. I was
> > wondering if you know how to regrid the RTMA file to match an NMMB
> > grid by using UPP (as in previous MET versions)?
> >
> > Any help will be greatly appreciated, Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Thursday, November 05, 2015 6:31 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Just wanted to follow up. Do you still see a shifting problem in
> > the regridding? And is this a showstopper for you guys? Is
adding
> > support for differing dx/dy as a patch for 5.1 a high priority or
> > can we delay until the next release?
> >
> > Thanks,
> > John
> >
> > On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway
<johnhg at ucar.edu>
> > wrote:
> >
> > > Jose,
> > >
> > > Are you able to send me the GRIB files as well by posting it to
> > > our anonymous ftp site?
> > >
> > > I took a look at the images you sent, but the skew you mention
> > > isn't jumping out at me. Before I saw that the values that
should
> > > be over the great lakes were shifted north by a large extent.
> > > What should I be looking for now?
> > >
> > > Rather than diffing the NMM and RTMA values, wouldn't it make
more
> > > sense to...
> > >
> > > (1) Run the RTMA through regrid_data_plane to regrid it
> > > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > > (3) Run plot_data_plane to plot the output from
regrid_data_plane
> > >
> > > Compare images from the input/output to make sure the values are
> > > all in the right places.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >>
> > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >>
> > >> Dear John,
> > >>
> > >> I just ran a simulation in which DX=DY, and the skewing was
much
> > >> smaller, but it's still observed.
> > >>
> > >> Best Regards,
> > >> Jose
> > >>
> > >>
> > >> -----Original Message-----
> > >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> Sent: Wednesday, November 04, 2015 4:28 PM
> > >> To: Jose Garcia-Rivera
> > >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >>
> > >> Jose,
> > >>
> > >> We're still trying to decide how to address this issue. We
> > >> definitely plan on adding support for grids where dx != dy, but
> > >> the question is whether we should add it as a patch for met-5.1
> > >> or just include it in the next release of the code.
> > >>
> > >> The enhancement will touch several files, and we worry about
the
> > >> potential for unintended consequences of the changes.
> > >>
> > >> Is supporting lambert conformal and polar stereographic grids
> > >> where dx != dy a requirement for your work? Is this
shortcoming
> > >> standing in your way of using MET?
> > >>
> > >> Thanks,
> > >> John
> > >>
> > >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> > >> <johnhg at ucar.edu>
> > >> wrote:
> > >>
> > >> > Jose, yes, running with dx = dy would solve the problem.
> > >> >
> > >> > John
> > >> >
> > >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT <
> > >> > met_help at ucar.edu> wrote:
> > >> >
> > >> >>
> > >> >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> >
> > >> >>
> > >> >> Dear John,
> > >> >>
> > >> >> Thank you for your insights. One last question: If we run
the
> > >> >> model with dx = dy, should that solve the problem?
> > >> >>
> > >> >> Best Regards,
> > >> >> Jose
> > >> >>
> > >> >> -----Original Message-----
> > >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> > >> >> To: Jose Garcia-Rivera
> > >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >>
> > >> >> Jose,
> > >> >>
> > >> >> Yes, that's correct. Not supporting differing Lambert
> > >> >> Conformal Dx and Dy values is a short-coming in MET.
However,
> > >> >> we're investigating the impact/cost of adding support for
it.
> > >> >> We may be able to add support by patching the build.
> > >> >>
> > >> >> I'll write you back once we figure out how to proceed.
> > >> >>
> > >> >> Thanks,
> > >> >> John
> > >> >>
> > >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT <
> > >> >> met_help at ucar.edu> wrote:
> > >> >>
> > >> >> >
> > >> >> > <URL:
> > >> >> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > >
> > >> >> >
> > >> >> > Hi John,
> > >> >> >
> > >> >> > Thanks so much for the information. So, that means that
I'll
> > >> >> > have the same problem when running any other analyses with
> > >> >> > this particular grid
> > >> >> using MET?
> > >> >> >
> > >> >> > Best Regards,
> > >> >> > Jose
> > >> >> >
> > >> >> > -----Original Message-----
> > >> >> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> > >> >> > To: Jose Garcia-Rivera
> > >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >> >
> > >> >> > Jose,
> > >> >> >
> > >> >> > I believe I've figured out the problem. The grid code in
> > >> >> > MET currently only supports Lambert Conformal grids where
> > >> >> > the grid spacing is the same in the x and y directions.
> > >> >> >
> > >> >> > Your grid is spaced 3.317km in the x-direction and 2.97km
in
> > >> >> > the y-direction. The MET code is reading the Dx value and
> > >> >> > using it for both x and y. This causes the skewing of the
> > >> >> > data to the north that
> > >> >> you're seeing.
> > >> >> >
> > >> >> > Really, there are 2 problems here.
> > >> >> >
> > >> >> > First, if we're requiring that Dx and Dy be the same for
> > >> >> > Lambert Conformal grids, MET should error out with a
> > >> >> > descriptive error message when they are not.
> > >> >> >
> > >> >> > Second, obviously it would be preferable for MET to
support
> > >> >> > Lambert Conformal grids where Dx and Dy are not the same!
> > >> >> > I'll ask around here for more information on that.
> > >> >> >
> > >> >> > Thanks for paying close attention to your data and
> > >> >> > uncovering this shortcoming in MET.
> > >> >> >
> > >> >> > John
> > >> >> >
> > >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via RT
<
> > >> >> > met_help at ucar.edu
> > >> >> > > wrote:
> > >> >> >
> > >> >> > >
> > >> >> > > <URL:
> > >> >> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > > >
> > >> >> > >
> > >> >> > > Hi John,
> > >> >> > >
> > >> >> > > Thanks for all the help. I want to add this pdf, showing
> > >> >> > > something odd with the regridding. The domains did not
> > >> >> > > match perfectly on the north side, it's like a skewing
> > >> >> > > effect. The south sides match perfectly, but note the
> > >> >> > > differences in the regridded RTMA file and the model
file
> > >> >> > > in Michigan, for
> example.
> > >> >> > > Maybe that has something to
> > >> >> do with it.
> > >> >> > >
> > >> >> > > Best Regards,
> > >> >> > > Jose
> > >> >> > >
> > >> >> > > -----Original Message-----
> > >> >> > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > >> >> > > To: Jose Garcia-Rivera
> > >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >> >> > >
> > >> >> > > Jose,
> > >> >> > >
> > >> >> > > I believe the problem is in our handling the domain on
> > >> >> > > which your model data resides. I ran the following
> > >> >> > > commands to plot temperature from your model data file:
> > >> >> > >
> > >> >> > > met-5.1/bin/plot_data_plane \
> > >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > >> >> > > 'name="TMP"; level="Z0";'
> > >> >> > >
> > >> >> > > convert -rotate 90 -background white -flatten \
> > >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > >> >> > >
> > >> >> > > The resulting image is attached and shows the shift you
> > described.
> > >> >> > >
> > >> >> > > So it isn't the regridding step that's the problem...
its
> > >> >> > > our handling of that domain. I plotted the same data
> > >> >> > > using IDV and it puts things in the correct place (also
attached).
> > >> >> > >
> > >> >> > > One issue is that the 2-m temperature data in GRIB
record
> > >> >> > > 296 contains a constant value of 273. The surface
> > >> >> > > temperature data in GRIB record
> > >> >> > > 276 contains more reasonable values. So that's the
field
> > >> >> > > I'm
> > >> >> > investigating.
> > >> >> > > But you might want to check why GRIB record 296 looks so
odd.
> > >> >> > >
> > >> >> > > I need to head out of the office, but will investigate
> > >> >> > > further
> > >> >> tomorrow.
> > >> >> > >
> > >> >> > > Thanks,
> > >> >> > > John
> > >> >> > >
> > >> >> > >
> > >> >> > >
> > >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via
RT
> > >> >> > > < met_help at ucar.edu
> > >> >> > > > wrote:
> > >> >> > >
> > >> >> > > >
> > >> >> > > > <URL:
> > >> >> > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > >> >> > > > >
> > >> >> > > >
> > >> >> > > > Dear John,
> > >> >> > > >
> > >> >> > > > Thank you so much for the quick reply.
> > >> >> > > >
> > >> >> > > > I followed the instructions, and ftp'ed the files to
> > >> >> > > "garcia-rivera_data".
> > >> >> > > >
> > >> >> > > > The nmbprs file is the .grb NMMB output file.
> > >> >> > > > The rtma file is the observational.
> > >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
> > >> >> > > > outputs in the same grid, and produces the plot.
> > >> >> > > >
> > >> >> > > > Any help will be greatly appreciated.
> > >> >> > > >
> > >> >> > > > Best Regards,
> > >> >> > > > Jose
> > >> >> > > >
> > >> >> > > > -----Original Message-----
> > >> >> > > > From: John Halley Gotway via RT
> > >> >> > > > [mailto:met_help at ucar.edu]
> > >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > >> >> > > > To: Jose Garcia-Rivera
> > >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > >> >> > > >
> > >> >> > > > Jose,
> > >> >> > > >
> > >> >> > > > Thanks for letting us know about the issue you're
seeing.
> > >> >> > > > I'd like to replicate the behavior your describe using
> > >> >> > > > the exact data
> > >> >> you used.
> > >> >> > > > Once I can replicate it, I'll be able to investigate
further.
> > >> >> > > >
> > >> >> > > > Are you able to send me the following 2 files you
used:
> > >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > >> >> > > >
> > >> >> > > > You can post them to our anonymous ftp site, following
> > >> >> > > > these
> > >> >> > > instructions:
> > >> >> > > >
> > >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#f
> > >> >> > > > tp
> > >> >> > > >
> > >> >> > > > Please write back to let me know when you've posted
the
> > >> >> > > > data and I'll go grab it.
> > >> >> > > >
> > >> >> > > > Thanks,
> > >> >> > > > John Halley Gotway
> > >> >> > > > met_help at ucar.edu
> > >> >> > > >
> > >> >> > > >
> > >> >> > > >
> > >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera via
> > >> >> > > > RT < met_help at ucar.edu
> > >> >> > > > > wrote:
> > >> >> > > >
> > >> >> > > > >
> > >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> > >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > >> >> > > > > Queue: met_help
> > >> >> > > > > Subject: Regridding problem
> > >> >> > > > > Owner: Nobody
> > >> >> > > > > Requestors: GarciaJ at imsg.com
> > >> >> > > > > Status: new
> > >> >> > > > > Ticket <URL:
> > >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=7397
> > >> >> > > > > 5
> > >> >> > > > > >
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > Hello,
> > >> >> > > > >
> > >> >> > > > > I'm a new MET user, and I'm attempting to regrid
RTMA
> > >> >> > > > > data to an NMMB output grid. I used the following
tool:
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > > > ./regrid_data_plane rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
> level="Z2";'
> > >> >> > > > > -method BILIN -width
> > >> >> > > > > 2 -v 1
> > >> >> > > > >
> > >> >> > > > > which worked. But, when I subtracted the fields from
> > >> >> > > > > one dataset to the other, there was a clear offset
on
> > >> >> > > > > the
> values.
> > >> >> > > > > So, the regridding was done indeed, but the
gridpoint
> > >> >> > > > > values did not seem to be assigned correctly with
> > >> >> > > > > respect to their
> > >> counterparts.
> > >> >> > > > >
> > >> >> > > > > I attached the plot, so you can see the offset
problem.
> > >> >> > > > > Any help will be greatly appreciated.
> > >> >> > > > >
> > >> >> > > > > Best Regards,
> > >> >> > > > > Jose
> > >> >> > > > >
> > >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M.
Systems
> > >> >> > > > > Group This e-mail, including attachments, may
include
> > >> >> > > > > confidential and/or proprietary information, and may
> > >> >> > > > > be used only by the person or entity to which it is
> > >> >> > > > > addressed. If the reader of this e-mail is not the
> > >> >> > > > > intended recipient or his or her authorized agent,
the
> > >> >> > > > > reader is hereby notified that any dissemination,
> > >> >> > > > > distribution or copying of this e-mail is
prohibited.
> > >> >> > > > > If you have received this e-mail in error, please
> > >> >> > > > > notify the sender by replying to this message and
> > >> >> > > > > delete this
> > >> >> > > > e-mail immediately.
> > >> >> > > > >
> > >> >> > > > >
> > >> >> > > >
> > >> >> > > > This e-mail, including attachments, may include
> > >> >> > > > confidential and/or proprietary information, and may
be
> > >> >> > > > used only by the person or entity to which it is
> > >> >> > > > addressed. If the reader of this e-mail is not the
> > >> >> > > > intended recipient or his or her authorized agent, the
> > >> >> > > > reader is hereby notified that any dissemination,
> > >> >> > > > distribution or copying of this e-mail is prohibited.
If
> > >> >> > > > you have received this e-mail in error, please notify
> > >> >> > > > the sender by replying to this message and delete this
> > >> >> > > e-mail immediately.
> > >> >> > > >
> > >> >> > > >
> > >> >> > > >
> > >> >> > >
> > >> >> > > This e-mail, including attachments, may include
> > >> >> > > confidential and/or proprietary information, and may be
> > >> >> > > used only by the person or entity to which it is
> > >> >> > > addressed. If the reader of this e-mail is not the
> > >> >> > > intended recipient or his or her authorized agent, the
> > >> >> > > reader is hereby notified that any dissemination,
> > >> >> > > distribution or copying of this e-mail is prohibited. If
> > >> >> > > you have received this e-mail in error, please notify
the
> > >> >> > > sender by replying to this message and delete this
> > >> >> > e-mail immediately.
> > >> >> > >
> > >> >> > >
> > >> >> >
> > >> >> > This e-mail, including attachments, may include
confidential
> > >> >> > and/or proprietary information, and may be used only by
the
> > >> >> > person or entity to which it is addressed. If the reader
of
> > >> >> > this e-mail is not the intended recipient or his or her
> > >> >> > authorized agent, the reader is hereby notified that any
> > >> >> > dissemination, distribution or copying of this e-mail is
> > >> >> > prohibited. If you have received this e-mail in error,
> > >> >> > please notify the sender by replying to this message and
> > >> >> > delete this
> > >> >> e-mail immediately.
> > >> >> >
> > >> >> >
> > >> >> >
> > >> >>
> > >> >> This e-mail, including attachments, may include confidential
> > >> >> and/or proprietary information, and may be used only by the
> > >> >> person or entity to which it is addressed. If the reader of
> > >> >> this e-mail is not the intended recipient or his or her
> > >> >> authorized agent, the reader is hereby notified that any
> > >> >> dissemination, distribution or copying of this e-mail is
> > >> >> prohibited. If you have received this e-mail in error,
please
> > >> >> notify the sender by replying to this message and delete
this
> > >> e-mail immediately.
> > >> >>
> > >> >>
> > >> >>
> > >> >
> > >>
> > >> This e-mail, including attachments, may include confidential
> > >> and/or proprietary information, and may be used only by the
> > >> person or entity to which it is addressed. If the reader of
this
> > >> e-mail is not the intended recipient or his or her authorized
> > >> agent, the reader is hereby notified that any dissemination,
> > >> distribution or copying of this e-mail is prohibited. If you
have
> > >> received this e-mail in error, please notify the sender by
> > >> replying to this message and delete this
> > e-mail immediately.
> > >>
> > >>
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Fri Nov 06 13:43:48 2015
Jose,
The cnvgrib utility can be downloaded here:
http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/
We talked internally about how to proceed with this dx != dy issue for
Lambert Confromal and Polar Stereographic grids. Unfortunately, our
funding for supporting MET is very tight for the coming months, and we
just
don't have the resources to add support as a patch for MET version
5.1.
However, we have added it to the list of development tasks for the
next
release. And as a temporary solution, I'll post a bugfix that checks
for
this dx != dy condition and prints a warning message about it.
Typically,
I would make this an error, but that would prevent you from being able
to
run your data through MET at all.
Thanks for bringing this short-coming in MET to our attention, and
sorry we
don't have a more timely solution for you.
Thanks,
John
On Fri, Nov 6, 2015 at 11:05 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Where is the cnvgrib utility?
>
> Thanks,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, November 06, 2015 12:15 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> The regridding tool is named copygb, but it only works on GRIB1
files.
> The RTMA data is in GRIB2. Here's one way of doing it:
>
> # Convert RTMA from GRIB2 to GRIB1
> cnvgrib -g21
> rtma2p5.t00z.2dvaranl_ndfd.grb2rtma2p5.t00z.2dvaranl_ndfd.grb
>
> # Run copygb to regrid GRIB1 file
> copygb \
> -xg'255 3 400 400 33302 -84951 8 -77459 3317 2970 0 64 38947
38947' \
> rtma2p5.t00z.2dvaranl_ndfd.grb
rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
>
> FYI, to do this I just followed an example from the MET online
tutorial:
>
>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run3.php
>
> Run wgrib on the output to verify that it's on the intended grid:
>
> wgrib -d 1 -V rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
> ...
> Lambert Conf: Lat1 33.302000 Lon1 -84.951000 Lov -77.459000
> Latin1 38.947000 Latin2 38.947000 LatSP 0.000000 LonSP
0.000000
> North Pole (400 x 400) Dx 3.317000 Dy 2.970000 scan 64 mode 8
...
>
> This data can now be run through MET. For example, running it
through
> grid-stat will give you statistics computed point-by-point and those
should
> be correct. But the underlying projection is still wrong because dx
!=
> dy. So any plots created by MET will show that skew to the north.
>
> It may also be possible to do the regridding using the "wgrib2" tool
> directly on the GRIB2 file. But I haven't done that as much. You
could
> look into that if you want to skip the GRIB2 to GRIB1 conversion
step.
>
> Thanks,
> John
>
>
>
>
> On Fri, Nov 6, 2015 at 9:51 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > Thanks for all the help. The target grid would be the original
grid I
> > sent you.
> >
> > This is the same grid in copy_gridnav.txt:
> >
> > 255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64 38947
38947
> > 44113 -69980
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Friday, November 06, 2015 11:49 AM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Sure, I can help with you that, but I'll need to know the "target"
> > grid definition. You said you reran with an NMMB grid with dx =
dy.
> > And you sent images of the output, but I don't think you sent me a
data
> file.
> >
> > If you send me that GRIB file, I can see the desired grid info and
> > send you a copygb command that'll put the RTMA data onto that
grid.
> >
> > Thanks,
> > John
> >
> > On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Dear John,
> > >
> > > We want to use MET, but we need to get the regridding right. I
was
> > > wondering if you know how to regrid the RTMA file to match an
NMMB
> > > grid by using UPP (as in previous MET versions)?
> > >
> > > Any help will be greatly appreciated, Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Thursday, November 05, 2015 6:31 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > Just wanted to follow up. Do you still see a shifting problem
in
> > > the regridding? And is this a showstopper for you guys? Is
adding
> > > support for differing dx/dy as a patch for 5.1 a high priority
or
> > > can we delay until the next release?
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway
<johnhg at ucar.edu>
> > > wrote:
> > >
> > > > Jose,
> > > >
> > > > Are you able to send me the GRIB files as well by posting it
to
> > > > our anonymous ftp site?
> > > >
> > > > I took a look at the images you sent, but the skew you mention
> > > > isn't jumping out at me. Before I saw that the values that
should
> > > > be over the great lakes were shifted north by a large extent.
> > > > What should I be looking for now?
> > > >
> > > > Rather than diffing the NMM and RTMA values, wouldn't it make
more
> > > > sense to...
> > > >
> > > > (1) Run the RTMA through regrid_data_plane to regrid it
> > > > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > > > (3) Run plot_data_plane to plot the output from
regrid_data_plane
> > > >
> > > > Compare images from the input/output to make sure the values
are
> > > > all in the right places.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >>
> > > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> > > >>
> > > >> Dear John,
> > > >>
> > > >> I just ran a simulation in which DX=DY, and the skewing was
much
> > > >> smaller, but it's still observed.
> > > >>
> > > >> Best Regards,
> > > >> Jose
> > > >>
> > > >>
> > > >> -----Original Message-----
> > > >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > >> Sent: Wednesday, November 04, 2015 4:28 PM
> > > >> To: Jose Garcia-Rivera
> > > >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >>
> > > >> Jose,
> > > >>
> > > >> We're still trying to decide how to address this issue. We
> > > >> definitely plan on adding support for grids where dx != dy,
but
> > > >> the question is whether we should add it as a patch for met-
5.1
> > > >> or just include it in the next release of the code.
> > > >>
> > > >> The enhancement will touch several files, and we worry about
the
> > > >> potential for unintended consequences of the changes.
> > > >>
> > > >> Is supporting lambert conformal and polar stereographic grids
> > > >> where dx != dy a requirement for your work? Is this
shortcoming
> > > >> standing in your way of using MET?
> > > >>
> > > >> Thanks,
> > > >> John
> > > >>
> > > >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> > > >> <johnhg at ucar.edu>
> > > >> wrote:
> > > >>
> > > >> > Jose, yes, running with dx = dy would solve the problem.
> > > >> >
> > > >> > John
> > > >> >
> > > >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT
<
> > > >> > met_help at ucar.edu> wrote:
> > > >> >
> > > >> >>
> > > >> >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> >
> > > >> >>
> > > >> >> Dear John,
> > > >> >>
> > > >> >> Thank you for your insights. One last question: If we run
the
> > > >> >> model with dx = dy, should that solve the problem?
> > > >> >>
> > > >> >> Best Regards,
> > > >> >> Jose
> > > >> >>
> > > >> >> -----Original Message-----
> > > >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> > > >> >> To: Jose Garcia-Rivera
> > > >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >> >>
> > > >> >> Jose,
> > > >> >>
> > > >> >> Yes, that's correct. Not supporting differing Lambert
> > > >> >> Conformal Dx and Dy values is a short-coming in MET.
However,
> > > >> >> we're investigating the impact/cost of adding support for
it.
> > > >> >> We may be able to add support by patching the build.
> > > >> >>
> > > >> >> I'll write you back once we figure out how to proceed.
> > > >> >>
> > > >> >> Thanks,
> > > >> >> John
> > > >> >>
> > > >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT
<
> > > >> >> met_help at ucar.edu> wrote:
> > > >> >>
> > > >> >> >
> > > >> >> > <URL:
> > > >> >> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> > >
> > > >> >> >
> > > >> >> > Hi John,
> > > >> >> >
> > > >> >> > Thanks so much for the information. So, that means that
I'll
> > > >> >> > have the same problem when running any other analyses
with
> > > >> >> > this particular grid
> > > >> >> using MET?
> > > >> >> >
> > > >> >> > Best Regards,
> > > >> >> > Jose
> > > >> >> >
> > > >> >> > -----Original Message-----
> > > >> >> > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > > >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> > > >> >> > To: Jose Garcia-Rivera
> > > >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >> >> >
> > > >> >> > Jose,
> > > >> >> >
> > > >> >> > I believe I've figured out the problem. The grid code
in
> > > >> >> > MET currently only supports Lambert Conformal grids
where
> > > >> >> > the grid spacing is the same in the x and y directions.
> > > >> >> >
> > > >> >> > Your grid is spaced 3.317km in the x-direction and
2.97km in
> > > >> >> > the y-direction. The MET code is reading the Dx value
and
> > > >> >> > using it for both x and y. This causes the skewing of
the
> > > >> >> > data to the north that
> > > >> >> you're seeing.
> > > >> >> >
> > > >> >> > Really, there are 2 problems here.
> > > >> >> >
> > > >> >> > First, if we're requiring that Dx and Dy be the same for
> > > >> >> > Lambert Conformal grids, MET should error out with a
> > > >> >> > descriptive error message when they are not.
> > > >> >> >
> > > >> >> > Second, obviously it would be preferable for MET to
support
> > > >> >> > Lambert Conformal grids where Dx and Dy are not the
same!
> > > >> >> > I'll ask around here for more information on that.
> > > >> >> >
> > > >> >> > Thanks for paying close attention to your data and
> > > >> >> > uncovering this shortcoming in MET.
> > > >> >> >
> > > >> >> > John
> > > >> >> >
> > > >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via
RT <
> > > >> >> > met_help at ucar.edu
> > > >> >> > > wrote:
> > > >> >> >
> > > >> >> > >
> > > >> >> > > <URL:
> > > >> >> > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> > > >
> > > >> >> > >
> > > >> >> > > Hi John,
> > > >> >> > >
> > > >> >> > > Thanks for all the help. I want to add this pdf,
showing
> > > >> >> > > something odd with the regridding. The domains did not
> > > >> >> > > match perfectly on the north side, it's like a skewing
> > > >> >> > > effect. The south sides match perfectly, but note the
> > > >> >> > > differences in the regridded RTMA file and the model
file
> > > >> >> > > in Michigan, for
> > example.
> > > >> >> > > Maybe that has something to
> > > >> >> do with it.
> > > >> >> > >
> > > >> >> > > Best Regards,
> > > >> >> > > Jose
> > > >> >> > >
> > > >> >> > > -----Original Message-----
> > > >> >> > > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > > >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > > >> >> > > To: Jose Garcia-Rivera
> > > >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > > >> >> > >
> > > >> >> > > Jose,
> > > >> >> > >
> > > >> >> > > I believe the problem is in our handling the domain on
> > > >> >> > > which your model data resides. I ran the following
> > > >> >> > > commands to plot temperature from your model data
file:
> > > >> >> > >
> > > >> >> > > met-5.1/bin/plot_data_plane \
> > > >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > > >> >> > > 'name="TMP"; level="Z0";'
> > > >> >> > >
> > > >> >> > > convert -rotate 90 -background white -flatten \
> > > >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > > >> >> > >
> > > >> >> > > The resulting image is attached and shows the shift
you
> > > described.
> > > >> >> > >
> > > >> >> > > So it isn't the regridding step that's the problem...
its
> > > >> >> > > our handling of that domain. I plotted the same data
> > > >> >> > > using IDV and it puts things in the correct place
(also
> attached).
> > > >> >> > >
> > > >> >> > > One issue is that the 2-m temperature data in GRIB
record
> > > >> >> > > 296 contains a constant value of 273. The surface
> > > >> >> > > temperature data in GRIB record
> > > >> >> > > 276 contains more reasonable values. So that's the
field
> > > >> >> > > I'm
> > > >> >> > investigating.
> > > >> >> > > But you might want to check why GRIB record 296 looks
so odd.
> > > >> >> > >
> > > >> >> > > I need to head out of the office, but will investigate
> > > >> >> > > further
> > > >> >> tomorrow.
> > > >> >> > >
> > > >> >> > > Thanks,
> > > >> >> > > John
> > > >> >> > >
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via
RT
> > > >> >> > > < met_help at ucar.edu
> > > >> >> > > > wrote:
> > > >> >> > >
> > > >> >> > > >
> > > >> >> > > > <URL:
> > > >> >> > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> > > > >
> > > >> >> > > >
> > > >> >> > > > Dear John,
> > > >> >> > > >
> > > >> >> > > > Thank you so much for the quick reply.
> > > >> >> > > >
> > > >> >> > > > I followed the instructions, and ftp'ed the files to
> > > >> >> > > "garcia-rivera_data".
> > > >> >> > > >
> > > >> >> > > > The nmbprs file is the .grb NMMB output file.
> > > >> >> > > > The rtma file is the observational.
> > > >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
> > > >> >> > > > outputs in the same grid, and produces the plot.
> > > >> >> > > >
> > > >> >> > > > Any help will be greatly appreciated.
> > > >> >> > > >
> > > >> >> > > > Best Regards,
> > > >> >> > > > Jose
> > > >> >> > > >
> > > >> >> > > > -----Original Message-----
> > > >> >> > > > From: John Halley Gotway via RT
> > > >> >> > > > [mailto:met_help at ucar.edu]
> > > >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > >> >> > > > To: Jose Garcia-Rivera
> > > >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > > >> >> > > >
> > > >> >> > > > Jose,
> > > >> >> > > >
> > > >> >> > > > Thanks for letting us know about the issue you're
seeing.
> > > >> >> > > > I'd like to replicate the behavior your describe
using
> > > >> >> > > > the exact data
> > > >> >> you used.
> > > >> >> > > > Once I can replicate it, I'll be able to investigate
> further.
> > > >> >> > > >
> > > >> >> > > > Are you able to send me the following 2 files you
used:
> > > >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > > >> >> > > >
> > > >> >> > > > You can post them to our anonymous ftp site,
following
> > > >> >> > > > these
> > > >> >> > > instructions:
> > > >> >> > > >
> > > >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php#f
> > > >> >> > > > tp
> > > >> >> > > >
> > > >> >> > > > Please write back to let me know when you've posted
the
> > > >> >> > > > data and I'll go grab it.
> > > >> >> > > >
> > > >> >> > > > Thanks,
> > > >> >> > > > John Halley Gotway
> > > >> >> > > > met_help at ucar.edu
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera
via
> > > >> >> > > > RT < met_help at ucar.edu
> > > >> >> > > > > wrote:
> > > >> >> > > >
> > > >> >> > > > >
> > > >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> > > >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > >> >> > > > > Queue: met_help
> > > >> >> > > > > Subject: Regridding problem
> > > >> >> > > > > Owner: Nobody
> > > >> >> > > > > Requestors: GarciaJ at imsg.com
> > > >> >> > > > > Status: new
> > > >> >> > > > > Ticket <URL:
> > > >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=7397
> > > >> >> > > > > 5
> > > >> >> > > > > >
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > Hello,
> > > >> >> > > > >
> > > >> >> > > > > I'm a new MET user, and I'm attempting to regrid
RTMA
> > > >> >> > > > > data to an NMMB output grid. I used the following
tool:
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > ./regrid_data_plane
rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > > >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
> > level="Z2";'
> > > >> >> > > > > -method BILIN -width
> > > >> >> > > > > 2 -v 1
> > > >> >> > > > >
> > > >> >> > > > > which worked. But, when I subtracted the fields
from
> > > >> >> > > > > one dataset to the other, there was a clear offset
on
> > > >> >> > > > > the
> > values.
> > > >> >> > > > > So, the regridding was done indeed, but the
gridpoint
> > > >> >> > > > > values did not seem to be assigned correctly with
> > > >> >> > > > > respect to their
> > > >> counterparts.
> > > >> >> > > > >
> > > >> >> > > > > I attached the plot, so you can see the offset
problem.
> > > >> >> > > > > Any help will be greatly appreciated.
> > > >> >> > > > >
> > > >> >> > > > > Best Regards,
> > > >> >> > > > > Jose
> > > >> >> > > > >
> > > >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M.
Systems
> > > >> >> > > > > Group This e-mail, including attachments, may
include
> > > >> >> > > > > confidential and/or proprietary information, and
may
> > > >> >> > > > > be used only by the person or entity to which it
is
> > > >> >> > > > > addressed. If the reader of this e-mail is not the
> > > >> >> > > > > intended recipient or his or her authorized agent,
the
> > > >> >> > > > > reader is hereby notified that any dissemination,
> > > >> >> > > > > distribution or copying of this e-mail is
prohibited.
> > > >> >> > > > > If you have received this e-mail in error, please
> > > >> >> > > > > notify the sender by replying to this message and
> > > >> >> > > > > delete this
> > > >> >> > > > e-mail immediately.
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > >
> > > >> >> > > > This e-mail, including attachments, may include
> > > >> >> > > > confidential and/or proprietary information, and may
be
> > > >> >> > > > used only by the person or entity to which it is
> > > >> >> > > > addressed. If the reader of this e-mail is not the
> > > >> >> > > > intended recipient or his or her authorized agent,
the
> > > >> >> > > > reader is hereby notified that any dissemination,
> > > >> >> > > > distribution or copying of this e-mail is
prohibited. If
> > > >> >> > > > you have received this e-mail in error, please
notify
> > > >> >> > > > the sender by replying to this message and delete
this
> > > >> >> > > e-mail immediately.
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > >
> > > >> >> > > This e-mail, including attachments, may include
> > > >> >> > > confidential and/or proprietary information, and may
be
> > > >> >> > > used only by the person or entity to which it is
> > > >> >> > > addressed. If the reader of this e-mail is not the
> > > >> >> > > intended recipient or his or her authorized agent, the
> > > >> >> > > reader is hereby notified that any dissemination,
> > > >> >> > > distribution or copying of this e-mail is prohibited.
If
> > > >> >> > > you have received this e-mail in error, please notify
the
> > > >> >> > > sender by replying to this message and delete this
> > > >> >> > e-mail immediately.
> > > >> >> > >
> > > >> >> > >
> > > >> >> >
> > > >> >> > This e-mail, including attachments, may include
confidential
> > > >> >> > and/or proprietary information, and may be used only by
the
> > > >> >> > person or entity to which it is addressed. If the reader
of
> > > >> >> > this e-mail is not the intended recipient or his or her
> > > >> >> > authorized agent, the reader is hereby notified that any
> > > >> >> > dissemination, distribution or copying of this e-mail is
> > > >> >> > prohibited. If you have received this e-mail in error,
> > > >> >> > please notify the sender by replying to this message and
> > > >> >> > delete this
> > > >> >> e-mail immediately.
> > > >> >> >
> > > >> >> >
> > > >> >> >
> > > >> >>
> > > >> >> This e-mail, including attachments, may include
confidential
> > > >> >> and/or proprietary information, and may be used only by
the
> > > >> >> person or entity to which it is addressed. If the reader
of
> > > >> >> this e-mail is not the intended recipient or his or her
> > > >> >> authorized agent, the reader is hereby notified that any
> > > >> >> dissemination, distribution or copying of this e-mail is
> > > >> >> prohibited. If you have received this e-mail in error,
please
> > > >> >> notify the sender by replying to this message and delete
this
> > > >> e-mail immediately.
> > > >> >>
> > > >> >>
> > > >> >>
> > > >> >
> > > >>
> > > >> This e-mail, including attachments, may include confidential
> > > >> and/or proprietary information, and may be used only by the
> > > >> person or entity to which it is addressed. If the reader of
this
> > > >> e-mail is not the intended recipient or his or her authorized
> > > >> agent, the reader is hereby notified that any dissemination,
> > > >> distribution or copying of this e-mail is prohibited. If you
have
> > > >> received this e-mail in error, please notify the sender by
> > > >> replying to this message and delete this
> > > e-mail immediately.
> > > >>
> > > >>
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #73975] Regridding problem
From: Jose Garcia-Rivera
Time: Fri Nov 06 15:28:41 2015
Dear John,
Thank you so much for your help. We'll use copygb for regridding, then
work with MET from there on. I'll keep you posted if any other
problems show up.
Have a great weekend!
Jose
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Friday, November 06, 2015 3:44 PM
To: Jose Garcia-Rivera
Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
Jose,
The cnvgrib utility can be downloaded here:
http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/
We talked internally about how to proceed with this dx != dy issue for
Lambert Confromal and Polar Stereographic grids. Unfortunately, our
funding for supporting MET is very tight for the coming months, and we
just don't have the resources to add support as a patch for MET
version 5.1.
However, we have added it to the list of development tasks for the
next release. And as a temporary solution, I'll post a bugfix that
checks for this dx != dy condition and prints a warning message about
it. Typically, I would make this an error, but that would prevent you
from being able to run your data through MET at all.
Thanks for bringing this short-coming in MET to our attention, and
sorry we don't have a more timely solution for you.
Thanks,
John
On Fri, Nov 6, 2015 at 11:05 AM, Jose Garcia-Rivera via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Hi John,
>
> Where is the cnvgrib utility?
>
> Thanks,
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, November 06, 2015 12:15 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> The regridding tool is named copygb, but it only works on GRIB1
files.
> The RTMA data is in GRIB2. Here's one way of doing it:
>
> # Convert RTMA from GRIB2 to GRIB1
> cnvgrib -g21
> rtma2p5.t00z.2dvaranl_ndfd.grb2rtma2p5.t00z.2dvaranl_ndfd.grb
>
> # Run copygb to regrid GRIB1 file
> copygb \
> -xg'255 3 400 400 33302 -84951 8 -77459 3317 2970 0 64 38947
38947' \
> rtma2p5.t00z.2dvaranl_ndfd.grb
> rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
>
> FYI, to do this I just followed an example from the MET online
tutorial:
>
>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copy
> gb/run3.php
>
> Run wgrib on the output to verify that it's on the intended grid:
>
> wgrib -d 1 -V rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
> ...
> Lambert Conf: Lat1 33.302000 Lon1 -84.951000 Lov -77.459000
> Latin1 38.947000 Latin2 38.947000 LatSP 0.000000 LonSP
0.000000
> North Pole (400 x 400) Dx 3.317000 Dy 2.970000 scan 64 mode 8
...
>
> This data can now be run through MET. For example, running it
through
> grid-stat will give you statistics computed point-by-point and those
> should be correct. But the underlying projection is still wrong
> because dx != dy. So any plots created by MET will show that skew
to the north.
>
> It may also be possible to do the regridding using the "wgrib2" tool
> directly on the GRIB2 file. But I haven't done that as much. You
> could look into that if you want to skip the GRIB2 to GRIB1
conversion step.
>
> Thanks,
> John
>
>
>
>
> On Fri, Nov 6, 2015 at 9:51 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Dear John,
> >
> > Thanks for all the help. The target grid would be the original
grid
> > I sent you.
> >
> > This is the same grid in copy_gridnav.txt:
> >
> > 255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64 38947
38947
> > 44113 -69980
> >
> > Best Regards,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Friday, November 06, 2015 11:49 AM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > Sure, I can help with you that, but I'll need to know the "target"
> > grid definition. You said you reran with an NMMB grid with dx =
dy.
> > And you sent images of the output, but I don't think you sent me a
> > data
> file.
> >
> > If you send me that GRIB file, I can see the desired grid info and
> > send you a copygb command that'll put the RTMA data onto that
grid.
> >
> > Thanks,
> > John
> >
> > On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Dear John,
> > >
> > > We want to use MET, but we need to get the regridding right. I
was
> > > wondering if you know how to regrid the RTMA file to match an
NMMB
> > > grid by using UPP (as in previous MET versions)?
> > >
> > > Any help will be greatly appreciated, Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Thursday, November 05, 2015 6:31 PM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > Just wanted to follow up. Do you still see a shifting problem
in
> > > the regridding? And is this a showstopper for you guys? Is
> > > adding support for differing dx/dy as a patch for 5.1 a high
> > > priority or can we delay until the next release?
> > >
> > > Thanks,
> > > John
> > >
> > > On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway
> > > <johnhg at ucar.edu>
> > > wrote:
> > >
> > > > Jose,
> > > >
> > > > Are you able to send me the GRIB files as well by posting it
to
> > > > our anonymous ftp site?
> > > >
> > > > I took a look at the images you sent, but the skew you mention
> > > > isn't jumping out at me. Before I saw that the values that
> > > > should be over the great lakes were shifted north by a large
extent.
> > > > What should I be looking for now?
> > > >
> > > > Rather than diffing the NMM and RTMA values, wouldn't it make
> > > > more sense to...
> > > >
> > > > (1) Run the RTMA through regrid_data_plane to regrid it
> > > > (2) Run plot_data_plane to plot the input to regrid_data_plane
> > > > (3) Run plot_data_plane to plot the output from
> > > > regrid_data_plane
> > > >
> > > > Compare images from the input/output to make sure the values
are
> > > > all in the right places.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >>
> > > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> > > >>
> > > >> Dear John,
> > > >>
> > > >> I just ran a simulation in which DX=DY, and the skewing was
> > > >> much smaller, but it's still observed.
> > > >>
> > > >> Best Regards,
> > > >> Jose
> > > >>
> > > >>
> > > >> -----Original Message-----
> > > >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > >> Sent: Wednesday, November 04, 2015 4:28 PM
> > > >> To: Jose Garcia-Rivera
> > > >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >>
> > > >> Jose,
> > > >>
> > > >> We're still trying to decide how to address this issue. We
> > > >> definitely plan on adding support for grids where dx != dy,
but
> > > >> the question is whether we should add it as a patch for met-
5.1
> > > >> or just include it in the next release of the code.
> > > >>
> > > >> The enhancement will touch several files, and we worry about
> > > >> the potential for unintended consequences of the changes.
> > > >>
> > > >> Is supporting lambert conformal and polar stereographic grids
> > > >> where dx != dy a requirement for your work? Is this
> > > >> shortcoming standing in your way of using MET?
> > > >>
> > > >> Thanks,
> > > >> John
> > > >>
> > > >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> > > >> <johnhg at ucar.edu>
> > > >> wrote:
> > > >>
> > > >> > Jose, yes, running with dx = dy would solve the problem.
> > > >> >
> > > >> > John
> > > >> >
> > > >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via RT
<
> > > >> > met_help at ucar.edu> wrote:
> > > >> >
> > > >> >>
> > > >> >> <URL:
> > > >> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> >
> > > >> >>
> > > >> >> Dear John,
> > > >> >>
> > > >> >> Thank you for your insights. One last question: If we run
> > > >> >> the model with dx = dy, should that solve the problem?
> > > >> >>
> > > >> >> Best Regards,
> > > >> >> Jose
> > > >> >>
> > > >> >> -----Original Message-----
> > > >> >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> > > >> >> To: Jose Garcia-Rivera
> > > >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >> >>
> > > >> >> Jose,
> > > >> >>
> > > >> >> Yes, that's correct. Not supporting differing Lambert
> > > >> >> Conformal Dx and Dy values is a short-coming in MET.
> > > >> >> However, we're investigating the impact/cost of adding
support for it.
> > > >> >> We may be able to add support by patching the build.
> > > >> >>
> > > >> >> I'll write you back once we figure out how to proceed.
> > > >> >>
> > > >> >> Thanks,
> > > >> >> John
> > > >> >>
> > > >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via RT
<
> > > >> >> met_help at ucar.edu> wrote:
> > > >> >>
> > > >> >> >
> > > >> >> > <URL:
> > > >> >> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> > >
> > > >> >> >
> > > >> >> > Hi John,
> > > >> >> >
> > > >> >> > Thanks so much for the information. So, that means that
> > > >> >> > I'll have the same problem when running any other
analyses
> > > >> >> > with this particular grid
> > > >> >> using MET?
> > > >> >> >
> > > >> >> > Best Regards,
> > > >> >> > Jose
> > > >> >> >
> > > >> >> > -----Original Message-----
> > > >> >> > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > > >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> > > >> >> > To: Jose Garcia-Rivera
> > > >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >> >> >
> > > >> >> > Jose,
> > > >> >> >
> > > >> >> > I believe I've figured out the problem. The grid code
in
> > > >> >> > MET currently only supports Lambert Conformal grids
where
> > > >> >> > the grid spacing is the same in the x and y directions.
> > > >> >> >
> > > >> >> > Your grid is spaced 3.317km in the x-direction and
2.97km
> > > >> >> > in the y-direction. The MET code is reading the Dx
value
> > > >> >> > and using it for both x and y. This causes the skewing
of
> > > >> >> > the data to the north that
> > > >> >> you're seeing.
> > > >> >> >
> > > >> >> > Really, there are 2 problems here.
> > > >> >> >
> > > >> >> > First, if we're requiring that Dx and Dy be the same for
> > > >> >> > Lambert Conformal grids, MET should error out with a
> > > >> >> > descriptive error message when they are not.
> > > >> >> >
> > > >> >> > Second, obviously it would be preferable for MET to
> > > >> >> > support Lambert Conformal grids where Dx and Dy are not
the same!
> > > >> >> > I'll ask around here for more information on that.
> > > >> >> >
> > > >> >> > Thanks for paying close attention to your data and
> > > >> >> > uncovering this shortcoming in MET.
> > > >> >> >
> > > >> >> > John
> > > >> >> >
> > > >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via
RT
> > > >> >> > < met_help at ucar.edu
> > > >> >> > > wrote:
> > > >> >> >
> > > >> >> > >
> > > >> >> > > <URL:
> > > >> >> > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > >> >> > > >
> > > >> >> > >
> > > >> >> > > Hi John,
> > > >> >> > >
> > > >> >> > > Thanks for all the help. I want to add this pdf,
showing
> > > >> >> > > something odd with the regridding. The domains did not
> > > >> >> > > match perfectly on the north side, it's like a skewing
> > > >> >> > > effect. The south sides match perfectly, but note the
> > > >> >> > > differences in the regridded RTMA file and the model
> > > >> >> > > file in Michigan, for
> > example.
> > > >> >> > > Maybe that has something to
> > > >> >> do with it.
> > > >> >> > >
> > > >> >> > > Best Regards,
> > > >> >> > > Jose
> > > >> >> > >
> > > >> >> > > -----Original Message-----
> > > >> >> > > From: John Halley Gotway via RT
> > > >> >> > > [mailto:met_help at ucar.edu]
> > > >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > > >> >> > > To: Jose Garcia-Rivera
> > > >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > > >> >> > >
> > > >> >> > > Jose,
> > > >> >> > >
> > > >> >> > > I believe the problem is in our handling the domain on
> > > >> >> > > which your model data resides. I ran the following
> > > >> >> > > commands to plot temperature from your model data
file:
> > > >> >> > >
> > > >> >> > > met-5.1/bin/plot_data_plane \
> > > >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > > >> >> > > 'name="TMP"; level="Z0";'
> > > >> >> > >
> > > >> >> > > convert -rotate 90 -background white -flatten \
> > > >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > > >> >> > >
> > > >> >> > > The resulting image is attached and shows the shift
you
> > > described.
> > > >> >> > >
> > > >> >> > > So it isn't the regridding step that's the problem...
> > > >> >> > > its our handling of that domain. I plotted the same
> > > >> >> > > data using IDV and it puts things in the correct place
> > > >> >> > > (also
> attached).
> > > >> >> > >
> > > >> >> > > One issue is that the 2-m temperature data in GRIB
> > > >> >> > > record
> > > >> >> > > 296 contains a constant value of 273. The surface
> > > >> >> > > temperature data in GRIB record
> > > >> >> > > 276 contains more reasonable values. So that's the
> > > >> >> > > field I'm
> > > >> >> > investigating.
> > > >> >> > > But you might want to check why GRIB record 296 looks
so odd.
> > > >> >> > >
> > > >> >> > > I need to head out of the office, but will investigate
> > > >> >> > > further
> > > >> >> tomorrow.
> > > >> >> > >
> > > >> >> > > Thanks,
> > > >> >> > > John
> > > >> >> > >
> > > >> >> > >
> > > >> >> > >
> > > >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera via
> > > >> >> > > RT < met_help at ucar.edu
> > > >> >> > > > wrote:
> > > >> >> > >
> > > >> >> > > >
> > > >> >> > > > <URL:
> > > >> >> > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=7397
> > > >> >> > > > 5
> > > >> >> > > > >
> > > >> >> > > >
> > > >> >> > > > Dear John,
> > > >> >> > > >
> > > >> >> > > > Thank you so much for the quick reply.
> > > >> >> > > >
> > > >> >> > > > I followed the instructions, and ftp'ed the files to
> > > >> >> > > "garcia-rivera_data".
> > > >> >> > > >
> > > >> >> > > > The nmbprs file is the .grb NMMB output file.
> > > >> >> > > > The rtma file is the observational.
> > > >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
> > > >> >> > > > outputs in the same grid, and produces the plot.
> > > >> >> > > >
> > > >> >> > > > Any help will be greatly appreciated.
> > > >> >> > > >
> > > >> >> > > > Best Regards,
> > > >> >> > > > Jose
> > > >> >> > > >
> > > >> >> > > > -----Original Message-----
> > > >> >> > > > From: John Halley Gotway via RT
> > > >> >> > > > [mailto:met_help at ucar.edu]
> > > >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > >> >> > > > To: Jose Garcia-Rivera
> > > >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
> > > >> >> > > > problem
> > > >> >> > > >
> > > >> >> > > > Jose,
> > > >> >> > > >
> > > >> >> > > > Thanks for letting us know about the issue you're
seeing.
> > > >> >> > > > I'd like to replicate the behavior your describe
using
> > > >> >> > > > the exact data
> > > >> >> you used.
> > > >> >> > > > Once I can replicate it, I'll be able to investigate
> further.
> > > >> >> > > >
> > > >> >> > > > Are you able to send me the following 2 files you
used:
> > > >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and nmbprs_d01.12
> > > >> >> > > >
> > > >> >> > > > You can post them to our anonymous ftp site,
following
> > > >> >> > > > these
> > > >> >> > > instructions:
> > > >> >> > > >
> > > >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php
> > > >> >> > > > #f
> > > >> >> > > > tp
> > > >> >> > > >
> > > >> >> > > > Please write back to let me know when you've posted
> > > >> >> > > > the data and I'll go grab it.
> > > >> >> > > >
> > > >> >> > > > Thanks,
> > > >> >> > > > John Halley Gotway
> > > >> >> > > > met_help at ucar.edu
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera
via
> > > >> >> > > > RT < met_help at ucar.edu
> > > >> >> > > > > wrote:
> > > >> >> > > >
> > > >> >> > > > >
> > > >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was acted
upon.
> > > >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > >> >> > > > > Queue: met_help
> > > >> >> > > > > Subject: Regridding problem
> > > >> >> > > > > Owner: Nobody
> > > >> >> > > > > Requestors: GarciaJ at imsg.com
> > > >> >> > > > > Status: new
> > > >> >> > > > > Ticket <URL:
> > > >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73
> > > >> >> > > > > 97
> > > >> >> > > > > 5
> > > >> >> > > > > >
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > Hello,
> > > >> >> > > > >
> > > >> >> > > > > I'm a new MET user, and I'm attempting to regrid
> > > >> >> > > > > RTMA data to an NMMB output grid. I used the
following tool:
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > > > ./regrid_data_plane
rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > > >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
> > level="Z2";'
> > > >> >> > > > > -method BILIN -width
> > > >> >> > > > > 2 -v 1
> > > >> >> > > > >
> > > >> >> > > > > which worked. But, when I subtracted the fields
from
> > > >> >> > > > > one dataset to the other, there was a clear offset
> > > >> >> > > > > on the
> > values.
> > > >> >> > > > > So, the regridding was done indeed, but the
> > > >> >> > > > > gridpoint values did not seem to be assigned
> > > >> >> > > > > correctly with respect to their
> > > >> counterparts.
> > > >> >> > > > >
> > > >> >> > > > > I attached the plot, so you can see the offset
problem.
> > > >> >> > > > > Any help will be greatly appreciated.
> > > >> >> > > > >
> > > >> >> > > > > Best Regards,
> > > >> >> > > > > Jose
> > > >> >> > > > >
> > > >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M.
> > > >> >> > > > > Systems Group This e-mail, including attachments,
> > > >> >> > > > > may include confidential and/or proprietary
> > > >> >> > > > > information, and may be used only by the person or
> > > >> >> > > > > entity to which it is addressed. If the reader of
> > > >> >> > > > > this e-mail is not the intended recipient or his
or
> > > >> >> > > > > her authorized agent, the reader is hereby
notified
> > > >> >> > > > > that any dissemination, distribution or copying of
this e-mail is prohibited.
> > > >> >> > > > > If you have received this e-mail in error, please
> > > >> >> > > > > notify the sender by replying to this message and
> > > >> >> > > > > delete this
> > > >> >> > > > e-mail immediately.
> > > >> >> > > > >
> > > >> >> > > > >
> > > >> >> > > >
> > > >> >> > > > This e-mail, including attachments, may include
> > > >> >> > > > confidential and/or proprietary information, and may
> > > >> >> > > > be used only by the person or entity to which it is
> > > >> >> > > > addressed. If the reader of this e-mail is not the
> > > >> >> > > > intended recipient or his or her authorized agent,
the
> > > >> >> > > > reader is hereby notified that any dissemination,
> > > >> >> > > > distribution or copying of this e-mail is
prohibited.
> > > >> >> > > > If you have received this e-mail in error, please
> > > >> >> > > > notify the sender by replying to this message and
> > > >> >> > > > delete this
> > > >> >> > > e-mail immediately.
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > > >
> > > >> >> > >
> > > >> >> > > This e-mail, including attachments, may include
> > > >> >> > > confidential and/or proprietary information, and may
be
> > > >> >> > > used only by the person or entity to which it is
> > > >> >> > > addressed. If the reader of this e-mail is not the
> > > >> >> > > intended recipient or his or her authorized agent, the
> > > >> >> > > reader is hereby notified that any dissemination,
> > > >> >> > > distribution or copying of this e-mail is prohibited.
If
> > > >> >> > > you have received this e-mail in error, please notify
> > > >> >> > > the sender by replying to this message and delete this
> > > >> >> > e-mail immediately.
> > > >> >> > >
> > > >> >> > >
> > > >> >> >
> > > >> >> > This e-mail, including attachments, may include
> > > >> >> > confidential and/or proprietary information, and may be
> > > >> >> > used only by the person or entity to which it is
> > > >> >> > addressed. If the reader of this e-mail is not the
> > > >> >> > intended recipient or his or her authorized agent, the
> > > >> >> > reader is hereby notified that any dissemination,
> > > >> >> > distribution or copying of this e-mail is prohibited. If
> > > >> >> > you have received this e-mail in error, please notify
the
> > > >> >> > sender by replying to this message and delete this
> > > >> >> e-mail immediately.
> > > >> >> >
> > > >> >> >
> > > >> >> >
> > > >> >>
> > > >> >> This e-mail, including attachments, may include
confidential
> > > >> >> and/or proprietary information, and may be used only by
the
> > > >> >> person or entity to which it is addressed. If the reader
of
> > > >> >> this e-mail is not the intended recipient or his or her
> > > >> >> authorized agent, the reader is hereby notified that any
> > > >> >> dissemination, distribution or copying of this e-mail is
> > > >> >> prohibited. If you have received this e-mail in error,
> > > >> >> please notify the sender by replying to this message and
> > > >> >> delete this
> > > >> e-mail immediately.
> > > >> >>
> > > >> >>
> > > >> >>
> > > >> >
> > > >>
> > > >> This e-mail, including attachments, may include confidential
> > > >> and/or proprietary information, and may be used only by the
> > > >> person or entity to which it is addressed. If the reader of
> > > >> this e-mail is not the intended recipient or his or her
> > > >> authorized agent, the reader is hereby notified that any
> > > >> dissemination, distribution or copying of this e-mail is
> > > >> prohibited. If you have received this e-mail in error, please
> > > >> notify the sender by replying to this message and delete this
> > > e-mail immediately.
> > > >>
> > > >>
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
> > > and/or proprietary information, and may be used only by the
person
> > > or entity to which it is addressed. If the reader of this e-mail
> > > is not the intended recipient or his or her authorized agent,
the
> > > reader is hereby notified that any dissemination, distribution
or
> > > copying of this e-mail is prohibited. If you have received this
> > > e-mail in error, please notify the sender by replying to this
> > > message and delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
> > entity to which it is addressed. If the reader of this e-mail is
not
> > the intended recipient or his or her authorized agent, the reader
is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
> > error, please notify the sender by replying to this message and
> > delete this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity
> to which it is addressed. If the reader of this e-mail is not the
> intended recipient or his or her authorized agent, the reader is
> hereby notified that any dissemination, distribution or copying of
> this e-mail is prohibited. If you have received this e-mail in
error,
> please notify the sender by replying to this message and delete this
e-mail immediately.
>
>
>
This e-mail, including attachments, may include confidential and/or
proprietary information, and may be used only by the person or entity
to which it is addressed. If the reader of this e-mail is not the
intended recipient or his or her authorized agent, the reader is
hereby notified that any dissemination, distribution or copying of
this e-mail is prohibited. If you have received this e-mail in error,
please notify the sender by replying to this message and delete this
e-mail immediately.
------------------------------------------------
Subject: Regridding problem
From: John Halley Gotway
Time: Fri Nov 06 17:23:20 2015
Jose,
I just posted a patch for met-5.1 about this issue, but all it does is
print a warning message when dx != dy:
http://www.dtcenter.org/met/users/support/known_issues/METv5.1/index.php
I'll go ahead and resolve this issue now.
Thanks,
John
On Fri, Nov 6, 2015 at 3:28 PM, Jose Garcia-Rivera via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
>
> Dear John,
>
> Thank you so much for your help. We'll use copygb for regridding,
then
> work with MET from there on. I'll keep you posted if any other
problems
> show up.
>
> Have a great weekend!
> Jose
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, November 06, 2015 3:44 PM
> To: Jose Garcia-Rivera
> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
>
> Jose,
>
> The cnvgrib utility can be downloaded here:
> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/
>
> We talked internally about how to proceed with this dx != dy issue
for
> Lambert Confromal and Polar Stereographic grids. Unfortunately, our
> funding for supporting MET is very tight for the coming months, and
we just
> don't have the resources to add support as a patch for MET version
5.1.
>
> However, we have added it to the list of development tasks for the
next
> release. And as a temporary solution, I'll post a bugfix that
checks for
> this dx != dy condition and prints a warning message about it.
Typically,
> I would make this an error, but that would prevent you from being
able to
> run your data through MET at all.
>
> Thanks for bringing this short-coming in MET to our attention, and
sorry
> we don't have a more timely solution for you.
>
> Thanks,
> John
>
> On Fri, Nov 6, 2015 at 11:05 AM, Jose Garcia-Rivera via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> >
> > Hi John,
> >
> > Where is the cnvgrib utility?
> >
> > Thanks,
> > Jose
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Sent: Friday, November 06, 2015 12:15 PM
> > To: Jose Garcia-Rivera
> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> >
> > Jose,
> >
> > The regridding tool is named copygb, but it only works on GRIB1
files.
> > The RTMA data is in GRIB2. Here's one way of doing it:
> >
> > # Convert RTMA from GRIB2 to GRIB1
> > cnvgrib -g21
> > rtma2p5.t00z.2dvaranl_ndfd.grb2rtma2p5.t00z.2dvaranl_ndfd.grb
> >
> > # Run copygb to regrid GRIB1 file
> > copygb \
> > -xg'255 3 400 400 33302 -84951 8 -77459 3317 2970 0 64 38947
38947' \
> > rtma2p5.t00z.2dvaranl_ndfd.grb
> > rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
> >
> > FYI, to do this I just followed an example from the MET online
tutorial:
> >
> >
> >
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copy
> > gb/run3.php
> >
> > Run wgrib on the output to verify that it's on the intended grid:
> >
> > wgrib -d 1 -V rtma2p5.t00z.2dvaranl_ndfd_regrid.grb
> > ...
> > Lambert Conf: Lat1 33.302000 Lon1 -84.951000 Lov -77.459000
> > Latin1 38.947000 Latin2 38.947000 LatSP 0.000000 LonSP
0.000000
> > North Pole (400 x 400) Dx 3.317000 Dy 2.970000 scan 64 mode
8 ...
> >
> > This data can now be run through MET. For example, running it
through
> > grid-stat will give you statistics computed point-by-point and
those
> > should be correct. But the underlying projection is still wrong
> > because dx != dy. So any plots created by MET will show that skew
to
> the north.
> >
> > It may also be possible to do the regridding using the "wgrib2"
tool
> > directly on the GRIB2 file. But I haven't done that as much. You
> > could look into that if you want to skip the GRIB2 to GRIB1
conversion
> step.
> >
> > Thanks,
> > John
> >
> >
> >
> >
> > On Fri, Nov 6, 2015 at 9:51 AM, Jose Garcia-Rivera via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > >
> > > Dear John,
> > >
> > > Thanks for all the help. The target grid would be the original
grid
> > > I sent you.
> > >
> > > This is the same grid in copy_gridnav.txt:
> > >
> > > 255 3 400 400 33371 -83884 8 -77459 2889 2970 0 64
38947
> 38947
> > > 44113 -69980
> > >
> > > Best Regards,
> > > Jose
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Sent: Friday, November 06, 2015 11:49 AM
> > > To: Jose Garcia-Rivera
> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > >
> > > Jose,
> > >
> > > Sure, I can help with you that, but I'll need to know the
"target"
> > > grid definition. You said you reran with an NMMB grid with dx =
dy.
> > > And you sent images of the output, but I don't think you sent me
a
> > > data
> > file.
> > >
> > > If you send me that GRIB file, I can see the desired grid info
and
> > > send you a copygb command that'll put the RTMA data onto that
grid.
> > >
> > > Thanks,
> > > John
> > >
> > > On Fri, Nov 6, 2015 at 8:54 AM, Jose Garcia-Rivera via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
>
> > > >
> > > > Dear John,
> > > >
> > > > We want to use MET, but we need to get the regridding right. I
was
> > > > wondering if you know how to regrid the RTMA file to match an
NMMB
> > > > grid by using UPP (as in previous MET versions)?
> > > >
> > > > Any help will be greatly appreciated, Jose
> > > >
> > > > -----Original Message-----
> > > > From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Sent: Thursday, November 05, 2015 6:31 PM
> > > > To: Jose Garcia-Rivera
> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > >
> > > > Jose,
> > > >
> > > > Just wanted to follow up. Do you still see a shifting problem
in
> > > > the regridding? And is this a showstopper for you guys? Is
> > > > adding support for differing dx/dy as a patch for 5.1 a high
> > > > priority or can we delay until the next release?
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Wed, Nov 4, 2015 at 3:03 PM, John Halley Gotway
> > > > <johnhg at ucar.edu>
> > > > wrote:
> > > >
> > > > > Jose,
> > > > >
> > > > > Are you able to send me the GRIB files as well by posting it
to
> > > > > our anonymous ftp site?
> > > > >
> > > > > I took a look at the images you sent, but the skew you
mention
> > > > > isn't jumping out at me. Before I saw that the values that
> > > > > should be over the great lakes were shifted north by a large
> extent.
> > > > > What should I be looking for now?
> > > > >
> > > > > Rather than diffing the NMM and RTMA values, wouldn't it
make
> > > > > more sense to...
> > > > >
> > > > > (1) Run the RTMA through regrid_data_plane to regrid it
> > > > > (2) Run plot_data_plane to plot the input to
regrid_data_plane
> > > > > (3) Run plot_data_plane to plot the output from
> > > > > regrid_data_plane
> > > > >
> > > > > Compare images from the input/output to make sure the values
are
> > > > > all in the right places.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > >
> > > > > On Wed, Nov 4, 2015 at 2:35 PM, Jose Garcia-Rivera via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > >>
> > > > >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975 >
> > > > >>
> > > > >> Dear John,
> > > > >>
> > > > >> I just ran a simulation in which DX=DY, and the skewing was
> > > > >> much smaller, but it's still observed.
> > > > >>
> > > > >> Best Regards,
> > > > >> Jose
> > > > >>
> > > > >>
> > > > >> -----Original Message-----
> > > > >> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > >> Sent: Wednesday, November 04, 2015 4:28 PM
> > > > >> To: Jose Garcia-Rivera
> > > > >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > > >>
> > > > >> Jose,
> > > > >>
> > > > >> We're still trying to decide how to address this issue. We
> > > > >> definitely plan on adding support for grids where dx != dy,
but
> > > > >> the question is whether we should add it as a patch for
met-5.1
> > > > >> or just include it in the next release of the code.
> > > > >>
> > > > >> The enhancement will touch several files, and we worry
about
> > > > >> the potential for unintended consequences of the changes.
> > > > >>
> > > > >> Is supporting lambert conformal and polar stereographic
grids
> > > > >> where dx != dy a requirement for your work? Is this
> > > > >> shortcoming standing in your way of using MET?
> > > > >>
> > > > >> Thanks,
> > > > >> John
> > > > >>
> > > > >> On Wed, Nov 4, 2015 at 11:50 AM, John Halley Gotway
> > > > >> <johnhg at ucar.edu>
> > > > >> wrote:
> > > > >>
> > > > >> > Jose, yes, running with dx = dy would solve the problem.
> > > > >> >
> > > > >> > John
> > > > >> >
> > > > >> > On Wed, Nov 4, 2015 at 11:44 AM, Jose Garcia-Rivera via
RT <
> > > > >> > met_help at ucar.edu> wrote:
> > > > >> >
> > > > >> >>
> > > > >> >> <URL:
> > > > >> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > >> >> >
> > > > >> >>
> > > > >> >> Dear John,
> > > > >> >>
> > > > >> >> Thank you for your insights. One last question: If we
run
> > > > >> >> the model with dx = dy, should that solve the problem?
> > > > >> >>
> > > > >> >> Best Regards,
> > > > >> >> Jose
> > > > >> >>
> > > > >> >> -----Original Message-----
> > > > >> >> From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > > > >> >> Sent: Wednesday, November 04, 2015 1:02 PM
> > > > >> >> To: Jose Garcia-Rivera
> > > > >> >> Subject: Re: [rt.rap.ucar.edu #73975] Regridding problem
> > > > >> >>
> > > > >> >> Jose,
> > > > >> >>
> > > > >> >> Yes, that's correct. Not supporting differing Lambert
> > > > >> >> Conformal Dx and Dy values is a short-coming in MET.
> > > > >> >> However, we're investigating the impact/cost of adding
support
> for it.
> > > > >> >> We may be able to add support by patching the build.
> > > > >> >>
> > > > >> >> I'll write you back once we figure out how to proceed.
> > > > >> >>
> > > > >> >> Thanks,
> > > > >> >> John
> > > > >> >>
> > > > >> >> On Wed, Nov 4, 2015 at 10:34 AM, Jose Garcia-Rivera via
RT <
> > > > >> >> met_help at ucar.edu> wrote:
> > > > >> >>
> > > > >> >> >
> > > > >> >> > <URL:
> > > > >> >> >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > >> >> > >
> > > > >> >> >
> > > > >> >> > Hi John,
> > > > >> >> >
> > > > >> >> > Thanks so much for the information. So, that means
that
> > > > >> >> > I'll have the same problem when running any other
analyses
> > > > >> >> > with this particular grid
> > > > >> >> using MET?
> > > > >> >> >
> > > > >> >> > Best Regards,
> > > > >> >> > Jose
> > > > >> >> >
> > > > >> >> > -----Original Message-----
> > > > >> >> > From: John Halley Gotway via RT
[mailto:met_help at ucar.edu]
> > > > >> >> > Sent: Wednesday, November 04, 2015 12:29 PM
> > > > >> >> > To: Jose Garcia-Rivera
> > > > >> >> > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > > > >> >> >
> > > > >> >> > Jose,
> > > > >> >> >
> > > > >> >> > I believe I've figured out the problem. The grid code
in
> > > > >> >> > MET currently only supports Lambert Conformal grids
where
> > > > >> >> > the grid spacing is the same in the x and y
directions.
> > > > >> >> >
> > > > >> >> > Your grid is spaced 3.317km in the x-direction and
2.97km
> > > > >> >> > in the y-direction. The MET code is reading the Dx
value
> > > > >> >> > and using it for both x and y. This causes the
skewing of
> > > > >> >> > the data to the north that
> > > > >> >> you're seeing.
> > > > >> >> >
> > > > >> >> > Really, there are 2 problems here.
> > > > >> >> >
> > > > >> >> > First, if we're requiring that Dx and Dy be the same
for
> > > > >> >> > Lambert Conformal grids, MET should error out with a
> > > > >> >> > descriptive error message when they are not.
> > > > >> >> >
> > > > >> >> > Second, obviously it would be preferable for MET to
> > > > >> >> > support Lambert Conformal grids where Dx and Dy are
not the
> same!
> > > > >> >> > I'll ask around here for more information on that.
> > > > >> >> >
> > > > >> >> > Thanks for paying close attention to your data and
> > > > >> >> > uncovering this shortcoming in MET.
> > > > >> >> >
> > > > >> >> > John
> > > > >> >> >
> > > > >> >> > On Wed, Nov 4, 2015 at 8:48 AM, Jose Garcia-Rivera via
RT
> > > > >> >> > < met_help at ucar.edu
> > > > >> >> > > wrote:
> > > > >> >> >
> > > > >> >> > >
> > > > >> >> > > <URL:
> > > > >> >> > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73975
> > > > >> >> > > >
> > > > >> >> > >
> > > > >> >> > > Hi John,
> > > > >> >> > >
> > > > >> >> > > Thanks for all the help. I want to add this pdf,
showing
> > > > >> >> > > something odd with the regridding. The domains did
not
> > > > >> >> > > match perfectly on the north side, it's like a
skewing
> > > > >> >> > > effect. The south sides match perfectly, but note
the
> > > > >> >> > > differences in the regridded RTMA file and the model
> > > > >> >> > > file in Michigan, for
> > > example.
> > > > >> >> > > Maybe that has something to
> > > > >> >> do with it.
> > > > >> >> > >
> > > > >> >> > > Best Regards,
> > > > >> >> > > Jose
> > > > >> >> > >
> > > > >> >> > > -----Original Message-----
> > > > >> >> > > From: John Halley Gotway via RT
> > > > >> >> > > [mailto:met_help at ucar.edu]
> > > > >> >> > > Sent: Tuesday, November 03, 2015 5:29 PM
> > > > >> >> > > To: Jose Garcia-Rivera
> > > > >> >> > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
problem
> > > > >> >> > >
> > > > >> >> > > Jose,
> > > > >> >> > >
> > > > >> >> > > I believe the problem is in our handling the domain
on
> > > > >> >> > > which your model data resides. I ran the following
> > > > >> >> > > commands to plot temperature from your model data
file:
> > > > >> >> > >
> > > > >> >> > > met-5.1/bin/plot_data_plane \
> > > > >> >> > > nmbprs_d01.00 nmbprs_d01.00_TMP_Z0.ps \
> > > > >> >> > > 'name="TMP"; level="Z0";'
> > > > >> >> > >
> > > > >> >> > > convert -rotate 90 -background white -flatten \
> > > > >> >> > > nmbprs_d01.00_TMP_Z0.ps nmbprs_d01.00_TMP_Z0.png
> > > > >> >> > >
> > > > >> >> > > The resulting image is attached and shows the shift
you
> > > > described.
> > > > >> >> > >
> > > > >> >> > > So it isn't the regridding step that's the
problem...
> > > > >> >> > > its our handling of that domain. I plotted the same
> > > > >> >> > > data using IDV and it puts things in the correct
place
> > > > >> >> > > (also
> > attached).
> > > > >> >> > >
> > > > >> >> > > One issue is that the 2-m temperature data in GRIB
> > > > >> >> > > record
> > > > >> >> > > 296 contains a constant value of 273. The surface
> > > > >> >> > > temperature data in GRIB record
> > > > >> >> > > 276 contains more reasonable values. So that's the
> > > > >> >> > > field I'm
> > > > >> >> > investigating.
> > > > >> >> > > But you might want to check why GRIB record 296
looks so
> odd.
> > > > >> >> > >
> > > > >> >> > > I need to head out of the office, but will
investigate
> > > > >> >> > > further
> > > > >> >> tomorrow.
> > > > >> >> > >
> > > > >> >> > > Thanks,
> > > > >> >> > > John
> > > > >> >> > >
> > > > >> >> > >
> > > > >> >> > >
> > > > >> >> > > On Tue, Nov 3, 2015 at 3:07 PM, Jose Garcia-Rivera
via
> > > > >> >> > > RT < met_help at ucar.edu
> > > > >> >> > > > wrote:
> > > > >> >> > >
> > > > >> >> > > >
> > > > >> >> > > > <URL:
> > > > >> >> > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=7397
> > > > >> >> > > > 5
> > > > >> >> > > > >
> > > > >> >> > > >
> > > > >> >> > > > Dear John,
> > > > >> >> > > >
> > > > >> >> > > > Thank you so much for the quick reply.
> > > > >> >> > > >
> > > > >> >> > > > I followed the instructions, and ftp'ed the files
to
> > > > >> >> > > "garcia-rivera_data".
> > > > >> >> > > >
> > > > >> >> > > > The nmbprs file is the .grb NMMB output file.
> > > > >> >> > > > The rtma file is the observational.
> > > > >> >> > > > The nmm_T2_bias.ncl subtracts the fields from both
> > > > >> >> > > > outputs in the same grid, and produces the plot.
> > > > >> >> > > >
> > > > >> >> > > > Any help will be greatly appreciated.
> > > > >> >> > > >
> > > > >> >> > > > Best Regards,
> > > > >> >> > > > Jose
> > > > >> >> > > >
> > > > >> >> > > > -----Original Message-----
> > > > >> >> > > > From: John Halley Gotway via RT
> > > > >> >> > > > [mailto:met_help at ucar.edu]
> > > > >> >> > > > Sent: Tuesday, November 03, 2015 4:04 PM
> > > > >> >> > > > To: Jose Garcia-Rivera
> > > > >> >> > > > Subject: Re: [rt.rap.ucar.edu #73975] Regridding
> > > > >> >> > > > problem
> > > > >> >> > > >
> > > > >> >> > > > Jose,
> > > > >> >> > > >
> > > > >> >> > > > Thanks for letting us know about the issue you're
seeing.
> > > > >> >> > > > I'd like to replicate the behavior your describe
using
> > > > >> >> > > > the exact data
> > > > >> >> you used.
> > > > >> >> > > > Once I can replicate it, I'll be able to
investigate
> > further.
> > > > >> >> > > >
> > > > >> >> > > > Are you able to send me the following 2 files you
used:
> > > > >> >> > > > rtma2p5.t23z.2dvaranl_nwrfc.grb2 and
nmbprs_d01.12
> > > > >> >> > > >
> > > > >> >> > > > You can post them to our anonymous ftp site,
following
> > > > >> >> > > > these
> > > > >> >> > > instructions:
> > > > >> >> > > >
> > > > >> >> > > >
http://www.dtcenter.org/met/users/support/met_help.php
> > > > >> >> > > > #f
> > > > >> >> > > > tp
> > > > >> >> > > >
> > > > >> >> > > > Please write back to let me know when you've
posted
> > > > >> >> > > > the data and I'll go grab it.
> > > > >> >> > > >
> > > > >> >> > > > Thanks,
> > > > >> >> > > > John Halley Gotway
> > > > >> >> > > > met_help at ucar.edu
> > > > >> >> > > >
> > > > >> >> > > >
> > > > >> >> > > >
> > > > >> >> > > > On Tue, Nov 3, 2015 at 1:53 PM, Jose Garcia-Rivera
via
> > > > >> >> > > > RT < met_help at ucar.edu
> > > > >> >> > > > > wrote:
> > > > >> >> > > >
> > > > >> >> > > > >
> > > > >> >> > > > > Tue Nov 03 13:53:37 2015: Request 73975 was
acted upon.
> > > > >> >> > > > > Transaction: Ticket created by GarciaJ at imsg.com
> > > > >> >> > > > > Queue: met_help
> > > > >> >> > > > > Subject: Regridding problem
> > > > >> >> > > > > Owner: Nobody
> > > > >> >> > > > > Requestors: GarciaJ at imsg.com
> > > > >> >> > > > > Status: new
> > > > >> >> > > > > Ticket <URL:
> > > > >> >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=73
> > > > >> >> > > > > 97
> > > > >> >> > > > > 5
> > > > >> >> > > > > >
> > > > >> >> > > > >
> > > > >> >> > > > >
> > > > >> >> > > > > Hello,
> > > > >> >> > > > >
> > > > >> >> > > > > I'm a new MET user, and I'm attempting to regrid
> > > > >> >> > > > > RTMA data to an NMMB output grid. I used the
following
> tool:
> > > > >> >> > > > >
> > > > >> >> > > > >
> > > > >> >> > > > >
> > > > >> >> > > > > ./regrid_data_plane
rtma2p5.t23z.2dvaranl_nwrfc.grb2
> > > > >> >> > > > > nmbprs_d01.12 regridded.nc -field 'name="TMP";
> > > level="Z2";'
> > > > >> >> > > > > -method BILIN -width
> > > > >> >> > > > > 2 -v 1
> > > > >> >> > > > >
> > > > >> >> > > > > which worked. But, when I subtracted the fields
from
> > > > >> >> > > > > one dataset to the other, there was a clear
offset
> > > > >> >> > > > > on the
> > > values.
> > > > >> >> > > > > So, the regridding was done indeed, but the
> > > > >> >> > > > > gridpoint values did not seem to be assigned
> > > > >> >> > > > > correctly with respect to their
> > > > >> counterparts.
> > > > >> >> > > > >
> > > > >> >> > > > > I attached the plot, so you can see the offset
problem.
> > > > >> >> > > > > Any help will be greatly appreciated.
> > > > >> >> > > > >
> > > > >> >> > > > > Best Regards,
> > > > >> >> > > > > Jose
> > > > >> >> > > > >
> > > > >> >> > > > > Jose M. Garcia-Rivera, PhD NWP Scientist I.M.
> > > > >> >> > > > > Systems Group This e-mail, including
attachments,
> > > > >> >> > > > > may include confidential and/or proprietary
> > > > >> >> > > > > information, and may be used only by the person
or
> > > > >> >> > > > > entity to which it is addressed. If the reader
of
> > > > >> >> > > > > this e-mail is not the intended recipient or his
or
> > > > >> >> > > > > her authorized agent, the reader is hereby
notified
> > > > >> >> > > > > that any dissemination, distribution or copying
of
> this e-mail is prohibited.
> > > > >> >> > > > > If you have received this e-mail in error,
please
> > > > >> >> > > > > notify the sender by replying to this message
and
> > > > >> >> > > > > delete this
> > > > >> >> > > > e-mail immediately.
> > > > >> >> > > > >
> > > > >> >> > > > >
> > > > >> >> > > >
> > > > >> >> > > > This e-mail, including attachments, may include
> > > > >> >> > > > confidential and/or proprietary information, and
may
> > > > >> >> > > > be used only by the person or entity to which it
is
> > > > >> >> > > > addressed. If the reader of this e-mail is not the
> > > > >> >> > > > intended recipient or his or her authorized agent,
the
> > > > >> >> > > > reader is hereby notified that any dissemination,
> > > > >> >> > > > distribution or copying of this e-mail is
prohibited.
> > > > >> >> > > > If you have received this e-mail in error, please
> > > > >> >> > > > notify the sender by replying to this message and
> > > > >> >> > > > delete this
> > > > >> >> > > e-mail immediately.
> > > > >> >> > > >
> > > > >> >> > > >
> > > > >> >> > > >
> > > > >> >> > >
> > > > >> >> > > This e-mail, including attachments, may include
> > > > >> >> > > confidential and/or proprietary information, and may
be
> > > > >> >> > > used only by the person or entity to which it is
> > > > >> >> > > addressed. If the reader of this e-mail is not the
> > > > >> >> > > intended recipient or his or her authorized agent,
the
> > > > >> >> > > reader is hereby notified that any dissemination,
> > > > >> >> > > distribution or copying of this e-mail is
prohibited. If
> > > > >> >> > > you have received this e-mail in error, please
notify
> > > > >> >> > > the sender by replying to this message and delete
this
> > > > >> >> > e-mail immediately.
> > > > >> >> > >
> > > > >> >> > >
> > > > >> >> >
> > > > >> >> > This e-mail, including attachments, may include
> > > > >> >> > confidential and/or proprietary information, and may
be
> > > > >> >> > used only by the person or entity to which it is
> > > > >> >> > addressed. If the reader of this e-mail is not the
> > > > >> >> > intended recipient or his or her authorized agent, the
> > > > >> >> > reader is hereby notified that any dissemination,
> > > > >> >> > distribution or copying of this e-mail is prohibited.
If
> > > > >> >> > you have received this e-mail in error, please notify
the
> > > > >> >> > sender by replying to this message and delete this
> > > > >> >> e-mail immediately.
> > > > >> >> >
> > > > >> >> >
> > > > >> >> >
> > > > >> >>
> > > > >> >> This e-mail, including attachments, may include
confidential
> > > > >> >> and/or proprietary information, and may be used only by
the
> > > > >> >> person or entity to which it is addressed. If the reader
of
> > > > >> >> this e-mail is not the intended recipient or his or her
> > > > >> >> authorized agent, the reader is hereby notified that any
> > > > >> >> dissemination, distribution or copying of this e-mail is
> > > > >> >> prohibited. If you have received this e-mail in error,
> > > > >> >> please notify the sender by replying to this message and
> > > > >> >> delete this
> > > > >> e-mail immediately.
> > > > >> >>
> > > > >> >>
> > > > >> >>
> > > > >> >
> > > > >>
> > > > >> This e-mail, including attachments, may include
confidential
> > > > >> and/or proprietary information, and may be used only by the
> > > > >> person or entity to which it is addressed. If the reader of
> > > > >> this e-mail is not the intended recipient or his or her
> > > > >> authorized agent, the reader is hereby notified that any
> > > > >> dissemination, distribution or copying of this e-mail is
> > > > >> prohibited. If you have received this e-mail in error,
please
> > > > >> notify the sender by replying to this message and delete
this
> > > > e-mail immediately.
> > > > >>
> > > > >>
> > > > >
> > > >
> > > > This e-mail, including attachments, may include confidential
> > > > and/or proprietary information, and may be used only by the
person
> > > > or entity to which it is addressed. If the reader of this e-
mail
> > > > is not the intended recipient or his or her authorized agent,
the
> > > > reader is hereby notified that any dissemination, distribution
or
> > > > copying of this e-mail is prohibited. If you have received
this
> > > > e-mail in error, please notify the sender by replying to this
> > > > message and delete this
> > > e-mail immediately.
> > > >
> > > >
> > > >
> > >
> > > This e-mail, including attachments, may include confidential
and/or
> > > proprietary information, and may be used only by the person or
> > > entity to which it is addressed. If the reader of this e-mail is
not
> > > the intended recipient or his or her authorized agent, the
reader is
> > > hereby notified that any dissemination, distribution or copying
of
> > > this e-mail is prohibited. If you have received this e-mail in
> > > error, please notify the sender by replying to this message and
> > > delete this
> > e-mail immediately.
> > >
> > >
> > >
> >
> > This e-mail, including attachments, may include confidential
and/or
> > proprietary information, and may be used only by the person or
entity
> > to which it is addressed. If the reader of this e-mail is not the
> > intended recipient or his or her authorized agent, the reader is
> > hereby notified that any dissemination, distribution or copying of
> > this e-mail is prohibited. If you have received this e-mail in
error,
> > please notify the sender by replying to this message and delete
this
> e-mail immediately.
> >
> >
> >
>
> This e-mail, including attachments, may include confidential and/or
> proprietary information, and may be used only by the person or
entity to
> which it is addressed. If the reader of this e-mail is not the
intended
> recipient or his or her authorized agent, the reader is hereby
notified
> that any dissemination, distribution or copying of this e-mail is
> prohibited. If you have received this e-mail in error, please notify
the
> sender by replying to this message and delete this e-mail
immediately.
>
>
>
------------------------------------------------
More information about the Met_help
mailing list