5/13/2012 by Adam

1. Added options to run simulations with NMR-type distance restraints;
introduced array dhpb1 to store right borders of distance restraints
(left borders stored in dhpb, if dhpb1(i).eq.0 ith restraint is harmonic type,
> 0, quartic (NMR-like) type).

2. Added option to imposed either Calpha or SC (Cbeta) restraints; the pertinent
array is ibecacarb.

3. The above changes have been introduced in readrtns.F and energy_p_new_barrier.F
Cleaned excessive output in MREMD.F (energies at exchange)

4. prng.f replaces prng_32.F in Intel Makefiles; program hangs up with prng_32