Re: Re: [ncl-talk] how to get the domain average of cape through wrf_SkewT2?

Dennis ,
   
   Thank you very much! That's what I need.

Yuhua



----您于 Tue, 21 Oct 2008 17:13:25 -0600 的来信中写到: ----


> 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
> >
> > 中国气象局上海台风研究所
> > -----------------------------------------
>

中国气象局上海台风研究所
-----------------------------------------
Shanghai Typhoon Instutite of China Meteorological Administration

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 21 2008 - 19:32:27 MDT