anchor.diff:     call f90sac_covar2(N,e,c)
      eig = sw_eig2x2(c)
      misfit = (misfit + minval(eig)/maxval(eig))/2.
  -   call f90sac_deletetrace(E)
  -   call f90sac_deletetrace(N)
  -   call f90sac_deletetrace(Z)
   
   contains
   
  
positive.code:       endif
        sum_dt = sum_dt + dt
     enddo
  
     !  Set window markers for SHEBA by moving end window along by the applied tshift
     t1%f = t1%f + sum_dt - modulo(sum_dt,t1%delta)
  
negative.code:    real :: dt, theta, sum_dt
     character(len=8) :: kcmpnm1, kcmpnm2
     integer :: i
     integer :: f90sac_suppress_warnings_previous
  
     ! See if we're being quiet or not: default to silet
     silent = .true.
  