anchor.diff:  !
   !  Splitting operators are applied in order
   
  +   use f90sac
  +
      implicit none
   
      type(SACtrace), intent(inout) :: t1, t2, t3
  
positive.code:    integer, parameter, private :: lu_stdin = 5, &
                                    lu_stdout = 6, &
                                    lu_stderr = 0
  
     ! Saved module variables which are private, but kept to prevent continuous
     ! allocation and deallocation (module variables are saved by default)
     type(SACtrace), private :: E_misfit, N_misfit, Z_misfit
  
     ! Hide internal routines
     private :: gaussian1, &
  
negative.code: !     freq     : Dominant frequency of wave
  !     noise    : Fractional amount of random noise
  !     wavetype : 'g'aussian (first-derivative) or 'r'icker
     use f90sac
  
     implicit none
  
     real :: misfit
     type(SACtrace) :: E,N,Z
     real, intent(in) :: phi1, dt1, phi2, dt2, spol
     real, intent(in), optional :: freq, delta, noise
     character(len=*), intent(in), optional :: wavetype
  