C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 18 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is k11 in component kinetic_constants (second_order_rate_constant).
C CONSTS(2) is k11_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(3) is k71 in component kinetic_constants (first_order_rate_constant).
C CONSTS(4) is k72 in component kinetic_constants (first_order_rate_constant).
C CONSTS(5) is k12 in component kinetic_constants (second_order_rate_constant).
C CONSTS(6) is k12_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(7) is k51 in component kinetic_constants (first_order_rate_constant).
C CONSTS(8) is k52 in component kinetic_constants (first_order_rate_constant).
C CONSTS(9) is k41 in component kinetic_constants (first_order_rate_constant).
C CONSTS(10) is k31 in component kinetic_constants (first_order_rate_constant).
C CONSTS(11) is k32 in component kinetic_constants (first_order_rate_constant).
C CONSTS(12) is k42 in component kinetic_constants (first_order_rate_constant).
C CONSTS(13) is k22 in component kinetic_constants (first_order_rate_constant).
C CONSTS(14) is k22_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(15) is k21 in component kinetic_constants (first_order_rate_constant).
C CONSTS(16) is k21_ in component kinetic_constants (first_order_rate_constant).
C STATES(1) is PB in component PB (micromolar).
C CONSTS(17) is kreset1 in component kreset1 (first_order_rate_constant).
C STATES(2) is A1 in component A1 (micromolar).
C STATES(3) is A2 in component A2 (micromolar).
C STATES(4) is PC1 in component PC1 (micromolar).
C STATES(5) is PF1 in component PF1 (micromolar).
C STATES(6) is PG1 in component PG1 (micromolar).
C STATES(7) is PC2 in component PC2 (micromolar).
C STATES(8) is PF2 in component PF2 (micromolar).
C STATES(9) is PD1 in component PD1 (micromolar).
C STATES(10) is PE1 in component PE1 (micromolar).
C STATES(11) is PD2 in component PD2 (micromolar).
C CONSTS(18) is PE2 in component PE2 (micromolar).
C RATES(1) is d/dt PB in component PB (micromolar).
C RATES(2) is d/dt A1 in component A1 (micromolar).
C RATES(3) is d/dt A2 in component A2 (micromolar).
C RATES(4) is d/dt PC1 in component PC1 (micromolar).
C RATES(9) is d/dt PD1 in component PD1 (micromolar).
C RATES(10) is d/dt PE1 in component PE1 (micromolar).
C RATES(5) is d/dt PF1 in component PF1 (micromolar).
C RATES(6) is d/dt PG1 in component PG1 (micromolar).
C RATES(7) is d/dt PC2 in component PC2 (micromolar).
C RATES(11) is d/dt PD2 in component PD2 (micromolar).
C RATES(8) is d/dt PF2 in component PF2 (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 110.0
      CONSTS(2) = 25.0
      CONSTS(3) = 0.1
      CONSTS(4) = 60.0
      CONSTS(5) = 110.0
      CONSTS(6) = 0.2
      CONSTS(7) = 7.0
      CONSTS(8) = 0.0
      CONSTS(9) = 60.0
      CONSTS(10) = 500.0
      CONSTS(11) = 50.0
      CONSTS(12) = 70.0
      CONSTS(13) = 100.0
      CONSTS(14) = 17.0
      CONSTS(15) = 100.0
      CONSTS(16) = 0.2
      STATES(1) = 0.1
      CONSTS(17) = 35.0
      STATES(2) = 0.1
      STATES(3) = 0.1
      STATES(4) = 0.1
      STATES(5) = 0.1
      STATES(6) = 0.1
      STATES(7) = 0.1
      STATES(8) = 0.1
      STATES(9) = 0.1
      STATES(10) = 0.1
      STATES(11) = 0.1
      CONSTS(18) = 0.1
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(2)*STATES(4)+ CONSTS(3)*STATES(5)+ CONSTS(17)*STATES(6)+ CONSTS(6)*STATES(7)+ CONSTS(4)*STATES(8)) - ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(5)*STATES(3)*STATES(1))
      RATES(2) =  CONSTS(2)*STATES(4) -  CONSTS(1)*STATES(2)*STATES(1)
      RATES(3) =  CONSTS(6)*STATES(7) -  CONSTS(5)*STATES(3)*STATES(1)
      RATES(4) = ( CONSTS(1)*STATES(2)*STATES(1) -  (CONSTS(2)+CONSTS(15))*STATES(4))+ CONSTS(16)*STATES(9)
      RATES(9) =  CONSTS(15)*STATES(4) -  (CONSTS(16)+CONSTS(10))*STATES(9)
      RATES(10) =  CONSTS(10)*STATES(9) -  CONSTS(9)*STATES(10)
      RATES(5) =  CONSTS(9)*STATES(10) -  (CONSTS(7)+CONSTS(3))*STATES(5)
      RATES(6) =  CONSTS(7)*STATES(5) -  CONSTS(17)*STATES(6)
      RATES(7) = ( CONSTS(5)*STATES(3)*STATES(1) -  (CONSTS(6)+CONSTS(13))*STATES(7))+ CONSTS(14)*STATES(11)
      RATES(11) =  CONSTS(13)*STATES(7) -  (CONSTS(14)+CONSTS(11))*STATES(11)
      RATES(8) =  CONSTS(12)*CONSTS(18) -  (CONSTS(8)+CONSTS(4))*STATES(8)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END