anchor.diff:  
   
   !  reference values for self test
  -  allocate (res_reference(nr_variables))
  +  allocate (res_reference(num_variables))
     res_reference(1) = -0.97182634532897016d0
     res_reference(2) = -8.9803044914016916d-3
     res_reference(3) = -1.3470456737102541d-2
  
positive.code:   !  and contract them with the first order densities
  !  first set up xcfun for first derivatives
    order = 1
    nr_variables = 4
    ierr = xc_eval_setup(id, XC_N_NX_NY_NZ, XC_CONTRACTED, order)
    if (ierr /= 0) then
      print *, 'xc_eval_setup failed with error ', ierr
  
negative.code:   deallocate (res_reference)
  
  !  reference values for self test
    allocate (res_reference(num_variables))
    res_reference(1) = -2.0795476275887927d0
    res_reference(2) = 2.1489169824673575d-2
    res_reference(3) = 4.1214059228412030d-2
  