<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns:mv="http://macVmlSchemaUri" xmlns="http://www.w3.org/TR/REC-html40"><head><meta name=Title content=""><meta name=Keywords content=""><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:Arial;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
span.msoIns
{mso-style-type:export-only;
mso-style-name:"";
text-decoration:underline;
color:teal;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:595.0pt 842.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style></head><body bgcolor=white lang=EN-GB link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='mso-fareast-language:EN-US'>I think you can do some more debugging yourself. <o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'>That error message comes from some internal code. You have to isolate what part of your script is causing this problem.<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'>Some good programming practice would help a lot here. Comment out parts of your code beginning from the end, while keeping the script functional, until you find out where the problem comes from.<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'>You also mention that your script works for one hemisphere but not the other. Look at your input data to see what has changed or what is it you assumed when you wrote this script that does not apply when you change hemisphere.<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'>This is not something for NCL users to do for you.<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'>/M<o:p></o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'><o:p> </o:p></span></p><div><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>Marston S. Ward, PhD<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>Department of Earth Sciences<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>University of Gothenburg, Sweden<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>Email: <a href="mailto:marston.johnston@gu.se"><span style='color:#0070C0'>marston.johnston@gu.se</span></a><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>SkypeID: marston.johnston <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>Phone: +46-31-7864901 <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>Only the fruitful thing is true!<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070C0'>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<o:p></o:p></span></p></div><p class=MsoNormal><span style='mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span style='mso-fareast-language:EN-US'><o:p> </o:p></span></p><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:12.0pt;color:black'>From: </span></b><span style='font-size:12.0pt;color:black'>ncl-talk <ncl-talk-bounces@ucar.edu> on behalf of Dinuka Silva <dinu.kera@gmail.com><br><b>Date: </b>Friday, 22 September 2017 at 10:51<br><b>To: </b><ncl-talk@ucar.edu><br><b>Subject: </b>[ncl-talk] plot Walker circulation<o:p></o:p></span></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><div><p class=MsoNormal>Hello All,<o:p></o:p></p></div><div><p class=MsoNormal>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. <o:p></o:p></p></div><div><p class=MsoNormal>the error is<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><div><p class=MsoNormal>fatal:["NclFile.c":2099]:Subscript out of range, error in subscript #3<o:p></o:p></p></div><div><p class=MsoNormal>fatal:Memory allocation failure:[errno=12]<o:p></o:p></p></div></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>I would appreciate any help. </span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>Thank you.</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>Dinuka.</span><o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>My ncl code is<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><div><p class=MsoNormal>begin<o:p></o:p></p></div><div><p class=MsoNormal>;----------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>ECMWF = " /home/dinu/ecmwf/monthly/"<o:p></o:p></p></div><div><p class=MsoNormal>files = systemfunc("ls "+ECMWF+"uwnd/u.mon.*.nc")<o:p></o:p></p></div><div><p class=MsoNormal>fi = addfiles(files, "r")<o:p></o:p></p></div><div><p class=MsoNormal>T = fi[:]->time<o:p></o:p></p></div><div><p class=MsoNormal>YMD = cd_calendar(T,2)<o:p></o:p></p></div><div><p class=MsoNormal>iTim1 = ind (YMD.ge.19790101.and.YMD.lt.20130101)<o:p></o:p></p></div><div><p class=MsoNormal>u = fi[:]->u (iTim1,{1000:200},{5:-10},{-40:-110}) ;Zonal wind<o:p></o:p></p></div><div><p class=MsoNormal>printVarSummary(u)<o:p></o:p></p></div><div><p class=MsoNormal>;------------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>ECMWF = " /home/dinu/ecmwf/monthly/"<o:p></o:p></p></div><div><p class=MsoNormal>fs = systemfunc("ls "+ECMWF+"wwnd/w.mon.*.nc")<o:p></o:p></p></div><div><p class=MsoNormal>fin1 = addfiles(fs, "r")<o:p></o:p></p></div><div><p class=MsoNormal>time = fin1[:]->time<o:p></o:p></p></div><div><p class=MsoNormal>ymd = cd_calendar(time,2)<o:p></o:p></p></div><div><p class=MsoNormal>iTim2 = ind (ymd.ge.19790101.and.ymd.lt.20130101)<o:p></o:p></p></div><div><p class=MsoNormal>w = fin1[:]->w (iTim2,{1000:200},{5:-10},{-40:-110}) ;Vertical wind<o:p></o:p></p></div><div><p class=MsoNormal>printVarSummary(w)<o:p></o:p></p></div><div><p class=MsoNormal>;-----------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>U = month_to_season(u,"JJA")<o:p></o:p></p></div><div><p class=MsoNormal>W = month_to_season(w,"JJA")<o:p></o:p></p></div><div><p class=MsoNormal>;------------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>scale=abs(U)/abs(W)<o:p></o:p></p></div><div><p class=MsoNormal>;printVarSummary(scale)<o:p></o:p></p></div><div><p class=MsoNormal>W = scale*W<o:p></o:p></p></div><div><p class=MsoNormal>;printVarSummary(W)<o:p></o:p></p></div><div><p class=MsoNormal>;----------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>Uavg = dim_avg_Wrap(dim_avg_Wrap(U(lev|:,lon|:,lat|:,time|:)))<o:p></o:p></p></div><div><p class=MsoNormal>Wavg = dim_avg_Wrap(dim_avg_Wrap(W(lev|:,lon|:,lat|:,time|:)))<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>;printVarSummary(Uavg)<o:p></o:p></p></div><div><p class=MsoNormal>;printVarSummary(Wavg)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>wks=gsn_open_wks("x11","")<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>res=True<o:p></o:p></p></div><div><p class=MsoNormal>res@cnLinesOn=False<o:p></o:p></p></div><div><p class=MsoNormal>res@cnLineLabelsOn=False<o:p></o:p></p></div><div><p class=MsoNormal>res@cnInfoLabelOn=False<o:p></o:p></p></div><div><p class=MsoNormal>res@gsnLeftString = ""<o:p></o:p></p></div><div><p class=MsoNormal>res@gsnRightString = ""<o:p></o:p></p></div><div><p class=MsoNormal>res@vcRefMagnitudeF = 3.0 ; define vector ref mag<o:p></o:p></p></div><div><p class=MsoNormal>res@vcRefLengthF = 0.045 ; define length of vec ref<o:p></o:p></p></div><div><p class=MsoNormal>res@vcGlyphStyle = "CurlyVector" ; turn on curly vectors<o:p></o:p></p></div><div><p class=MsoNormal>res@vcMinDistanceF = 0.01 ; thin out vectors<o:p></o:p></p></div><div><p class=MsoNormal>res@vcMapDirection = False<o:p></o:p></p></div><div><p class=MsoNormal>res@tiMainString = "Walker circulation"<o:p></o:p></p></div><div><p class=MsoNormal>;-------------------------------------------------------------------<o:p></o:p></p></div><div><p class=MsoNormal>plot=gsn_csm_pres_hgt_vector(wks,Uavg,Uavg,Wavg,res)<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>end<o:p></o:p></p></div></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div><div id=DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2><p class=MsoNormal><o:p> </o:p></p><table class=MsoNormalTable border=1 cellpadding=0 style='border:none;border-top:solid #D3D4DE 1.0pt'><tr><td width=55 style='width:41.25pt;border:none;padding:9.75pt .75pt .75pt .75pt'><p class=MsoNormal><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"><span style='text-decoration:none'><img border=0 width=46 height=29 id="_x0000_i1025" src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif"></span></a><o:p></o:p></p></td><td width=470 style='width:352.5pt;border:none;padding:9.0pt .75pt .75pt .75pt'><p class=MsoNormal style='line-height:13.5pt'><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#41424E'>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"><span style='color:#4453EA'>www.avast.com</span></a> <o:p></o:p></span></p></td></tr></table></div><p class=MsoNormal>_______________________________________________ ncl-talk mailing list ncl-talk@ucar.edu List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talk <o:p></o:p></p></div></body></html>