[ncl-talk] Problem of drawing Walker circulation

[Dennis Shea]

The Walker Circulation:   Generally: U (Zonal wind component)   >> W

Hence the U dominates the plot. You should **scale** the vertical W the
'match' the zonal numbers (or vice versa)

Maybe       scale = abs(Uavg)/abs(Wavg)

   W = scale*W    ; keep meta data
   W at long_name = "W scaled by: "+scale

On Tue, May 30, 2017 at 10:18 AM, Priyanka Silva <
priyankakeraminiya6 at gmail.com> wrote:

> 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]
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170530/7da2f54c/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/7da2f54c/attachment.png