anchor.diff:    real(dp), intent(out) :: f(nfun)
   
   ! Globar variables
  -  common /gvars/ Pi, t, mu, x0, epsf, Gkinv, Gkqinv
  +  common /gvars/ Pi, t, mu, x0, epsf, V, Uc, Uf, Gkinv, Gkqinv
     common /freq/  omega, qx, qy
  -  real(dp)    :: Pi, t, mu, x0, epsf, qx, qy
  +  real(dp)    :: Pi, t, mu, x0, epsf, V, Uc, Uf, qx, qy
     complex(dp) :: omega, Gkinv(16), Gkqinv(16)
   
   ! Local variables
  
positive.code: ! Global variables
  !     omega  complex.
  !            Frequency (energy) used for G
    common /gvars/ Pi, t, mu, x0, epsf, Gkinv, Gkqinv
    common /freq/  omega, qx, qy
    real(dp)    :: Pi, t, mu, x0, epsf, qx, qy
    complex(dp) :: omega, Gkinv(16), Gkqinv(16)
  
  ! Local variables
  
negative.code:   call cpu_time(start)
    call date_and_time(date,time)
  
  ! frequency
    omega = dcmplx(om,del)
  
  ! Save results to file
    open(log, file="susy_qpi.log", position="append", status="unknown")
    write(somega, '(f0.2,"+",f0.2,"i")') om, del
    filename = date//"_"//time//"_w="//trim(somega)//"_susy_qpi.dat"
    write(log,*)
    write(log,*) "======================================================="
    write(log,*)
    write(log,*) "START: "//filename
    write(log, *) "t:     ", t
    write(log, *) "mu:    ", mu
    write(log, *) "x0:    ", x0
    write(log, *) "epsf:  ", epsf
    write(log, *) "V:     ", V
    write(log, *) "Uc:    ", Uc
    write(log, *) "Uf:    ", Uf
    write(log, *) "omega: ", omega
    close(log)
    open(dat, file=filename, status="unknown")
  
  