use io_units
use MD_data, only:D_ban,IP
use geometry_data
+! use prng ! prng.f90 or prng_32.f90
implicit none
!
!-----------------------------------------------------------------------------
integer :: L,JL,JK,J,JM1,K
real(kind=8) :: A(*)
real(kind=8) :: D(L)
+!el real(kind=8) :: E(L)
real(kind=8) :: E(*)!el E(L)
real(kind=8) :: F
real(kind=8) :: G