added two input parameters nsp and nrxn

code/m_parameters.f90

@@ -31,6 +31,7 @@ module m_parameters
      real :: pflame,pflold,oldu
 
      ! Solver Control
+     integer :: nsp, nrxn
      integer :: ncyc=0,int_pr
      real :: absolute_tolerence=1e-8
      real :: dt,tf,t_now,t_uf,dt_uf
@@ -109,6 +110,8 @@ contains
         CALL READ_PARAMETER (hrp)
         CALL READ_PARAMETER (le_a)
         CALL READ_PARAMETER (le_x)
+        CALL READ_PARAMETER (nsp)
+        CALL READ_PARAMETER (nrxn)
 
         ! READ(itape,*) cdum,table_size
         ! WRITE(otape,*) cdum,table_size

code/sample/itape.py

@@ -37,6 +37,8 @@ beta1     5
 hrp       0.75
 le_a      1.
 le_x      0.20
+nsp       1
+nrxn      1
 '''
 ).split()