title(max 50 char.)= HOD test 
**** GENERAL *************************************************************************************************************
401 500	Number of first and last trajectory to be generated
6.      Initial distance of center of mass from the surface (Ang)
**** VIBRATION ***********************************************************************************************************
.false. 450.0 1.5       Use Boltzmann sampling of vibrational levels,  Nozzle temperature (K) and maximum relative vibrational energy considered (eV) ; 
Below (column by column): [i] Normal Mode is quasi-classical (if false, Q and VQ for that normal mode are set to Qmin and 0) / [ii] if true, number of quanta in vibration:
.TRUE. 1       ! v1    OH stretch 
.TRUE. 0       ! v2    OD stretch 
.TRUE. 0       ! v3    bending 
**** TRANSLATION *********************************************************************************************************
.true.           Use a Velocity distribution for normal translational energy. Parameters for velocitty distribution on the following line.
3026.00 246.00 3.0  Stream Velocity (m/s), velocity width (m/s), Maximum translational energy considered (eV)
.true. 0.025      Shift velocity distribution. If true, by how much (eV)
0.0	          If do not use Velocity distribution, initial normal translational energy of the molecule (eV)
**** ROTATION ************************************************************************************************************
2 -2 0 2   Rotational state of the molecule: (J, M, Ka, Kc) 
.true.     if J=0, sample orientation of molecule isotropically (if .false. next line for euler angles that specify orientation ) 
0.0 0.0 0.0	  Euler angles phi, theta, phi (Deg) [see definition of angles in Rotation_ZYZ function in RotationalState.f90]
**** IMPACT SITE *********************************************************************************************************
.true.  Assign Random impact parameter
0.0	A1 impact parameter (direct coordinate)
0.0	A2 impact parameter (direct coordinate)
7.601166280316	Surface lattice constant (Ang) 
1 10 1 999     Nsurfa Nsurfb Nframea Nframeb - Choose initial conditions of the surface by choosing one random number in [Nsurfa;Nsurfb], which represents the surface number, and one in [Nframea;Nframeb], which represents the frame number.
**************************************************************************************************************************