<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
<font size="-1">hi all, <br>
<br>
I had problems merging grids with the nek tools.<br>
When the final merged grid was having more than 100000 elements
genmap was failing. <br>
<br>
I think the problem was due to an inconsistency between the
writing formats for the BC section in n2to3,nekmerge and genmap
(at least in my version, see below). In my case i have solved the
problem by changing the formats in </font><font size="-1">nekmerge.f.</font><font
size="-1"> <br>
<br>
just my 2 cents<br>
francesco<br>
<br>
<br>
<br>
------------------------------<br>
n2to3.f lines 441-... <br>
------------------------------<br>
<br>
if (neln.lt.1000) then<br>
write(11,20) cbc(k,e),id,k,(bc(j,k,e),j=1,5)<br>
elseif (neln.lt.100000) then<br>
write(11,21) cbc(k,e),id,k,(bc(j,k,e),j=1,5)<br>
else<br>
write(11,22) cbc(k,e),id,k,(bc(j,k,e),j=1,5)<br>
endif<br>
....<br>
20 FORMAT(1x,A3,2I3,5G14.7)<br>
21 FORMAT(1x,A3,i5,i1,5G14.7)<br>
22 FORMAT(1x,A3,i6,i1,5G14.7)<br>
<br>
<br>
<br>
</font><font size="-1">----------------------------------<br>
nekmerge.f lines 664-... <br>
----------------------------------</font><font size="-1"><br>
<br>
if (nel.lt.1000) then<br>
write(11,20) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5)<br>
elseif (nel.lt.100000) then<br>
write(11,21) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5)<br>
else<br>
write(11,22) cbc(f,e,fld),e,(bc(j,f,e,fld),j=1,5)<br>
endif<br>
...<br>
20 format(1x,a3,2i3,5g14.6)<br>
21 format(1x,a3,i5,i1,5g14.6)<br>
22 format(1x,a3,i10,5g14.6)<br>
<br>
<br>
<br>
<br>
</font><font size="-1">--------------------------------<br>
genmap.f lines 507-... <br>
--------------------------------</font><br>
<font size="-1"><br>
if (nel.lt.1000) then<br>
read(io,50,err=510,end=600) <br>
$ chtemp,<br>
$ cbc(f,e),id1,id2,<br>
$ (bc(ii,f,e),ii=1,nbcrea)<br>
50 format(a1,a3,2i3,5g14.7)<br>
elseif (nel.lt.100 000) then<br>
read(io,51,err=520,end=600) <br>
$ chtemp,<br>
$ cbc(f,e),id1,id2,<br>
$ (bc(ii,f,e),ii=1,nbcrea)<br>
51 format(a1,a3,i5,i1,5g14.7)<br>
elseif (nel.lt.1 000 000) then<br>
read(io,52,err=530,end=600) <br>
$ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)<br>
52 format(1x,a3,i7,5g14.6)<br>
elseif (nel.lt.10 000 000) then<br>
read(io,53,err=540,end=600) <br>
$ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)<br>
53 format(1x,a3,i7,5e20.12)<br>
else<br>
read(io,*,err=550,end=600) <br>
$ cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)<br>
endif<br>
<br>
<br>
<br>
</font>
</body>
</html>