C
C There are a total of 4 entries in the algebraic variable array.
C There are a total of 6 entries in each of the rate and state variable arrays.
C There are a total of 28 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is H_1 in component parameters (dimensionless).
C CONSTS(2) is H_2 in component parameters (dimensionless).
C CONSTS(3) is H_3 in component parameters (dimensionless).
C CONSTS(4) is H_4 in component parameters (dimensionless).
C CONSTS(5) is K_1 in component parameters (dimensionless).
C CONSTS(6) is K_2 in component parameters (dimensionless).
C CONSTS(7) is K_3 in component parameters (dimensionless).
C CONSTS(8) is K_4 in component parameters (dimensionless).
C CONSTS(9) is V_M1 in component parameters (first_order_rate_constant).
C CONSTS(10) is V_M3 in component parameters (first_order_rate_constant).
C CONSTS(11) is U_M1 in component parameters (first_order_rate_constant).
C CONSTS(12) is U_M3 in component parameters (first_order_rate_constant).
C CONSTS(13) is V_2 in component parameters (first_order_rate_constant).
C CONSTS(14) is V_4 in component parameters (first_order_rate_constant).
C CONSTS(15) is U_2 in component parameters (first_order_rate_constant).
C CONSTS(16) is U_4 in component parameters (first_order_rate_constant).
C CONSTS(17) is K_c1 in component parameters (micromolar).
C CONSTS(18) is K_c2 in component parameters (micromolar).
C CONSTS(19) is K_d1 in component parameters (micromolar).
C CONSTS(20) is K_d2 in component parameters (micromolar).
C CONSTS(21) is v_d1 in component parameters (flux).
C CONSTS(22) is v_d2 in component parameters (flux).
C CONSTS(23) is v_i1 in component parameters (flux).
C CONSTS(24) is v_i2 in component parameters (flux).
C CONSTS(25) is k_d1 in component parameters (first_order_rate_constant).
C CONSTS(26) is k_d2 in component parameters (first_order_rate_constant).
C CONSTS(27) is K_im1 in component parameters (dimensionless).
C CONSTS(28) is K_im2 in component parameters (dimensionless).
C STATES(1) is C_1 in component C_1 (micromolar).
C STATES(2) is M_2 in component M_2 (dimensionless).
C STATES(3) is X_1 in component X_1 (dimensionless).
C STATES(4) is M_1 in component M_1 (dimensionless).
C ALGBRC(1) is V_1 in component V_1 (first_order_rate_constant).
C ALGBRC(2) is V_3 in component V_3 (first_order_rate_constant).
C STATES(5) is C_2 in component C_2 (micromolar).
C STATES(6) is X_2 in component X_2 (dimensionless).
C ALGBRC(3) is U_1 in component U_1 (first_order_rate_constant).
C ALGBRC(4) is U_3 in component U_3 (first_order_rate_constant).
C RATES(1) is d/dt C_1 in component C_1 (micromolar).
C RATES(4) is d/dt M_1 in component M_1 (dimensionless).
C RATES(3) is d/dt X_1 in component X_1 (dimensionless).
C RATES(5) is d/dt C_2 in component C_2 (micromolar).
C RATES(2) is d/dt M_2 in component M_2 (dimensionless).
C RATES(6) is d/dt X_2 in component X_2 (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.01
      CONSTS(2) = 0.01
      CONSTS(3) = 0.01
      CONSTS(4) = 0.01
      CONSTS(5) = 0.01
      CONSTS(6) = 0.01
      CONSTS(7) = 0.01
      CONSTS(8) = 0.01
      CONSTS(9) = 0.3
      CONSTS(10) = 0.1
      CONSTS(11) = 0.3
      CONSTS(12) = 0.1
      CONSTS(13) = 0.15
      CONSTS(14) = 0.05
      CONSTS(15) = 0.15
      CONSTS(16) = 0.05
      CONSTS(17) = 0.5
      CONSTS(18) = 0.5
      CONSTS(19) = 0.02
      CONSTS(20) = 0.02
      CONSTS(21) = 0.025
      CONSTS(22) = 0.025
      CONSTS(23) = 0.05
      CONSTS(24) = 0.05
      CONSTS(25) = 0.001
      CONSTS(26) = 0.001
      CONSTS(27) = 0.03
      CONSTS(28) = 0.03
      STATES(1) = 2
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 1
      STATES(5) = 0
      STATES(6) = 0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (( CONSTS(23)*CONSTS(27))/(CONSTS(27)+STATES(2)) - ( CONSTS(21)*STATES(3)*STATES(1))/(CONSTS(19)+STATES(1))) -  CONSTS(25)*STATES(1)
      RATES(5) = (( CONSTS(24)*CONSTS(28))/(CONSTS(28)+STATES(4)) - ( CONSTS(22)*STATES(6)*STATES(5))/(CONSTS(20)+STATES(5))) -  CONSTS(26)*STATES(5)
      ALGBRC(1) =  (STATES(1)/(CONSTS(17)+STATES(1)))*CONSTS(9)
      RATES(4) = ( ALGBRC(1)*(1.00000 - STATES(4)))/(CONSTS(5)+(1.00000 - STATES(4))) - ( CONSTS(13)*STATES(4))/(CONSTS(6)+STATES(4))
      ALGBRC(2) =  STATES(4)*CONSTS(10)
      RATES(3) = ( ALGBRC(2)*(1.00000 - STATES(3)))/(CONSTS(7)+(1.00000 - STATES(3))) - ( CONSTS(14)*STATES(3))/(CONSTS(8)+STATES(3))
      ALGBRC(3) =  (STATES(5)/(CONSTS(18)+STATES(5)))*CONSTS(11)
      RATES(2) = ( ALGBRC(3)*(1.00000 - STATES(2)))/(CONSTS(1)+(1.00000 - STATES(2))) - ( CONSTS(15)*STATES(2))/(CONSTS(2)+STATES(2))
      ALGBRC(4) =  STATES(2)*CONSTS(12)
      RATES(6) = ( ALGBRC(4)*(1.00000 - STATES(6)))/(CONSTS(3)+(1.00000 - STATES(6))) - ( CONSTS(16)*STATES(6))/(CONSTS(4)+STATES(6))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  (STATES(1)/(CONSTS(17)+STATES(1)))*CONSTS(9)
      ALGBRC(2) =  STATES(4)*CONSTS(10)
      ALGBRC(3) =  (STATES(5)/(CONSTS(18)+STATES(5)))*CONSTS(11)
      ALGBRC(4) =  STATES(2)*CONSTS(12)
      RETURN
      END