Hi Micha,<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="direction:ltr" bgcolor="#FFFFFF" text="#000000">
    The next addition might be how feed these crosspoints X and Y
    straight into the easting and northing of r.basins...<br>
     <br></div></blockquote></div>you mean something like (untested):<br><br>#!/usr/bin/env python<br><br>fin = open("crosspoints","r")<br>fout = open("commands","w")<br>linesout = (line.rstrip().split() for line in fin)<br>
cmd = 'r.basin map=dem prefix=bas%s easting=%s northing=%s threshold=10000\n'<br>fout.writelines(cmd % (n, easting, northing) for n,(easting,northing) in enumerate(linesout))<br><br>?<br><br>ciao madi<br><br clear="all">
<br>-- <br>Dr. Margherita Di Leo<br>