anchor.diff:    logical :: ok
     logical, external :: detEq
     integer, external :: omp_get_thread_num
  -
  +  double precision :: contrib(N_states)
  +  
  +  contrib = 0d0
     delta_ij_mrcc = 0d0
     delta_ii_mrcc = 0d0
     delta_ij_s2_mrcc = 0d0
  
positive.code:       n = n - 1
        if(n /= 0) then
          call mrcc_part_dress(delta_ij_mrcc_sto, delta_ii_mrcc_sto, delta_ij_s2_mrcc_sto, &
              delta_ii_s2_mrcc_sto, gen,n,buf,N_int,omask,myCoef)
        endif
      end do
      deallocate(buf)
  
negative.code:   integer, external :: omp_get_thread_num
    double precision :: coefs(N_det_non_ref), myCoef
    integer :: n_in_teeth
    double precision :: contrib(N_states), curn, in_teeth_step, curlim, curnorm
    
    contrib = 0d0
    read(*,*) n_in_teeth
    !n_in_teeth = 2
    in_teeth_step = 1d0 / dfloat(n_in_teeth)
  