anchor.diff:    Uf  = 1e-3_dp
   
   ! frequency
  -  omega = dcmplx(20.0_dp,0.5_dp)
  +  omega = dcmplx(om,del)
   
   ! Set up matrix as an array of size 16
     Gkinv = (/ complex(dp) :: &
  
positive.code:   real             :: start, end
    character(len=8) :: date
    character(len=4) :: time
  
  ! Variables for dcuhre
  !     ndim   integer.
  !            number of variables. 1 < ndim <=  15.
  
negative.code: 
  ! 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: "//filename
    close(log)
    open(dat, file=filename, status="unknown")
  
  ! Set constants
    Pi=4.0_dp*datan(1.0_dp)
  