<div dir="ltr"><div class="gmail_default" style="font-size:small">See example shapefile_13.ncl at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml">http://www.ncl.ucar.edu/Applications/shapefiles.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">You will want to download the USA shapefiles from <a href="http://gadm.org/country">gadm.org/country</a>, and use the USA_adm1.shp shapefile, which contains the outlines for each state.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 1, 2016 at 1:49 PM, Adv <span dir="ltr"><<a href="mailto:advita6@gmail.com" target="_blank">advita6@gmail.com</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>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 = "<a href="http://air.mon.mean.nc" target="_blank">air.mon.mean.nc</a>"<br> f = addfile(diri+fili,"r") <br><br> YYYYMM = cd_calendar( f->time, -1)<br><br> iStrt = ind(YYYYMM.eq.ymStrt)<br> iLast = ind(YYYYMM.eq.ymLast)<br> xt = f->air(:,::-1,:)<br> printVarSummary(xt) <wbr> ; [time| 720]x[lat| 91]x[lon| 180]<br></div>; Minnesota., Montana......<br><div><div> x = xt(iStrt:iLast,{40:49},{104:<wbr>97})<br> x1 = xt(iStrt:iLast,{45:49},{116:<wbr>104})<br> x2 = xt(iStrt:iLast,{43.5:49},{97:<wbr>93})<br> x3 = xt(iStrt:iLast,{43.5:49},{97.<wbr>5:90})<br> delete(xt)<br> printVarSummary(x) <wbr> ; [time| 720]x[lat| 91]x[lon| 180]<br>;*****************************<wbr>*******************<br>; Calculate the Annual Mean (slopes) <br>;*****************************<wbr>*******************<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>
<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>