anchor.diff:  c...  parameters
   c
         double precision zero
  -      integer ione
  -      parameter(zero=0.0d0,ione=1)
  +      parameter(zero=0.0d0)
   c
   c...  local variables
   c
         double precision str(6)
         integer i
   c
  -      call dscal(nstr,zero,str,ione)
  +      call fill(str,zero,nstr)
         if(idout.gt.0) write(kw,700) iel,igauss,(str(i),i=1,nstr)
         if(idsk.eq.0) write(kp,710) (str(i),i=1,nstr)
         if(idsk.eq.1) write(kp) str
  
positive.code: c       
  
  c version
  c $Id: mat1.f,v 1.2 2004/05/20 21:13:21 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  
negative.code:       integer i,j
        double precision e,pr,pr1,pr2,pr3,fac,dd,od,ss
  c
        e=prop(2)
        pr=prop(3)
        pr1=1.0d0-pr
  