anchor.diff:    implicit none
     include "inc/susy_qpi.f90"
   
  +  open(10, file="susy_qpi.log", position="append", status="unknown")
  +  write(10,*) "******************START: susy_qpi********************"
  +  close(10)
   
   
   ! Set constants
  
positive.code:   omega = dcmplx(om,del)
  
  ! Save results to file
    open(log, file="susy_qpi.log", position="append", status="unknown")
    write(somega, '(f0.2,"+",f0.2,"i")') om, del
    filename = date//"_"//time//"_w="//trim(somega)//"_susy_qpi.dat"
    write(log,*)
    write(log,*) "======================================================="
    write(log,*)
  
negative.code: !     Log results
        if (ifail>0) then
          open(log, file="susy_qpi.log", position="append", status="old")
          write(log,*) "ERROR: DCUHRE exit code ", ifail
          close(log)
        end if
        write(dat,*)  qx,  qy, result(1)
  