<div dir="ltr"><div><div><div><div><div><div>Hi Saharwardi,<br><br></div>The error is on line 22 and if I'm not mistaken your line 22 might be this one spi(nr,:) = dim_spi_n(prc, runlen(nr), opt, 0)<br><br></div>Your SPI must be more than 2 dimensions.<br></div><br></div>See example 2 from here <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml</a><br><br></div>hope it helps<br></div>Izidine<br><div><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Izidine Pinto<br><br><br><br><br></div></div></div></div>
<br><div class="gmail_quote">On 3 May 2017 at 07:58, Md. Saquib Saharwardi 1610303 <span dir="ltr"><<a href="mailto:saquib16@iiserb.ac.in" target="_blank">saquib16@iiserb.ac.in</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"><div><div><div>Dear NCL-user,<br></div> I am trying to plot SPI (standardized precipitation Index) for 30-years data. It is always showing the problem while execution of scripts.<br></div><br>please provide me suggestions on that.<br><br></div>I haver attached the script with error below.<br><br>fili = "<a href="http://precip.1981.2010.nc" target="_blank">precip.1981.2010.nc</a>"<br> f = addfile(fili, "r")<br><br> yyyymm = f->time<br><br> ntim = dimsizes(yyyymm)<br> if ((ntim%12).ne.0) then<br> print("ntim must be divisable by 12")<br> exit<br> end if<br> print(ntim)<br><br> prc = f->precip<br> prc@_FillValue = -999.0 ; change from the default<br> pmsg = prc@_FillValue ; convenience<br><br> printVarSummary(prc)<br> printMinMax(prc,0)<br><br> runlen = (/ 3, 6, 12, 24, 36, 48 /)<br> nrun = dimsizes(runlen)<br><br> opt = False<br> ;opt@spi_type = 3 ; 6.3.0 onward; only if opt=True<br><br> spi = new( (/nrun,ntim/), typeof(prc), pmsg)<br> do nr=0,nrun-1<br> spi(nr,:) = dim_spi_n(prc, runlen(nr), opt, 0)<br> end do<br><br>fatal:Dimension sizes on right hand side of assignment do not match dimension sizes of left hand side<br>fatal:["Execute.c":8573]:<wbr>Execute: Error occurred at or near line 22<br><br clear="all"><div><div><div><div><div><div class="m_-4285432753145418797gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div><br><br><br></div><div>Thankyou in Advance<br><br></div><div>With regards<br></div><div><span style="font-size:12.8px">******************************<wbr>******************************<wbr>*****************</span><br></div>Md. Saquib Saharwardi<div>Research Scholar</div><div>Department of Earth and Environmental Sciences</div><div>Indian Institute of Science Education and Research Bhopal (IISERB)</div><div>Bhopal- 462066, Madhya Pradesh</div><div>Email: <a href="mailto:saquib16@iiserb.ac.in" target="_blank">saquib16@iiserb.ac.in</a></div><div>Mob: <a href="tel:+91%2074711%2012404" value="+917471112404" target="_blank">+91-7471112404</a><br></div><div><a href="http://home.iiserb.ac.in/%7Ekumarp/index.html" style="font-size:12.8px" target="_blank">CGMLI (Home Page)</a><br></div></div></div></div></div></div></div></div></div>
</div></div></div></div></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>