anchor.diff:  
   ! Save results to file
     open(log, file="susy_qpi.log", position="append", status="unknown")
  +  write(somega, '(f6.2,"+",f6.2,"i")') om, del
  +  filename = date//"_"//time//"_w="//trim(somega)//"_susy_qpi.dat"
     write(log,*)
     write(log,*) "======================================================="
     write(log,*)
  -  write(log,*) "START: "//date//"_"//time//"_susy_qpi.dat"
  +  write(log,*) "START: "//filename
     close(log)
  -  open(dat, file=date//"_"//time//"_susy_qpi.dat", status="unknown")
  +  open(dat, file=filename, status="unknown")
   
   ! Set constants
     Pi=4.0_dp*datan(1.0_dp)
  
positive.code:   Uf  = 1e-3_dp
  
  ! frequency
    omega = dcmplx(20.0_dp,0.5_dp)
  
  ! Set up matrix as an array of size 16
    Gkinv = (/ complex(dp) :: &
  
negative.code:   write(log,*) "TOTAL TIME: ", end-start
    close(log)
  
  end subroutine write_data
  
  !==============================================================
  
  