anchor.diff:          integer, intent(out) :: nspawn
           type(excit), intent(out) :: connection
   
  -        real(p),parameter      :: P__=0.2, Po_=(1.0-P__)*0.5, P_o=Po_
  +        real(p),parameter      :: P__=0.05, Po_=(1.0-P__)*0.5, P_o=Po_
           real(p),parameter, dimension(3) :: P_double_elt_type =(/P__, Po_, P_o /)
   
           real(p)          :: choose_double_elt_type
  
positive.code:         ! **however, the values of P__ etc will ideally be reference-determinant dependent, what is the best way to order this sequence?**
  
  
  elttype:if(choose_double_elt_type <= P__ ) then
              
              !      __   __ 
              !    ./  \./  \.
  
negative.code:             call create_excited_det(cdet_excit%f, connection_jk, f_excit_2)
              ! (ii)  calculate the connection to this excited determinant
              connection = get_excitation(cdet%f,f_excit_2)
              ! 6. deallocate the cdet_excit (allocated in create_excited_det_complete)
              call dealloc_det_info(cdet_excit)
  
  
           !******************* this code does not work in the case of looping back on itself *********************
  