<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Feb 23, 2011, at 2:37 AM, Marcelino II Villafuerte wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hi all,<div><br></div><div>Aside from reading and writing a shapefile in grads (<a href="http://www.iges.org/grads/gadoc/gadocindex.html">http://www.iges.org/grads/gadoc/gadocindex.html</a>)</div><div>is it also possible to produce an areal average [aave() function] of a specific parameter e.g. rainfall, </div> <div>in a specific province or region prescribed in a shapefile?</div></blockquote>No. The contents of shapefiles are treated in GrADS as graphical elements (points, lines, or polygons) to be drawn in a geo-registered context. </div><div><br></div><div>To do what you want, you would have to draw the polygon that corresponds to the province or region of interest, then overlay a 'gxout grid' plot of your data, and then use the 'q pos' interface to manually flag grid points that you want to include in your area average, creating a customized mask. Take a look at 'set defval'. This is not a straightforward task, and would require a bit of script writing. </div><div>--Jennifer</div><div><br></div><div><br></div><div><br><blockquote type="cite"><div><br></div><div><br clear="all">Thanks in advance!<br>-- <br>---------------------------------------<br>Marcelino Q. Villafuerte II<br>Weather Specialist I<br> IAAS, CAD<br>PAGASA-DOST<br><br> </div> _______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>http://gradsusr.org/mailman/listinfo/gradsusr<br></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div>--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></div></span> </div><br></body></html>