anchor.diff:              call usage('calc_deformation: error in progress_indicator (usage:none)')
           endif
       endif
  +
  +
  +    ! End the ordered directive, so threads can be FREE!!!!
  +    !$OMP END ORDERED
  +
  +
   enddo
  +
  +
  +! End the parallel do loop
  +!$OMP END DO
  +
  +! Stop parallelizing - the rest of the subroutine is serial
  +!$OMP END PARALLEL
  +
   
   ! Close files that were opened for writing
   inquire(file=displacement_file,number=file_unit,opened=isThisUnitOpen)
  
positive.code:     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
  
negative.code:                     tshrmx, tnor, coul
  double precision :: nvec(3), svec(3), tstr, tupd
  
  ! integer :: OMP_GET_NUM_THREADS, OMP_GET_THREAD_NUM
  
  
  if (verbosity.ge.1) then
      write(stdout,*) 'calc_deformation: starting'
  