diff --git a/code/m_parameters.f90 b/code/m_parameters.f90
index 5f17f29..a096fa2 100644
--- a/code/m_parameters.f90
+++ b/code/m_parameters.f90
@@ -11,6 +11,7 @@ module m_parameters
 
      ! Transport Properties
      real :: vis,sc,diff
+     real :: d_turb
      real :: scp,prp,lep,vis0p,rod
      real :: lewis, le_a, le_x
 
@@ -86,33 +87,41 @@ contains
        END IF
 
        OPEN(otape,FILE='otape')
-        CALL READ_PARAMETER (nx)
-        CALL READ_PARAMETER (l_0)
-        CALL READ_PARAMETER (int_pr)
-        CALL READ_PARAMETER (tar_lo)
-        CALL READ_PARAMETER (dt)
-        CALL READ_PARAMETER (sc)
-        CALL READ_PARAMETER (vis)
-        CALL READ_PARAMETER (pre)
-        CALL READ_PARAMETER (ac)
-        CALL READ_PARAMETER (bc)
-        CALL READ_PARAMETER (u0)
-        CALL READ_PARAMETER (tf)
-        CALL READ_PARAMETER (dt_uf)
-        CALL READ_PARAMETER (ctmp)
-        CALL READ_PARAMETER (c_cut)
-        CALL READ_PARAMETER (c_ref)
-        CALL READ_PARAMETER (min_wr)
-        CALL READ_PARAMETER (prof_wr)
-        CALL READ_PARAMETER (lewis)
-        CALL READ_PARAMETER (lambda1)
-        CALL READ_PARAMETER (lambda2)
-        CALL READ_PARAMETER (beta1)
-        CALL READ_PARAMETER (hrp)
-        CALL READ_PARAMETER (le_a)
-        CALL READ_PARAMETER (le_x)
-        CALL READ_PARAMETER (nsp)
-        CALL READ_PARAMETER (nrxn)
+        CALL READ_PARAMETER (nx)       ! n grid points
+        CALL READ_PARAMETER (l_0)      ! domain length
+        CALL READ_PARAMETER (int_pr)   ! print interval
+        CALL READ_PARAMETER (tar_lo)   ! target flame location
+        CALL READ_PARAMETER (dt)       ! time step size
+        CALL READ_PARAMETER (sc)       ! Schmidt Number
+        CALL READ_PARAMETER (vis)      ! viscosity
+
+        CALL READ_PARAMETER (pre)      ! 1-step chemistry
+        CALL READ_PARAMETER (ac)       ! 1-step chemistry
+        CALL READ_PARAMETER (bc)       ! 1-step chemistry
+
+        CALL READ_PARAMETER (u0)       ! inlet velocity
+        CALL READ_PARAMETER (tf)       ! end time
+        CALL READ_PARAMETER (dt_uf)    ! flame control time step
+        CALL READ_PARAMETER (ctmp)     ! unburned c value
+
+        CALL READ_PARAMETER (c_cut)    ! reaction rate profile parameter
+        CALL READ_PARAMETER (c_ref)    ! reaction rate profile parameter
+        CALL READ_PARAMETER (min_wr)   ! reaction rate profile parameter
+        CALL READ_PARAMETER (prof_wr)  ! reaction rate profile parameter
+
+        CALL READ_PARAMETER (lewis)    ! Lewis number
+
+        CALL READ_PARAMETER (lambda1)  ! 2-step chemistry
+        CALL READ_PARAMETER (lambda2)  ! 2-step chemistry
+        CALL READ_PARAMETER (beta1)    ! 2-step chemistry
+        CALL READ_PARAMETER (hrp)      ! 2-step chemistry
+        CALL READ_PARAMETER (le_a)     ! 2-step chemistry
+        CALL READ_PARAMETER (le_x)     ! 2-step chemistry
+
+        CALL READ_PARAMETER (nsp)      ! number of species
+        CALL READ_PARAMETER (nrxn)     ! number of reactions
+
+        CALL READ_PARAMETER (d_turb)   ! turbulent diffusivity
 
         ! READ(itape,*) cdum,table_size
         ! WRITE(otape,*) cdum,table_size
diff --git a/code/sample/itape.py b/code/sample/itape.py
index 994dd4f..ea3c0bf 100755
--- a/code/sample/itape.py
+++ b/code/sample/itape.py
@@ -39,6 +39,7 @@ le_a      1.
 le_x      0.20
 nsp       2
 nrxn      1
+d_turb    0
 '''
 ).split()