anchor.diff:    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
  +    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
  +    write(log,*) "ERROR: Singular matrix "
  +    write(log,*) Gkqinv
       close(log)
     end if
   
  
positive.code:   implicit none
    include "inc/susy_qpi.f90"
  
  
  
  ! Set constants
  
negative.code:   omega = dcmplx(om,del)
  
  ! Save results to file
    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 = trim(somega)//"_susy_qpi.dat"
    write(log,*)
    write(log,*) "======================================================="
    write(log,*)
  