A problem about drawing a Hadley Cell by the NCEP/NCAR Reanalysis Data

Hello,
 I'm a starter of NCL, I have a problem about drawing a Hadley Cell along 120E by the NCEP/NCAR Reanalysis Data. The problem is that the graphic seems something wrong when I use "gsn_csm_pres_hgt_streamline". Although I want to use the function "vinth2p", I can not find the Var "hyam" and "hybm" (such as the script: http://www.ncl.ucar.edu/Applications/Scripts/stream_8.ncl) in the NCEP/NCAR Reanalysis Data. Could you please tell me how to solve the problem? Thanks for your help!
The program is as below:
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"
begin
dir="/media/work/"
vf=addfile(dir+"vwnd.2009.nc","r")
wf=addfile(dir+"omega.2009.nc","r")
vv=vf->vwnd
v=short2flt(vf->vwnd)
copy_VarMeta(vv,v)
ww=wf->omega
w=short2flt(wf->omega)
copy_VarMeta(ww,w)

wAve=avg(w(49,:,:,{120}))
vAve=avg(v(49,:,:,{120}))
scale=fabs(vAve/wAve)
ws=w*scale
copy_VarCoords(w,ws)

wks=gsn_open_wks("x11","hadley")
res=True

res_at_cnLinesOn=False
res_at_cnLineLabelsOn=False
res_at_cnInfoLabelOn=False
res_at_stMinArrowSpacingF=0.008
res_at_stArrowLengthF=0.008
plot=gsn_csm_pres_hgt_streamline(wks,v(49,0:11,{-40:40},{120}),v(49,0:11,{-40:40},{120}),-ws(49,0:11,{-40:40},{120}),res)
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Apr 13 2009 - 08:15:53 MDT