<br><font size=2 face="sans-serif">Margherita,</font>
<br><font size=2 face="sans-serif"> You
may be running into the open file limit for linux. </font><a href="http://www.cyberciti.biz/faq/linux-increase-the-maximum-number-of-open-files/"><font size=2 face="sans-serif">http://www.cyberciti.biz/faq/linux-increase-the-maximum-number-of-open-files/,
although 1573 files seems a bit low. I generally batch my
r.patch'ing in groups of 120 or so with a python script , then go back
and r.patch the results. It works well if the file names are related
to the file locations. </font></a>
<br>
<br><font size=2 face="sans-serif">See ugly python script below where elevation
grid data with the "temp_"grid tiles were named by the x coordinate
were patched together ( yes, I know you can do this better now with direct
python calls to grass:-)):</font>
<br>
<br><font size=2 face="sans-serif"># grass_tile_merge.py</font>
<br><font size=2 face="sans-serif">import os,string,re,glob</font>
<br><font size=2 face="sans-serif">grdarr=['D2370','D2371','d2371','d2372','D2372','D2373']</font>
<br>
<br><font size=2 face="sans-serif">for grdarrlist in grdarr:</font>
<br><font size=2 face="sans-serif"> gtarr=['0','1','2','3','4','5','6','7','8','9']</font>
<br><font size=2 face="sans-serif"> for
gtarrl in gtarr:</font>
<br><font size=2 face="sans-serif">
gtarrb=['0','1','2','3','4','5','6','7','8','9']</font>
<br><font size=2 face="sans-serif">
for gtarrl2 in gtarrb:</font>
<br><font size=2 face="sans-serif">
gtarrc=['0','1','2','3','4','5','6','7','8','9']</font>
<br><font size=2 face="sans-serif">
for gtarrl3
in gtarrc:</font>
<br><font size=2 face="sans-serif">
grdstr=""</font>
<br><font size=2 face="sans-serif">
gtemp="temp_"+grdarrlist+gtarrl+gtarrl2+gtarrl3</font>
<br><font size=2 face="sans-serif">
#print grpitem</font>
<br><font size=2 face="sans-serif">
greplist="/imagery/grass/ncstpftlidar_83/newdem/fcell/"+grdarrlist+gtarrl+gtarrl2+gtarrl3+"*_ras"</font>
<br><font size=2 face="sans-serif">
grdlist=glob.glob(greplist) # creates a list
of those grid files</font>
<br><font size=2 face="sans-serif">
for t1 in grdlist:</font>
<br><font size=2 face="sans-serif">
t1base=os.path.basename(t1)</font>
<br><font size=2 face="sans-serif">
grdstr=grdstr+t1base+","</font>
<br><font size=2 face="sans-serif">
#grdarr.append(t1base)</font>
<br><font size=2 face="sans-serif">
grdstr=grdstr[:-1]</font>
<br><font size=2 face="sans-serif">
print grdstr</font>
<br><font size=2 face="sans-serif">
regtite="g.region rast=%s"%(grdstr)</font>
<br><font size=2 face="sans-serif">
os.system(regtite)</font>
<br><font size=2 face="sans-serif">
patchgrp="r.patch input=%s output=%s"%(grdstr,gtemp)</font>
<br><font size=2 face="sans-serif">
print patchgrp</font>
<br><font size=2 face="sans-serif">
os.system(patchgrp)</font>
<br><font size=2 face="sans-serif">#http://grass.fbk.eu/grass64/manuals/html64_user/r.patch.html</font>
<br>
<br><font size=2 face="sans-serif">Hope this helps.</font>
<br>
<br><font size=2 face="sans-serif">Doug</font>
<br>
<br>
<br>
<br><font size=2 face="sans-serif">Doug Newcomb
<br>
USFWS<br>
Raleigh, NC<br>
919-856-4520 ext. 14 doug_newcomb@fws.gov<br>
---------------------------------------------------------------------------------------------------------<br>
The opinions I express are my own and are not representative of the official
policy of the U.S.Fish and Wildlife Service or Dept. of the Interior.
Life is too short for undocumented, proprietary data formats.</font>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>Margherita Di Leo <dileomargherita@gmail.com></b>
</font>
<br><font size=1 face="sans-serif">Sent by: grass-dev-bounces@lists.osgeo.org</font>
<p><font size=1 face="sans-serif">07/13/2012 11:27 AM</font>
<td width=59%>
<table width=100%>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">To</font></div>
<td><font size=1 face="sans-serif">grass-dev@lists.osgeo.org</font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">Subject</font></div>
<td><font size=1 face="sans-serif">[GRASS-dev] limits of r.patch</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br><font size=3>Hi All,<br>
<br>
I probably am requiring too much, but this is just for sharing an experience.
I've just imported in GRASS 1573 ASTER GDEM tiles, using r.external, and
everything seemed to be OK. Then I ran:<br>
<br>
GRASS 6.4.3svn (ASTER_GDEM):~ > r.patch in=`g.mlist pat=AST* sep=,`
out=aster_gdem<br>
ERROR 1: TIFFOpen:/forest/ASTER-GDEM/Dem_lzw/ASTGTM_N45E010_dem_lzw.tif:
Too many open files<br>
WARNING: Unable to open raster map <ASTGTM_N45E011_dem_lzw.tif@PERMANENT><br>
ERROR: One or more input raster maps not found<br>
<br>
I'm on Red Hat Enterprise Linux Server Rel. 6.3, RAM 16 GB.<br>
<br>
Cheers,<br>
<br>
-- </font><font size=3 color=#5f5f5f><br>
Dr. Margherita Di Leo</font><font size=3><br>
<br>
<br>
</font><tt><font size=2>_______________________________________________<br>
grass-dev mailing list<br>
grass-dev@lists.osgeo.org<br>
</font></tt><a href="http://lists.osgeo.org/mailman/listinfo/grass-dev"><tt><font size=2>http://lists.osgeo.org/mailman/listinfo/grass-dev</font></tt></a>
<br>