[Met_help] [rt.rap.ucar.edu #89736] History for specific vertical levels in ensemble_stat

[John Halley Gotway via RT]

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



Hello there,
    I am going to calculate some ensemble 
fields (Ensemble Mean Field, Ensemble Standard Deviation Field, Ensemble Maximum Field...) with ensemble_stat. It works well with TMP, UGRD, VGRD on pressure surface (such as 149:6993632:d=19031500:TMP:kpds5=11:kpds6=100:kpds7=700:TR=0:P1=6:P2=0:TimeU=1:700 mb:6hr fcst:NAve=0). When dealing with CAPE that is not on pressure surface, I can't properly specify the corresponding level. 


   field = [
      {
         name    = "CAPE";
         level     = ["???"];
         cat_thresh = [ NA ];
      }
   ];
 The records in grib file are :
'''
>>>wgrib fcst_d01.2019-03-15_06_00_00 | grep CAPE
289:13518592:d=19031500:CAPE:kpds5=157:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr fcst:NAve=0
303:14111208:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=46080:TR=0:P1=6:P2=0:TimeU=1:180-0 mb above gnd:6hr fcst:NAve=0
305:14172146:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=23040:TR=0:P1=6:P2=0:TimeU=1:90-0 mb above gnd:6hr fcst:NAve=0
307:14245234:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=65280:TR=0:P1=6:P2=0:TimeU=1:255-0 mb above gnd:6hr fcst:NAve=0
'''
,where the file "fcst_d01.2019-03-15_06_00_00" is output from UPP. The corresponding wrf_cntrl.parm is :
'''
 (CNVCT AVBL POT ENRGY) SCAL=( 3.0)
 L=(11110 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000 00000)
'''




Thanks for any suggestions,
Jiao B. F.
    



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

Subject: specific vertical levels in ensemble_stat
From: John Halley Gotway
Time: Mon Apr 15 12:28:07 2019

Julio,

I see that you're trying to access CAPE.  For surface variables, the
level
information is usually encoded as a height of 0.  So I'd suggest using
"*Z0*".
You can also use "*L0*" (for level value of 0).  Whereas Z0 matches
only
heights above ground, L0 matches any level type with a value of 0.

A couple of suggestions...
- Whenever evaluating new fields in MET, I always run them through
plot_data_plane first.  That way, I can see what the data looks like
and
make sure I have the name/level settings set correctly.
- When using Z0 or L0 for other variable types (in particular HGT),
you may
get multiple matching records.  If that happens, look in
met-8.0/share/met/config/READ for "GRIB1_..." for a list of config
file
options to filter down to the record you actually want to evaluate.

Thanks,
John Halley Gotway


On Sun, Apr 14, 2019 at 9:27 AM 焦宝峰 via RT <met_help at ucar.edu> wrote:

>
> Sun Apr 14 09:27:07 2019: Request 89736 was acted upon.
> Transaction: Ticket created by nuistjiao at 163.com
>        Queue: met_help
>      Subject: specific vertical levels in ensemble_stat
>        Owner: Nobody
>   Requestors: nuistjiao at 163.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=89736 >
>
>
>
>
> Hello there,
>     I am going to calculate some ensemble  fields (Ensemble Mean
Field,
> Ensemble Standard Deviation Field, Ensemble Maximum Field...) with
> ensemble_stat. It works well with TMP, UGRD, VGRD on pressure
surface (such
> as
>
149:6993632:d=19031500:TMP:kpds5=11:kpds6=100:kpds7=700:TR=0:P1=6:P2=0:TimeU=1:700
> mb:6hr fcst:NAve=0). When dealing with CAPE that is not on pressure
> surface, I can't properly specify the corresponding level.
>
>
>    field = [
>       {
>          name    = "CAPE";
>          level     = ["???"];
>          cat_thresh = [ NA ];
>       }
>    ];
>  The records in grib file are :
> '''
> >>>wgrib fcst_d01.2019-03-15_06_00_00 | grep CAPE
>
289:13518592:d=19031500:CAPE:kpds5=157:kpds6=1:kpds7=0:TR=0:P1=6:P2=0:TimeU=1:sfc:6hr
> fcst:NAve=0
>
303:14111208:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=46080:TR=0:P1=6:P2=0:TimeU=1:180-
0
> mb above gnd:6hr fcst:NAve=0
>
305:14172146:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=23040:TR=0:P1=6:P2=0:TimeU=1:90-
0
> mb above gnd:6hr fcst:NAve=0
>
307:14245234:d=19031500:CAPE:kpds5=157:kpds6=116:kpds7=65280:TR=0:P1=6:P2=0:TimeU=1:255-
0
> mb above gnd:6hr fcst:NAve=0
> '''
> ,where the file "fcst_d01.2019-03-15_06_00_00" is output from UPP.
The
> corresponding wrf_cntrl.parm is :
> '''
>  (CNVCT AVBL POT ENRGY) SCAL=( 3.0)
>  L=(11110 00000 00000 00000 00000 00000 00000 00000 00000 00000
00000
> 00000 00000 00000)
> '''
>
>
>
>
> Thanks for any suggestions,
> Jiao B. F.
>
>
>
>

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