anchor.diff:  call test_value(temp,0.13533300490703207d0,'driver2 (anneal_dp_array): final temp')
   call test_value(obj,-5.8331140982675264d-002,'driver2 (anneal_dp_array): final obj')
   close(12,status='delete')
  +write(stdout,*) 'subroutine anneal_dp_array() - driver2() passed unit test'
  +write(stdout,*)
   
   
   ! driver3 (anneal_dp_array; module)
  
positive.code: 
  implicit none
  
  integer :: j, ios, it, model_best_int(2)
  character(len=512) :: line
  character(len=1) :: ch
  double precision :: temp, obj, model_best_dp(2)
  
  
  ! driver1 (anneal_int_array)
  
negative.code: call test_value(temp,0.13533300490703207d0,'driver1 (anneal_int_array): final temp')
  call test_value(obj,-0.34999999999999998d0,'driver1 (anneal_int_array): final obj')
  close(11,status='delete')
  write(stdout,*) 'subroutine anneal_int_array() - driver1() passed unit test'
  write(stdout,*)
  
  
  ! driver2 (anneal_dp_array)
  