anchor.diff:  !     Log results
         if (ifail>0) then
           open(log, file="susy_qpi.log", position="append", status="old")
  -        write (log,*) "ERROR: DCUHRE exit code ", ifail
  +        write(log,*) "ERROR: DCUHRE exit code ", ifail
           close(log)
         end if
         write(dat,*)  qx,  qy, result(1)
  
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:     end do
    end do
  
    open(10, file="susy_qpi.log", position="append", status="old")
    write(10,*) "********************END: susy_qpi********************"
    close(10)
  end program susy_qpi
  
  !==============================================================
  ! write_data  Writes data to OMEGA_susy_qpi.dat
  !   om        real part of frequency
  !   del       imaginary part of frequency
  !==============================================================
  