anchor.diff:  double precision, allocatable :: faults(:,:)      ! Fault parameter array
   double precision, allocatable :: tensile(:,:)     ! Tensile source parameter array
   double precision, allocatable :: stations(:,:)    ! Station location array
  +character(len=1), allocatable :: sta_char(:)      ! Station character (optional)
   double precision, allocatable :: targets(:,:)     ! Target/receiver geometry array
   double precision :: centroid(3)                   ! Moment weighted mean position of EQ
   
  
positive.code:     !!$OMP DO
  
  
      ! Calculate deformation from each fault and tensile source at the station
      do iFlt = 1,nfaults+ntensile
  
  
negative.code: use o92util, only: station_file, &
                     nstations, &
                     stations, &
                     sta_char, &
                     isStationFileDefined, &
                     coord_type, &
                     auto_mode, &
  