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-talkReceived on Wed Apr 15 2009 - 16:24:14 MDT
This archive was generated by hypermail 2.2.0 : Thu Apr 16 2009 - 10:05:49 MDT