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
This archive was generated by hypermail 2.1.8 : Fri Nov 06 2009 - 09:15:29 MST