Re: Query about NCL cape_3d function's result?

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Tue Feb 25 2014 - 21:54:46 MST

It should be in the call to wrf_user_getvar. When you use “-1” for the third argument, it calculates that quantity for all time steps:

>>> cape3d = wrf_user_getvar(a,"cape_3d”,-1)

—Mary

On Feb 24, 2014, at 7:17 PM, mmkamal@uwaterloo.ca wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 25 21:54:59 2014

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