anchor.diff:  
       elseif (tag.eq.'-az') then
           disp_output_mode = 'amz'
  +
  +    elseif (tag.eq.'-gmt') then
  +        i = i + 1
  +        call get_command_argument(i,gmt_fault_file)
   
       elseif (trim(tag).eq.'-v'.or.trim(tag).eq.'-verbose'.or.trim(tag).eq.'-verbosity') then
           i = i + 1
  
positive.code: character(len=512) :: mag_file                    ! source faults: ... mag format
  character(len=512) :: flt_file                    ! source faults: ... slip wid len format
  character(len=512) :: tns_file                    ! volume source: ... opening wid len
  logical :: isFaultFileDefined                     ! input fault tag
  logical :: isTensileFileDefined                   ! input tensile source tag
  character(len=16) :: fault_type                   ! point or finite source tag
  
negative.code:                    fsp_file, &
                     mag_file, &
                     flt_file, &
                     gmt_fault_file, &
                     isFaultFileDefined, &
                     empirical_relation, &
                     slip_threshold, &
  