[Met_help] [rt.rap.ucar.edu #51653] History for mask_grid and mask_poly

Paul Oldenburg via RT met_help at ucar.edu
Mon Dec 5 15:42:02 MST 2011


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

Hi,

We tried to use MET to validate the simulation results of wrf against the observation from ncep gdas prepburf files.
We have used CONUS.poly for mask_poly (since our simulation is in CONUS) and choose "mask_grid[] = [ ];". Is 
this correct?

Thanks,


Shujia



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

Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: Paul Oldenburg
Time: Wed Nov 23 11:28:20 2011

Shujia,

You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
observations over conus:

mask_grid[] = [ ];

mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];

Please let us know if you have any other questions.

Thanks,

Paul


On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>
> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>        Queue: met_help
>      Subject: mask_grid and mask_poly
>        Owner: Nobody
>   Requestors: shujia.zhou-1 at nasa.gov
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>
>
> Hi,
>
> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
> this correct?
>
> Thanks,
>
>
> Shujia
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: zhining.tao at nasa.gov
Time: Wed Nov 23 12:28:13 2011

Hello Paul,

Thank you for the quick reply. We do have several other questions. In
order to compare our WRF results with the observations using the MET
program, we:

1) use "p_interp" to put the WRF output into the pressure level, the
namelist.pinterp looks:

--------------------------
path_to_output = '.'
process        = 'all'
fields         = 'U, V, T2, SMOIS, QVAPOR'
/
met_em_output  = .TRUE.
split_output   = .TRUE.
debug          = .TRUE.
bit64          = .TRUE.

! process can be "all", "list", "all_list"
!       "all" - process all fields in file
!               diagnostics PRES, TT, GHT, RH are also calculated
!       "list" - list fields required in "fields"
!               available diagnostics are PRES, TT, GHT, RH

&interp_in
interp_levels  = 1000.,987.5,975.,962.5,950.,937.5,925.,
                 900.,875.,850.,825.,800.,750.,700.,650.,
                 600.,550.,500.,450.,400.,350.,300.,250.,
                 225.,200.,175.,150.,137.5,125.,112.5,100.,
                 87.5,75.,62.5,50.,37.5,25.,12.5,
/
extrapolate    = 1,
interp_method  = 2,
unstagger_grid = .TRUE.
---------------------------------

Here we use the default 38 pressure levels for interpolation. Is it
properly written?


2) use "pb2nc" to convert the GRIB observation file into a netcdf
file. In the configure file, what should we put in the message_type[
]? We want to compare the model results at surface, 850 mb, 500 mb,
and 300 mb to the corresponding observations. Will the "ADPSFC"
suffice or both "ADPSFC" and "ADPUPA" be included? In other words, do
the "ADPSFC" data include upper level measurements or just the surface
observation?


3) use "point_stat" to do the point comparison. If we want to compare
wind at 850 mb, should we do:
   fcst_field[ ] = ["UU(0,9,*,*)", "VV(0,9,*,*)"]
   obs_field[ ] = ["WIND/P850", "WDIR/P850"]

Note that 850 mb is the 10th element in the interp_levels array. Does
the array index start at "0" or at "1"? Does the "0" in the UU array
represents the time step 1?
In addition, does "TMP/Z2" represent the temperature at 2 m?

The threshold set in fcst_thresh seems not to work properly. E.g., we
set the threshold of TMP/Z2 as "gt300", but the mean tmp is around 295
K from the point_stat output. Is there anything wrong?

Again, for message_tpe[ ], should we put both "ADPSFC" and "ADPUPA" in
order to compares data at both surface and upper air (e.g., 850 mb)?

We highly appreciate your quick response.

Zhining Tao


On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:

