<div>Dear Folks,</div>
<div>&nbsp;</div>
<div>I am trying to extimate linear trend in an interannual series of one-monthly data from a multy-year monthly data file.</div>
<div>I can estimate properly the linear trend using the script below, if the&nbsp;data file has the one-monthly data series previously separately. The problem is when using the&nbsp;full monthly multy-year data file.</div>
<div>&nbsp;</div>
<div>How can I just select one given month from the data , for instance June, and estimate the linear trend for the June-series alone, having a multy-year data file with all the months of the year?<br clear="all"></div>
<div>I have tried using the function monmask (with argument 6), but I cannot solve it properly.</div>
<div>&nbsp;</div>
<div>Any suggestion I will much appreciate.</div>
<div>&nbsp;</div>
<div>Eduardo.<br></div>
<div>&nbsp;</div>
<div>reinit&#39;<br>&#39;sdfopen d:\anterior\era\<a href="http://mes.nc">mes.nc</a>&#39;</div>
<div>&#39;set clopts -1 -1 0.16&#39;<br>&#39;set mproj sps&#39;</div>
<div>k=70</div>
<div><br>&#39;set grads off&#39;<br>&#39;set lev &#39;k&#39;&#39;<br>&#39;set mproj sps&#39;<br>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;</div>
<div><br>&#39;set lat -60&#39;<br>&#39;set lon&nbsp; -180&#39;<br>&#39;set t 1 132&#39;</div>
<div>* contruyo variable lineal 1, 2, ..., 276 to para calcular m y c de y = mx + c<br>&#39;xvar = 100*t*0 + 1&#39;<br>&#39;xvar = const(xvar,1,-u)&#39;<br>&#39;xvar = 133 - sum(xvar,t+0,t=132)&#39;</div>
<div>&nbsp;</div>
<div>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;<br>&#39;set t 1&#39;</div>
<div><br>&#39;define coeff = tregr(xvar,100*t,t=1, t=132)&#39;<br>&#39;define tave = ave(100*t, t=1, t=132)&#39;<br>&#39;define xvarave = ave(xvar, t=1, t=132)&#39;<br>&#39;define trend = coeff*(xvar - xvarave) + tave&#39;</div>