<div dir="ltr"><div><div><div>The Walker Circulation: Generally: U (Zonal wind component) >> W <br><br>Hence the U dominates the plot. You should **scale** the vertical W the 'match' the zonal numbers (or vice versa)<br><br></div>Maybe scale = abs(Uavg)/abs(Wavg)<br></div><div><br></div> W = scale*W ; keep meta data<br></div> W@long_name = "W scaled by: "+scale<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 30, 2017 at 10:18 AM, Priyanka Silva <span dir="ltr"><<a href="mailto:priyankakeraminiya6@gmail.com" target="_blank">priyankakeraminiya6@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Ncl Talk,<div><br></div><div>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</div><div><br></div><div>Thank you</div><div><br></div><div>Priyanka</div><div><br></div><div><br></div><div><br></div><div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/diagnostics_<wbr>cam.ncl"</div><div><br></div><div>begin</div><div>;-----------------------------<wbr>------------------------------<wbr>-----------</div><div>ECMWF = " /home/priya/ecmwf/monthly/"</div><div>files = systemfunc("ls "+ECMWF+"uwnd/u.mon.*.nc")</div><div>fi = addfiles(files, "r")</div><div>T = fi[:]->time</div><div>YMD = cd_calendar(T,2)</div><div>iTim1 = ind (<a href="http://YMD.ge.19790101.and.YMD.lt">YMD.ge.19790101.and.YMD.lt</a>.<wbr>20130101)</div><div>u = fi[:]->u (iTim1,{1000:200},{-5:10},{40:<wbr>130}) ;Zonal wind</div><div>printVarSummary(u)</div><div>;-----------------------------<wbr>------------------------------<wbr>-------------</div><div>ECMWF = " /home/priya/ecmwf/monthly/"</div><div>fs = systemfunc("ls "+ECMWF+"wwnd/w.mon.*.nc")</div><div>fin1 = addfiles(fs, "r")</div><div>time = fin1[:]->time</div><div>ymd = cd_calendar(time,2)</div><div>iTim2 = ind (<a href="http://ymd.ge.19790101.and.ymd.lt">ymd.ge.19790101.and.ymd.lt</a>.<wbr>20130101)</div><div>w = fin1[:]->w (iTim2,{1000:200},{-5:10},{40:<wbr>130}) ;Vertical wind</div><div>printVarSummary(w)</div><div>;-----------------------------<wbr>------------------------------<wbr>------------</div><div>U = month_to_season(u,"JJA")</div><div>W = month_to_season(w,"JJA")</div><div>;-----------------------------<wbr>------------------------------<wbr>-----------</div><div>;printVarSummary(U)</div><div>;printVarSummary(W)</div><div>;-----------------------------<wbr>------------------------------<wbr>----------</div><div>Uavg = dim_avg_Wrap(dim_avg_Wrap(U(<wbr>lev|:,lon|:,lat|:,time|:)))</div><div>Wavg = dim_avg_Wrap(dim_avg_Wrap(W(<wbr>lev|:,lon|:,lat|:,time|:)))</div><div><br></div><div>Uv = dim_avg_Wrap(U(lev|:,lon|:,<wbr>lat|:,time|:))</div><div>printVarSummary(Uv)</div><div><br></div><div>printVarSummary(Uavg)</div><div>printVarSummary(Wavg)</div><div><br></div><div>wks=gsn_open_wks("x11","")</div><div>res=True</div><div>res@cnLinesOn=False</div><div>res@cnLineLabelsOn=False</div><div>res@cnInfoLabelOn=False</div><div>res@stMinArrowSpacingF=0.008</div><div>res@stArrowLengthF=0.008</div><div>res@gsnLeftString = ""</div><div>res@gsnRightString = ""</div><div><br></div><div>plot=gsn_csm_pres_hgt_<wbr>streamline(wks,Uavg,Uavg,Wavg,<wbr>res)</div><div><br></div><div>end</div></div><div><br></div><div><img src="cid:ii_15c59b5f39de164a" alt="Inline image 1" width="489" height="394"><br></div><div><br></div><p class="MsoNormal"><span></span></p></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>