anchor.diff:  endif
   if (iWantTraction) then
       if (verbosity.ge.1) then
  -        write(0,*) 'calc_deformation: calculating resolved tractions'
  +        write(stderr,*) 'calc_deformation: calculating resolved tractions'
       endif
       if (normal_file.ne.'') then
           open(unit=131,file=normal_file,status='unknown')
  
positive.code: module triutil_module
  
  character(len=512) :: fault_file       ! source faults (x1 y1 z1 x2 y2 z2 x3 y3 z3 ss ds ts)
  character(len=512) :: station_file     ! station locations (stlo stla stdp)
  
negative.code: !--------------------------------------------------------------------------------------------------!
  
  subroutine read_stations()
  !----
  ! Read the input station file, in format:
  !     x1 y1 z1
  !----
  
  use io, only: stderr, line_count
  use triutil, only: station_file, nstations, stations, isStationFileDefined
  implicit none
  
  ! Local variables
  integer :: i, j, ios
  character(len=512) :: input_line
  