C
C There are a total of 13 entries in the algebraic variable array.
C There are a total of 3 entries in each of the rate and state variable arrays.
C There are a total of 19 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(17) is F in component main (dimensionless).
C CONSTS(1) is R_T in component main (dimensionless).
C ALGBRC(1) is R_off in component main (dimensionless).
C STATES(1) is D in component main (dimensionless).
C STATES(2) is A_1 in component main (dimensionless).
C STATES(3) is A_2 in component main (dimensionless).
C ALGBRC(5) is lambda_off in component main (dimensionless).
C ALGBRC(6) is lambda_on in component main (dimensionless).
C ALGBRC(2) is lambda_D in component main (dimensionless).
C ALGBRC(3) is lambda_A_1 in component main (dimensionless).
C ALGBRC(7) is lambda_A_2 in component main (dimensionless).
C ALGBRC(4) is lambda_A2_cyc in component main (dimensionless).
C CONSTS(16) is Ca in component main (dimensionless).
C CONSTS(15) is Ca_50 in component main (dimensionless).
C ALGBRC(12) is k_on in component XB_RU_interaction (per_second).
C ALGBRC(13) is k_off in component XB_RU_interaction (per_second).
C CONSTS(2) is k_0_on in component main (per_second).
C CONSTS(3) is k_0_off in component main (per_second).
C CONSTS(4) is k_Ca_on in component main (per_second).
C CONSTS(5) is k_Ca_off in component main (per_second).
C ALGBRC(10) is f in component XB_XB_interaction (per_second).
C ALGBRC(11) is f_prime in component XB_XB_interaction (per_second).
C CONSTS(6) is f_0 in component main (per_second).
C CONSTS(7) is f_prime_0 in component main (per_second).
C CONSTS(8) is h in component main (per_second).
C CONSTS(9) is h_prime in component main (per_second).
C CONSTS(10) is g in component main (per_second).
C CONSTS(11) is n_H in component main (dimensionless).
C CONSTS(12) is u in component main (dimensionless).
C CONSTS(13) is w in component main (dimensionless).
C CONSTS(14) is v in component main (dimensionless).
C CONSTS(18) is k_u_on in component RU_rate_constant (per_second).
C CONSTS(19) is k_u_off in component RU_rate_constant (per_second).
C ALGBRC(8) is k_w_on in component RU_RU_interaction (per_second).
C ALGBRC(9) is k_w_off in component RU_RU_interaction (per_second).
C RATES(1) is d/dt D in component main (dimensionless).
C RATES(2) is d/dt A_1 in component main (dimensionless).
C RATES(3) is d/dt A_2 in component main (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      STATES(1) = 0
      STATES(2) = 0
      STATES(3) = 0
      CONSTS(2) = 0
      CONSTS(3) = 100
      CONSTS(4) = 120
      CONSTS(5) = 50
      CONSTS(6) = 50
      CONSTS(7) = 400
      CONSTS(8) = 8
      CONSTS(9) = 6
      CONSTS(10) = 4
      CONSTS(11) = 1
      CONSTS(12) = 1
      CONSTS(13) = 1
      CONSTS(14) = 1
      CONSTS(15) = CONSTS(5)/CONSTS(4)
      CONSTS(16) =  CONSTS(15)*100.000
      CONSTS(17) = 1.00000/(1.00000+CONSTS(16)/CONSTS(15) ** - CONSTS(11))
      CONSTS(18) = CONSTS(2)+( (CONSTS(4) - CONSTS(2))*CONSTS(16))/(CONSTS(15)+CONSTS(16))
      CONSTS(19) = CONSTS(3)+( (CONSTS(5) - CONSTS(3))*CONSTS(16))/(CONSTS(15)+CONSTS(16))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(3) =  CONSTS(8)*STATES(2) -  (CONSTS(9)+CONSTS(10))*STATES(3)
      ALGBRC(7) = STATES(3)/CONSTS(1)
      ALGBRC(10) =  CONSTS(6)*1.00000+ ALGBRC(7)*(EXP(CONSTS(14) - 1.00000) - 1.00000) ** 2.00000
      ALGBRC(11) =  CONSTS(7)*1.00000+ ALGBRC(7)*(EXP(- (CONSTS(14) - 1.00000)) - 1.00000) ** 2.00000
      RATES(2) = ( ALGBRC(10)*STATES(1)+ CONSTS(9)*STATES(3)) -  (ALGBRC(11)+CONSTS(8))*STATES(2)
      ALGBRC(1) = CONSTS(1) - (STATES(1)+STATES(2)+STATES(3))
      ALGBRC(6) = (STATES(1)+STATES(2)+STATES(3))/CONSTS(1)
      ALGBRC(8) =  CONSTS(18)*1.00000+ ALGBRC(6)*(CONSTS(12) - 1.00000) ** 2.00000
      ALGBRC(12) =  ALGBRC(8)*1.00000+ ALGBRC(7)*(EXP(CONSTS(13) - 1.00000) - 1.00000) ** 2.00000
      ALGBRC(9) =  CONSTS(19)*CONSTS(12) -  ALGBRC(6)*(CONSTS(12) - 1.00000) ** 2.00000
      ALGBRC(13) =  ALGBRC(9)*1.00000+ ALGBRC(7)*(EXP(- (CONSTS(13) - 1.00000)) - 1.00000) ** 2.00000
      RATES(1) = ( ALGBRC(12)*ALGBRC(1)+ ALGBRC(11)*STATES(2)+ CONSTS(10)*STATES(3)) -  (ALGBRC(13)+ALGBRC(10))*STATES(1)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = STATES(3)/CONSTS(1)
      ALGBRC(10) =  CONSTS(6)*1.00000+ ALGBRC(7)*(EXP(CONSTS(14) - 1.00000) - 1.00000) ** 2.00000
      ALGBRC(11) =  CONSTS(7)*1.00000+ ALGBRC(7)*(EXP(- (CONSTS(14) - 1.00000)) - 1.00000) ** 2.00000
      ALGBRC(1) = CONSTS(1) - (STATES(1)+STATES(2)+STATES(3))
      ALGBRC(6) = (STATES(1)+STATES(2)+STATES(3))/CONSTS(1)
      ALGBRC(8) =  CONSTS(18)*1.00000+ ALGBRC(6)*(CONSTS(12) - 1.00000) ** 2.00000
      ALGBRC(12) =  ALGBRC(8)*1.00000+ ALGBRC(7)*(EXP(CONSTS(13) - 1.00000) - 1.00000) ** 2.00000
      ALGBRC(9) =  CONSTS(19)*CONSTS(12) -  ALGBRC(6)*(CONSTS(12) - 1.00000) ** 2.00000
      ALGBRC(13) =  ALGBRC(9)*1.00000+ ALGBRC(7)*(EXP(- (CONSTS(13) - 1.00000)) - 1.00000) ** 2.00000
      ALGBRC(2) = STATES(1)/CONSTS(1)
      ALGBRC(3) = STATES(2)/CONSTS(1)
      ALGBRC(4) = STATES(3)/(STATES(1)+STATES(2)+STATES(3))
      ALGBRC(5) = ALGBRC(1)/CONSTS(1)
      RETURN
      END