> Shujia,
>
> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
> observations over conus:
>
> mask_grid[] = [ ];
>
> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>
> Please let us know if you have any other questions.
>
> Thanks,
>
> Paul
>
>
> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>
>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>       Queue: met_help
>>     Subject: mask_grid and mask_poly
>>       Owner: Nobody
>>  Requestors: shujia.zhou-1 at nasa.gov
>>      Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>
>>
>> Hi,
>>
>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>> this correct?
>>
>> Thanks,
>>
>>
>> Shujia
>>
>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: Paul Oldenburg
Time: Wed Nov 23 12:32:15 2011

Zhining,

I see that you're not currently a registered user of the MET software.
Prior to providing direct support in the use of
MET, we require that users fill out the free online registration form.
That enables us to ensure that we remain in
compliance with the international export laws of the United States.

To fill out the registration, just go here:
   http://www.dtcenter.org/met/users/downloads/index.php
Enter you email address and hit enter to be taken to the form.

Please write back to let me know when you've registered.

Thanks,

Paul



On 11/23/2011 12:28 PM, zhining.tao at nasa.gov via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>
> Hello Paul,
>
> Thank you for the quick reply. We do have several other questions.
In order to compare our WRF results with the observations using the
MET program, we:
>
> 1) use "p_interp" to put the WRF output into the pressure level, the
namelist.pinterp looks:
>
> --------------------------
> path_to_output = '.'
> process        = 'all'
> fields         = 'U, V, T2, SMOIS, QVAPOR'
> /
> met_em_output  = .TRUE.
> split_output   = .TRUE.
> debug          = .TRUE.
> bit64          = .TRUE.
>
> ! process can be "all", "list", "all_list"
> !       "all" - process all fields in file
> !               diagnostics PRES, TT, GHT, RH are also calculated
> !       "list" - list fields required in "fields"
> !               available diagnostics are PRES, TT, GHT, RH
>
> &interp_in
> interp_levels  = 1000.,987.5,975.,962.5,950.,937.5,925.,
>                  900.,875.,850.,825.,800.,750.,700.,650.,
>                  600.,550.,500.,450.,400.,350.,300.,250.,
>                  225.,200.,175.,150.,137.5,125.,112.5,100.,
>                  87.5,75.,62.5,50.,37.5,25.,12.5,
> /
> extrapolate    = 1,
> interp_method  = 2,
> unstagger_grid = .TRUE.
> ---------------------------------
>
> Here we use the default 38 pressure levels for interpolation. Is it
properly written?
>
>
> 2) use "pb2nc" to convert the GRIB observation file into a netcdf
file. In the configure file, what should we put in the message_type[
]? We want to compare the model results at surface, 850 mb, 500 mb,
and 300 mb to the corresponding observations. Will the "ADPSFC"
suffice or both "ADPSFC" and "ADPUPA" be included? In other words, do
the "ADPSFC" data include upper level measurements or just the surface
observation?
>
>
> 3) use "point_stat" to do the point comparison. If we want to
compare wind at 850 mb, should we do:
>    fcst_field[ ] = ["UU(0,9,*,*)", "VV(0,9,*,*)"]
>    obs_field[ ] = ["WIND/P850", "WDIR/P850"]
>
> Note that 850 mb is the 10th element in the interp_levels array.
Does the array index start at "0" or at "1"? Does the "0" in the UU
array represents the time step 1?
> In addition, does "TMP/Z2" represent the temperature at 2 m?
>
> The threshold set in fcst_thresh seems not to work properly. E.g.,
we set the threshold of TMP/Z2 as "gt300", but the mean tmp is around
295 K from the point_stat output. Is there anything wrong?
>
> Again, for message_tpe[ ], should we put both "ADPSFC" and "ADPUPA"
in order to compares data at both surface and upper air (e.g., 850
mb)?
>
> We highly appreciate your quick response.
>
> Zhining Tao
>
>
> On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:
>
>> Shujia,
>>
>> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
>> observations over conus:
>>
>> mask_grid[] = [ ];
>>
>> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>>
>> Please let us know if you have any other questions.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>>
>>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>>       Queue: met_help
>>>     Subject: mask_grid and mask_poly
>>>       Owner: Nobody
>>>  Requestors: shujia.zhou-1 at nasa.gov
>>>      Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>>
>>>
>>> Hi,
>>>
>>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>>> this correct?
>>>
>>> Thanks,
>>>
>>>
>>> Shujia
>>>
>>
>>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: zhining.tao at nasa.gov
Time: Wed Nov 23 12:38:46 2011

