anchor.diff:  
   i = 1
   do while (i.le.narg)
  +
       call get_command_argument(i,tag)
  +
       if (trim(tag).eq.'-flt') then
           i = i + 1
           call get_command_argument(i,fault_file)
  +
           isFaultFileDefined = .true.
  +
       elseif (trim(tag).eq.'-sta') then
           i = i + 1
           call get_command_argument(i,station_file)
  +
           isStationFileDefined = .true.
  +
       elseif (trim(tag).eq.'-trg') then
           i = i + 1
           call get_command_argument(i,target_file)
  +
           isTargetFileDefined = .true.
  +
       elseif (trim(tag).eq.'-haf') then
           i = i + 1
           call get_command_argument(i,halfspace_file)
  +
       elseif (trim(tag).eq.'-disp') then
           i = i + 1
           call get_command_argument(i,displacement_file)
  +
           isOutputFileDefined = .true.
           iWantDisp = .true.
  +
       elseif (trim(tag).eq.'-strain') then
           i = i + 1
           call get_command_argument(i,strain_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
  +
       elseif (trim(tag).eq.'-stress') then
           i = i + 1
           call get_command_argument(i,stress_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
           iWantStress = .true.
  +
       elseif (trim(tag).eq.'-estress') then
           i = i + 1
           call get_command_argument(i,estress_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
           iWantStress = .true.
  +
       elseif (trim(tag).eq.'-normal') then
           i = i + 1
           call get_command_argument(i,normal_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
           iWantStress = .true.
           iWantTraction = .true.
  +
       elseif (trim(tag).eq.'-shear'.or.trim(tag).eq.'-shearmax') then
           i = i + 1
           call get_command_argument(i,shear_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
           iWantStress = .true.
           iWantTraction = .true.
  +
       elseif (trim(tag).eq.'-coul') then
           i = i + 1
           call get_command_argument(i,coulomb_file)
  +
           isOutputFileDefined = .true.
           iWantStrain = .true.
           iWantStress = .true.
           iWantTraction = .true.
  +
       elseif (trim(tag).eq.'-geographic'.or.trim(tag).eq.'-geo') then
           coord_type = 'geographic'
  +
       elseif (trim(tag).eq.'-cartesian'.or.trim(tag).eq.'-xy') then
           coord_type = 'cartesian'
  -    elseif (trim(tag).eq.'-poisson') then
  -        i = i + 1
  -        call get_command_argument(i,tag)
  -        read(tag,*,iostat=ios,err=9001) poisson
  -    elseif (trim(tag).eq.'-lame') then
  -        i = i + 1
  -        call get_command_argument(i,tag)
  -        read(tag,*,iostat=ios,err=9001) lame
  -    elseif (trim(tag).eq.'-mu') then
  -        i = i + 1
  -        call get_command_argument(i,tag)
  -        read(tag,*,iostat=ios,err=9001) mu
  +
       elseif (trim(tag).eq.'-v'.or.trim(tag).eq.'-verbose'.or.trim(tag).eq.'-verbosity') then
           i = i + 1
           call get_command_argument(i,tag)
           read(tag,*,iostat=ios,err=9001) verbosity
  +
       else
           call usage('!! gcmdln: no option '//trim(tag))
  -    endif
  +
  +    endif
  +
       i = i + 1
  +
   enddo
   
   if (verbosity.ge.1) then
  -    write(0,*) 'gcmdln read:'
  -    write(0,*) 'fault_file:          ',fault_file
  -    write(0,*) 'station_file:        ',station_file
  -    write(0,*) 'target_file:         ',target_file
  -    write(0,*) 'halfspace_file:      ',halfspace_file
  -    write(0,*) 'displacement_file:   ',displacement_file
  -    write(0,*) 'strain_file:         ',strain_file
  -    write(0,*) 'stress_file:         ',stress_file
  -    write(0,*) 'estress_file:        ',estress_file
  -    write(0,*) 'normal_file:         ',normal_file
  -    write(0,*) 'shear_file:          ',shear_file
  -    write(0,*) 'coulomb_file:        ',coulomb_file
  -    write(0,*) 'isFaultFileDefined:  ',isFaultFileDefined
  -    write(0,*) 'isStationFileDefined:',isStationFileDefined
  -    write(0,*) 'isTargetFileDefined: ',isTargetFileDefined
  -    write(0,*) 'isOutputFileDefined: ',isOutputFileDefined
  -    write(0,*) 'iWantDisp:           ',iWantDisp
  -    write(0,*) 'iWantStrain:         ',iWantStrain
  -    write(0,*) 'iWantStress:         ',iWantStress
  -    write(0,*) 'iWantTraction:       ',iWantTraction
  -    write(0,*) 'poisson:             ',poisson
  -    write(0,*) 'lame:                ',lame
  -    write(0,*) 'mu:                  ',mu
  -    write(0,*) 'coord_type:          ',coord_type
  -    write(0,*) 'verbosity:           ',verbosity
  +    write(stderr,*) 'gcmdln read:'
  +    write(stderr,*) 'fault_file:          ',fault_file
  +    write(stderr,*) 'station_file:        ',station_file
  +    write(stderr,*) 'target_file:         ',target_file
  +    write(stderr,*) 'halfspace_file:      ',halfspace_file
  +    write(stderr,*) 'displacement_file:   ',displacement_file
  +    write(stderr,*) 'strain_file:         ',strain_file
  +    write(stderr,*) 'stress_file:         ',stress_file
  +    write(stderr,*) 'estress_file:        ',estress_file
  +    write(stderr,*) 'normal_file:         ',normal_file
  +    write(stderr,*) 'shear_file:          ',shear_file
  +    write(stderr,*) 'coulomb_file:        ',coulomb_file
  +    write(stderr,*) 'isFaultFileDefined:  ',isFaultFileDefined
  +    write(stderr,*) 'isStationFileDefined:',isStationFileDefined
  +    write(stderr,*) 'isTargetFileDefined: ',isTargetFileDefined
  +    write(stderr,*) 'isOutputFileDefined: ',isOutputFileDefined
  +    write(stderr,*) 'iWantDisp:           ',iWantDisp
  +    write(stderr,*) 'iWantStrain:         ',iWantStrain
  +    write(stderr,*) 'iWantStress:         ',iWantStress
  +    write(stderr,*) 'iWantTraction:       ',iWantTraction
  +    write(stderr,*) 'coord_type:          ',coord_type
  +    write(stderr,*) 'verbosity:           ',verbosity
   endif
   
   9001 if (ios.ne.0) then
  -    write(0,*) 'gcmdln: error reading string ',trim(tag)
  +    write(stderr,*) 'gcmdln: error reading string ',trim(tag)
       call usage('exiting triutil from gcmdln')
   endif
   
  
positive.code: !--------------------------------------------------------------------------------------------------!
  
  subroutine calc_displacement(disp,x,y,z)
  use tri_disloc, only: tri_disloc_disp, tri_center
  use triutil_module, only: nfaults, faults, poisson, coord_type
  implicit none
  ! Arguments
  double precision :: disp(3), x, y, z
  ! Local variables
  double precision :: sta_coord(3), tri_coord(3,4), slip(3), tri_coord_new(3,4)
  double precision :: center(3), dist, az
  
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')
  