! units and constants
       double precision tempAr,sigmaAr,timeAr
       double precision u2vel,u2velsq,u2len,u2mass,u2partmass
       double precision u2temp, u2E,u2F,uvsq2temp
       double precision massAr, eps, boltz, iboltz, avogadro
       double precision kgparticle

       parameter (pi=3.14159D0)
   
       parameter (boltz = 1.38065D-23)  ! boltzman constant 1.3806503 x 10-23 m^2 kg s^-2 K^-1 
       parameter (iboltz = 7.242965D22) ! inverse kb 
       parameter (eps = 1.654D-21)  ! joules
       parameter (massAr = 39.948)      ! amu per mole
       parameter (tempAr   = 119.8)
       parameter (sigmaAr  = 3.405D-10)
       parameter (timeAr   = 2.156D-12)
       parameter (avogadro   = 6.022D23)

       parameter (u2vel   = sigmaAr / timeAr)
       parameter (u2velsq = u2vel ** 2)
       parameter (u2len   = sigmaAr)
       parameter (u2temp  = tempAr)
       parameter (u2mass  = 0.03994)   ! kg / mole
       parameter (u2partmass  = 6.632D-26)
       parameter (u2E = eps)
       parameter (u2F = eps / sigmaAr)
       parameter (uvsq2temp = u2velsq * u2partmass * iboltz)
       parameter ( kgparticle = u2mass / avogadro )