Paul, now I registered.

Thanks,

Zhining


On Nov 23, 2011, at 2:32 PM, Paul Oldenburg via RT wrote:

> Zhining,
>
> I see that you're not currently a registered user of the MET
software.  Prior to providing direct support in the use of
> MET, we require that users fill out the free online registration
form.  That enables us to ensure that we remain in
> compliance with the international export laws of the United States.
>
> To fill out the registration, just go here:
>   http://www.dtcenter.org/met/users/downloads/index.php
> Enter you email address and hit enter to be taken to the form.
>
> Please write back to let me know when you've registered.
>
> Thanks,
>
> Paul
>
>
>
> On 11/23/2011 12:28 PM, zhining.tao at nasa.gov via RT wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>
>> Hello Paul,
>>
>> Thank you for the quick reply. We do have several other questions.
In order to compare our WRF results with the observations using the
MET program, we:
>>
>> 1) use "p_interp" to put the WRF output into the pressure level,
the namelist.pinterp looks:
>>
>> --------------------------
>> path_to_output = '.'
>> process        = 'all'
>> fields         = 'U, V, T2, SMOIS, QVAPOR'
>> /
>> met_em_output  = .TRUE.
>> split_output   = .TRUE.
>> debug          = .TRUE.
>> bit64          = .TRUE.
>>
>> ! process can be "all", "list", "all_list"
>> !       "all" - process all fields in file
>> !               diagnostics PRES, TT, GHT, RH are also calculated
>> !       "list" - list fields required in "fields"
>> !               available diagnostics are PRES, TT, GHT, RH
>>
>> &interp_in
>> interp_levels  = 1000.,987.5,975.,962.5,950.,937.5,925.,
>>                 900.,875.,850.,825.,800.,750.,700.,650.,
>>                 600.,550.,500.,450.,400.,350.,300.,250.,
>>                 225.,200.,175.,150.,137.5,125.,112.5,100.,
>>                 87.5,75.,62.5,50.,37.5,25.,12.5,
>> /
>> extrapolate    = 1,
>> interp_method  = 2,
>> unstagger_grid = .TRUE.
>> ---------------------------------
>>
>> Here we use the default 38 pressure levels for interpolation. Is it
properly written?
>>
>>
>> 2) use "pb2nc" to convert the GRIB observation file into a netcdf
file. In the configure file, what should we put in the message_type[
]? We want to compare the model results at surface, 850 mb, 500 mb,
and 300 mb to the corresponding observations. Will the "ADPSFC"
suffice or both "ADPSFC" and "ADPUPA" be included? In other words, do
the "ADPSFC" data include upper level measurements or just the surface
observation?
>>
>>
>> 3) use "point_stat" to do the point comparison. If we want to
compare wind at 850 mb, should we do:
>>   fcst_field[ ] = ["UU(0,9,*,*)", "VV(0,9,*,*)"]
>>   obs_field[ ] = ["WIND/P850", "WDIR/P850"]
>>
>> Note that 850 mb is the 10th element in the interp_levels array.
Does the array index start at "0" or at "1"? Does the "0" in the UU
array represents the time step 1?
>> In addition, does "TMP/Z2" represent the temperature at 2 m?
>>
>> The threshold set in fcst_thresh seems not to work properly. E.g.,
we set the threshold of TMP/Z2 as "gt300", but the mean tmp is around
295 K from the point_stat output. Is there anything wrong?
>>
>> Again, for message_tpe[ ], should we put both "ADPSFC" and "ADPUPA"
in order to compares data at both surface and upper air (e.g., 850
mb)?
>>
>> We highly appreciate your quick response.
>>
>> Zhining Tao
>>
>>
>> On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:
>>
>>> Shujia,
>>>
>>> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
>>> observations over conus:
>>>
>>> mask_grid[] = [ ];
>>>
>>> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>>>
>>> Please let us know if you have any other questions.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>>>
>>>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>>>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>>>      Queue: met_help
>>>>    Subject: mask_grid and mask_poly
>>>>      Owner: Nobody
>>>> Requestors: shujia.zhou-1 at nasa.gov
>>>>     Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>>>
>>>>
>>>> Hi,
>>>>
>>>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>>>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>>>> this correct?
>>>>
>>>> Thanks,
>>>>
>>>>
>>>> Shujia
>>>>
>>>
>>>
>>
>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: Paul Oldenburg
Time: Wed Nov 23 12:45:32 2011

