anchor.diff:  
       ! Displacement: ux, uy, uz  OR  az, uhor, uz
       if (displacement_file.ne.'') then
  -        if (sta_char(iSta).eq.'#'.or.sta_char(iSta).eq.'>') then
  +        if (sta_char(iSta)(1:1).eq.'#'.or.sta_char(iSta)(1:1).eq.'>') then
               write(101,'(A)') sta_char(iSta)
           elseif (disp_output_mode.eq.'enz') then
               write(101,*) stations(iSta,:),disp
  
positive.code:     !!$OMP DO
  
  
      if (sta_char(iSta).ne.'#'.and.sta_char(iSta).ne.'>') then
  
      ! Calculate deformation from each fault and tensile source at the station
      do iFlt = 1,nfaults+ntensile
  
negative.code: do i = 1,nstations
      read(13,'(A)') input_line
      input_line = adjustl(input_line)
      sta_char(i) = input_line
      if (sta_char(i)(1:1).eq.'#'.or.sta_char(i)(1:1).eq.'>') then
          ! This line is a comment (#) or a segment header (>)
          cycle
      endif
  