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

Hello,

I modified your script slightly with some monthly NCEP Reanalysis data
that I have here. I added a few streamline resources, resulting in mixed
results. Here is the graphics portion of my script:

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

res_at_cnLinesOn=False
res_at_cnLineLabelsOn=False
res_at_cnInfoLabelOn=False

  res_at_stArrowLengthF = 0.008 ; default is dynamic
  res_at_stLengthCheckCount = 90 ; default is 35
  res_at_stArrowStride = 1
  res_at_stLineStartStride = 1 ; default is 2
  res_at_stMinArrowSpacingF = 0.0 ; default is 0.0
; res_at_stStepSizeF = 0.001 ; default is dynamic
plot=gsn_csm_pres_hgt_streamline(wks,v(0:11,:,{120}),v(0:11,:,{170}),ws(0:11,:,{120}),res)
end

I would recommend varying/playing with the 6 st resources I used above,
and see if you can create a plot that works for you.
Good luck,
Adam

Áõ²®Ææ wrote:
> 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 <mailto:res_at_cnLinesOn=False>
> res_at_cnLineLabelsOn=False <mailto:res_at_cnLineLabelsOn=False>
> res_at_cnInfoLabelOn=False <mailto:res_at_cnInfoLabelOn=False>
> res_at_stMinArrowSpacingF=0.008 <mailto:res_at_stMinArrowSpacingF=0.008>
> res_at_stArrowLengthF=0.008 <mailto: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

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk