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: 
      elseif (tag.eq.'-az') then
          disp_output_mode = 'amz'
  
      elseif (trim(tag).eq.'-v'.or.trim(tag).eq.'-verbose'.or.trim(tag).eq.'-verbosity') then
          i = i + 1
  
negative.code: flt_file = ''
  mag_file = ''
  tns_file = ''
  gmt_fault_file = ''
  isFaultFileDefined = .false.
  isTensileFileDefined = .false.
  fault_type = 'rect'
  