anchor.diff:    Kmat = 0.0d0
   
     call AOMat_to_MO()
  +
  +  Coef = NOCoef
  +
  +  call Fock_to_MO()
   
     Erdmft = 0.0d0
     E1e  = 0.0d0
  
positive.code: !#                 Run 1-RDMFT
  !#############################################
  subroutine run_rdmft()
    use module_data, only      : dim_1e,Hcore,NOEps
    implicit none
    integer                   :: i,j
    double precision          :: E1e,E2ec,E2ex,E2e
  
negative.code:   do i=1,dim_1e
        E1e = E1e + 1.0d0*NOEps(i,1)*Hcore(i,i)
      do j=1,dim_1e
  !      E2ec = E2ec + 1.0d0*NOEps(i,1)*NOEps(j,1)*Jmat(i,j,1)
  !      E2ex = E2ex - 0.5d0*NOEps(i,1)*NOEps(j,1)*Kmat(i,j,1)
        E2e = NOEps(i,1)*NOEps(j,1)*(Fock(i,j,1)-Hcore(i,j))      
      enddo
    enddo
  
  !  E2e = E2ex + E2ec
    Erdmft = E1e + E2e!+ E2ec + E2ex
  
    write(*,*) "         E1e = ", E1e   
    write(*,*) "         E2ex = ", E2ex   
  