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...  parameters
  c
        double precision zero
        integer ione
        parameter(zero=0.0d0,ione=1)
  c
  c...  local variables
  c
        double precision str(6)
        integer i
  c
        call dscal(nstr,zero,str,ione)
        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
  
negative.code: c...  parameters
  c
        character labelp(3)*15,modelname*24
        data labelp/"Density",
       &            "Young's modulus",
       &            "Poisson's ratio"/
        parameter(modelname="Isotropic Linear Elastic")
        integer mattype
        parameter(mattype=1)
  c
  c...  local variables
  c
        integer i
  c
  c...  open output files
  c
  