anchor.diff:  c
   c                             Charles A. Williams
   c                       Rensselaer Polytechnic Institute
  -c                        (C) 2004  All Rights Reserved
  +c                        (C) 2005  All Rights Reserved
   c
  -c  Copyright 2004 Rensselaer Polytechnic Institute.
   c  All worldwide rights reserved.  A license to use, copy, modify and
   c  distribute this software for non-commercial research purposes only
   c  is hereby granted, provided that this copyright notice and
  
positive.code: c
  c...  write element connectivities
  c
        do i=1,numelt
          indien=infiel(1,i)
          imat=infiel(2,i)
          ietype=infiel(3,i)
          nen=infetype(2,ietype)
          ngauss=infetype(1,ietype)
          ngpts=ngpts+ngauss
          indtype=ione
          if(ietype.eq.28) indtype=itwo
          if(ietype.eq.29) indtype=ithree
          if(ietype.eq.30) indtype=ifour
          if(ietype.eq.31) indtype=ifive
          if(ietype.gt.31) indtype=isix
          if(ietype.eq.59) indtype=iseven
          if(ietype.eq.60) indtype=ieight
          if(ietype.eq.61) indtype=inine
          if(ietype.eq.62) indtype=10
          do j=1,nen
            itmp(j)=ien(indien+j-1)
          end do
          write(kucd,"(2i7,2x,a4,20i7)") i,imat,eltype(indtype),
       &   (itmp(inducd(j,indtype)),j=1,nen)
        end do
        close(kucd)
  c
  c...  write Gauss point info
  c
        filenm=ucdroot(i1:i2)//".gmesh.inp"
        open(kucd,file=filenm,status="new")
        write(kucd,"(5i7)") ngpts,numelg,nngattr,negattr,nmattr
  c
  c...  write Gauss coordinates
  c
        igpt=izero
        do i=1,numelt
          indien=infiel(1,i)
          ietype=infiel(3,i)
          nen=infetype(2,ietype)
          ngauss=infetype(1,ietype)
          call lcoord(x,xl,ien(indien),nen,numnp)
          do l=1,ngauss
            igpt=igpt+ione
            xg(1)=zero
            xg(2)=zero
            xg(3)=zero
            do j=1,nen
              xg(1)=xg(1)+xl(1,j)*sh(4,j,l,ietype)
              xg(2)=xg(2)+xl(2,j)*sh(4,j,l,ietype)
              xg(3)=xg(3)+xl(3,j)*sh(4,j,l,ietype)
            end do
            write(kucd,"(i7,3(2x,1pe15.8))") igpt,(xg(j),j=1,nsd)
          end do
  
negative.code:       end
  c
  c version
  c $Id: write_ucd_mesh.f,v 1.6 2005/03/26 00:41:01 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  