[Met_help] [rt.rap.ucar.edu #65680] History for Stage IV to model grid
John Halley Gotway via RT
met_help at ucar.edu
Mon Jun 2 14:18:57 MDT 2014
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Greetings:
Is there a way to regrid Stage IV data onto a local WRF grid using copygb
from UPP ? I also have wgrib and wgrib 2 installed (but not MET).
Cheers and best,
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies* (OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 11:23:36 2014
OK - I just figured it out - using:
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
and:
http://mailman.ucar.edu/pipermail/met_help/2009-December/001043.html
Cheers and many thanks,
On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT
<met_help at ucar.edu
> wrote:
> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
> "Stage IV to model grid",
> a summary of which appears below.
>
> There is no need to reply to this message right now. Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #65680].
>
> Please include the string:
>
> [rt.rap.ucar.edu #65680]
>
> in the subject line of all future correspondence about this issue.
To do
> so,
> you may reply to this message.
>
> Thank you,
> met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> Greetings:
>
> Is there a way to regrid Stage IV data onto a local WRF grid using
copygb
> from UPP ? I also have wgrib and wgrib 2 installed (but not MET).
>
> Cheers and best,
>
> --
> -------------------------------------------------------------
> Alexandre Fierro, PhD
> Research Scientist-
> National Severe Storms Laboratory (NSSL/NOAA)
> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>
> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
> why it is called the Present"
>
> "There are only 10 types of people in the world:
> Those who understand binary, and those who don't"
>
> "My opinions are my own and not representative of OU, NSSL,
> AOML, HRD, LANL or any affiliates."
> ^.^
> (o o)
> /( V )\
> ---m---m----
>
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65680] AutoReply: Stage IV to model grid
From: John Halley Gotway
Time: Tue Mar 04 11:27:03 2014
Alex,
Great, I'm glad you were able to figure it out. We have several
examples of running copygb for different projections in the METv4.1
online tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
The big limitation is that copygb only works on GRIB1 data. If you
have GRIB2 data, you can...
(1) Use cnvgrib to convert to GRIB1 and run copygb.
(2) Or wgrib2 does have the ability to regrid data but it'd pretty
limited.
Hope that helps.
Thanks,
John Halley Gotway
met_help at ucar.edu
On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>
> OK - I just figured it out - using:
>
>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>
> and:
>
> http://mailman.ucar.edu/pipermail/met_help/2009-December/001043.html
>
> Cheers and many thanks,
>
>
> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT
<met_help at ucar.edu
>> wrote:
>
>> Greetings,
>>
>> This message has been automatically generated in response to the
>> creation of a trouble ticket regarding:
>> "Stage IV to model grid",
>> a summary of which appears below.
>>
>> There is no need to reply to this message right now. Your ticket
has been
>> assigned an ID of [rt.rap.ucar.edu #65680].
>>
>> Please include the string:
>>
>> [rt.rap.ucar.edu #65680]
>>
>> in the subject line of all future correspondence about this issue.
To do
>> so,
>> you may reply to this message.
>>
>> Thank you,
>> met_help at ucar.edu
>>
>>
-------------------------------------------------------------------------
>> Greetings:
>>
>> Is there a way to regrid Stage IV data onto a local WRF grid using
copygb
>> from UPP ? I also have wgrib and wgrib 2 installed (but not MET).
>>
>> Cheers and best,
>>
>> --
>> -------------------------------------------------------------
>> Alexandre Fierro, PhD
>> Research Scientist-
>> National Severe Storms Laboratory (NSSL/NOAA)
>> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>
>> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
>> why it is called the Present"
>>
>> "There are only 10 types of people in the world:
>> Those who understand binary, and those who don't"
>>
>> "My opinions are my own and not representative of OU, NSSL,
>> AOML, HRD, LANL or any affiliates."
>> ^.^
>> (o o)
>> /( V )\
>> ---m---m----
>>
>>
>
>
------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 13:06:11 2014
Dear John:
Many thanks for your prompt reply and for the useful tips.
While I have your here, I'd have the following query: Is there a way
to
run copygb but on *only one single record* (i.e., 1 variable) within
the
WRF Netcdf file? I am currently only interested in APCP and wanted to
know
if copybg had the ability to only focus its interpolation magic on one
record similar to wgrib e.g. -which would save a lot of cpu time:
wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
where 279 is record# 279=APCP in the WRF GRIB1 converted file.
Cheers and best,
Alex
On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Alex,
>
> Great, I'm glad you were able to figure it out. We have several
examples
> of running copygb for different projections in the METv4.1 online
tutorial:
>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>
> The big limitation is that copygb only works on GRIB1 data. If you
have
> GRIB2 data, you can...
> (1) Use cnvgrib to convert to GRIB1 and run copygb.
> (2) Or wgrib2 does have the ability to regrid data but it'd pretty
limited.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >
> > OK - I just figured it out - using:
> >
> >
>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
> >
> > and:
> >
> > http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
> >
> > Cheers and many thanks,
> >
> >
> > On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
> met_help at ucar.edu
> >> wrote:
> >
> >> Greetings,
> >>
> >> This message has been automatically generated in response to the
> >> creation of a trouble ticket regarding:
> >> "Stage IV to model grid",
> >> a summary of which appears below.
> >>
> >> There is no need to reply to this message right now. Your ticket
has
> been
> >> assigned an ID of [rt.rap.ucar.edu #65680].
> >>
> >> Please include the string:
> >>
> >> [rt.rap.ucar.edu #65680]
> >>
> >> in the subject line of all future correspondence about this
issue. To do
> >> so,
> >> you may reply to this message.
> >>
> >> Thank you,
> >> met_help at ucar.edu
> >>
> >>
>
-------------------------------------------------------------------------
> >> Greetings:
> >>
> >> Is there a way to regrid Stage IV data onto a local WRF grid
using
> copygb
> >> from UPP ? I also have wgrib and wgrib 2 installed (but not MET).
> >>
> >> Cheers and best,
> >>
> >> --
> >> -------------------------------------------------------------
> >> Alexandre Fierro, PhD
> >> Research Scientist-
> >> National Severe Storms Laboratory (NSSL/NOAA)
> >> *The Cooperative Institute for Mesoscale Meteorological Studies*
> (OU/NOAA)
> >> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>
> >> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That
> is
> >> why it is called the Present"
> >>
> >> "There are only 10 types of people in the world:
> >> Those who understand binary, and those who don't"
> >>
> >> "My opinions are my own and not representative of OU, NSSL,
> >> AOML, HRD, LANL or any affiliates."
> >> ^.^
> >> (o o)
> >> /( V )\
> >> ---m---m----
> >>
> >>
> >
> >
>
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 13:19:44 2014
Seems that again I found the answer to my own question:
I just changed the default wrf_cntrl.parm to only contain:
1 KGTYPE******I5*******:(00255)********START OF THIS OUTPUT
GRID******
2 IMDLTY *I5* :(00125)
3 DATSET *A6* :(WRFPRS)
4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000
00000 00000 00000)
correct ?
On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro
<alexandre.o.fierro at gmail.com>wrote:
> Dear John:
>
> Many thanks for your prompt reply and for the useful tips.
>
> While I have your here, I'd have the following query: Is there a
way to
> run copygb but on *only one single record* (i.e., 1 variable) within
the
> WRF Netcdf file? I am currently only interested in APCP and wanted
to know
> if copybg had the ability to only focus its interpolation magic on
one
> record similar to wgrib e.g. -which would save a lot of cpu time:
>
> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
>
> where 279 is record# 279=APCP in the WRF GRIB1 converted file.
>
> Cheers and best,
>
> Alex
>
>
> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Alex,
>>
>> Great, I'm glad you were able to figure it out. We have several
examples
>> of running copygb for different projections in the METv4.1 online
tutorial:
>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>>
>> The big limitation is that copygb only works on GRIB1 data. If you
have
>> GRIB2 data, you can...
>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
>> (2) Or wgrib2 does have the ability to regrid data but it'd pretty
>> limited.
>>
>> Hope that helps.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>> >
>> > OK - I just figured it out - using:
>> >
>> >
>>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>> >
>> > and:
>> >
>> > http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
>> >
>> > Cheers and many thanks,
>> >
>> >
>> > On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
>> met_help at ucar.edu
>> >> wrote:
>> >
>> >> Greetings,
>> >>
>> >> This message has been automatically generated in response to the
>> >> creation of a trouble ticket regarding:
>> >> "Stage IV to model grid",
>> >> a summary of which appears below.
>> >>
>> >> There is no need to reply to this message right now. Your
ticket has
>> been
>> >> assigned an ID of [rt.rap.ucar.edu #65680].
>> >>
>> >> Please include the string:
>> >>
>> >> [rt.rap.ucar.edu #65680]
>> >>
>> >> in the subject line of all future correspondence about this
issue. To
>> do
>> >> so,
>> >> you may reply to this message.
>> >>
>> >> Thank you,
>> >> met_help at ucar.edu
>> >>
>> >>
>>
-------------------------------------------------------------------------
>> >> Greetings:
>> >>
>> >> Is there a way to regrid Stage IV data onto a local WRF grid
using
>> copygb
>> >> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
>> >>
>> >> Cheers and best,
>> >>
>> >> --
>> >> -------------------------------------------------------------
>> >> Alexandre Fierro, PhD
>> >> Research Scientist-
>> >> National Severe Storms Laboratory (NSSL/NOAA)
>> >> *The Cooperative Institute for Mesoscale Meteorological Studies*
>> (OU/NOAA)
>> >> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>> >>
>> >> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift; That
>> is
>> >> why it is called the Present"
>> >>
>> >> "There are only 10 types of people in the world:
>> >> Those who understand binary, and those who don't"
>> >>
>> >> "My opinions are my own and not representative of OU, NSSL,
>> >> AOML, HRD, LANL or any affiliates."
>> >> ^.^
>> >> (o o)
>> >> /( V )\
>> >> ---m---m----
>> >>
>> >>
>> >
>> >
>>
>>
>
>
> --
> -------------------------------------------------------------
> Alexandre Fierro, PhD
> Research Scientist-
> National Severe Storms Laboratory (NSSL/NOAA)
> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
>
> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>
> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
> why it is called the Present"
>
> "There are only 10 types of people in the world:
> Those who understand binary, and those who don't"
>
> "My opinions are my own and not representative of OU, NSSL,
> AOML, HRD, LANL or any affiliates."
> ^.^
> (o o)
> /( V )\
> ---m---m----
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65680] AutoReply: Stage IV to model grid
From: John Halley Gotway
Time: Tue Mar 04 13:29:05 2014
Alex,
The wrfcntrl.parm file tells the unified post-processor what fields
you'd like post-processed. You could configure it to only dump out a
single GRIB record, but typically users want more than one
output field. The unified post-processor (UPP) and copygb are not the
same thing. With wrfcntrl.parm, you're controlling UPP. But you
could also run copygb to regrid only a subset of records
available in the input GRIB file. It's up to you to decide what makes
the most sense for your project.
You can configure copygb to run on a subset of the GRIB records using
the "-k" option. You can read copygb.doc for details
(http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt)
and see the examples at the bottom. Example (2) uses "-k'4*-1 61'" to
select out precipitation whose GRIB number
is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out 250mb
height.
But I think the most powerful way of doing it is using wgrib and
copygb together as shown in example (8). The "-kw" option tells
copygb to read the GRIB record numbers from standard input. For
example:
wgrib sample.grb | grep RH | copygb -xg212 -kw sample_g212.grb
This would filter out the GRIB records for RH, regrid them to NCEP
grid 212, and write the output in the file named sample_g212.grb.
wgrib is also pretty useful. You can use it to easily select one or
more records from an input GRIB file. But if you need to interpolate
to a new grid, you should use copygb.
Hope that helps.
John
On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>
> Seems that again I found the answer to my own question:
>
> I just changed the default wrf_cntrl.parm to only contain:
>
> 1 KGTYPE******I5*******:(00255)********START OF THIS OUTPUT
> GRID******
> 2 IMDLTY *I5* :(00125)
> 3 DATSET *A6* :(WRFPRS)
> 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
> 5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000
> 00000 00000 00000)
>
> correct ?
>
>
>
>
> On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro
<alexandre.o.fierro at gmail.com>wrote:
>
>> Dear John:
>>
>> Many thanks for your prompt reply and for the useful tips.
>>
>> While I have your here, I'd have the following query: Is there a
way to
>> run copygb but on *only one single record* (i.e., 1 variable)
within the
>> WRF Netcdf file? I am currently only interested in APCP and wanted
to know
>> if copybg had the ability to only focus its interpolation magic on
one
>> record similar to wgrib e.g. -which would save a lot of cpu time:
>>
>> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
>>
>> where 279 is record# 279=APCP in the WRF GRIB1 converted file.
>>
>> Cheers and best,
>>
>> Alex
>>
>>
>> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
>> met_help at ucar.edu> wrote:
>>
>>> Alex,
>>>
>>> Great, I'm glad you were able to figure it out. We have several
examples
>>> of running copygb for different projections in the METv4.1 online
tutorial:
>>>
>>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>>>
>>> The big limitation is that copygb only works on GRIB1 data. If
you have
>>> GRIB2 data, you can...
>>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
>>> (2) Or wgrib2 does have the ability to regrid data but it'd pretty
>>> limited.
>>>
>>> Hope that helps.
>>>
>>> Thanks,
>>> John Halley Gotway
>>> met_help at ucar.edu
>>>
>>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>>>>
>>>> OK - I just figured it out - using:
>>>>
>>>>
>>>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>>>>
>>>> and:
>>>>
>>>> http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
>>>>
>>>> Cheers and many thanks,
>>>>
>>>>
>>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
>>> met_help at ucar.edu
>>>>> wrote:
>>>>
>>>>> Greetings,
>>>>>
>>>>> This message has been automatically generated in response to the
>>>>> creation of a trouble ticket regarding:
>>>>> "Stage IV to model grid",
>>>>> a summary of which appears below.
>>>>>
>>>>> There is no need to reply to this message right now. Your
ticket has
>>> been
>>>>> assigned an ID of [rt.rap.ucar.edu #65680].
>>>>>
>>>>> Please include the string:
>>>>>
>>>>> [rt.rap.ucar.edu #65680]
>>>>>
>>>>> in the subject line of all future correspondence about this
issue. To
>>> do
>>>>> so,
>>>>> you may reply to this message.
>>>>>
>>>>> Thank you,
>>>>> met_help at ucar.edu
>>>>>
>>>>>
>>>
-------------------------------------------------------------------------
>>>>> Greetings:
>>>>>
>>>>> Is there a way to regrid Stage IV data onto a local WRF grid
using
>>> copygb
>>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
>>>>>
>>>>> Cheers and best,
>>>>>
>>>>> --
>>>>> -------------------------------------------------------------
>>>>> Alexandre Fierro, PhD
>>>>> Research Scientist-
>>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>>> *The Cooperative Institute for Mesoscale Meteorological Studies*
>>> (OU/NOAA)
>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>>
>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift; That
>>> is
>>>>> why it is called the Present"
>>>>>
>>>>> "There are only 10 types of people in the world:
>>>>> Those who understand binary, and those who don't"
>>>>>
>>>>> "My opinions are my own and not representative of OU, NSSL,
>>>>> AOML, HRD, LANL or any affiliates."
>>>>> ^.^
>>>>> (o o)
>>>>> /( V )\
>>>>> ---m---m----
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>> --
>> -------------------------------------------------------------
>> Alexandre Fierro, PhD
>> Research Scientist-
>> National Severe Storms Laboratory (NSSL/NOAA)
>> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
>>
>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>
>> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
>> why it is called the Present"
>>
>> "There are only 10 types of people in the world:
>> Those who understand binary, and those who don't"
>>
>> "My opinions are my own and not representative of OU, NSSL,
>> AOML, HRD, LANL or any affiliates."
>> ^.^
>> (o o)
>> /( V )\
>> ---m---m----
>>
>
>
>
------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 13:44:23 2014
Hi John:
Once again, many thanks for your prompt reply and useful tips !
I did not specify clearly my goal, my mistake.
I just realized that I might not need to run unipost if the following
is
possible: Is it possible to convert the WRF NETCDF data for one single
record (APCP=rec#279) to a GRIB1 file ?
Cheers and many thanks for your help,
Alex-
On Tue, Mar 4, 2014 at 2:29 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Alex,
>
> The wrfcntrl.parm file tells the unified post-processor what fields
you'd
> like post-processed. You could configure it to only dump out a
single GRIB
> record, but typically users want more than one
> output field. The unified post-processor (UPP) and copygb are not
the
> same thing. With wrfcntrl.parm, you're controlling UPP. But you
could
> also run copygb to regrid only a subset of records
> available in the input GRIB file. It's up to you to decide what
makes the
> most sense for your project.
>
> You can configure copygb to run on a subset of the GRIB records
using the
> "-k" option. You can read copygb.doc for details
> (
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt)
> and see the examples at the bottom. Example (2) uses "-k'4*-1 61'"
to
> select out precipitation whose GRIB number
> is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out 250mb
height.
>
> But I think the most powerful way of doing it is using wgrib and
copygb
> together as shown in example (8). The "-kw" option tells copygb to
read
> the GRIB record numbers from standard input. For example:
> wgrib sample.grb | grep RH | copygb -xg212 -kw sample_g212.grb
>
> This would filter out the GRIB records for RH, regrid them to NCEP
grid
> 212, and write the output in the file named sample_g212.grb.
>
> wgrib is also pretty useful. You can use it to easily select one or
more
> records from an input GRIB file. But if you need to interpolate to
a new
> grid, you should use copygb.
>
> Hope that helps.
>
> John
>
> On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >
> > Seems that again I found the answer to my own question:
> >
> > I just changed the default wrf_cntrl.parm to only contain:
> >
> > 1 KGTYPE******I5*******:(00255)********START OF THIS OUTPUT
> > GRID******
> > 2 IMDLTY *I5* :(00125)
> > 3 DATSET *A6* :(WRFPRS)
> > 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
> > 5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000
00000
> 00000
> > 00000 00000 00000)
> >
> > correct ?
> >
> >
> >
> >
> > On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro <
> alexandre.o.fierro at gmail.com>wrote:
> >
> >> Dear John:
> >>
> >> Many thanks for your prompt reply and for the useful tips.
> >>
> >> While I have your here, I'd have the following query: Is there a
way to
> >> run copygb but on *only one single record* (i.e., 1 variable)
within the
> >> WRF Netcdf file? I am currently only interested in APCP and
wanted to
> know
> >> if copybg had the ability to only focus its interpolation magic
on one
> >> record similar to wgrib e.g. -which would save a lot of cpu
time:
> >>
> >> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
> >>
> >> where 279 is record# 279=APCP in the WRF GRIB1 converted file.
> >>
> >> Cheers and best,
> >>
> >> Alex
> >>
> >>
> >> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >>> Alex,
> >>>
> >>> Great, I'm glad you were able to figure it out. We have several
> examples
> >>> of running copygb for different projections in the METv4.1
online
> tutorial:
> >>>
> >>>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
> >>>
> >>> The big limitation is that copygb only works on GRIB1 data. If
you
> have
> >>> GRIB2 data, you can...
> >>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
> >>> (2) Or wgrib2 does have the ability to regrid data but it'd
pretty
> >>> limited.
> >>>
> >>> Hope that helps.
> >>>
> >>> Thanks,
> >>> John Halley Gotway
> >>> met_help at ucar.edu
> >>>
> >>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
> >>>>
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >>>>
> >>>> OK - I just figured it out - using:
> >>>>
> >>>>
> >>>
>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
> >>>>
> >>>> and:
> >>>>
> >>>> http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
> >>>>
> >>>> Cheers and many thanks,
> >>>>
> >>>>
> >>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
> >>> met_help at ucar.edu
> >>>>> wrote:
> >>>>
> >>>>> Greetings,
> >>>>>
> >>>>> This message has been automatically generated in response to
the
> >>>>> creation of a trouble ticket regarding:
> >>>>> "Stage IV to model grid",
> >>>>> a summary of which appears below.
> >>>>>
> >>>>> There is no need to reply to this message right now. Your
ticket has
> >>> been
> >>>>> assigned an ID of [rt.rap.ucar.edu #65680].
> >>>>>
> >>>>> Please include the string:
> >>>>>
> >>>>> [rt.rap.ucar.edu #65680]
> >>>>>
> >>>>> in the subject line of all future correspondence about this
issue. To
> >>> do
> >>>>> so,
> >>>>> you may reply to this message.
> >>>>>
> >>>>> Thank you,
> >>>>> met_help at ucar.edu
> >>>>>
> >>>>>
> >>>
>
-------------------------------------------------------------------------
> >>>>> Greetings:
> >>>>>
> >>>>> Is there a way to regrid Stage IV data onto a local WRF grid
using
> >>> copygb
> >>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
> >>>>>
> >>>>> Cheers and best,
> >>>>>
> >>>>> --
> >>>>> -------------------------------------------------------------
> >>>>> Alexandre Fierro, PhD
> >>>>> Research Scientist-
> >>>>> National Severe Storms Laboratory (NSSL/NOAA)
> >>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
> >>> (OU/NOAA)
> >>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>>>>
> >>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
> That
> >>> is
> >>>>> why it is called the Present"
> >>>>>
> >>>>> "There are only 10 types of people in the world:
> >>>>> Those who understand binary, and those who don't"
> >>>>>
> >>>>> "My opinions are my own and not representative of OU, NSSL,
> >>>>> AOML, HRD, LANL or any affiliates."
> >>>>> ^.^
> >>>>> (o o)
> >>>>> /( V )\
> >>>>> ---m---m----
> >>>>>
> >>>>>
> >>>>
> >>>>
> >>>
> >>>
> >>
> >>
> >> --
> >> -------------------------------------------------------------
> >> Alexandre Fierro, PhD
> >> Research Scientist-
> >> National Severe Storms Laboratory (NSSL/NOAA)
> >> *The Cooperative Institute for Mesoscale Meteorological Studies*
> (OU/NOAA)
> >>
> >> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>
> >> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That
> is
> >> why it is called the Present"
> >>
> >> "There are only 10 types of people in the world:
> >> Those who understand binary, and those who don't"
> >>
> >> "My opinions are my own and not representative of OU, NSSL,
> >> AOML, HRD, LANL or any affiliates."
> >> ^.^
> >> (o o)
> >> /( V )\
> >> ---m---m----
> >>
> >
> >
> >
>
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 13:53:12 2014
Hi John:
I just checked and what I initially did was correct: Changing the
../parm/wrf_cntrl.parm and running unipost.exe just does that. The
reason
why I am surprised is that the process (unipost.exe on 8 procs)
remains
painfully slow - Perhaps this is due to the fact that unipost has to
read
the entire netcdf record prior to dumping the data? Is there a way to
accelerate the process ?
i) convert wrf.netcdf into a GRIB1 file for only one record (APCP) -
no
interp of the WRF data is necessary as I actually interp the stage IV
data
onto the WRF grid using copygb later on.
Cheers and many thanks again for your patience,
Alex-
On Tue, Mar 4, 2014 at 2:43 PM, Alex Fierro
<alexandre.o.fierro at gmail.com>wrote:
> Hi John:
>
> Once again, many thanks for your prompt reply and useful tips !
>
> I did not specify clearly my goal, my mistake.
>
> I just realized that I might not need to run unipost if the
following is
> possible: Is it possible to convert the WRF NETCDF data for one
single
> record (APCP=rec#279) to a GRIB1 file ?
>
> Cheers and many thanks for your help,
>
> Alex-
>
>
>
>
> On Tue, Mar 4, 2014 at 2:29 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Alex,
>>
>> The wrfcntrl.parm file tells the unified post-processor what fields
you'd
>> like post-processed. You could configure it to only dump out a
single GRIB
>> record, but typically users want more than one
>> output field. The unified post-processor (UPP) and copygb are not
the
>> same thing. With wrfcntrl.parm, you're controlling UPP. But you
could
>> also run copygb to regrid only a subset of records
>> available in the input GRIB file. It's up to you to decide what
makes
>> the most sense for your project.
>>
>> You can configure copygb to run on a subset of the GRIB records
using the
>> "-k" option. You can read copygb.doc for details
>> (
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt)
>> and see the examples at the bottom. Example (2) uses "-k'4*-1 61'"
to
>> select out precipitation whose GRIB number
>> is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out 250mb
height.
>>
>> But I think the most powerful way of doing it is using wgrib and
copygb
>> together as shown in example (8). The "-kw" option tells copygb to
read
>> the GRIB record numbers from standard input. For example:
>> wgrib sample.grb | grep RH | copygb -xg212 -kw sample_g212.grb
>>
>> This would filter out the GRIB records for RH, regrid them to NCEP
grid
>> 212, and write the output in the file named sample_g212.grb.
>>
>> wgrib is also pretty useful. You can use it to easily select one
or more
>> records from an input GRIB file. But if you need to interpolate to
a new
>> grid, you should use copygb.
>>
>> Hope that helps.
>>
>> John
>>
>> On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>> >
>> > Seems that again I found the answer to my own question:
>> >
>> > I just changed the default wrf_cntrl.parm to only contain:
>> >
>> > 1 KGTYPE******I5*******:(00255)********START OF THIS
OUTPUT
>> > GRID******
>> > 2 IMDLTY *I5* :(00125)
>> > 3 DATSET *A6* :(WRFPRS)
>> > 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
>> > 5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000
00000
>> 00000
>> > 00000 00000 00000)
>> >
>> > correct ?
>> >
>> >
>> >
>> >
>> > On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro <
>> alexandre.o.fierro at gmail.com>wrote:
>> >
>> >> Dear John:
>> >>
>> >> Many thanks for your prompt reply and for the useful tips.
>> >>
>> >> While I have your here, I'd have the following query: Is there
a way
>> to
>> >> run copygb but on *only one single record* (i.e., 1 variable)
within
>> the
>> >> WRF Netcdf file? I am currently only interested in APCP and
wanted to
>> know
>> >> if copybg had the ability to only focus its interpolation magic
on one
>> >> record similar to wgrib e.g. -which would save a lot of cpu
time:
>> >>
>> >> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
>> >>
>> >> where 279 is record# 279=APCP in the WRF GRIB1 converted file.
>> >>
>> >> Cheers and best,
>> >>
>> >> Alex
>> >>
>> >>
>> >> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
>> >> met_help at ucar.edu> wrote:
>> >>
>> >>> Alex,
>> >>>
>> >>> Great, I'm glad you were able to figure it out. We have
several
>> examples
>> >>> of running copygb for different projections in the METv4.1
online
>> tutorial:
>> >>>
>> >>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>> >>>
>> >>> The big limitation is that copygb only works on GRIB1 data. If
you
>> have
>> >>> GRIB2 data, you can...
>> >>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
>> >>> (2) Or wgrib2 does have the ability to regrid data but it'd
pretty
>> >>> limited.
>> >>>
>> >>> Hope that helps.
>> >>>
>> >>> Thanks,
>> >>> John Halley Gotway
>> >>> met_help at ucar.edu
>> >>>
>> >>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>> >>>>
>> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680
>
>> >>>>
>> >>>> OK - I just figured it out - using:
>> >>>>
>> >>>>
>> >>>
>>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>> >>>>
>> >>>> and:
>> >>>>
>> >>>> http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
>> >>>>
>> >>>> Cheers and many thanks,
>> >>>>
>> >>>>
>> >>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
>> >>> met_help at ucar.edu
>> >>>>> wrote:
>> >>>>
>> >>>>> Greetings,
>> >>>>>
>> >>>>> This message has been automatically generated in response to
the
>> >>>>> creation of a trouble ticket regarding:
>> >>>>> "Stage IV to model grid",
>> >>>>> a summary of which appears below.
>> >>>>>
>> >>>>> There is no need to reply to this message right now. Your
ticket
>> has
>> >>> been
>> >>>>> assigned an ID of [rt.rap.ucar.edu #65680].
>> >>>>>
>> >>>>> Please include the string:
>> >>>>>
>> >>>>> [rt.rap.ucar.edu #65680]
>> >>>>>
>> >>>>> in the subject line of all future correspondence about this
issue.
>> To
>> >>> do
>> >>>>> so,
>> >>>>> you may reply to this message.
>> >>>>>
>> >>>>> Thank you,
>> >>>>> met_help at ucar.edu
>> >>>>>
>> >>>>>
>> >>>
>>
-------------------------------------------------------------------------
>> >>>>> Greetings:
>> >>>>>
>> >>>>> Is there a way to regrid Stage IV data onto a local WRF grid
using
>> >>> copygb
>> >>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
>> >>>>>
>> >>>>> Cheers and best,
>> >>>>>
>> >>>>> --
>> >>>>> -------------------------------------------------------------
>> >>>>> Alexandre Fierro, PhD
>> >>>>> Research Scientist-
>> >>>>> National Severe Storms Laboratory (NSSL/NOAA)
>> >>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
>> >>> (OU/NOAA)
>> >>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>> >>>>>
>> >>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
>> That
>> >>> is
>> >>>>> why it is called the Present"
>> >>>>>
>> >>>>> "There are only 10 types of people in the world:
>> >>>>> Those who understand binary, and those who don't"
>> >>>>>
>> >>>>> "My opinions are my own and not representative of OU, NSSL,
>> >>>>> AOML, HRD, LANL or any affiliates."
>> >>>>> ^.^
>> >>>>> (o o)
>> >>>>> /( V )\
>> >>>>> ---m---m----
>> >>>>>
>> >>>>>
>> >>>>
>> >>>>
>> >>>
>> >>>
>> >>
>> >>
>> >> --
>> >> -------------------------------------------------------------
>> >> Alexandre Fierro, PhD
>> >> Research Scientist-
>> >> National Severe Storms Laboratory (NSSL/NOAA)
>> >> *The Cooperative Institute for Mesoscale Meteorological Studies*
>> (OU/NOAA)
>> >>
>> >> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>> >>
>> >> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift; That
>> is
>> >> why it is called the Present"
>> >>
>> >> "There are only 10 types of people in the world:
>> >> Those who understand binary, and those who don't"
>> >>
>> >> "My opinions are my own and not representative of OU, NSSL,
>> >> AOML, HRD, LANL or any affiliates."
>> >> ^.^
>> >> (o o)
>> >> /( V )\
>> >> ---m---m----
>> >>
>> >
>> >
>> >
>>
>>
>
>
> --
> -------------------------------------------------------------
> Alexandre Fierro, PhD
> Research Scientist-
> National Severe Storms Laboratory (NSSL/NOAA)
> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
>
> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>
> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
> why it is called the Present"
>
> "There are only 10 types of people in the world:
> Those who understand binary, and those who don't"
>
> "My opinions are my own and not representative of OU, NSSL,
> AOML, HRD, LANL or any affiliates."
> ^.^
> (o o)
> /( V )\
> ---m---m----
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65680] AutoReply: Stage IV to model grid
From: John Halley Gotway
Time: Tue Mar 04 13:58:38 2014
Alex,
I concur - configuring wrf_cntrl.parm to only dump out a single record
will do what you want.
As to why it's so slow, I have no idea. And I'm really not the right
person to answer this question. Support for UPP is provided through
wrfhelp at ucar.edu. You could ask them if there's any way to
expedite the postprocessing of a single record.
Assuming you're running WRF-ARW, another option would be running the
pinterp utility
(http://www.mmm.ucar.edu/wrf/OnLineTutorial/Tools/p_interp.htm) to
post-process your output. That may (or may not)
be faster. MET can read the pinterp output, but I don't know if
there's a way to only post-process a single record.
Thanks,
John
On 03/04/2014 01:53 PM, Alex Fierro via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>
> Hi John:
>
> I just checked and what I initially did was correct: Changing the
> ../parm/wrf_cntrl.parm and running unipost.exe just does that. The
reason
> why I am surprised is that the process (unipost.exe on 8 procs)
remains
> painfully slow - Perhaps this is due to the fact that unipost has to
read
> the entire netcdf record prior to dumping the data? Is there a way
to
> accelerate the process ?
>
> i) convert wrf.netcdf into a GRIB1 file for only one record (APCP) -
no
> interp of the WRF data is necessary as I actually interp the stage
IV data
> onto the WRF grid using copygb later on.
>
> Cheers and many thanks again for your patience,
>
> Alex-
>
>
> On Tue, Mar 4, 2014 at 2:43 PM, Alex Fierro
<alexandre.o.fierro at gmail.com>wrote:
>
>> Hi John:
>>
>> Once again, many thanks for your prompt reply and useful tips !
>>
>> I did not specify clearly my goal, my mistake.
>>
>> I just realized that I might not need to run unipost if the
following is
>> possible: Is it possible to convert the WRF NETCDF data for one
single
>> record (APCP=rec#279) to a GRIB1 file ?
>>
>> Cheers and many thanks for your help,
>>
>> Alex-
>>
>>
>>
>>
>> On Tue, Mar 4, 2014 at 2:29 PM, John Halley Gotway via RT <
>> met_help at ucar.edu> wrote:
>>
>>> Alex,
>>>
>>> The wrfcntrl.parm file tells the unified post-processor what
fields you'd
>>> like post-processed. You could configure it to only dump out a
single GRIB
>>> record, but typically users want more than one
>>> output field. The unified post-processor (UPP) and copygb are not
the
>>> same thing. With wrfcntrl.parm, you're controlling UPP. But you
could
>>> also run copygb to regrid only a subset of records
>>> available in the input GRIB file. It's up to you to decide what
makes
>>> the most sense for your project.
>>>
>>> You can configure copygb to run on a subset of the GRIB records
using the
>>> "-k" option. You can read copygb.doc for details
>>> (
>>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt)
>>> and see the examples at the bottom. Example (2) uses "-k'4*-1
61'" to
>>> select out precipitation whose GRIB number
>>> is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out 250mb
height.
>>>
>>> But I think the most powerful way of doing it is using wgrib and
copygb
>>> together as shown in example (8). The "-kw" option tells copygb
to read
>>> the GRIB record numbers from standard input. For example:
>>> wgrib sample.grb | grep RH | copygb -xg212 -kw sample_g212.grb
>>>
>>> This would filter out the GRIB records for RH, regrid them to NCEP
grid
>>> 212, and write the output in the file named sample_g212.grb.
>>>
>>> wgrib is also pretty useful. You can use it to easily select one
or more
>>> records from an input GRIB file. But if you need to interpolate
to a new
>>> grid, you should use copygb.
>>>
>>> Hope that helps.
>>>
>>> John
>>>
>>> On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>>>>
>>>> Seems that again I found the answer to my own question:
>>>>
>>>> I just changed the default wrf_cntrl.parm to only contain:
>>>>
>>>> 1 KGTYPE******I5*******:(00255)********START OF THIS
OUTPUT
>>>> GRID******
>>>> 2 IMDLTY *I5* :(00125)
>>>> 3 DATSET *A6* :(WRFPRS)
>>>> 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
>>>> 5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000
00000
>>> 00000
>>>> 00000 00000 00000)
>>>>
>>>> correct ?
>>>>
>>>>
>>>>
>>>>
>>>> On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro <
>>> alexandre.o.fierro at gmail.com>wrote:
>>>>
>>>>> Dear John:
>>>>>
>>>>> Many thanks for your prompt reply and for the useful tips.
>>>>>
>>>>> While I have your here, I'd have the following query: Is there
a way
>>> to
>>>>> run copygb but on *only one single record* (i.e., 1 variable)
within
>>> the
>>>>> WRF Netcdf file? I am currently only interested in APCP and
wanted to
>>> know
>>>>> if copybg had the ability to only focus its interpolation magic
on one
>>>>> record similar to wgrib e.g. -which would save a lot of cpu
time:
>>>>>
>>>>> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
>>>>>
>>>>> where 279 is record# 279=APCP in the WRF GRIB1 converted file.
>>>>>
>>>>> Cheers and best,
>>>>>
>>>>> Alex
>>>>>
>>>>>
>>>>> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Alex,
>>>>>>
>>>>>> Great, I'm glad you were able to figure it out. We have
several
>>> examples
>>>>>> of running copygb for different projections in the METv4.1
online
>>> tutorial:
>>>>>>
>>>>>>
>>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>>>>>>
>>>>>> The big limitation is that copygb only works on GRIB1 data. If
you
>>> have
>>>>>> GRIB2 data, you can...
>>>>>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
>>>>>> (2) Or wgrib2 does have the ability to regrid data but it'd
pretty
>>>>>> limited.
>>>>>>
>>>>>> Hope that helps.
>>>>>>
>>>>>> Thanks,
>>>>>> John Halley Gotway
>>>>>> met_help at ucar.edu
>>>>>>
>>>>>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680
>
>>>>>>>
>>>>>>> OK - I just figured it out - using:
>>>>>>>
>>>>>>>
>>>>>>
>>>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>>>>>>>
>>>>>>> and:
>>>>>>>
>>>>>>> http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
>>>>>>>
>>>>>>> Cheers and many thanks,
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
>>>>>> met_help at ucar.edu
>>>>>>>> wrote:
>>>>>>>
>>>>>>>> Greetings,
>>>>>>>>
>>>>>>>> This message has been automatically generated in response to
the
>>>>>>>> creation of a trouble ticket regarding:
>>>>>>>> "Stage IV to model grid",
>>>>>>>> a summary of which appears below.
>>>>>>>>
>>>>>>>> There is no need to reply to this message right now. Your
ticket
>>> has
>>>>>> been
>>>>>>>> assigned an ID of [rt.rap.ucar.edu #65680].
>>>>>>>>
>>>>>>>> Please include the string:
>>>>>>>>
>>>>>>>> [rt.rap.ucar.edu #65680]
>>>>>>>>
>>>>>>>> in the subject line of all future correspondence about this
issue.
>>> To
>>>>>> do
>>>>>>>> so,
>>>>>>>> you may reply to this message.
>>>>>>>>
>>>>>>>> Thank you,
>>>>>>>> met_help at ucar.edu
>>>>>>>>
>>>>>>>>
>>>>>>
>>>
-------------------------------------------------------------------------
>>>>>>>> Greetings:
>>>>>>>>
>>>>>>>> Is there a way to regrid Stage IV data onto a local WRF grid
using
>>>>>> copygb
>>>>>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
>>>>>>>>
>>>>>>>> Cheers and best,
>>>>>>>>
>>>>>>>> --
>>>>>>>> -------------------------------------------------------------
>>>>>>>> Alexandre Fierro, PhD
>>>>>>>> Research Scientist-
>>>>>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>>>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
>>>>>> (OU/NOAA)
>>>>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>>>>>
>>>>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
>>> That
>>>>>> is
>>>>>>>> why it is called the Present"
>>>>>>>>
>>>>>>>> "There are only 10 types of people in the world:
>>>>>>>> Those who understand binary, and those who don't"
>>>>>>>>
>>>>>>>> "My opinions are my own and not representative of OU, NSSL,
>>>>>>>> AOML, HRD, LANL or any affiliates."
>>>>>>>> ^.^
>>>>>>>> (o o)
>>>>>>>> /( V )\
>>>>>>>> ---m---m----
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> -------------------------------------------------------------
>>>>> Alexandre Fierro, PhD
>>>>> Research Scientist-
>>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>>> *The Cooperative Institute for Mesoscale Meteorological Studies*
>>> (OU/NOAA)
>>>>>
>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>>
>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift; That
>>> is
>>>>> why it is called the Present"
>>>>>
>>>>> "There are only 10 types of people in the world:
>>>>> Those who understand binary, and those who don't"
>>>>>
>>>>> "My opinions are my own and not representative of OU, NSSL,
>>>>> AOML, HRD, LANL or any affiliates."
>>>>> ^.^
>>>>> (o o)
>>>>> /( V )\
>>>>> ---m---m----
>>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>> --
>> -------------------------------------------------------------
>> Alexandre Fierro, PhD
>> Research Scientist-
>> National Severe Storms Laboratory (NSSL/NOAA)
>> *The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
>>
>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>
>> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That is
>> why it is called the Present"
>>
>> "There are only 10 types of people in the world:
>> Those who understand binary, and those who don't"
>>
>> "My opinions are my own and not representative of OU, NSSL,
>> AOML, HRD, LANL or any affiliates."
>> ^.^
>> (o o)
>> /( V )\
>> ---m---m----
>>
>
>
>
------------------------------------------------
Subject: Stage IV to model grid
From: Alex Fierro
Time: Tue Mar 04 15:22:53 2014
Dear John:
OK - Many thanks once again for your prompt reply. Much appreciated. I
guess there is no way to make the exc run faster. Reasoning that
unipost
would only extract a 2D field and then do the conversion to GRIB1, I
was
convinced the latter operation would run in a flash, but it appears
that
unipost must read the entire WRF (2GB) file beforehand.
Another question I'd have is related to this field called "mode",
which in
the case of stage IV data is defaulted to 8 (as in the copygb
tutorial),
whilst the regridded data unipost dumps has mode=136 (below). Does it
matter if this field matches or not (i.e., between the re gridded
Stage IV
GRIB1 data and the GRIB1 converted WRF data)? I just changed the 8 in
the
copygb command with a 136 (when running copygb on the Stage IV data)
and
that solved the problem but wanted to know for sure.
wgrib -V WRFPRS06.tm00 | grep -A 6 -in APCP
2503:rec 279:417837088:date 2013053100 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-6hr acc:
2504: APCP=Total precipitation [kg/m^2]
2505- timerange 4 P1 0 P2 6 TimeU 1 nx 1199 ny 799 GDS grid 3
num_in_ave
0 missing 0
2506- center 7 subcenter 0 process 125 Table 2 scan: WE:SN
winds(grid)
2507- Lambert Conf: Lat1 21.641000 Lon1 -120.450000 Lov -98.000000
2508- Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP
0.000000
2509- North Pole (1199 x 799) Dx 4.000000 Dy 4.000000 scan 64
mode 136
2510- min/max data 0 164.26 num bits 15 BDS_Ref 0 DecScale 2
BinScale 0
On Tue, Mar 4, 2014 at 2:58 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Alex,
>
> I concur - configuring wrf_cntrl.parm to only dump out a single
record
> will do what you want.
>
> As to why it's so slow, I have no idea. And I'm really not the
right
> person to answer this question. Support for UPP is provided through
> wrfhelp at ucar.edu. You could ask them if there's any way to
> expedite the postprocessing of a single record.
>
> Assuming you're running WRF-ARW, another option would be running the
> pinterp utility (
> http://www.mmm.ucar.edu/wrf/OnLineTutorial/Tools/p_interp.htm) to
> post-process your output. That may (or may not)
> be faster. MET can read the pinterp output, but I don't know if
there's a
> way to only post-process a single record.
>
> Thanks,
> John
>
> On 03/04/2014 01:53 PM, Alex Fierro via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >
> > Hi John:
> >
> > I just checked and what I initially did was correct: Changing the
> > ../parm/wrf_cntrl.parm and running unipost.exe just does that. The
reason
> > why I am surprised is that the process (unipost.exe on 8 procs)
remains
> > painfully slow - Perhaps this is due to the fact that unipost has
to read
> > the entire netcdf record prior to dumping the data? Is there a way
to
> > accelerate the process ?
> >
> > i) convert wrf.netcdf into a GRIB1 file for only one record (APCP)
- no
> > interp of the WRF data is necessary as I actually interp the stage
IV
> data
> > onto the WRF grid using copygb later on.
> >
> > Cheers and many thanks again for your patience,
> >
> > Alex-
> >
> >
> > On Tue, Mar 4, 2014 at 2:43 PM, Alex Fierro <
> alexandre.o.fierro at gmail.com>wrote:
> >
> >> Hi John:
> >>
> >> Once again, many thanks for your prompt reply and useful tips !
> >>
> >> I did not specify clearly my goal, my mistake.
> >>
> >> I just realized that I might not need to run unipost if the
following is
> >> possible: Is it possible to convert the WRF NETCDF data for one
single
> >> record (APCP=rec#279) to a GRIB1 file ?
> >>
> >> Cheers and many thanks for your help,
> >>
> >> Alex-
> >>
> >>
> >>
> >>
> >> On Tue, Mar 4, 2014 at 2:29 PM, John Halley Gotway via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >>> Alex,
> >>>
> >>> The wrfcntrl.parm file tells the unified post-processor what
fields
> you'd
> >>> like post-processed. You could configure it to only dump out a
single
> GRIB
> >>> record, but typically users want more than one
> >>> output field. The unified post-processor (UPP) and copygb are
not the
> >>> same thing. With wrfcntrl.parm, you're controlling UPP. But
you could
> >>> also run copygb to regrid only a subset of records
> >>> available in the input GRIB file. It's up to you to decide what
makes
> >>> the most sense for your project.
> >>>
> >>> You can configure copygb to run on a subset of the GRIB records
using
> the
> >>> "-k" option. You can read copygb.doc for details
> >>> (
> >>>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt
> )
> >>> and see the examples at the bottom. Example (2) uses "-k'4*-1
61'" to
> >>> select out precipitation whose GRIB number
> >>> is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out
250mb
> height.
> >>>
> >>> But I think the most powerful way of doing it is using wgrib and
copygb
> >>> together as shown in example (8). The "-kw" option tells copygb
to
> read
> >>> the GRIB record numbers from standard input. For example:
> >>> wgrib sample.grb | grep RH | copygb -xg212 -kw
sample_g212.grb
> >>>
> >>> This would filter out the GRIB records for RH, regrid them to
NCEP grid
> >>> 212, and write the output in the file named sample_g212.grb.
> >>>
> >>> wgrib is also pretty useful. You can use it to easily select
one or
> more
> >>> records from an input GRIB file. But if you need to interpolate
to a
> new
> >>> grid, you should use copygb.
> >>>
> >>> Hope that helps.
> >>>
> >>> John
> >>>
> >>> On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
> >>>>
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >>>>
> >>>> Seems that again I found the answer to my own question:
> >>>>
> >>>> I just changed the default wrf_cntrl.parm to only contain:
> >>>>
> >>>> 1 KGTYPE******I5*******:(00255)********START OF THIS
OUTPUT
> >>>> GRID******
> >>>> 2 IMDLTY *I5* :(00125)
> >>>> 3 DATSET *A6* :(WRFPRS)
> >>>> 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
> >>>> 5 L=(10000 00000 00000 00000 00000 00000 00000 00000 00000
00000
> >>> 00000
> >>>> 00000 00000 00000)
> >>>>
> >>>> correct ?
> >>>>
> >>>>
> >>>>
> >>>>
> >>>> On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro <
> >>> alexandre.o.fierro at gmail.com>wrote:
> >>>>
> >>>>> Dear John:
> >>>>>
> >>>>> Many thanks for your prompt reply and for the useful tips.
> >>>>>
> >>>>> While I have your here, I'd have the following query: Is
there a way
> >>> to
> >>>>> run copygb but on *only one single record* (i.e., 1 variable)
within
> >>> the
> >>>>> WRF Netcdf file? I am currently only interested in APCP and
wanted to
> >>> know
> >>>>> if copybg had the ability to only focus its interpolation
magic on
> one
> >>>>> record similar to wgrib e.g. -which would save a lot of cpu
time:
> >>>>>
> >>>>> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
> >>>>>
> >>>>> where 279 is record# 279=APCP in the WRF GRIB1 converted
file.
> >>>>>
> >>>>> Cheers and best,
> >>>>>
> >>>>> Alex
> >>>>>
> >>>>>
> >>>>> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
> >>>>> met_help at ucar.edu> wrote:
> >>>>>
> >>>>>> Alex,
> >>>>>>
> >>>>>> Great, I'm glad you were able to figure it out. We have
several
> >>> examples
> >>>>>> of running copygb for different projections in the METv4.1
online
> >>> tutorial:
> >>>>>>
> >>>>>>
> >>>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
> >>>>>>
> >>>>>> The big limitation is that copygb only works on GRIB1 data.
If you
> >>> have
> >>>>>> GRIB2 data, you can...
> >>>>>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
> >>>>>> (2) Or wgrib2 does have the ability to regrid data but it'd
pretty
> >>>>>> limited.
> >>>>>>
> >>>>>> Hope that helps.
> >>>>>>
> >>>>>> Thanks,
> >>>>>> John Halley Gotway
> >>>>>> met_help at ucar.edu
> >>>>>>
> >>>>>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
> >>>>>>>
> >>>>>>> OK - I just figured it out - using:
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>
>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
> >>>>>>>
> >>>>>>> and:
> >>>>>>>
> >>>>>>>
> http://mailman.ucar.edu/pipermail/met_help/2009-December/001043.html
> >>>>>>>
> >>>>>>> Cheers and many thanks,
> >>>>>>>
> >>>>>>>
> >>>>>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
> >>>>>> met_help at ucar.edu
> >>>>>>>> wrote:
> >>>>>>>
> >>>>>>>> Greetings,
> >>>>>>>>
> >>>>>>>> This message has been automatically generated in response
to the
> >>>>>>>> creation of a trouble ticket regarding:
> >>>>>>>> "Stage IV to model grid",
> >>>>>>>> a summary of which appears below.
> >>>>>>>>
> >>>>>>>> There is no need to reply to this message right now. Your
ticket
> >>> has
> >>>>>> been
> >>>>>>>> assigned an ID of [rt.rap.ucar.edu #65680].
> >>>>>>>>
> >>>>>>>> Please include the string:
> >>>>>>>>
> >>>>>>>> [rt.rap.ucar.edu #65680]
> >>>>>>>>
> >>>>>>>> in the subject line of all future correspondence about this
issue.
> >>> To
> >>>>>> do
> >>>>>>>> so,
> >>>>>>>> you may reply to this message.
> >>>>>>>>
> >>>>>>>> Thank you,
> >>>>>>>> met_help at ucar.edu
> >>>>>>>>
> >>>>>>>>
> >>>>>>
> >>>
>
-------------------------------------------------------------------------
> >>>>>>>> Greetings:
> >>>>>>>>
> >>>>>>>> Is there a way to regrid Stage IV data onto a local WRF
grid using
> >>>>>> copygb
> >>>>>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
> >>>>>>>>
> >>>>>>>> Cheers and best,
> >>>>>>>>
> >>>>>>>> --
> >>>>>>>>
-------------------------------------------------------------
> >>>>>>>> Alexandre Fierro, PhD
> >>>>>>>> Research Scientist-
> >>>>>>>> National Severe Storms Laboratory (NSSL/NOAA)
> >>>>>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
> >>>>>> (OU/NOAA)
> >>>>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>>>>>>>
> >>>>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
> >>> That
> >>>>>> is
> >>>>>>>> why it is called the Present"
> >>>>>>>>
> >>>>>>>> "There are only 10 types of people in the world:
> >>>>>>>> Those who understand binary, and those who don't"
> >>>>>>>>
> >>>>>>>> "My opinions are my own and not representative of OU, NSSL,
> >>>>>>>> AOML, HRD, LANL or any affiliates."
> >>>>>>>> ^.^
> >>>>>>>> (o o)
> >>>>>>>> /( V )\
> >>>>>>>> ---m---m----
> >>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>>
> >>>>>
> >>>>>
> >>>>> --
> >>>>> -------------------------------------------------------------
> >>>>> Alexandre Fierro, PhD
> >>>>> Research Scientist-
> >>>>> National Severe Storms Laboratory (NSSL/NOAA)
> >>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
> >>> (OU/NOAA)
> >>>>>
> >>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>>>>
> >>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
> That
> >>> is
> >>>>> why it is called the Present"
> >>>>>
> >>>>> "There are only 10 types of people in the world:
> >>>>> Those who understand binary, and those who don't"
> >>>>>
> >>>>> "My opinions are my own and not representative of OU, NSSL,
> >>>>> AOML, HRD, LANL or any affiliates."
> >>>>> ^.^
> >>>>> (o o)
> >>>>> /( V )\
> >>>>> ---m---m----
> >>>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>>
> >>
> >>
> >> --
> >> -------------------------------------------------------------
> >> Alexandre Fierro, PhD
> >> Research Scientist-
> >> National Severe Storms Laboratory (NSSL/NOAA)
> >> *The Cooperative Institute for Mesoscale Meteorological Studies*
> (OU/NOAA)
> >>
> >> Los Alamos National Laboratory, Los Alamos, NM (LANL)
> >>
> >> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That
> is
> >> why it is called the Present"
> >>
> >> "There are only 10 types of people in the world:
> >> Those who understand binary, and those who don't"
> >>
> >> "My opinions are my own and not representative of OU, NSSL,
> >> AOML, HRD, LANL or any affiliates."
> >> ^.^
> >> (o o)
> >> /( V )\
> >> ---m---m----
> >>
> >
> >
> >
>
>
--
-------------------------------------------------------------
Alexandre Fierro, PhD
Research Scientist-
National Severe Storms Laboratory (NSSL/NOAA)
*The Cooperative Institute for Mesoscale Meteorological Studies*
(OU/NOAA)
Los Alamos National Laboratory, Los Alamos, NM (LANL)
"Yesterday is History, Tomorrow is a Mystery and Today is a Gift; That
is
why it is called the Present"
"There are only 10 types of people in the world:
Those who understand binary, and those who don't"
"My opinions are my own and not representative of OU, NSSL,
AOML, HRD, LANL or any affiliates."
^.^
(o o)
/( V )\
---m---m----
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65680] AutoReply: Stage IV to model grid
From: John Halley Gotway
Time: Tue Mar 04 15:40:41 2014
Alex,
I don't think you need to worry about it. Running "wgrib -V -GDS10"
to dump the grid description section, I think that the 8 refers to
octet number 17:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table7.html
That specifies various things about the data.
A value of 8 is stored in 8 bits as 00001000, from most to least
significant bit. So bit number 5 is 1 and the rest are 0. Looking at
that table, that means the winds are defined as grid-relative.
A values of 136 is stored in 8 bits as 1001000. So bit numbers 1 and
5 are 1 and the rest are 0. That means the "Direction increments are
not given". I'm not exactly sure what that means, but I
don't think you need to worry about it.
John
On 03/04/2014 03:22 PM, Alex Fierro via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>
> Dear John:
>
> OK - Many thanks once again for your prompt reply. Much appreciated.
I
> guess there is no way to make the exc run faster. Reasoning that
unipost
> would only extract a 2D field and then do the conversion to GRIB1, I
was
> convinced the latter operation would run in a flash, but it appears
that
> unipost must read the entire WRF (2GB) file beforehand.
>
> Another question I'd have is related to this field called "mode",
which in
> the case of stage IV data is defaulted to 8 (as in the copygb
tutorial),
> whilst the regridded data unipost dumps has mode=136 (below). Does
it
> matter if this field matches or not (i.e., between the re gridded
Stage IV
> GRIB1 data and the GRIB1 converted WRF data)? I just changed the 8
in the
> copygb command with a 136 (when running copygb on the Stage IV data)
and
> that solved the problem but wanted to know for sure.
>
> wgrib -V WRFPRS06.tm00 | grep -A 6 -in APCP
> 2503:rec 279:417837088:date 2013053100 APCP kpds5=61 kpds6=1 kpds7=0
> levels=(0,0) grid=255 sfc 0-6hr acc:
> 2504: APCP=Total precipitation [kg/m^2]
> 2505- timerange 4 P1 0 P2 6 TimeU 1 nx 1199 ny 799 GDS grid 3
num_in_ave
> 0 missing 0
> 2506- center 7 subcenter 0 process 125 Table 2 scan: WE:SN
winds(grid)
> 2507- Lambert Conf: Lat1 21.641000 Lon1 -120.450000 Lov -98.000000
> 2508- Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP
0.000000
> 2509- North Pole (1199 x 799) Dx 4.000000 Dy 4.000000 scan 64
mode 136
> 2510- min/max data 0 164.26 num bits 15 BDS_Ref 0 DecScale 2
BinScale 0
>
>
> On Tue, Mar 4, 2014 at 2:58 PM, John Halley Gotway via RT
<met_help at ucar.edu
>> wrote:
>
>> Alex,
>>
>> I concur - configuring wrf_cntrl.parm to only dump out a single
record
>> will do what you want.
>>
>> As to why it's so slow, I have no idea. And I'm really not the
right
>> person to answer this question. Support for UPP is provided
through
>> wrfhelp at ucar.edu. You could ask them if there's any way to
>> expedite the postprocessing of a single record.
>>
>> Assuming you're running WRF-ARW, another option would be running
the
>> pinterp utility (
>> http://www.mmm.ucar.edu/wrf/OnLineTutorial/Tools/p_interp.htm) to
>> post-process your output. That may (or may not)
>> be faster. MET can read the pinterp output, but I don't know if
there's a
>> way to only post-process a single record.
>>
>> Thanks,
>> John
>>
>> On 03/04/2014 01:53 PM, Alex Fierro via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>>>
>>> Hi John:
>>>
>>> I just checked and what I initially did was correct: Changing the
>>> ../parm/wrf_cntrl.parm and running unipost.exe just does that. The
reason
>>> why I am surprised is that the process (unipost.exe on 8 procs)
remains
>>> painfully slow - Perhaps this is due to the fact that unipost has
to read
>>> the entire netcdf record prior to dumping the data? Is there a way
to
>>> accelerate the process ?
>>>
>>> i) convert wrf.netcdf into a GRIB1 file for only one record (APCP)
- no
>>> interp of the WRF data is necessary as I actually interp the stage
IV
>> data
>>> onto the WRF grid using copygb later on.
>>>
>>> Cheers and many thanks again for your patience,
>>>
>>> Alex-
>>>
>>>
>>> On Tue, Mar 4, 2014 at 2:43 PM, Alex Fierro <
>> alexandre.o.fierro at gmail.com>wrote:
>>>
>>>> Hi John:
>>>>
>>>> Once again, many thanks for your prompt reply and useful tips !
>>>>
>>>> I did not specify clearly my goal, my mistake.
>>>>
>>>> I just realized that I might not need to run unipost if the
following is
>>>> possible: Is it possible to convert the WRF NETCDF data for one
single
>>>> record (APCP=rec#279) to a GRIB1 file ?
>>>>
>>>> Cheers and many thanks for your help,
>>>>
>>>> Alex-
>>>>
>>>>
>>>>
>>>>
>>>> On Tue, Mar 4, 2014 at 2:29 PM, John Halley Gotway via RT <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>> Alex,
>>>>>
>>>>> The wrfcntrl.parm file tells the unified post-processor what
fields
>> you'd
>>>>> like post-processed. You could configure it to only dump out a
single
>> GRIB
>>>>> record, but typically users want more than one
>>>>> output field. The unified post-processor (UPP) and copygb are
not the
>>>>> same thing. With wrfcntrl.parm, you're controlling UPP. But
you could
>>>>> also run copygb to regrid only a subset of records
>>>>> available in the input GRIB file. It's up to you to decide what
makes
>>>>> the most sense for your project.
>>>>>
>>>>> You can configure copygb to run on a subset of the GRIB records
using
>> the
>>>>> "-k" option. You can read copygb.doc for details
>>>>> (
>>>>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/copygb.txt
>> )
>>>>> and see the examples at the bottom. Example (2) uses "-k'4*-1
61'" to
>>>>> select out precipitation whose GRIB number
>>>>> is 61. Example (6) uses "-k'4*-1 7 100 250'" to select out
250mb
>> height.
>>>>>
>>>>> But I think the most powerful way of doing it is using wgrib and
copygb
>>>>> together as shown in example (8). The "-kw" option tells copygb
to
>> read
>>>>> the GRIB record numbers from standard input. For example:
>>>>> wgrib sample.grb | grep RH | copygb -xg212 -kw
sample_g212.grb
>>>>>
>>>>> This would filter out the GRIB records for RH, regrid them to
NCEP grid
>>>>> 212, and write the output in the file named sample_g212.grb.
>>>>>
>>>>> wgrib is also pretty useful. You can use it to easily select
one or
>> more
>>>>> records from an input GRIB file. But if you need to interpolate
to a
>> new
>>>>> grid, you should use copygb.
>>>>>
>>>>> Hope that helps.
>>>>>
>>>>> John
>>>>>
>>>>> On 03/04/2014 01:19 PM, Alex Fierro via RT wrote:
>>>>>>
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>>>>>>
>>>>>> Seems that again I found the answer to my own question:
>>>>>>
>>>>>> I just changed the default wrf_cntrl.parm to only contain:
>>>>>>
>>>>>> 1 KGTYPE******I5*******:(00255)********START OF THIS
OUTPUT
>>>>>> GRID******
>>>>>> 2 IMDLTY *I5* :(00125)
>>>>>> 3 DATSET *A6* :(WRFPRS)
>>>>>> 4 (ACM TOTAL PRECIP ) SCAL=( 5.0)
>>>>>> 5 L=(10000 00000 00000 00000 00000 00000 00000 00000
00000 00000
>>>>> 00000
>>>>>> 00000 00000 00000)
>>>>>>
>>>>>> correct ?
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Tue, Mar 4, 2014 at 2:05 PM, Alex Fierro <
>>>>> alexandre.o.fierro at gmail.com>wrote:
>>>>>>
>>>>>>> Dear John:
>>>>>>>
>>>>>>> Many thanks for your prompt reply and for the useful tips.
>>>>>>>
>>>>>>> While I have your here, I'd have the following query: Is
there a way
>>>>> to
>>>>>>> run copygb but on *only one single record* (i.e., 1 variable)
within
>>>>> the
>>>>>>> WRF Netcdf file? I am currently only interested in APCP and
wanted to
>>>>> know
>>>>>>> if copybg had the ability to only focus its interpolation
magic on
>> one
>>>>>>> record similar to wgrib e.g. -which would save a lot of cpu
time:
>>>>>>>
>>>>>>> wgrib -d 279 WRFPRS01.tm00 -o WRFPRS01.binary
>>>>>>>
>>>>>>> where 279 is record# 279=APCP in the WRF GRIB1 converted
file.
>>>>>>>
>>>>>>> Cheers and best,
>>>>>>>
>>>>>>> Alex
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Mar 4, 2014 at 12:27 PM, John Halley Gotway via RT <
>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Alex,
>>>>>>>>
>>>>>>>> Great, I'm glad you were able to figure it out. We have
several
>>>>> examples
>>>>>>>> of running copygb for different projections in the METv4.1
online
>>>>> tutorial:
>>>>>>>>
>>>>>>>>
>>>>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/copygb/index.php
>>>>>>>>
>>>>>>>> The big limitation is that copygb only works on GRIB1 data.
If you
>>>>> have
>>>>>>>> GRIB2 data, you can...
>>>>>>>> (1) Use cnvgrib to convert to GRIB1 and run copygb.
>>>>>>>> (2) Or wgrib2 does have the ability to regrid data but it'd
pretty
>>>>>>>> limited.
>>>>>>>>
>>>>>>>> Hope that helps.
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>> John Halley Gotway
>>>>>>>> met_help at ucar.edu
>>>>>>>>
>>>>>>>> On 03/04/2014 11:23 AM, Alex Fierro via RT wrote:
>>>>>>>>>
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65680 >
>>>>>>>>>
>>>>>>>>> OK - I just figured it out - using:
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>
>>
http://www.mmm.ucar.edu/wrf/users/tutorial/200907/Basic/2009jul_wpp_joint_tutorial.pdf
>>>>>>>>>
>>>>>>>>> and:
>>>>>>>>>
>>>>>>>>>
>> http://mailman.ucar.edu/pipermail/met_help/2009-
December/001043.html
>>>>>>>>>
>>>>>>>>> Cheers and many thanks,
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Tue, Mar 4, 2014 at 11:29 AM, met_help at ucar.edu via RT <
>>>>>>>> met_help at ucar.edu
>>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> Greetings,
>>>>>>>>>>
>>>>>>>>>> This message has been automatically generated in response
to the
>>>>>>>>>> creation of a trouble ticket regarding:
>>>>>>>>>> "Stage IV to model grid",
>>>>>>>>>> a summary of which appears below.
>>>>>>>>>>
>>>>>>>>>> There is no need to reply to this message right now. Your
ticket
>>>>> has
>>>>>>>> been
>>>>>>>>>> assigned an ID of [rt.rap.ucar.edu #65680].
>>>>>>>>>>
>>>>>>>>>> Please include the string:
>>>>>>>>>>
>>>>>>>>>> [rt.rap.ucar.edu #65680]
>>>>>>>>>>
>>>>>>>>>> in the subject line of all future correspondence about this
issue.
>>>>> To
>>>>>>>> do
>>>>>>>>>> so,
>>>>>>>>>> you may reply to this message.
>>>>>>>>>>
>>>>>>>>>> Thank you,
>>>>>>>>>> met_help at ucar.edu
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>
>>>>>
>>
-------------------------------------------------------------------------
>>>>>>>>>> Greetings:
>>>>>>>>>>
>>>>>>>>>> Is there a way to regrid Stage IV data onto a local WRF
grid using
>>>>>>>> copygb
>>>>>>>>>> from UPP ? I also have wgrib and wgrib 2 installed (but not
MET).
>>>>>>>>>>
>>>>>>>>>> Cheers and best,
>>>>>>>>>>
>>>>>>>>>> --
>>>>>>>>>>
-------------------------------------------------------------
>>>>>>>>>> Alexandre Fierro, PhD
>>>>>>>>>> Research Scientist-
>>>>>>>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>>>>>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
>>>>>>>> (OU/NOAA)
>>>>>>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>>>>>>>
>>>>>>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
>>>>> That
>>>>>>>> is
>>>>>>>>>> why it is called the Present"
>>>>>>>>>>
>>>>>>>>>> "There are only 10 types of people in the world:
>>>>>>>>>> Those who understand binary, and those who don't"
>>>>>>>>>>
>>>>>>>>>> "My opinions are my own and not representative of OU, NSSL,
>>>>>>>>>> AOML, HRD, LANL or any affiliates."
>>>>>>>>>> ^.^
>>>>>>>>>> (o o)
>>>>>>>>>> /( V )\
>>>>>>>>>> ---m---m----
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> --
>>>>>>> -------------------------------------------------------------
>>>>>>> Alexandre Fierro, PhD
>>>>>>> Research Scientist-
>>>>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>>>>> *The Cooperative Institute for Mesoscale Meteorological
Studies*
>>>>> (OU/NOAA)
>>>>>>>
>>>>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>>>>
>>>>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a
Gift;
>> That
>>>>> is
>>>>>>> why it is called the Present"
>>>>>>>
>>>>>>> "There are only 10 types of people in the world:
>>>>>>> Those who understand binary, and those who don't"
>>>>>>>
>>>>>>> "My opinions are my own and not representative of OU, NSSL,
>>>>>>> AOML, HRD, LANL or any affiliates."
>>>>>>> ^.^
>>>>>>> (o o)
>>>>>>> /( V )\
>>>>>>> ---m---m----
>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> -------------------------------------------------------------
>>>> Alexandre Fierro, PhD
>>>> Research Scientist-
>>>> National Severe Storms Laboratory (NSSL/NOAA)
>>>> *The Cooperative Institute for Mesoscale Meteorological Studies*
>> (OU/NOAA)
>>>>
>>>> Los Alamos National Laboratory, Los Alamos, NM (LANL)
>>>>
>>>> "Yesterday is History, Tomorrow is a Mystery and Today is a Gift;
That
>> is
>>>> why it is called the Present"
>>>>
>>>> "There are only 10 types of people in the world:
>>>> Those who understand binary, and those who don't"
>>>>
>>>> "My opinions are my own and not representative of OU, NSSL,
>>>> AOML, HRD, LANL or any affiliates."
>>>> ^.^
>>>> (o o)
>>>> /( V )\
>>>> ---m---m----
>>>>
>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
More information about the Met_help
mailing list