anchor.diff:      write(stdout,*) 'iWantStrain:              ',iWantStrain
       write(stdout,*) 'iWantStress:              ',iWantStress
       write(stdout,*) 'iWantTraction:            ',iWantTraction
  +    write(stdout,*) 'isParallel:               ',isParallel
  +    write(stdout,*) 'nthreads:                 ',nthreads
       write(stdout,*) 'coord_type:               ',trim(coord_type)
       write(stdout,*) 'prog:                     ',iWantProg
       write(stdout,*) 'debug_mode:               ',trim(debug_mode)
  
positive.code:                     tshrmx, tnor, coul
  double precision :: nvec(3), svec(3), tstr, tupd
  
  ! integer :: OMP_GET_NUM_THREADS, OMP_GET_THREAD_NUM
  
  
  if (verbosity.ge.1) then
  
negative.code: double precision, allocatable :: targets(:,:)     ! Target/receiver geometry array
  double precision :: centroid(3)                   ! Moment weighted mean position of EQ
  
  ! Parallel options
  logical :: isParallel
  integer :: nthreads
  
  ! Debugging
  character(len=32) :: debug_mode                   ! Verbose output for debugging
  
  end module
  