      subroutine lhmono(ifil)
      use module
      !
      ! lecture des integrales monoelectoniques
      ! stockage dans une matice carree
      !
      IMPLICIT REAL(8)(A-H,O-X,Z),logical(1)(y)
      double precision,allocatable,dimension (:) :: ajk
      allocatable hmono1(:,:)
      character(80) aa
!      yw=.false.
!      hmono=0.d0
!      if(ycl) then
!         allocate (hmono1(norb,norb))
!         hmono1=0.d0
!         rewind 26
!         do
!            read(26,*) i,j,k,l,a
!            if(i==0) exit
!            if(k==0) then
!               hmono1(i,j)=a
!               hmono1(j,i)=a
!            endif
!         enddo
!         if(yw)print*,'hmono1 dans lhmono'
!         if(yw)call script(hmono1,norb,norb,norb,norb)
!         do is=1,nsym
!            ii=sum(isym(1:is-1))
!            do i=1,norb
!               if(itsym(i)/=is) cycle
!               ii=ii+1
!               jj=sum(isym(1:is-1))
!               do j=1,norb
!                  if(itsym(j)/=is) cycle
!                  jj=jj+1
!      !            print*,ii,jj,i,j
!                  hmono(ii,jj)=hmono1(i,j)
!               enddo
!            enddo
!         enddo
!         deallocate (hmono1)
!      else
!         rewind ifil
!         read(ifil) aa
!         ii=0
!         do while (ii==0)
!            read(ifil) aa
!            read(ifil)
!            read(ifil)
!            ii=index(aa,'ONEINT MATRIX (SYM)')
!         enddo
!         backspace (ifil)
!         backspace (ifil)
!         read(ifil) nn,nsym
!         if(yw)print*,'nn,nsym',nn,nsym
!         allocate (ajk(nn))
!         call lir8(ajk,nn,ifil)
!         ij=0
!         do is=1,nsym
!            do i=1,isym(is)
!               i1=i+sum(isym(1:is-1))
!               do j=1,i
!                  j1=j+sum(isym(1:is-1))
!                  ij=ij+1
!                  hmono(i1,j1)=ajk(ij)
!                  hmono(j1,i1)=ajk(ij)
!               enddo
!            enddo
!         enddo
!         deallocate (ajk)
!      endif
!      ! Phase des orbitales
!      do k=1,norb
!         if(iphas(ngel+k)==-1) then
!            hmono(:,k)=-hmono(:,k)
!            hmono(k,:)=-hmono(k,:)
!         endif
!      enddo
!      if(yw) then
!         print*,'hmono dans lhmono (dans l''ordre des symetries'
!         call script(hmono,norb,norb,norb,norb)
!      endif
      return
      end
