anchor.diff:    omega = dcmplx(om,del)
   
   ! Save results to OMEGA_susy_qpi.dat
  -  write(somega, '("w=",f0.2,"+",f0.2,"i")') om, del
  -  filename = date//"_"//time//"/"//trim(somega)//"_susy_qpi.dat"
  +  write(somega, '("v=", f0.2, "w=", f0.2, "+", f0.2, "i")') v, om, del
  +  filename = date//"_"//time//"/"//trim(somega)//".dat"
     open(dat, file=filename, status="new")
   
   ! Start log
  
positive.code:   write(10, *) "Uf:    ", Uf
    close(10)
  
  ! Set up matrix as an array of size 16
    Gkinv = (/ complex(dp) :: &
         0,  V, -Uc,  0, &
  
negative.code:      -Uc,  0,   0,  V, &
         0, Uf,   V,  0 /)
  
      do j=0,20
        call write_data(j*0.1_dp,0.1_dp)
      end do
    end do
  
  