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

From: <lbqnj_at_nyahnyahspammersnyahnyah>
Date: Tue, 14 Apr 2009 20:47:11 +0800

Hello, Adam
    Thanks for your kindness! But I'm afraid what you said in the last E-mail can not solve my problem. Now I've sent the figures to you, the file 'right.jpg' is made by Grads, while the one 'wrong.jpg' is by NCL. The script is the same as what you've sent to me. Could you please show me what' wrong? Thanks!
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

right.jpg wrong.jpg
Received on Tue Apr 14 2009 - 06:47:11 MDT

This archive was generated by hypermail 2.2.0 : Thu Apr 16 2009 - 10:05:49 MDT