Zhining,

1.  Regarding p_interp, I do not know very much about this software.
I think the best approach would be to test your
model output using point_stat and then make adjustments to
namelist.pinterp if necessary.  If you need help with
p_interp, please contact WRF help at wrfhelp at ucar.edu.

2.  If you are interested in using observations both from the surface
and upper air, you must use both message types
"ADPSFC" and "ADPUPA".  Try:

message_type[] = [ "ADPUPA", "ADPSFC" ];

or, you could also use:

message_type[] = [ "ANYSFC" ];

See the PointsStatConfig message_type comments for more details.

3.
> Does the array index start at "0" or at "1"?
0

> Does the "0" in the UU array represents the time step 1?
It depends on how p_interp output is set up, but probably.

> In addition, does "TMP/Z2" represent the temperature at 2 m?
Yes

The fcst_thresh setting is used to build contingency tables for the
data, not to perform conditional verification.
Conditional verification is not supported in MET at this time.

Please let me know if you have any other questions.

Paul


On 11/23/2011 12:28 PM, zhining.tao at nasa.gov via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>
> Hello Paul,
>
> Thank you for the quick reply. We do have several other questions.
In order to compare our WRF results with the observations using the
MET program, we:
>
> 1) use "p_interp" to put the WRF output into the pressure level, the
namelist.pinterp looks:
>
> --------------------------
> path_to_output = '.'
> process        = 'all'
> fields         = 'U, V, T2, SMOIS, QVAPOR'
> /
> met_em_output  = .TRUE.
> split_output   = .TRUE.
> debug          = .TRUE.
> bit64          = .TRUE.
>
> ! process can be "all", "list", "all_list"
> !       "all" - process all fields in file
> !               diagnostics PRES, TT, GHT, RH are also calculated
> !       "list" - list fields required in "fields"
> !               available diagnostics are PRES, TT, GHT, RH
>
> &interp_in
> interp_levels  = 1000.,987.5,975.,962.5,950.,937.5,925.,
>                  900.,875.,850.,825.,800.,750.,700.,650.,
>                  600.,550.,500.,450.,400.,350.,300.,250.,
>                  225.,200.,175.,150.,137.5,125.,112.5,100.,
>                  87.5,75.,62.5,50.,37.5,25.,12.5,
> /
> extrapolate    = 1,
> interp_method  = 2,
> unstagger_grid = .TRUE.
> ---------------------------------
>
> Here we use the default 38 pressure levels for interpolation. Is it
properly written?
>
>
> 2) use "pb2nc" to convert the GRIB observation file into a netcdf
file. In the configure file, what should we put in the message_type[
]? We want to compare the model results at surface, 850 mb, 500 mb,
and 300 mb to the corresponding observations. Will the "ADPSFC"
suffice or both "ADPSFC" and "ADPUPA" be included? In other words, do
the "ADPSFC" data include upper level measurements or just the surface
observation?
>
>
> 3) use "point_stat" to do the point comparison. If we want to
compare wind at 850 mb, should we do:
>    fcst_field[ ] = ["UU(0,9,*,*)", "VV(0,9,*,*)"]
>    obs_field[ ] = ["WIND/P850", "WDIR/P850"]
>
> Note that 850 mb is the 10th element in the interp_levels array.
Does the array index start at "0" or at "1"? Does the "0" in the UU
array represents the time step 1?
> In addition, does "TMP/Z2" represent the temperature at 2 m?
>
> The threshold set in fcst_thresh seems not to work properly. E.g.,
we set the threshold of TMP/Z2 as "gt300", but the mean tmp is around
295 K from the point_stat output. Is there anything wrong?
>
> Again, for message_tpe[ ], should we put both "ADPSFC" and "ADPUPA"
in order to compares data at both surface and upper air (e.g., 850
mb)?
>
> We highly appreciate your quick response.
>
> Zhining Tao
>
>
> On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:
>
>> Shujia,
>>
>> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
>> observations over conus:
>>
>> mask_grid[] = [ ];
>>
>> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>>
>> Please let us know if you have any other questions.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>>
>>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>>       Queue: met_help
>>>     Subject: mask_grid and mask_poly
>>>       Owner: Nobody
>>>  Requestors: shujia.zhou-1 at nasa.gov
>>>      Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>>
>>>
>>> Hi,
>>>
>>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>>> this correct?
>>>
>>> Thanks,
>>>
>>>
>>> Shujia
>>>
>>
>>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: zhining.tao at nasa.gov
Time: Wed Nov 23 13:31:52 2011

