anchor.diff:  
           ! Stress tensor: sxx, syy, szz, sxy, sxz, syz
           if (stress_file.ne.'') then
  -            write(121,*) stations(iSta,:),sts(1,1),sts(2,2),sts(3,3),sts(1,2),sts(1,3),sts(2,3)
  +            if (sta_char(iSta).eq.'#'.or.sta_char(iSta).eq.'>') then
  +                write(121,'(A)') sta_char(iSta)
  +            else
  +                write(121,*) stations(iSta,:),sts(1,1),sts(2,2),sts(3,3),sts(1,2),sts(1,3),sts(2,3)
  +            endif
           endif
   
           ! Maximum (effective) shear stress (aka, second invariant of deviatoric stress tensor): ests
           if (estress_file.ne.'') then
               call max_shear_stress(sts,ests)
  -            write(122,*) stations(iSta,:),ests
  +            if (sta_char(iSta).eq.'#'.or.sta_char(iSta).eq.'>') then
  +                write(122,'(A)') sta_char(iSta)
  +            else
  +                write(122,*) stations(iSta,:),ests
  +            endif
           endif
       endif
   
  
positive.code:     !!$OMP DO
  
  
      ! Calculate deformation from each fault and tensile source at the station
      do iFlt = 1,nfaults+ntensile
  
  
negative.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
          elseif (disp_output_mode.eq.'amz') then
              write(101,*) stations(iSta,:), &
  