<div dir="ltr"><div>Hi Reza,</div><div>I am not an expert on MJO. However, in this case the problem is that you changed time selection only at the top of the script (which is for calculating climatology) and you did not change time selection for plotting anomalies at these lines:</div><div>
<pre>        ntBegin   = ind(yyyymmdd.eq.<b><span style="color:rgb(255,0,0)">19961001</span></b>)
        ntEnd     = ind(yyyymmdd.eq.<b><span style="color:rgb(255,0,0)">19970401</span></b>)</pre>

</div><div>So for the selected times of 20170101 to 20181231, the above will result in missing values and hence failure of the script. In the attached, I changed the above lines to:</div><div></div><div>           ntBegin   = ind(yyyymmdd.eq.ymdStrt)<br>           ntEnd     = ind(yyyymmdd.eq.ymdLast)</div><div><br></div><div>And that works just fine. </div><div><br></div><div><b>NOTE that in the attached I just chose times from 20150101 to 20181231, it is up to you to decide which periods you want for climatology and which periods you want for plotting anomalies. For anomalies it is always recommend to use several years to get appropriate climatology, 2 years would in my opinion is too short for climatology. <br></b></div><div><b><span style="color:rgb(0,0,255)">Do not forget to change the path for "diri" in the attached.</span><br></b></div><div><b><br></b></div><div><b>hope that helps,</b></div><div><b>Rashed</b><br></div><div><br></div><div><br></div><div><br></div></div>