[Nek5000-users] nek tools bug ?

[nek5000-users at lists.mcs.anl.gov]

hi all,

I had problems merging grids with the nek tools.
When the final merged grid was having more than 100000 elements genmap 
was failing.

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 nekmerge.f.

just my 2 cents
francesco



------------------------------
n2to3.f lines 441-...
------------------------------

                   if (neln.lt.1000) then
                      write(11,20) cbc(k,e),id,k,(bc(j,k,e),j=1,5)
                   elseif (neln.lt.100000) then
                      write(11,21) cbc(k,e),id,k,(bc(j,k,e),j=1,5)
                   else
                      write(11,22) cbc(k,e),id,k,(bc(j,k,e),j=1,5)
                   endif
....
    20 FORMAT(1x,A3,2I3,5G14.7)
    21 FORMAT(1x,A3,i5,i1,5G14.7)
    22 FORMAT(1x,A3,i6,i1,5G14.7)



----------------------------------
nekmerge.f lines 664-...
----------------------------------

             if (nel.lt.1000) then
                write(11,20) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5)
             elseif (nel.lt.100000) then
                write(11,21) cbc(f,e,fld),e,f,(bc(j,f,e,fld),j=1,5)
             else
                write(11,22) cbc(f,e,fld),e,(bc(j,f,e,fld),j=1,5)
             endif
...
    20       format(1x,a3,2i3,5g14.6)
    21       format(1x,a3,i5,i1,5g14.6)
    22       format(1x,a3,i10,5g14.6)




--------------------------------
genmap.f lines 507-...
--------------------------------

          if (nel.lt.1000) then
             read(io,50,err=510,end=600)
      $      chtemp,
      $      cbc(f,e),id1,id2,
      $      (bc(ii,f,e),ii=1,nbcrea)
    50       format(a1,a3,2i3,5g14.7)
          elseif (nel.lt.100 000) then
             read(io,51,err=520,end=600)
      $      chtemp,
      $      cbc(f,e),id1,id2,
      $      (bc(ii,f,e),ii=1,nbcrea)
    51       format(a1,a3,i5,i1,5g14.7)
          elseif (nel.lt.1 000 000) then
             read(io,52,err=530,end=600)
      $      cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)
    52       format(1x,a3,i7,5g14.6)
          elseif (nel.lt.10 000 000) then
             read(io,53,err=540,end=600)
      $      cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)
    53       format(1x,a3,i7,5e20.12)
          else
             read(io,*,err=550,end=600)
      $      cbc(f,e),id1,(bc(ii,f,e),ii=1,nbcrea)
          endif



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110930/5e162ec1/attachment.html>