anchor.diff:  c                                          Young's modulus
   c                                          Poisson's ratio
   c
  -      subroutine mat_prt_1(prop,nprop,matnum,idout,idsk,kw,kp,ierr,
  -     & ofile,pfile)
  +      subroutine mat_prt_1(prop,nprop,matnum,idout,idsk,kw,kp,
  +     & ofile,pfile,ierr,errstrng)
   c
   c...  subroutine to output material properties for material model 1.
   c
   c     Error codes:
  -c         0:  No error
  -c         1:  Error opening output file
  -c         2:  Not used
  -c         3:  Write error
  -c
  -      include "implicit.inc"
  +c         2:  Error opening output file
  +c         4:  Write error
  +c
  +      include "implicit.inc"
  +c
  +c...  parameter definitions
  +c
  +      include "nconsts.inc"
   c
   c...  subroutine arguments
   c
         integer nprop,matnum,idout,idsk,kw,kp,ierr
         double precision prop(nprop)
  -      character ofile*(*),pfile*(*)
  +      character ofile*(*),pfile*(*),errstrng*(*)
   c
   c...  parameters
   c
  
positive.code: c
  c...  open output files
  c
        ierr=0
        if(idout.gt.0) open(kw,file=ofile,status="old",access="append")
        if(idsk.eq.0) open(kp,file=pfile,status="old",access="append")
        if(idsk.eq.1) open(kp,file=pfile,status="old",access="append",
       & form="unformatted")
  c
  c...  output plot results
  c
        if(idsk.eq.0) then
  	write(kp,"(3i7)") matnum,mattype,nprop
  	write(kp,"(1pe15.8,20(2x,1pe15.8))") (prop(i),i=1,nprop)
        else if(idsk.eq.1) then
  	write(kp) matnum,mattype,nprop
  	write(kp) prop
        end if
  c
  c...  output ascii results, if desired
  c
        if(idout.gt.0) then
  	write(kw,700) matnum,modelname,nprop
  	do i=1,nprop
  	  write(kw,710) labelp(i),prop(i)
          end do
        end if
  c
   700  format("Material number:   ",i7,/,10x,a80,/,
       &       "Number of properties:  ",i7,/)
   710  format(15x,a15,3x,1pe15.8)
  c
        return
        end
  c
  c
  
negative.code: c     solution.
  c     Note that only the upper triangle is used (or available), as
  c     dmat is assumed to always be symmetric.
  c
        include "implicit.inc"
  c
  