C
C There are a total of 4 entries in the algebraic variable array.
C There are a total of 20 entries in each of the rate and state variable arrays.
C There are a total of 26 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is S_0 in component S_0 (dimensionless).
C CONSTS(16) is R_0 in component R (first_order_rate_constant).
C CONSTS(17) is u_s in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(1) is phi in component phi (first_order_rate_constant).
C STATES(2) is S_1 in component S_1 (dimensionless).
C CONSTS(18) is R_1 in component R (first_order_rate_constant).
C CONSTS(1) is epsilon_s in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is alpha in component kinetic_parameters (dimensionless).
C CONSTS(3) is u in component kinetic_parameters (first_order_rate_constant).
C STATES(3) is S_2 in component S_2 (dimensionless).
C CONSTS(19) is R_2 in component R (first_order_rate_constant).
C STATES(4) is S_3 in component S_3 (dimensionless).
C CONSTS(20) is R_3 in component R (first_order_rate_constant).
C STATES(5) is S_4 in component S_4 (dimensionless).
C CONSTS(21) is R_4 in component R (first_order_rate_constant).
C STATES(6) is S_5 in component S_5 (dimensionless).
C CONSTS(22) is R_5 in component R (first_order_rate_constant).
C STATES(7) is S_6 in component S_6 (dimensionless).
C CONSTS(23) is R_6 in component R (first_order_rate_constant).
C STATES(8) is S_7 in component S_7 (dimensionless).
C CONSTS(24) is R_7 in component R (first_order_rate_constant).
C STATES(9) is S_8 in component S_8 (dimensionless).
C CONSTS(25) is R_8 in component R (first_order_rate_constant).
C STATES(10) is M_0 in component M_0 (dimensionless).
C CONSTS(4) is epsilon_m in component kinetic_parameters (first_order_rate_constant).
C CONSTS(26) is u_m in component kinetic_parameters (first_order_rate_constant).
C STATES(11) is M_1 in component M_1 (dimensionless).
C STATES(12) is M_1_2 in component M_1 (dimensionless).
C STATES(13) is M_2 in component M_2 (dimensionless).
C STATES(14) is M_3 in component M_3 (dimensionless).
C STATES(15) is M_4 in component M_4 (dimensionless).
C STATES(16) is M_5 in component M_5 (dimensionless).
C STATES(17) is M_6 in component M_6 (dimensionless).
C STATES(18) is M_7 in component M_7 (dimensionless).
C STATES(19) is M_8 in component M_8 (dimensionless).
C CONSTS(5) is r_0 in component R (first_order_rate_constant).
C CONSTS(6) is r_1 in component R (first_order_rate_constant).
C CONSTS(7) is r_2 in component R (first_order_rate_constant).
C CONSTS(8) is r_3 in component R (first_order_rate_constant).
C CONSTS(9) is r_4 in component R (first_order_rate_constant).
C CONSTS(10) is r_5 in component R (first_order_rate_constant).
C CONSTS(11) is r_6 in component R (first_order_rate_constant).
C CONSTS(12) is r_7 in component R (first_order_rate_constant).
C CONSTS(13) is r_8 in component R (first_order_rate_constant).
C CONSTS(14) is a in component kinetic_parameters (dimensionless).
C STATES(20) is w in component w (dimensionless).
C ALGBRC(2) is total_cells in component total_cells (dimensionless).
C ALGBRC(3) is stable_total in component total_cells (dimensionless).
C ALGBRC(4) is mutant_total in component total_cells (dimensionless).
C CONSTS(15) is beta in component kinetic_parameters (dimensionless).
C RATES(1) is d/dt S_0 in component S_0 (dimensionless).
C RATES(2) is d/dt S_1 in component S_1 (dimensionless).
C RATES(3) is d/dt S_2 in component S_2 (dimensionless).
C RATES(4) is d/dt S_3 in component S_3 (dimensionless).
C RATES(5) is d/dt S_4 in component S_4 (dimensionless).
C RATES(6) is d/dt S_5 in component S_5 (dimensionless).
C RATES(7) is d/dt S_6 in component S_6 (dimensionless).
C RATES(8) is d/dt S_7 in component S_7 (dimensionless).
C RATES(9) is d/dt S_8 in component S_8 (dimensionless).
C RATES(10) is d/dt M_0 in component M_0 (dimensionless).
C RATES(11) is d/dt M_1 in component M_1 (dimensionless).
C RATES(12) is d/dt M_1_2 in component M_1 (dimensionless).
C RATES(13) is d/dt M_2 in component M_2 (dimensionless).
C RATES(14) is d/dt M_3 in component M_3 (dimensionless).
C RATES(15) is d/dt M_4 in component M_4 (dimensionless).
C RATES(16) is d/dt M_5 in component M_5 (dimensionless).
C RATES(17) is d/dt M_6 in component M_6 (dimensionless).
C RATES(18) is d/dt M_7 in component M_7 (dimensionless).
C RATES(19) is d/dt M_8 in component M_8 (dimensionless).
C RATES(20) is d/dt w in component w (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.5
      STATES(2) = 0
      CONSTS(1) = 0.99
      CONSTS(2) = 0.6
      CONSTS(3) = 0.07
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 0.5
      CONSTS(4) = 0.1
      STATES(11) = 0
      STATES(12) = 0
      STATES(13) = 0
      STATES(14) = 0
      STATES(15) = 0
      STATES(16) = 0
      STATES(17) = 0
      STATES(18) = 0
      STATES(19) = 0
      CONSTS(5) = 0.5
      CONSTS(6) = 0.6
      CONSTS(7) = 0.7
      CONSTS(8) = 0.8
      CONSTS(9) = 0.9
      CONSTS(10) = 1
      CONSTS(11) = 1.1
      CONSTS(12) = 1.2
      CONSTS(13) = 1.3
      CONSTS(14) = 0.5
      STATES(20) = 0
      CONSTS(15) = 0.2
      CONSTS(16) = CONSTS(5)
      CONSTS(17) =  CONSTS(3)*(1.00000 - ( CONSTS(15)*CONSTS(1))/1.00000)
      CONSTS(18) =  CONSTS(6)*(1.00000 - CONSTS(14))
      CONSTS(19) =  CONSTS(7)*(1.00000 - CONSTS(14))
      CONSTS(20) =  CONSTS(8)*(1.00000 - CONSTS(14))
      CONSTS(21) =  CONSTS(9)*(1.00000 - CONSTS(14))
      CONSTS(22) =  CONSTS(10)*(1.00000 - CONSTS(14))
      CONSTS(23) =  CONSTS(11)*(1.00000 - CONSTS(14))
      CONSTS(24) =  CONSTS(12)*(1.00000 - CONSTS(14))
      CONSTS(25) =  CONSTS(13)*(1.00000 - CONSTS(14))
      CONSTS(26) =  CONSTS(3)*(1.00000 - ( CONSTS(15)*CONSTS(4))/1.00000)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(11) =  STATES(10)*1.00000
      ALGBRC(1) =  ((1.00000 - CONSTS(17))/1.00000)*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5))+ ((1.00000 - CONSTS(26))/1.00000)*( CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15))
      RATES(1) =  (( CONSTS(16)*STATES(1))/1.00000)*(1.00000 - CONSTS(17)) -  ALGBRC(1)*STATES(1)
      RATES(2) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(16)*STATES(1))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(18)*STATES(2))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(2)
      RATES(3) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(18)*STATES(2))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(19)*STATES(3))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(3)
      RATES(4) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(19)*STATES(3))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(20)*STATES(4))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(4)
      RATES(5) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(20)*STATES(4))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(21)*STATES(5))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(5)
      RATES(6) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(21)*STATES(5))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(22)*STATES(6))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(6)
      RATES(7) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(22)*STATES(6))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(23)*STATES(7))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(7)
      RATES(8) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(23)*STATES(7))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(24)*STATES(8))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(8)
      RATES(9) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(24)*STATES(8))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(25)*STATES(9))/1.00000)*(1.00000 - CONSTS(17))) -  ALGBRC(1)*STATES(9)
      RATES(10) =  (( CONSTS(16)*STATES(10))/1.00000)*(1.00000 - CONSTS(26)) -  ALGBRC(1)*STATES(10)
      RATES(12) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(16)*STATES(10))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(18)*STATES(11))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(11)
      RATES(13) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(18)*STATES(11))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(19)*STATES(13))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(13)
      RATES(14) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(19)*STATES(13))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(20)*STATES(14))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(14)
      RATES(15) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(20)*STATES(15))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(21)*STATES(15))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(15)
      RATES(16) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(21)*STATES(15))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(22)*STATES(16))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(16)
      RATES(17) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(22)*STATES(16))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(23)*STATES(17))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(17)
      RATES(18) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(23)*STATES(17))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(24)*STATES(18))/1.00000)*(1.00000 - CONSTS(26))) -  ALGBRC(1)*STATES(18)
      RATES(19) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(24)*STATES(18))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(25)*STATES(19))/1.00000)*((1.00000 - CONSTS(26))+ (( CONSTS(2)*CONSTS(3))/1.00000)*(1.00000 - CONSTS(4)))) -  ALGBRC(1)*STATES(9)
      RATES(20) =  (( (1.00000 - CONSTS(2))*CONSTS(3))/1.00000)*( (1.00000 - CONSTS(1))*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5)+ CONSTS(22)*STATES(6)+ CONSTS(23)*STATES(7)+ CONSTS(24)*STATES(8)+ CONSTS(25)*STATES(9))+ (1.00000 - CONSTS(4))*( CONSTS(16)*STATES(10)+ CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15)+ CONSTS(22)*STATES(16)+ CONSTS(23)*STATES(17)+ CONSTS(24)*STATES(18)+ CONSTS(25)*STATES(19))) -  ALGBRC(1)*STATES(20)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  ((1.00000 - CONSTS(17))/1.00000)*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5))+ ((1.00000 - CONSTS(26))/1.00000)*( CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15))
      ALGBRC(2) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)
      ALGBRC(3) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)
      ALGBRC(4) = STATES(10)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)
      RETURN
      END