rm and Gaussian styles the momentum of the initialized atoms is
  also zeroed, but only if all atoms are being initialized 
for uniform and Gaussian styles, RN are generated with Park/Miller RNG
for velocity style in 2-d simulations, still specify z velocity component,
  even though it is ignored
</PRE>
<HR>
<H3>
<A NAME="_cch3_931276727">diagnostic</A></H3>
<UL>
    <LI>
    1st parameter = nametag of a user routine added to diagnostic.f file 
    <LI>
    2nd parameter = call this user routine every this 