anchor.diff:      !!$OMP DO
   
   
  -    if (sta_char(iSta).ne.'#'.and.sta_char(iSta).ne.'>') then
  +    if (sta_char(iSta)(1:1).ne.'#'.and.sta_char(iSta)(1:1).ne.'>') then
   
       ! Calculate deformation from each fault and tensile source at the station
       do iFlt = 1,nfaults+ntensile
  
positive.code: 
      ! Displacement: ux, uy, uz  OR  az, uhor, uz
      if (displacement_file.ne.'') then
          if (sta_char(iSta).eq.'#'.or.sta_char(iSta).eq.'>') then
              write(101,'(A)') sta_char(iSta)
          elseif (disp_output_mode.eq.'enz') then
              write(101,*) stations(iSta,:),disp
  
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
  