Thanks a lot, Paul.

Happy Thanksgiving.

Zhining

On Nov 23, 2011, at 2:45 PM, Paul Oldenburg via RT wrote:

> Zhining,
>
> 1.  Regarding p_interp, I do not know very much about this software.
I think the best approach would be to test your
> model output using point_stat and then make adjustments to
namelist.pinterp if necessary.  If you need help with
> p_interp, please contact WRF help at wrfhelp at ucar.edu.
>
> 2.  If you are interested in using observations both from the
surface and upper air, you must use both message types
> "ADPSFC" and "ADPUPA".  Try:
>
> message_type[] = [ "ADPUPA", "ADPSFC" ];
>
> or, you could also use:
>
> message_type[] = [ "ANYSFC" ];
>
> See the PointsStatConfig message_type comments for more details.
>
> 3.
>> Does the array index start at "0" or at "1"?
> 0
>
>> Does the "0" in the UU array represents the time step 1?
> It depends on how p_interp output is set up, but probably.
>
>> In addition, does "TMP/Z2" represent the temperature at 2 m?
> Yes
>
> The fcst_thresh setting is used to build contingency tables for the
data, not to perform conditional verification.
> Conditional verification is not supported in MET at this time.
>
> Please let me know if you have any other questions.
>
> Paul
>
>
> On 11/23/2011 12:28 PM, zhining.tao at nasa.gov via RT wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>
>> Hello Paul,
>>
>> Thank you for the quick reply. We do have several other questions.
In order to compare our WRF results with the observations using the
MET program, we:
>>
>> 1) use "p_interp" to put the WRF output into the pressure level,
the namelist.pinterp looks:
>>
>> --------------------------
>> path_to_output = '.'
>> process        = 'all'
>> fields         = 'U, V, T2, SMOIS, QVAPOR'
>> /
>> met_em_output  = .TRUE.
>> split_output   = .TRUE.
>> debug          = .TRUE.
>> bit64          = .TRUE.
>>
>> ! process can be "all", "list", "all_list"
>> !       "all" - process all fields in file
>> !               diagnostics PRES, TT, GHT, RH are also calculated
>> !       "list" - list fields required in "fields"
>> !               available diagnostics are PRES, TT, GHT, RH
>>
>> &interp_in
>> interp_levels  = 1000.,987.5,975.,962.5,950.,937.5,925.,
>>                 900.,875.,850.,825.,800.,750.,700.,650.,
>>                 600.,550.,500.,450.,400.,350.,300.,250.,
>>                 225.,200.,175.,150.,137.5,125.,112.5,100.,
>>                 87.5,75.,62.5,50.,37.5,25.,12.5,
>> /
>> extrapolate    = 1,
>> interp_method  = 2,
>> unstagger_grid = .TRUE.
>> ---------------------------------
>>
>> Here we use the default 38 pressure levels for interpolation. Is it
properly written?
>>
>>
>> 2) use "pb2nc" to convert the GRIB observation file into a netcdf
file. In the configure file, what should we put in the message_type[
]? We want to compare the model results at surface, 850 mb, 500 mb,
and 300 mb to the corresponding observations. Will the "ADPSFC"
suffice or both "ADPSFC" and "ADPUPA" be included? In other words, do
the "ADPSFC" data include upper level measurements or just the surface
observation?
>>
>>
>> 3) use "point_stat" to do the point comparison. If we want to
compare wind at 850 mb, should we do:
>>   fcst_field[ ] = ["UU(0,9,*,*)", "VV(0,9,*,*)"]
>>   obs_field[ ] = ["WIND/P850", "WDIR/P850"]
>>
>> Note that 850 mb is the 10th element in the interp_levels array.
Does the array index start at "0" or at "1"? Does the "0" in the UU
array represents the time step 1?
>> In addition, does "TMP/Z2" represent the temperature at 2 m?
>>
>> The threshold set in fcst_thresh seems not to work properly. E.g.,
we set the threshold of TMP/Z2 as "gt300", but the mean tmp is around
295 K from the point_stat output. Is there anything wrong?
>>
>> Again, for message_tpe[ ], should we put both "ADPSFC" and "ADPUPA"
in order to compares data at both surface and upper air (e.g., 850
mb)?
>>
>> We highly appreciate your quick response.
>>
>> Zhining Tao
>>
>>
>> On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:
>>
>>> Shujia,
>>>
>>> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
>>> observations over conus:
>>>
>>> mask_grid[] = [ ];
>>>
>>> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>>>
>>> Please let us know if you have any other questions.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>>>
>>>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>>>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>>>      Queue: met_help
>>>>    Subject: mask_grid and mask_poly
>>>>      Owner: Nobody
>>>> Requestors: shujia.zhou-1 at nasa.gov
>>>>     Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>>>
>>>>
>>>> Hi,
>>>>
>>>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>>>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>>>> this correct?
>>>>
>>>> Thanks,
>>>>
>>>>
>>>> Shujia
>>>>
>>>
>>>
>>
>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51653] mask_grid and mask_poly
From: Shujia Zhou
Time: Wed Nov 23 14:15:56 2011

