C
C There are a total of 19 entries in the algebraic variable array.
C There are a total of 9 entries in each of the rate and state variable arrays.
C There are a total of 30 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is NO in component NO (micromolar).
C STATES(2) is O_2m in component O_2m (micromolar).
C CONSTS(1) is O_2 in component NO (micromolar).
C STATES(3) is NO_2 in component NO (micromolar).
C STATES(4) is N2O3 in component N2O3 (micromolar).
C STATES(5) is GSNO in component GSNO (micromolar).
C STATES(6) is CcOX in component NO (micromolar).
C STATES(7) is FeLn in component FeLn (micromolar).
C CONSTS(29) is r_1NO in component NO (flux).
C ALGBRC(1) is r_4NO in component NO (flux).
C ALGBRC(2) is r_12aNO in component NO (flux).
C ALGBRC(3) is r_12bNOp in component NO (flux).
C ALGBRC(4) is r_12bNOm in component NO (flux).
C ALGBRC(5) is r_14NO in component NO (flux).
C ALGBRC(6) is r_15NO in component NO (flux).
C ALGBRC(7) is r_16NO in component NO (flux).
C CONSTS(2) is k_1NO in component model_constant (flux).
C CONSTS(3) is k_4NO in component model_constant (second_order_rate_constant).
C CONSTS(4) is k_12aNO in component model_constant (rate2).
C CONSTS(5) is k_12bNOp in component model_constant (second_order_rate_constant).
C CONSTS(6) is k_12bNOm in component model_constant (first_order_rate_constant).
C CONSTS(7) is k_14NO in component model_constant (first_order_rate_constant).
C CONSTS(8) is k_15NO in component model_constant (second_order_rate_constant).
C CONSTS(9) is k_16NO in component model_constant (second_order_rate_constant).
C CONSTS(10) is SOD in component O_2m (micromolar).
C CONSTS(30) is r_2NO in component O_2m (flux).
C ALGBRC(8) is r_5NO in component O_2m (flux).
C ALGBRC(9) is r_10NO in component O_2m (flux).
C CONSTS(11) is k_2NO in component model_constant (flux).
C CONSTS(12) is k_5NO in component model_constant (second_order_rate_constant).
C CONSTS(13) is k_10NO in component model_constant (rate2).
C STATES(8) is ONOO_m in component ONOO_m (micromolar).
C STATES(9) is GSH in component GSH (micromolar).
C CONSTS(14) is GPX in component ONOO_m (micromolar).
C CONSTS(15) is CO_2 in component ONOO_m (micromolar).
C CONSTS(16) is Cyt_c in component ONOO_m (micromolar).
C ALGBRC(10) is r_6NO in component ONOO_m (flux).
C ALGBRC(11) is r_7NO in component ONOO_m (flux).
C ALGBRC(13) is r_8NO in component ONOO_m (flux).
C ALGBRC(15) is r_9NO in component ONOO_m (flux).
C CONSTS(17) is k_6NO in component model_constant (second_order_rate_constant).
C CONSTS(18) is k_7NO in component model_constant (second_order_rate_constant).
C CONSTS(19) is k_8NO in component model_constant (second_order_rate_constant).
C CONSTS(20) is k_9NO in component model_constant (second_order_rate_constant).
C ALGBRC(14) is GSSG in component GSH (micromolar).
C ALGBRC(12) is FeLnNO in component GSH (micromolar).
C CONSTS(21) is FeLn_0 in component model_constant (micromolar).
C CONSTS(22) is GSH_0 in component model_constant (micromolar).
C ALGBRC(16) is r_11NO in component GSH (flux).
C ALGBRC(17) is r_m in component GSH (flux).
C ALGBRC(19) is r_17NO in component GSH (flux).
C CONSTS(23) is k_11NO in component model_constant (second_order_rate_constant).
C CONSTS(24) is v_m in component model_constant (flux).
C CONSTS(25) is k_m in component model_constant (micromolar).
C CONSTS(26) is k_17NO in component model_constant (second_order_rate_constant).
C ALGBRC(18) is r_13NO in component N2O3 (flux).
C CONSTS(27) is k_13NO in component model_constant (first_order_rate_constant).
C CONSTS(28) is k_17bNO in component model_constant (second_order_rate_constant).
C RATES(1) is d/dt NO in component NO (micromolar).
C RATES(6) is d/dt CcOX in component NO (micromolar).
C RATES(3) is d/dt NO_2 in component NO (micromolar).
C RATES(2) is d/dt O_2m in component O_2m (micromolar).
C RATES(8) is d/dt ONOO_m in component ONOO_m (micromolar).
C RATES(9) is d/dt GSH in component GSH (micromolar).
C RATES(5) is d/dt GSNO in component GSNO (micromolar).
C RATES(4) is d/dt N2O3 in component N2O3 (micromolar).
C RATES(7) is d/dt FeLn in component FeLn (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0
      STATES(2) = 0
      CONSTS(1) = 35
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0.1
      STATES(7) = 0.05
      CONSTS(2) = 1
      CONSTS(3) = 6700
      CONSTS(4) = 0.000006
      CONSTS(5) = 1100
      CONSTS(6) = 81000
      CONSTS(7) = 0.0002
      CONSTS(8) = 100
      CONSTS(9) = 1.21
      CONSTS(10) = 10
      CONSTS(11) = 0.1
      CONSTS(12) = 2400
      CONSTS(13) = 0.0006
      STATES(8) = 0
      STATES(9) = 10000
      CONSTS(14) = 5.8
      CONSTS(15) = 1000
      CONSTS(16) = 400
      CONSTS(17) = 0.00135
      CONSTS(18) = 2
      CONSTS(19) = 0.058
      CONSTS(20) = 0.025
      CONSTS(21) = 0.05
      CONSTS(22) = 10000
      CONSTS(23) = 66
      CONSTS(24) = 320
      CONSTS(25) = 50
      CONSTS(26) = 66
      CONSTS(27) = 1600
      CONSTS(28) = 0.0002
      CONSTS(29) = CONSTS(2)
      CONSTS(30) = CONSTS(11)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) =  CONSTS(4)*STATES(1) ** 2.00000*CONSTS(1)
      ALGBRC(3) =  CONSTS(5)*STATES(3)*STATES(1)
      ALGBRC(4) =  CONSTS(6)*STATES(4)
      RATES(3) = ( 2.00000*ALGBRC(2) - ALGBRC(3))+ALGBRC(4)
      ALGBRC(6) =  CONSTS(8)*STATES(6)*STATES(1)
      RATES(6) = - ALGBRC(6)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(2)
      ALGBRC(5) =  CONSTS(7)*STATES(5)
      ALGBRC(7) =  CONSTS(9)*STATES(7)*STATES(1)
      RATES(1) = (((((CONSTS(29) - ALGBRC(1)) -  2.00000*ALGBRC(2)) - ALGBRC(3))+ALGBRC(4)+ALGBRC(5)) - ALGBRC(6)) - ALGBRC(7)
      ALGBRC(8) =  CONSTS(12)*CONSTS(10)*STATES(2)
      ALGBRC(9) =  CONSTS(13)*STATES(5) ** 2.00000*STATES(2)
      RATES(2) = ((CONSTS(30) - ALGBRC(1)) - ALGBRC(8)) - ALGBRC(9)
      ALGBRC(10) =  CONSTS(17)*STATES(8)*STATES(9)
      ALGBRC(11) =  CONSTS(18)*STATES(8)*CONSTS(14)
      ALGBRC(13) =  CONSTS(19)*STATES(8)*CONSTS(15)
      ALGBRC(15) =  CONSTS(20)*STATES(8)*CONSTS(16)
      RATES(8) = (((ALGBRC(1) - ALGBRC(10)) - ALGBRC(11)) - ALGBRC(13)) - ALGBRC(15)
      ALGBRC(16) =  CONSTS(23)*STATES(4)*STATES(9)
      ALGBRC(18) =  CONSTS(27)*STATES(4)
      RATES(4) = ((- ALGBRC(16)+ALGBRC(3)) - ALGBRC(4)) - ALGBRC(18)
      ALGBRC(14) = ((CONSTS(22) - STATES(9)) - STATES(5))/2.00000
      ALGBRC(17) = ( CONSTS(24)*ALGBRC(14))/(CONSTS(25)+ALGBRC(14))
      ALGBRC(12) = CONSTS(21) - STATES(7)
      ALGBRC(19) =  CONSTS(26)*ALGBRC(12)*STATES(9)
      RATES(9) = ((- ALGBRC(10) - ALGBRC(16))+ 2.00000*ALGBRC(17)) - ALGBRC(19)
      RATES(5) = (((ALGBRC(10) -  2.00000*ALGBRC(9))+ALGBRC(16)) - ALGBRC(5))+ALGBRC(19)
      RATES(7) = - ALGBRC(7)+ALGBRC(19)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) =  CONSTS(4)*STATES(1) ** 2.00000*CONSTS(1)
      ALGBRC(3) =  CONSTS(5)*STATES(3)*STATES(1)
      ALGBRC(4) =  CONSTS(6)*STATES(4)
      ALGBRC(6) =  CONSTS(8)*STATES(6)*STATES(1)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(2)
      ALGBRC(5) =  CONSTS(7)*STATES(5)
      ALGBRC(7) =  CONSTS(9)*STATES(7)*STATES(1)
      ALGBRC(8) =  CONSTS(12)*CONSTS(10)*STATES(2)
      ALGBRC(9) =  CONSTS(13)*STATES(5) ** 2.00000*STATES(2)
      ALGBRC(10) =  CONSTS(17)*STATES(8)*STATES(9)
      ALGBRC(11) =  CONSTS(18)*STATES(8)*CONSTS(14)
      ALGBRC(13) =  CONSTS(19)*STATES(8)*CONSTS(15)
      ALGBRC(15) =  CONSTS(20)*STATES(8)*CONSTS(16)
      ALGBRC(16) =  CONSTS(23)*STATES(4)*STATES(9)
      ALGBRC(18) =  CONSTS(27)*STATES(4)
      ALGBRC(14) = ((CONSTS(22) - STATES(9)) - STATES(5))/2.00000
      ALGBRC(17) = ( CONSTS(24)*ALGBRC(14))/(CONSTS(25)+ALGBRC(14))
      ALGBRC(12) = CONSTS(21) - STATES(7)
      ALGBRC(19) =  CONSTS(26)*ALGBRC(12)*STATES(9)
      RETURN
      END