[ncl-talk] Problem of drawing Walker circulation

[Priyanka Silva]

Hello Ncl Talk,

I want to plot Walker circulation. I tried to plot that but I could not get
correct plot. Could you give me some advice on me solve this problem. My
code and plot are as follows.I use NCL version 6.3.0

Thank you

Priyanka




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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"

begin
;----------------------------------------------------------------------
ECMWF = " /home/priya/ecmwf/monthly/"
files = systemfunc("ls "+ECMWF+"uwnd/u.mon.*.nc")
fi = addfiles(files, "r")
T = fi[:]->time
YMD = cd_calendar(T,2)
iTim1 = ind (YMD.ge.19790101.and.YMD.lt.20130101)
u = fi[:]->u (iTim1,{1000:200},{-5:10},{40:130})  ;Zonal wind
printVarSummary(u)
;------------------------------------------------------------------------
ECMWF = " /home/priya/ecmwf/monthly/"
fs = systemfunc("ls "+ECMWF+"wwnd/w.mon.*.nc")
fin1 = addfiles(fs, "r")
time = fin1[:]->time
ymd = cd_calendar(time,2)
iTim2 = ind (ymd.ge.19790101.and.ymd.lt.20130101)
w = fin1[:]->w (iTim2,{1000:200},{-5:10},{40:130}) ;Vertical wind
printVarSummary(w)
;-----------------------------------------------------------------------
U = month_to_season(u,"JJA")
W = month_to_season(w,"JJA")
;----------------------------------------------------------------------
;printVarSummary(U)
;printVarSummary(W)
;---------------------------------------------------------------------
Uavg = dim_avg_Wrap(dim_avg_Wrap(U(lev|:,lon|:,lat|:,time|:)))
Wavg = dim_avg_Wrap(dim_avg_Wrap(W(lev|:,lon|:,lat|:,time|:)))

Uv = dim_avg_Wrap(U(lev|:,lon|:,lat|:,time|:))
printVarSummary(Uv)

printVarSummary(Uavg)
printVarSummary(Wavg)

wks=gsn_open_wks("x11","")
res=True
res at cnLinesOn=False
res at cnLineLabelsOn=False
res at cnInfoLabelOn=False
res at stMinArrowSpacingF=0.008
res at stArrowLengthF=0.008
res at gsnLeftString         = ""
res at gsnRightString         = ""

plot=gsn_csm_pres_hgt_streamline(wks,Uavg,Uavg,Wavg,res)

end

[image: Inline image 1]
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170530/fe7ea840/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 29371 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170530/fe7ea840/attachment.png