<div dir="ltr"><div>Hello All,<br></div><div>I'm trying to plot Walker circulation. I could plot it for the eastern hemisphere but when I use the same code for western hemisphere I could not plot. </div><div>the error is</div><div><br></div><div><div>fatal:["NclFile.c":2099]:Subscript out of range, error in subscript #3</div><div>fatal:Memory allocation failure:[errno=12]</div></div><div><br></div><div><span style="font-size:12.8px">I would appreciate any help. </span><br></div><div><span style="font-size:12.8px">Thank you.</span></div><div><span style="font-size:12.8px">Dinuka.</span></div><div><br></div><div>My ncl code is</div><div><br></div><div><br></div><div><div>begin</div><div>;----------------------------------------------------------------------</div><div>ECMWF = " /home/dinu/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 (YMD.ge.19790101.and.YMD.lt.20130101)</div><div>u = fi[:]->u (iTim1,{1000:200},{5:-10},{-40:-110})           ;Zonal wind</div><div>printVarSummary(u)</div><div>;------------------------------------------------------------------------</div><div>ECMWF = " /home/dinu/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 (ymd.ge.19790101.and.ymd.lt.20130101)</div><div>w = fin1[:]->w (iTim2,{1000:200},{5:-10},{-40:-110})           ;Vertical wind</div><div>printVarSummary(w)</div><div>;-----------------------------------------------------------------------</div><div>U = month_to_season(u,"JJA")</div><div>W = month_to_season(w,"JJA")</div><div>;------------------------------------------------------------------------</div><div>scale=abs(U)/abs(W)</div><div>;printVarSummary(scale)</div><div>W = scale*W</div><div>;printVarSummary(W)</div><div>;----------------------------------------------------------------------</div><div>Uavg = dim_avg_Wrap(dim_avg_Wrap(U(lev|:,lon|:,lat|:,time|:)))</div><div>Wavg = dim_avg_Wrap(dim_avg_Wrap(W(lev|:,lon|:,lat|:,time|:)))</div><div><br></div><div>;printVarSummary(Uavg)</div><div>;printVarSummary(Wavg)</div><div><br></div><div>wks=gsn_open_wks("x11","")</div><div><br></div><div>res=True</div><div>res@cnLinesOn=False</div><div>res@cnLineLabelsOn=False</div><div>res@cnInfoLabelOn=False</div><div>res@gsnLeftString         = ""</div><div>res@gsnRightString         = ""</div><div>res@vcRefMagnitudeF = 3.0                ; define vector ref mag</div><div>res@vcRefLengthF    = 0.045              ; define length of vec ref</div><div>res@vcGlyphStyle    = "CurlyVector"      ; turn on curly vectors</div><div>res@vcMinDistanceF  = 0.01               ; thin out vectors</div><div>res@vcMapDirection  = False</div><div>res@tiMainString         = "Walker circulation"</div><div>;-------------------------------------------------------------------</div><div>plot=gsn_csm_pres_hgt_vector(wks,Uavg,Uavg,Wavg,res)</div><div><br></div><div>end</div></div><div><br></div><div><br></div></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br>
<table style="border-top:1px solid #d3d4de">
        <tr>
        <td style="width:55px;padding-top:13px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:12px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Virus-free. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link" target="_blank" style="color:#4453ea">www.avast.com</a>
                </td>
        </tr>
</table><a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div>