anchor.diff:          endif
       enddo
   
  +    endif ! (sta_char(iSta).ne.'#'.and.sta_char(iSta).ne.'>')
  +
   
       !! End the parallel do loop
       !!$OMP END DO
  
positive.code: 
          ! 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)
          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
          endif
      endif
  
  
negative.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
  
  