anchor.diff:  
   !--------------------------------------------------------------------------------------------------!
   
  -subroutine line_count(nlines,filename)
  -implicit none
  -integer nlines, ios
  -character(len=*) :: filename
  -open(unit=11,file=filename,status='old')
  -nlines = 0
  -do
  -    read(11,*,iostat=ios)
  -    if (ios.eq.0) then
  -        nlines = nlines + 1
  -    else
  -        exit
  -    endif
  -enddo
  -close(11)
  -return
  -end subroutine line_count
  -
  -!--------------------------------------------------------------------------------------------------!
  -
   subroutine gcmdln()
  -use triutil_module, only: fault_file, station_file, target_file, halfspace_file, &
  -                          displacement_file, &
  -                          strain_file, &
  -                          stress_file, estress_file, &
  -                          normal_file, shear_file, coulomb_file, &
  -                          coord_type, &
  -             isFaultFileDefined, isStationFileDefined, isTargetFileDefined, isOutputFileDefined, &
  -             iWantDisp, iWantStrain, iWantStress, iWantTraction, poisson, lame, mu, verbosity
  -implicit none
  +
  +use io, only: stderr, verbosity
  +use triutil, only: fault_file, &
  +                   station_file, &
  +                   target_file, &
  +                   halfspace_file, &
  +                   displacement_file, &
  +                   strain_file, &
  +                   stress_file, &
  +                   estress_file, &
  +                   normal_file, &
  +                   shear_file, &
  +                   coulomb_file, &
  +                   coord_type, &
  +                   isFaultFileDefined, &
  +                   isStationFileDefined, &
  +                   isTargetFileDefined, &
  +                   isOutputFileDefined, &
  +                   iWantDisp, &
  +                   iWantStrain, &
  +                   iWantStress, &
  +                   iWantTraction, &
  +                   poisson, &
  +                   lame, &
  +                   shearmod
  +implicit none
  +
   ! Local variables
   integer :: i, narg, ios
   character(len=512) :: tag
   
  +! Initialize control variables
   fault_file = ''
   station_file = ''
   target_file = ''
  
positive.code: endif
  if (iWantStress) then
      if (verbosity.ge.1) then
          write(0,*) 'calc_deformation: calculating stresses'
      endif
      if (stress_file.ne.'') then
          open(unit=121,file=stress_file,status='unknown')
  
negative.code: endif
  if (iWantTraction) then
      if (verbosity.ge.1) then
          write(stderr,*) 'calc_deformation: calculating resolved tractions'
      endif
      if (normal_file.ne.'') then
          open(unit=131,file=normal_file,status='unknown')
  