Hamish: ... > use 'g.region rast=tile1,tile2,tile3,tile5' first to set the computational > region to the outside bounds of all the tiles listed.
OR, given that you keep the "original" naming convention of the (elevation) tiles without the ".tif" extension (e.g. ASTGTM_N34E024_dem, ASTGTM_N35E023_dem, etc.), you can do: g.region rast=`g.mlist pat=AST*dem sep=,` -pa --- @Hamish: Something against using g.region's "-a" flag? Nikos _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user