Re: Query about NCL cape_3d function's result?

From: <mmkamal_at_nyahnyahspammersnyahnyah>
Date: Mon Feb 24 2014 - 19:17:33 MST

Hi Mary,

I am a bit confused about where should I "-1" instead of "it". Could
you please be more specific.

Thanks
Kamal

Quoting Mary Haley <haley@ucar.edu>:

> 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 19:17:52 2014

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