anchor.diff:    end do
   
   end program susy_qpi
  +
  +!==============================================================
  +! write_data  Writes data to YYYYMMDD_HHSS_w=OMEGA_susy_qpi.dat
  +!   om        real part of frequency
  +!   del       imaginary part of frequency
  +!==============================================================
   
   subroutine write_data(om, del)
     implicit none
  
positive.code: program susy_qpi
    implicit none
  
  
negative.code: end subroutine write_data
  
  !==============================================================
  ! sG0         Integrand for dcuhre
  !   ndim      number of variables
  !   z         variables for integration
  !   nfun      number of components of integral
  !   f         value of integrand
  !==============================================================
  
  subroutine sG0(ndim, z, nfun, f)
    implicit none
  