<div dir="ltr"><div><div><div>Hi,<br><br>I am interested in the following steps:-<br></div>1) I have a 2 D matrix where there are 133 rows and 12 columns called air<br></div>2) I want to c alculate trend between years which go row 1880:2012 and eachcolumn of air. <br>
</div><div>3) Error is unidentified z<br></div><div><br><br>;======================================<br>; create year<br>;======================================<br><br> years=ispan(1880,2012,1)<br> <br>;======================================<br>
; read data; add some meta data<br>;======================================<br><br> diri = "~/Documents/PhD_June_2015/Data_AIR/"<br> fili = "air.txt"<br> air = readAsciiTable(diri+fili, nmos, "float",2)<br>
air@_FillValue = msg<br><br> air@long_name = "All India Precip"<br> air@units = "mm"<br> ; <br> opt = True<br> opt@fout = "all_india_data"<br><br> write_matrix(air, nmos+"f7.2", opt)<br>
;=========================================<br>;Doing trend analysis<br>;=========================================<br><br>do nmo=0,11<br> z(nmo)=regline(years,air(:,nmo))<br>end do<br> <br></div>;======================================<br>
; create year<br>;======================================<br><br> years=ispan(1880,2012,1)<br> <br>;======================================<br>; read data; add some meta data<br>;======================================<br>
<br> diri = "~/Documents/PhD_June_2015/Data_AIR/"<br> fili = "air.txt"<br> air = readAsciiTable(diri+fili, nmos, "float",2)<br> air@_FillValue = msg<br><br> air@long_name = "All India Precip"<br>
air@units = "mm"<br> ; <br> opt = True<br> opt@fout = "all_india_data"<br><br> write_matrix(air, nmos+"f7.2", opt)<br>;=========================================<br>;Doing trend analysis<br>
;=========================================<br><br>do nmo=0,11<br> z(nmo)=regline(years,air(:,nmo))<br>end do<br> <br></div>