C
C There are a total of 4 entries in the algebraic variable array.
C There are a total of 2 entries in each of the rate and state variable arrays.
C There are a total of 13 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C ALGBRC(1) is B in component B (micromolar).
C CONSTS(1) is F in component B (micromolar).
C CONSTS(2) is n2 in component B (dimensionless).
C CONSTS(3) is K2 in component B (per_micromolar).
C STATES(1) is Ca in component Ca (micromolar).
C ALGBRC(4) is R in component R (flux).
C CONSTS(4) is Vmax in component R (flux).
C CONSTS(5) is Km in component R (micromolar).
C ALGBRC(3) is fu in component R (dimensionless).
C CONSTS(13) is ISF in component R (dimensionless).
C CONSTS(6) is age in component R (dimensionless).
C STATES(2) is C in component C (micromolar).
C CONSTS(7) is Va in component Ca (ml).
C CONSTS(8) is Vv in component Ca (ml).
C CONSTS(9) is Qc in component model_constants (flow).
C ALGBRC(2) is Cv in component Cv (micromolar).
C CONSTS(10) is Q in component model_constants (flow).
C CONSTS(11) is P in component model_constants (dimensionless).
C CONSTS(12) is V in component C (ml).
C RATES(1) is d/dt Ca in component Ca (micromolar).
C RATES(2) is d/dt C in component C (micromolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.48
      CONSTS(2) = 1
      CONSTS(3) = 0.8532
      STATES(1) = 6.6685
      CONSTS(4) = 9.433e-3
      CONSTS(5) = 198
      CONSTS(6) = 5
      STATES(2) = 0
      CONSTS(7) = 2148
      CONSTS(8) = 3431
      CONSTS(9) = 6445.65
      CONSTS(10) = 1221.34
      CONSTS(11) = 15.61
      CONSTS(12) = 1454
      CONSTS(13) = - 8.32120+ 2.04010*CONSTS(6)+ 4.19620*arbitrary_log( CONSTS(6)*365.000, 10)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = (( CONSTS(10)*STATES(2))/CONSTS(11))/CONSTS(9)
      RATES(1) = ( CONSTS(9)*(ALGBRC(2) - STATES(1)))/(CONSTS(7)+CONSTS(8))
      ALGBRC(1) = ( CONSTS(1)*CONSTS(2)*CONSTS(3)*STATES(1))/(1.00000+ CONSTS(3)*CONSTS(1))
      ALGBRC(3) = CONSTS(1)/(CONSTS(1)+ALGBRC(1))
      ALGBRC(4) = ( CONSTS(13)*CONSTS(4)*ALGBRC(3)*STATES(2))/(CONSTS(5)+ ALGBRC(3)*STATES(2))
      RATES(2) = ( CONSTS(10)*(STATES(1) - STATES(2)/CONSTS(11)) -  ALGBRC(4)*1.00000)/CONSTS(12)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = (( CONSTS(10)*STATES(2))/CONSTS(11))/CONSTS(9)
      ALGBRC(1) = ( CONSTS(1)*CONSTS(2)*CONSTS(3)*STATES(1))/(1.00000+ CONSTS(3)*CONSTS(1))
      ALGBRC(3) = CONSTS(1)/(CONSTS(1)+ALGBRC(1))
      ALGBRC(4) = ( CONSTS(13)*CONSTS(4)*ALGBRC(3)*STATES(2))/(CONSTS(5)+ ALGBRC(3)*STATES(2))
      RETURN
      END