C
C There are a total of 11 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 27 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C ALGBRC(3) is Rate_Ca_influx_across_the_SR in component Ca_influx_across_the_SR (flux).
C CONSTS(1) is Ca_e in component extracellular_calcium (molar).
C STATES(1) is Ca_f in component fuzzy_space_calcium (molar).
C CONSTS(2) is k1 in component Ca_influx_across_the_SR (first_order_rate_constant).
C STATES(2) is Ca_2_S1 in component Ca_bound_to_the_SRRC_fast_activating_binding_site (molar).
C STATES(3) is Ca_S2 in component Ca_bound_to_the_SRRC_slow_inactivating_binding_site (molar).
C STATES(4) is S1 in component SRRC_fast_activating_binding_site (molar).
C STATES(5) is S2 in component SRRC_slow_inactivating_binding_site (molar).
C STATES(6) is Ca_s in component SR_calcium (molar).
C ALGBRC(1) is dCa2_S1_dt in component Ca_movement_through_the_SRRC (flux).
C ALGBRC(2) is dCa_S2_dt in component Ca_movement_through_the_SRRC (flux).
C CONSTS(3) is k_on1 in component Ca_movement_through_the_SRRC (second_order_rate_constant).
C CONSTS(4) is k_off1 in component Ca_movement_through_the_SRRC (first_order_rate_constant).
C CONSTS(5) is k_on2 in component Ca_movement_through_the_SRRC (second_order_rate_constant).
C CONSTS(6) is k_off2 in component Ca_movement_through_the_SRRC (first_order_rate_constant).
C CONSTS(7) is k_s in component Ca_movement_through_the_SRRC (first_order_rate_constant).
C STATES(7) is cas1 in component Ca_movement_through_the_SRRC (dimensionless).
C STATES(8) is cas2 in component Ca_movement_through_the_SRRC (dimensionless).
C ALGBRC(5) is dcas1_dt in component Ca_movement_through_the_SRRC (first_order_rate_constant).
C ALGBRC(6) is dcas2_dt in component Ca_movement_through_the_SRRC (first_order_rate_constant).
C ALGBRC(4) is r_o in component Ca_movement_through_the_SRRC (dimensionless).
C ALGBRC(7) is Rate_Ca_movement_through_the_SRRC in component Ca_movement_through_the_SRRC (flux).
C CONSTS(8) is Km_NaCaX in component Ca_efflux_across_the_SR_by_NaCa_exchange (molar).
C CONSTS(9) is Vmax_NaCaX in component Ca_efflux_across_the_SR_by_NaCa_exchange (flux).
C ALGBRC(8) is Rate_Ca_efflux_across_the_SR_by_NaCa_exchange in component Ca_efflux_across_the_SR_by_NaCa_exchange (flux).
C STATES(9) is Ca_c in component cytosolic_calcium (molar).
C CONSTS(10) is kf in component Ca_movement_between_the_fuzzy_space_and_cytosol (first_order_rate_constant).
C ALGBRC(9) is Rate_Ca_movement_between_the_fuzzy_space_and_cytosol in component Ca_movement_between_the_fuzzy_space_and_cytosol (flux).
C CONSTS(11) is Km_s in component Ca_uptake_by_SR_Ca_ATPase (molar).
C CONSTS(12) is Vmax_s in component Ca_uptake_by_SR_Ca_ATPase (flux).
C ALGBRC(10) is Rate_Ca_uptake_by_SR_Ca_ATPase in component Ca_uptake_by_SR_Ca_ATPase (flux).
C STATES(10) is Ca_CSQ in component calsequestrin_bound_calcium (molar).
C CONSTS(13) is K_ons in component Ca_buffering_in_the_SR (second_order_rate_constant).
C CONSTS(14) is K_offs in component Ca_buffering_in_the_SR (first_order_rate_constant).
C CONSTS(15) is Bmax_s in component Ca_buffering_in_the_SR (molar).
C ALGBRC(11) is Rate_Ca_buffering_in_the_SR in component Ca_buffering_in_the_SR (flux).
C CONSTS(16) is Rt in component fuzzy_space_calcium (molar).
C CONSTS(17) is Bmax_f1 in component fuzzy_space_calcium (molar).
C CONSTS(18) is Bmax_f2 in component fuzzy_space_calcium (molar).
C CONSTS(19) is Kb_f1 in component fuzzy_space_calcium (molar).
C CONSTS(20) is Kb_f2 in component fuzzy_space_calcium (molar).
C CONSTS(21) is V_f in component fuzzy_space_calcium (dimensionless).
C CONSTS(22) is Bmax_c in component cytosolic_calcium (molar).
C CONSTS(23) is dye_c in component cytosolic_calcium (molar).
C CONSTS(24) is Kb_c in component cytosolic_calcium (molar).
C CONSTS(25) is Kb_dye in component cytosolic_calcium (molar).
C CONSTS(26) is V_c in component cytosolic_calcium (dimensionless).
C CONSTS(27) is V_s in component SR_calcium (dimensionless).
C RATES(7) is d/dt cas1 in component Ca_movement_through_the_SRRC (dimensionless).
C RATES(8) is d/dt cas2 in component Ca_movement_through_the_SRRC (dimensionless).
C RATES(1) is d/dt Ca_f in component fuzzy_space_calcium (molar).
C RATES(9) is d/dt Ca_c in component cytosolic_calcium (molar).
C RATES(6) is d/dt Ca_s in component SR_calcium (molar).
C RATES(10) is d/dt Ca_CSQ in component calsequestrin_bound_calcium (molar).
C RATES(4) is d/dt S1 in component SRRC_fast_activating_binding_site (molar).
C RATES(5) is d/dt S2 in component SRRC_slow_inactivating_binding_site (molar).
C RATES(2) is d/dt Ca_2_S1 in component Ca_bound_to_the_SRRC_fast_activating_binding_site (molar).
C RATES(3) is d/dt Ca_S2 in component Ca_bound_to_the_SRRC_slow_inactivating_binding_site (molar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.002
      STATES(1) = 0.12e-6
      CONSTS(2) = 0.2
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 201e-6
      CONSTS(3) = 2000000000
      CONSTS(4) = 1400
      CONSTS(5) = 13000000
      CONSTS(6) = 3.9
      CONSTS(7) = 9
      STATES(7) = 0
      STATES(8) = 0
      CONSTS(8) = 0.000036
      CONSTS(9) = 0.0012
      STATES(9) = 1e-7
      CONSTS(10) = 2500
      CONSTS(11) = 0.00000025
      CONSTS(12) = 0.000525
      STATES(10) = 0
      CONSTS(13) = 8772
      CONSTS(14) = 5.596536
      CONSTS(15) = 0.008
      CONSTS(16) = 0.00000015
      CONSTS(17) = 0.0002
      CONSTS(18) = 0.0011
      CONSTS(19) = 0.000017
      CONSTS(20) = 0.000013
      CONSTS(21) = 0.0013
      CONSTS(22) = 0.00012
      CONSTS(23) = 0
      CONSTS(24) = 0.00000096
      CONSTS(25) = 2e-7
      CONSTS(26) = 0.9287
      CONSTS(27) = 0.07
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(4) -  (CONSTS(4) ** 2.00000/( CONSTS(3)*STATES(1)))*STATES(2)
      RATES(4) = - ALGBRC(1)
      ALGBRC(2) =  CONSTS(5)*STATES(1)*STATES(5) -  CONSTS(6)*STATES(3)
      RATES(5) = - ALGBRC(2)
      RATES(2) = ALGBRC(1)
      RATES(3) = ALGBRC(2)
      ALGBRC(5) =  CONSTS(3)*STATES(1)*(1.00000 - STATES(7)) -  (CONSTS(4) ** 2.00000/( CONSTS(3)*STATES(1)))*STATES(7)
      RATES(7) = ALGBRC(5)
      ALGBRC(6) =  CONSTS(5)*STATES(1)*(1.00000 - STATES(8)) -  CONSTS(6)*STATES(8)
      RATES(8) = ALGBRC(6)
      ALGBRC(3) =  CONSTS(2)*(CONSTS(1) - STATES(1))
      ALGBRC(4) =  STATES(7)*(1.00000 - STATES(8))
      ALGBRC(7) =  CONSTS(7)*ALGBRC(4)*(STATES(6) - STATES(1))
      ALGBRC(8) = ( CONSTS(9)*STATES(1))/(CONSTS(8)+STATES(1))
      ALGBRC(9) =  CONSTS(10)*(STATES(1) - STATES(9))
      RATES(1) = ((ALGBRC(7) - ( CONSTS(16)*(ALGBRC(5)+ALGBRC(6))+ALGBRC(9)+ALGBRC(8)))+ALGBRC(3))/(( CONSTS(17)*CONSTS(19))/STATES(1)+CONSTS(19) ** 2.00000+( CONSTS(18)*CONSTS(20))/STATES(1)+CONSTS(20) ** 2.00000+CONSTS(21))
      ALGBRC(10) = ( CONSTS(12)*(STATES(9) ** 2.00000 - STATES(6) ** 2.00000/7000.00 ** 2.00000))/(CONSTS(11) ** 2.00000+STATES(9) ** 2.00000+STATES(6) ** 2.00000/7000.00 ** 2.00000)
      RATES(9) = (ALGBRC(9) - ALGBRC(10))/(( CONSTS(22)*CONSTS(24))/STATES(9)+CONSTS(24) ** 2.00000+( CONSTS(23)*CONSTS(25))/STATES(9)+CONSTS(25) ** 2.00000+CONSTS(26))
      ALGBRC(11) =  CONSTS(13)*STATES(6)*(CONSTS(15) - STATES(10)) -  CONSTS(14)*STATES(10)
      RATES(6) = (ALGBRC(10) - ALGBRC(7))/CONSTS(27) - ALGBRC(11)
      RATES(10) = ALGBRC(11)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(4) -  (CONSTS(4) ** 2.00000/( CONSTS(3)*STATES(1)))*STATES(2)
      ALGBRC(2) =  CONSTS(5)*STATES(1)*STATES(5) -  CONSTS(6)*STATES(3)
      ALGBRC(5) =  CONSTS(3)*STATES(1)*(1.00000 - STATES(7)) -  (CONSTS(4) ** 2.00000/( CONSTS(3)*STATES(1)))*STATES(7)
      ALGBRC(6) =  CONSTS(5)*STATES(1)*(1.00000 - STATES(8)) -  CONSTS(6)*STATES(8)
      ALGBRC(3) =  CONSTS(2)*(CONSTS(1) - STATES(1))
      ALGBRC(4) =  STATES(7)*(1.00000 - STATES(8))
      ALGBRC(7) =  CONSTS(7)*ALGBRC(4)*(STATES(6) - STATES(1))
      ALGBRC(8) = ( CONSTS(9)*STATES(1))/(CONSTS(8)+STATES(1))
      ALGBRC(9) =  CONSTS(10)*(STATES(1) - STATES(9))
      ALGBRC(10) = ( CONSTS(12)*(STATES(9) ** 2.00000 - STATES(6) ** 2.00000/7000.00 ** 2.00000))/(CONSTS(11) ** 2.00000+STATES(9) ** 2.00000+STATES(6) ** 2.00000/7000.00 ** 2.00000)
      ALGBRC(11) =  CONSTS(13)*STATES(6)*(CONSTS(15) - STATES(10)) -  CONSTS(14)*STATES(10)
      RETURN
      END