anchor.diff: -program e_curve
  +program truncate
  +  read_wf = .True.
  +  SOFT_TOUCH read_wf
  +  call run
  +end
  +
  +subroutine run
     use bitmasks
     implicit none
     integer :: i,j,k, kk, nab, m, l
  
positive.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))
  
    read(*,*) ci_threshold
  
    norm_sort(0) = 0.d0
    iorder(0) = 0
  
negative.code:   integer, allocatable :: iorder(:)
    double precision , allocatable :: norm_sort(:)
    double precision :: e_0(N_states)
  
    PROVIDE mo_bielec_integrals_in_map H_apply_buffer_allocated
  
  