anchor.diff:              write(stdout,*) stn
           endif
       endif
  -
   
   
       ! Displacement: ux, uy, uz  OR  az, uhor, uz
  
positive.code:         write(stdout,'(X,"(DEBUG)",X,A5,I6,3F12.3)') 'ista=',iSta,stations(iSta,:)
      endif
  
  !$OMP PARALLEL PRIVATE(iFlt,ierr), &
  !$OMP&         PRIVATE(evlo,evla,evdp,str,dip,rak,slip_mag,wid,len,dist,az), &
  !$OMP&         PRIVATE(slip,mom,disptmp,stntmp), &
  !$OMP&         FIRSTPRIVATE(sta_coord), &
  !$OMP&         REDUCTION(+:disp), REDUCTION(+:stn)
  
  !$OMP DO
  
      ! Calculate deformation from each fault and tensile source at the station
      do iFlt = 1,nfaults+ntensile
  
negative.code:         endif
      enddo
  
  
      ! End the parallel do loop
      !$OMP END DO
  
      ! Stop parallelizing - the rest of the calculation for this station is serial
      !$OMP END PARALLEL
  
  
  
      if (debug_mode.eq.'calc_deformation'.or.debug_mode.eq.'all') then
          if(iWantDisp) then
  