Query about NCL cape_3d function's result?

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
Received on Sat Feb 22 01:31:36 2014