SUBROUTINE DEEPREC(ENERGY,EECRT,IOUT)
  use sebomd_module, only : sebomd_obj
  implicit none
! compute machine precision for current energy
  double precision a, b, c, ENERGY, EECRT
  integer iout, iprec
  iprec = sebomd_obj%iprec
  a = ENERGY
  b = 0.1d0
  c = ENERGY
  do while ( (a+b).ne.(c) )
    b = b*0.1d0
  end do
! write(iout,'("DEEPREC",d20.10)') b
  EECRT = b*(10.0d0**iprec)
! write(0,'(3f30.20,i5)') energy, b, eecrt, iprec
  RETURN
END