Paul:


Thank you for your prompt reply.


Shujia
On Nov 23, 2011, at 1:28 PM, Paul Oldenburg via RT wrote:

> Shujia,
>
> You are correct.  You can use the following two settings in your
point_stat config file to mask out only PrepBUFR
> observations over conus:
>
> mask_grid[] = [ ];
>
> mask_poly[] = [ "MET_BASE/data/poly/CONUS.poly" ];
>
> Please let us know if you have any other questions.
>
> Thanks,
>
> Paul
>
>
> On 11/23/2011 10:37 AM, Shujia Zhou via RT wrote:
>>
>> Wed Nov 23 10:37:41 2011: Request 51653 was acted upon.
>> Transaction: Ticket created by shujia.zhou-1 at nasa.gov
>>       Queue: met_help
>>     Subject: mask_grid and mask_poly
>>       Owner: Nobody
>>  Requestors: shujia.zhou-1 at nasa.gov
>>      Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51653 >
>>
>>
>> Hi,
>>
>> We tried to use MET to validate the simulation results of wrf
against the observation from ncep gdas prepburf files.
>> We have used CONUS.poly for mask_poly (since our simulation is in
CONUS) and choose "mask_grid[] = [ ];". Is
>> this correct?
>>
>> Thanks,
>>
>>
>> Shujia
>>
>
>


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


More information about the Met_help mailing list