----- Original Message -----
From: B.Q. Liu
To: asphilli_at_cgd.ucar.edu
Cc: ncl-talk_at_ucar.edu
Sent: Tuesday, April 14, 2009 8:47 PM
Subject: Re: [ncl-talk] A problem about drawing a Hadley Cell by the NCEP/NCAR Reanalysis Data
Hello,
Thanks for Adam's help! But I'm afraid what he said in the last E-mail can not solve my problem. Now I've sent the figures to you, the right figure ('right.jpg') is made by GrADs, while the wrong one ('wrong.jpg') is made by NCL. Could you please show me what' wrong?
The script is as below (I have chosen all latitudes to draw a globle Hadley circulation along 120E this time.)
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")
v=short2flt(vf->vwnd)
w=short2flt(wf->omega)
wAve=avg(w(49,:,:,{120})) ; you don't need to scale the data
vAve=avg(v(49,:,:,{120})) ; if you don't want to.
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,:,{120}),v(49,0:11,{-40:40},{120}),ws(49,0:11,:,{120}),res)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.2.0 : Thu Apr 16 2009 - 10:05:49 MDT