C
C There are a total of 9 entries in the algebraic variable array.
C There are a total of 7 entries in each of the rate and state variable arrays.
C There are a total of 37 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is R_T in component Parameters (dimensionless).
C CONSTS(2) is K_1 in component Parameters (micromolar).
C CONSTS(3) is K_2 in component Parameters (micromolar).
C CONSTS(4) is k_r in component Parameters (per_second).
C CONSTS(5) is k_p in component Parameters (per_second).
C CONSTS(6) is k_e in component Parameters (per_second).
C CONSTS(7) is xi in component Parameters (dimensionless).
C CONSTS(8) is G_T in component Parameters (dimensionless).
C CONSTS(9) is k_deg in component Parameters (per_second).
C CONSTS(10) is k_a in component Parameters (per_second).
C CONSTS(11) is k_d in component Parameters (per_second).
C CONSTS(12) is PIP_2_T in component Parameters (dimensionless).
C CONSTS(13) is r_r in component Parameters (per_second).
C CONSTS(14) is delta in component Parameters (dimensionless).
C CONSTS(15) is K_c in component Parameters (micromolar).
C CONSTS(16) is alpha in component Parameters (per_second).
C CONSTS(17) is N_a in component Parameters (per_micromole).
C CONSTS(18) is v in component Parameters (litre).
C CONSTS(19) is epsilon_r in component Parameters (dimensionless).
C CONSTS(20) is d_1 in component Parameters (micromolar).
C CONSTS(21) is d_2 in component Parameters (micromolar).
C CONSTS(22) is d_3 in component Parameters (micromolar).
C CONSTS(23) is d_5 in component Parameters (micromolar).
C CONSTS(24) is a_2 in component Parameters (per_micromolar_per_second).
C CONSTS(25) is B_e in component Parameters (micromolar).
C CONSTS(26) is K_e in component Parameters (micromolar).
C CONSTS(27) is B_ER in component Parameters (micromolar).
C CONSTS(28) is K_ER in component Parameters (micromolar).
C CONSTS(29) is B_x in component Parameters (micromolar).
C CONSTS(30) is K_x in component Parameters (micromolar).
C CONSTS(31) is k_3 in component Parameters (micromolar).
C CONSTS(32) is eta_1 in component Parameters (per_second).
C CONSTS(33) is eta_2 in component Parameters (per_second).
C CONSTS(34) is eta_3 in component Parameters (flux).
C CONSTS(35) is C_T in component Parameters (micromolar).
C CONSTS(36) is L in component ligand (micromolar).
C STATES(1) is RS in component RS (dimensionless).
C STATES(2) is RS_p in component RS_p (dimensionless).
C STATES(3) is G in component G_GTP (dimensionless).
C ALGBRC(1) is rho_r in component rho_r (dimensionless).
C STATES(4) is IP_3 in component IP_3 (micromolar).
C STATES(5) is PIP_2 in component PIP_2 (dimensionless).
C ALGBRC(2) is r_h in component r_h (per_second).
C STATES(6) is C in component C (micromolar).
C ALGBRC(9) is C_ER in component C_ER (micromolar).
C ALGBRC(5) is beta in component beta (dimensionless).
C STATES(7) is h in component h (dimensionless).
C ALGBRC(3) is m_infinit in component m_infinit (dimensionless).
C ALGBRC(8) is h_infinit in component h_infinit (dimensionless).
C ALGBRC(6) is tau_h in component tau_h (second).
C ALGBRC(4) is zeta in component zeta (micromolar).
C ALGBRC(7) is gamma in component gamma (dimensionless).
C CONSTS(37) is RS_E in component RS_E (dimensionless).
C RATES(1) is d/dt RS in component RS (dimensionless).
C RATES(2) is d/dt RS_p in component RS_p (dimensionless).
C RATES(3) is d/dt G in component G_GTP (dimensionless).
C RATES(4) is d/dt IP_3 in component IP_3 (micromolar).
C RATES(5) is d/dt PIP_2 in component PIP_2 (dimensionless).
C RATES(6) is d/dt C in component C (micromolar).
C RATES(7) is d/dt h in component h (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 2e4
      CONSTS(2) = 5
      CONSTS(3) = 100
      CONSTS(4) = 1.75e-4
      CONSTS(5) = 0.03
      CONSTS(6) = 6e-3
      CONSTS(7) = 0.85
      CONSTS(8) = 1e5
      CONSTS(9) = 1.25
      CONSTS(10) = 0.017
      CONSTS(11) = 0.15
      CONSTS(12) = 5e7
      CONSTS(13) = 10
      CONSTS(14) = 1.238e-3
      CONSTS(15) = 0.4
      CONSTS(16) = 2.781e-5
      CONSTS(17) = 6.02252e17
      CONSTS(18) = 5e-13
      CONSTS(19) = 0.185
      CONSTS(20) = 0.13
      CONSTS(21) = 1.05
      CONSTS(22) = 0.943
      CONSTS(23) = 0.0823
      CONSTS(24) = 0.2
      CONSTS(25) = 150
      CONSTS(26) = 10
      CONSTS(27) = 120000
      CONSTS(28) = 1200
      CONSTS(29) = 50
      CONSTS(30) = 0.2
      CONSTS(31) = 0.4
      CONSTS(32) = 575
      CONSTS(33) = 5.2
      CONSTS(34) = 45
      CONSTS(35) = 67
      CONSTS(36) = 1000
      STATES(1) = 1.7e4
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0.01
      STATES(5) = 49997000
      STATES(6) = 0.0961
      STATES(7) = 0.6155
      CONSTS(37) =  (( CONSTS(4)*(1.00000+( (CONSTS(5)/CONSTS(6))*(CONSTS(3)+CONSTS(36)))/(CONSTS(2)+CONSTS(36))))/(CONSTS(4)+( CONSTS(5)*CONSTS(36))/(CONSTS(2)+CONSTS(36))+( (( CONSTS(4)*CONSTS(5))/CONSTS(6))*(CONSTS(3)+CONSTS(36)))/(CONSTS(2)+CONSTS(36))))*CONSTS(7)*CONSTS(1)+ (1.00000 - CONSTS(7))*CONSTS(1)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(4)*CONSTS(1) -  (CONSTS(4)+( CONSTS(5)*CONSTS(36))/(CONSTS(2)+CONSTS(36)))*STATES(1)) -  CONSTS(4)*STATES(2)
      RATES(2) =  CONSTS(36)*(( CONSTS(5)*STATES(1))/(CONSTS(2)+CONSTS(36)) - ( CONSTS(6)*STATES(2))/(CONSTS(3)+CONSTS(36)))
      ALGBRC(1) = ( CONSTS(36)*STATES(1))/( CONSTS(7)*CONSTS(1)*(CONSTS(2)+CONSTS(36)))
      RATES(3) =  CONSTS(10)*(CONSTS(14)+ALGBRC(1))*(CONSTS(8) - STATES(3)) -  CONSTS(11)*STATES(3)
      ALGBRC(2) =  (( CONSTS(16)*STATES(6))/(CONSTS(15)+STATES(6)))*STATES(3)
      RATES(4) = ( ALGBRC(2)*STATES(5))/( CONSTS(17)*CONSTS(18)) -  CONSTS(9)*STATES(4)
      RATES(5) = ( - (ALGBRC(2)+CONSTS(13))*STATES(5) -  CONSTS(13)*CONSTS(17)*CONSTS(18)*STATES(4))+ CONSTS(13)*CONSTS(12)
      ALGBRC(4) = ( CONSTS(21)*(STATES(4)+CONSTS(20)))/(STATES(4)+CONSTS(22))
      ALGBRC(8) = ALGBRC(4)/(ALGBRC(4)+STATES(6))
      ALGBRC(6) =  CONSTS(24)*(ALGBRC(4)+STATES(6)) ** - 1.00000
      RATES(7) = (ALGBRC(8) - STATES(7))/ALGBRC(6)
      ALGBRC(7) = 1.00000+CONSTS(25)/(CONSTS(26)+STATES(6))+CONSTS(29)/(CONSTS(30)+STATES(6)) ** - 1.00000
      ALGBRC(9) =  (CONSTS(28)/( CONSTS(27)*CONSTS(19)))*(CONSTS(35) - STATES(6)/ALGBRC(7))
      ALGBRC(5) = 1.00000+( CONSTS(26)*CONSTS(25))/CONSTS(26)+STATES(6) ** 2.00000+( CONSTS(30)*CONSTS(29))/CONSTS(30)+STATES(6) ** 2.00000 ** - 1.00000
      ALGBRC(3) = ( (STATES(4)/(CONSTS(20)+STATES(4)))*STATES(6))/(CONSTS(23)+STATES(6))
      RATES(6) =  ALGBRC(5)*( CONSTS(19)*( CONSTS(32)*ALGBRC(3) ** 3.00000*STATES(7) ** 3.00000+CONSTS(33))*(ALGBRC(9) - STATES(6)) - ( CONSTS(34)*STATES(6) ** 2.00000)/(CONSTS(31) ** 2.00000+STATES(6) ** 2.00000))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(36)*STATES(1))/( CONSTS(7)*CONSTS(1)*(CONSTS(2)+CONSTS(36)))
      ALGBRC(2) =  (( CONSTS(16)*STATES(6))/(CONSTS(15)+STATES(6)))*STATES(3)
      ALGBRC(4) = ( CONSTS(21)*(STATES(4)+CONSTS(20)))/(STATES(4)+CONSTS(22))
      ALGBRC(8) = ALGBRC(4)/(ALGBRC(4)+STATES(6))
      ALGBRC(6) =  CONSTS(24)*(ALGBRC(4)+STATES(6)) ** - 1.00000
      ALGBRC(7) = 1.00000+CONSTS(25)/(CONSTS(26)+STATES(6))+CONSTS(29)/(CONSTS(30)+STATES(6)) ** - 1.00000
      ALGBRC(9) =  (CONSTS(28)/( CONSTS(27)*CONSTS(19)))*(CONSTS(35) - STATES(6)/ALGBRC(7))
      ALGBRC(5) = 1.00000+( CONSTS(26)*CONSTS(25))/CONSTS(26)+STATES(6) ** 2.00000+( CONSTS(30)*CONSTS(29))/CONSTS(30)+STATES(6) ** 2.00000 ** - 1.00000
      ALGBRC(3) = ( (STATES(4)/(CONSTS(20)+STATES(4)))*STATES(6))/(CONSTS(23)+STATES(6))
      RETURN
      END