anchor.diff:      write(stderr,*) '-parallel [NTHREADS] Calculate deformation in parallel'
       write(stderr,*) '-geo|-xy             Use geographic (default) or cartesian coordinates'
       write(stderr,*) '-az                  Displacement vector outputs (AZ HMAG Z)'
  +    write(stderr,*) '-gmt FILE            GMT psxy -SJ file (lon lat str len proj_wid)'
       write(stderr,*) '-prog                Turn on progress indicator'
       write(stderr,*) '-v LVL               Turn on verbose mode'
       write(stderr,*) '-debug [ROUTINE]     Turn on debugging'
  
positive.code:     enddo
  
      nfaults = nfaults + mag_data%nflt
  
      if (debug_mode.eq.'read_faults'.or.debug_mode.eq.'all') then
          write(stdout,*) '(DEBUG) read_faults: finished read_mag'
  
negative.code: flt_file = ''
  mag_file = ''
  tns_file = ''
  gmt_fault_file = ''
  isFaultFileDefined = .false.
  isTensileFileDefined = .false.
  fault_type = 'rect'
  