[ncl-talk] Problem of drawing Walker circulation

[David Brown]

Or you might try setting vcMapDirection to False. See
http://www.ncl.ucar.edu/Document/Graphics/Resources/vc.shtml#vcMapDirection

 -dave

On Tue, May 30, 2017 at 3:40 PM, Dennis Shea <shea at ucar.edu> wrote:

> 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
>>
>>
>
> _______________________________________________
> 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/72b850f3/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/72b850f3/attachment.png