!---- shock parameters ---------------------------------------------------------
C                       ! shock type : 'C' or 'J', Steady state : 'S', PDR: 'P' or 'p' (without GRAIN reactions)
3                       ! Nfluids : 1, 2 ou 3
1.0                       ! Bbeta -> Bfield = Bbeta * sqrt(nH)
UUUUUU  	                ! Vs -> shock speed (km/s)
1.0e3                   ! Vn - Vi initial (cm s-1)
OOOOOO                     ! op_H2 -> initial H2 ortho/para ratio (999.9 -> ETL)
10.0                    ! T(n,i,e) -> initial gas temperature (K)
1e3                    ! nH_init -> initial value for n(H) + 2.0 n(H2) + n(H+) (cm-3)
15                      ! Tgrains -> initial grain temperature (K)
1                       ! Cool_KN -> 1: Kaufman & Neufeld cooling
!---- environment --------------------------------------------------------------
3.0D-17                 ! Zeta -> cosmic ray ionization rate (s-1)
1                     ! RAD  -> flux radiation (multiplicative factor)
0.1        !0.133352143216332                    ! Av0  -> initial extinction (magnitudes)
5.34D-22      !-5.34D-22               ! inv_Av_conv_factor   -> E(B-V)/R (if zero, Av is not integrated) 
1e20    ! 9.36e+20			! N_H2_0 -> column-density of H2 at surface
0e17 			! N_CO_0 -> column-density of CO at surface
N			! Av_of_NH2 -> whether we guess Av uniquely from NH2 (Y/N) otherwise, use Av(NH+2NH2)
!---- numerical parameters -----------------------------------------------------
100000                  ! Nstep_max   -> max number of integration steps
5                       ! Nstep_w -> number of steps between 2 outputs
149                      ! NH2_lev -> Number of H2 levels included
10                       ! NH2_lines_out -> Max number of H2 lines in output file
BOTH                     ! H_H2_flag -> H-H2 collisions : DRF, MM or BOTH
1                       ! iforH2 -> Formation on grain model (1, 2, 3, 4)
2                       ! ikinH2 -> Kinetic energy of H2 newly formed (1, 2)
0.333333D15                  ! XLL/nH -> caracteristic viscous length (cm)
1.00D-7                 ! Eps_V -> precision of computation
1.00D7                 ! timeJ -> shock age (years)
1.00D7                 ! duration_max -> max. shock duration (years)
1                       ! Force_I_C -> 1: Force Ion Conservation
!---- output specifications ----------------------------------------------------
FD                      ! species: 'AD' (cm-3), 'CD' (cm-2) or 'FD' (n(x)/nH)
CD                      ! H2 levels: 'AD' (cm-3), 'CD' (cm-2) or 'ln(N/g)'
integrated              ! H2 lines: 'local' (erg/s/cm3) or 'integrated' (erg/s/cm2/sr)
N			! flag_analysis -> Output chemical analysis (best reactions) (Y/N) 
!-------------------------------------------------------------------------------

 INTEGER            :: iforH2 = 1  ! Flag : H2 formation on grains
                                   !   0: 1/3 of 4.4781 eV in internal energy (=> 17249 K) (Allen, 1999)
                                   !   1: Proportional to Boltzman Distrib at 17249 K
                                   !   2: Dissociation limit : v = 14, J = 0,1 (4.4781 eV)
                                   !   3: v = 6, J = 0,1
                                   !   4: fraction = relative populations at t, initialised as H2_lev%density
                                   !                 and changed during integration
 INTEGER            :: ikinH2 = 2  ! Flag : H2 formation energy released as kinetic energy
                                   !   1: 0.5 * (4.4781 - internal)
                                   !   2: Inf(1.4927 eV, 4.4781 - internal)