<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>I am looking for suggestions on a query of worldclim data (<a href="http://worldclim.org/">http://worldclim.org/</a>).</div><div><br></div><div>I have used something similar to below for previous data extraction where I am looking for data associated with a particular site, or list of sites, where a single lon lat coordinate will return a single value. However, in this case, I need to pull an entire longitudinal band for a particular latitude. </div><div><br></div><div>For example instead of grabbing data from a single coordinates with "east_north = -92.50, 46.51”, I want to grab all the data in the longitudinal band at latitude 46.51. </div><div><br></div><div>If it were possible something like “east_north = *, 46.51” </div><div><br></div><div>Any suggestions on this point would be much appreciated. Thanks in advance - Kirk</div><div><br></div><div><br></div>==============================================<div><br></div><div><br><div apple-content-edited="true">
<div><div>if test "$GISBASE" = ""; then</div><div> echo "You must be in GRASS to run this program"</div><div> exit</div><div>fi</div><div><br></div><div># Use input text file for coords in the format: lon lat (easting northing) single </div><div># space between coords. Example: -98.42072 55.91481. Must use real coordinates, no blanks. </div><div># Script will create file for each $MAP in the list. g.mlist can also take the "pattern"</div><div># argument which takes regular expressions. For example: pattern="*" returns all </div><div># maps in the database. If script fails check for hidden characters in text file, best to </div><div># use vi or emacs to eliminate them</div><div><br></div><div>for MAP in `g.mlist type=rast pattern=“HIST_tmin*"`; do</div><div> r.what --verbose input="$MAP" < ~/Desktop/ai/coords.txt > ~/Desktop/ai/"$MAP".txt;</div><div><br></div><div>echo "$MAP"</div><div><br></div><div>done</div></div><div><br></div><br class="Apple-interchange-newline">
</div>
<br></div></body></html>