Re: Query about NCL cape_3d function's result?

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon Feb 24 2014 - 09:09:46 MST

Hi,

Note: if you set "it" to -1 instead of a time index, the CAPE will be done across all time steps on the file.
I don't believe it is giving you an average.

I'm forwarding this to wrfhelp@ucar.edu just in case.

--Mary

On Feb 22, 2014, at 1:31 AM, mmkamal@uwaterloo.ca wrote:

> Hi,
>
> I have used the NCL function to calculate CAPE from WRF output. But, I
> am not sure whether the following code will give me only 1st time step
> or temporal average. I wan to plot monthly CAPE. Any help would be
> greatly appreciated!
>
>
> Thanks
> Kamal
>
>
> =====================================================
>
> a = addfile("wrfout_d01_2002-07-01.nc","r")
>
> times = wrf_user_getvar(a,"times",-1)
> ntimes = dimsizes(times)
>
> it = ntimes-1
> res@TimeLabel = times(it)
>
> cape3d = wrf_user_getvar(a,"cape_3d",it)
>
> ; extract the data from the multi-dimensional fields
>
> cape = cape3d(0,:,:,:)
> cin = cape3d(1,:,:,:)
>
> contour = wrf_contour(a,wks,mcape,opts)
> plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
>
> ==================================================================
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 24 09:10:00 2014

This archive was generated by hypermail 2.1.8 : Mon Mar 03 2014 - 14:26:18 MST