anchor.diff:    omega = dcmplx(om,del)
   
   ! Save results to file
  -  open(log, file="susy_qpi.log", position="append", status="unknown")
  +  open(log, file="susy_qpi.log", position="append", status="old")
     write(somega, '(f0.2,"+",f0.2,"i")') om, del
  -  filename = date//"_"//time//"_w="//trim(somega)//"_susy_qpi.dat"
  +!   filename = date//"_"//time//"_w="//trim(somega)//"_susy_qpi.dat"
  +  filename = trim(somega)//"_susy_qpi.dat"
     write(log,*)
     write(log,*) "======================================================="
     write(log,*)
  
positive.code:   call invert4x4(Gkinv, Gk, status)
    if ( status>0 ) then
      open(log, file="susy_qpi.log", position="append", status="old")
      write (log,*) "ERROR: Singular matrix "
      write (log,*) Gkinv
      close(log)
    end if
    call invert4x4(Gkqinv, Gkq, status)
    if ( status>0 ) then
      open(log, file="susy_qpi.log", position="append", status="old")
      write (log,*) "ERROR: Singular matrix "
      write (log,*) Gkqinv
      close(log)
    end if
  
  
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)
  