anchor.diff:        end
   c
   c version
  -c $Id: read_connect.f,v 1.6 2005/02/24 00:17:58 willic3 Exp $
  +c $Id: read_connect.f,v 1.7 2005/03/28 19:53:38 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code:       if(idout.gt.izero) then
          open(kw,file=ofile,status="old",access="append")
          npage=50
          do i=1,numelt
            imat=infiel(2,i)
            ietype=infiel(3,i)
            ngauss=infetype(1,ietype)
            nen=infetype(2,ietype)
            i1=infiel(1,i)
            i2=infiel(1,i)+nen-1
            if(i.eq.ione.or.mod(n,npage).eq.izero) then
              write(kw,2000)(head(j),j,j=1,nen)
              write(kw,2500)
            end if
            write(kw,3000) i,imat,ietype,ngauss,(ien(j),j=i1,i2)
          end do
          write(kw,4000)
          close(kw)
  
negative.code:           errstrng="read_connect"
            return
          end if
  c
  c...  store element nodes in the proper family and update the pointer
  c     array for that family
  c
          iloc=ivftmp(imat)
          call iquate(ien(1,iloc),itmp,nen)
          ivftmp(imat)=ivftmp(imat)+ione
        end do
        close(kr)
  c
  c...  output plot info, if desired
  c
        if(idsk.eq.ione) then
          open(kp,file=pfile,err=40,status="old",access="append")
          write(kp,1000,err=50) numelv
          ielg=izero
          do ifam=1,nvfamilies
            nelfamily=ivfamily(1,ifam)
            matmod=ivfamily(2,ifam)
            do ielf=1,nelfamily
              ielg=ielg+ione
              write(kp,1000) ielg,matmod,ietypev,ngauss,
       &       (ien(j,ielg),j=1,nen)
            end do
          end do
          close(kp)
        else if(idsk.eq.itwo) then
          open(kp,file=pfile,err=40,status="old",access="append",
       &   form="unformatted")
          write(kp,err=50) numelv,nen,ngauss
          write(kp,err=50) ivfamily
          write(kp,err=50) ien
          close(kp)
        end if
  