anchor.diff:              if(nwinkx.ne.izero.and.unlckf) call unlock(zcg,btot,idx,
        &       iwinkx,idhist,ibond,diforc,histry,nstep,numnp,nwinkx,nhist,
        &       neq,numdif,lastep,itwo)
  +cdebug            write(6,*) "After unlock:"
   	    call daxpy(neq,one,zcg,ione,btot,ione)
   c
   c...  reform time-dependent material and stiffness matrices if
  
positive.code: c...  open output files for appending, if necessary
  c
  cdebug      write(6,*) "Hello from viscos_f!"
  cdebug      write(6,*) "From viscos_f, ngauss,nsd,gauss:",ngauss,nsd,
  cdebug     & (gauss(idb),idb=1,(nsd+1)*ngauss)
  c
        if(idout.gt.ione) open(kw,file=ofile,status="old",access="append")
        if(idsk.eq.izero) open(kp,file=pfile,status="old",access="append")
  
negative.code:       end
  c
  c version
  c $Id: viscos.f,v 1.4 2004/08/12 02:34:54 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  