%(x): for note      Energy input: [  Ea   |  H0   ]                        %#ok
 (1): H2 + 2#   <-> 2H#           [ 1.209  -0.757 ]                        %#ok
 (2): N2 + 2#   <-> 2N#           [ 0.658  -1.175 ]                        %#ok
 (3): N# + H#   <-> NH# + #       [ 1.429  -0.083 ]                        %#ok
 (4): NH# + H#  <-> NH2# + #      [ 1.592  -0.158 ]                        %#ok
 (5): NH2# + H# <-> NH3 + 2#      [ 1.986   1.435 ]                        %#ok
CalcDRC = 1;       % calculate the degree of rate control
T = 673;           % reaction temperature 
P_H2_FROZ = 75;    % H2 pressure
P_N2_FROZ = 25;    % N2 pressure
P_NH3_FROZ = 1;    % NH3 pressure, ~1% coversion
Q_v_INIT = 1;      % initial coverage of free site
ThermoMode = 4;    % only include entropy correction                       
BarrierMode = 1;   % deal the adsorption with collision theory            
Mr_NH3 = -2;       % the barrier of the NH3 desorption is used by given                                                                                                                      
npar = 8;          % the process number for parallel computing                                              

