anchor.diff:  program susy_qpi
     implicit none
  +  integer, parameter :: dp=kind(0.d0)
  +  integer :: i,j
  +
  +  do i=1,100
  +    do j=1,10
  +      call write_data(i*0.5_dp,j*0.2_dp)
  +    end do
  +  end do
  +
  +end program susy_qpi
  +
  +subroutine write_data(om, del)
  +  implicit none
     external sG0
   
   ! Double precision
     integer, parameter :: dp=kind(0.d0)
  +
  +! Parameters
  +!    omega = om + i * del
  +  real(dp), intent(in) :: om, del
   
   ! Global variables
   !   epsq = -2 * t ( cos(kx) + cos(ky) ) - mu
  
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:   write(log,*) "TOTAL TIME: ", end-start
    close(log)
  
  end subroutine write_data
  
  !==============================================================
  
  