<div dir="ltr"><div><div>Hi ,<br></div>Could someone help me to compute area average of specific US state. ?<br><br>Thank you<br>Adv<br><br>  ymStrt   = 196501     <br>  ymLast   = 200512<br><br>   fili   = &quot;<a href="http://air.mon.mean.nc">air.mon.mean.nc</a>&quot;<br>   f      = addfile(diri+fili,&quot;r&quot;) <br><br>   YYYYMM = cd_calendar( f-&gt;time, -1)<br><br>   iStrt  = ind(YYYYMM.eq.ymStrt)<br>   iLast  = ind(YYYYMM.eq.ymLast)<br>   xt      = f-&gt;air(:,::-1,:)<br>   printVarSummary(xt)                            ; [time| 720]x[lat| 91]x[lon| 180]<br></div>; Minnesota., Montana......<br><div><div>   x      = xt(iStrt:iLast,{40:49},{104:97})<br>   x1      = xt(iStrt:iLast,{45:49},{116:104})<br>   x2      = xt(iStrt:iLast,{43.5:49},{97:93})<br>   x3      = xt(iStrt:iLast,{43.5:49},{97.5:90})<br>   delete(xt)<br>   printVarSummary(x)                            ; [time| 720]x[lat| 91]x[lon| 180]<br>;************************************************<br>; Calculate the Annual Mean (slopes)    <br>;************************************************<br>rc=  dim_avg_n(x,(/1,2/))<br>rc1=  dim_avg_n(x1,(/1,2/))<br>rc2=  dim_avg_n(x2,(/1,2/))<br>rc3=  dim_avg_n(x3,(/1,2/))<br><br>   printVarSummary(rc)<br>t=(rc+rc1+rc2+rc3)/4<br>print(t)<br>return<br><br></div></div></div>