anchor.diff:    integer, allocatable :: iorder(:)
     double precision , allocatable :: norm_sort(:)
     double precision :: e_0(N_states)
  -  if (.not.read_wf) then
  -    stop 'Please set read_wf to true'
  -  endif
   
     PROVIDE mo_bielec_integrals_in_map H_apply_buffer_allocated
   
  
positive.code: program e_curve
    use bitmasks
    implicit none
    integer :: i,j,k, kk, nab, m, l
  
negative.code:   double precision, allocatable  :: u_0(:,:), v_0(:,:)
    allocate(u_t(N_states,N_det),v_t(N_states,N_det),s_t(N_states,N_det))
    allocate(u_0(N_det,N_states),v_0(N_det,N_states))
  
    norm_sort(0) = 0.d0
    iorder(0) = 0
  