vertical velocity(omega)

Dear NCL users,
I want to calculate the vertical velocity with the help of NCEP reanalysis data. So in this I have used the programme is given below but it is not running, so please any one could help me in this regard.
Thanks in advance
kalim
;******************
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
 
  begin
   a     = addfile ("/data/ncep_data/pressure/omega.2005.nc" , "r")
   omega1 = short2flt(w->omega(152,:,:,:))
   ;lev  = a ->levelist
   ;lat  = a ->latitude
   ;lon  = a ->longitude
   omega = dim_avg_Wrap( omega1(levelist|:,longitude|:,{latitude|15:20}))
 
   wks = gsn_open_wks("X11" ,"omega")
   gsn_define_colormap(wks,"BlWhRe")
 
   res=True
  res@gsnDraw         = False
  res@gsnFrame        = False
  res@cnInfoLabelOn = False
  res@gsnContourZeroLineThicknessF = 2.
  res@gsnContourNegLineDashPattern = 1
  plot  = gsn_csm_pres_hgt(wks,omega(:,{70:100}),res)
  draw(plot)
  frame (wks)
  end
 
 

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Nov 2 19:26:11 2009