Re: how to get the domain average of cape through wrf_SkewT2?

Hello,

I am not sure I understand the question.

You could manually edit the skewT code:
    $NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl

In the function "function skewT_PlotData", add the following
at the very end of the script

after:
   end if ; end "PlotWindH"

add the following 6 lines
   skewt_bkgd_at_cape = (/ cape /)
   skewt_bkgd_at_precip_water = pwat ; cm
   skewt_bkgd_at_showalter_index = shox
   skewt_bkgd_at_plcl = plcl
   skewt_bkgd_at_tlcl = tlcl

before
   return (skewt_bkgd)

=============================================

Then, you could do the following:
[snip]
   skewt_data = skewT_PlotData (wks, skewt_bkgd, p,tc,tdc,z \
                                     , wspd,wdir, dataOpts)
   draw (skewt_bkgd)
   draw (skewt_data)
   frame(wks)

printVarSummary(skewt_bkgd)
printVarSummary(skewt_data)

    CAPE = skwt_bkgd_at_cape ; joules

Yuhua wrote:
> Hi,
>
> I want to get the domain average of cape through wrf_SkewT2.ncl
> which called the skewt_func.ncl, and I can output the cape of each grid
> in the skewt_func.ncl, but how to accumulate all the grids? I want to
> return the cape value from skewt_func.ncl to wrf_SkewT2.ncl , but I failed.
>
> Thanks very much!
>
> Yuhua
>
> 中国气象局上海台风研究所
> -----------------------------------------
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 21 2008 - 17:13:25 MDT