<div dir="ltr"><div><div><div><div>according to this answer  <a href="https://www.ncl.ucar.edu/Support/talk_archives/2013/2968.html" target="_blank">https://www.ncl.ucar.edu/Suppo<wbr>rt/talk_archives/2013/2968.<wbr>html</a><br><br></div>I need the 10 positive and negative lag regression  of my data x(time), P1(lat,lon,time) but I dont loop trough lat lon .<br><br></div><div>This is my code :<br><br></div>dim_s   = dimsizes(P1)  ;;; <br> Nend = 5000 ; length of  x<br> nlagr= 21 ; +- 10 lag <br>  yLag = 10<br><br> Regdata              = new((/dim_s(0),dim_s(1),nlagr/),"float")<br><br>n=-1<br>do n=0,nlagr-1<br>      n= n+1<br>       Regdata(:,:,n) =regline(x(0:Nend-n-yLag), P1(:,:,yLag+n:Nend))<br>end do ;<br><br><br></div>This give me an error :<br>fatal:Subscript out of range, error in subscript #2<br>fatal:An error occurred reading P1<br><br></div><div>I try to replace regline by regCoef and get the same error message.<br></div><div><br></div>Thanks in advance!!<br><div><div><br></div></div></div>