C
C There are a total of 15 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 29 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is A_iso in component parameters (mN_per_mm2_per_mM).
C CONSTS(2) is k_off in component parameters (per_msec).
C CONSTS(3) is k_on in component parameters (per_mM_per_msec).
C CONSTS(4) is k_onI in component parameters (per_msec).
C CONSTS(5) is k_offi in component parameters (per_msec).
C CONSTS(6) is k_OFF in component parameters (dimensionless).
C CONSTS(7) is k_xboff in component parameters (dimensionless).
C CONSTS(8) is k_tmoff in component parameters (per_msec).
C CONSTS(9) is k_tmonc in component parameters (per_msec).
C CONSTS(10) is k_tmRU in component parameters (dimensionless).
C CONSTS(11) is k_tmxb in component parameters (dimensionless).
C CONSTS(12) is k_12 in component parameters (per_mM_per_msec).
C CONSTS(13) is k_minus12 in component parameters (per_msec).
C CONSTS(14) is k_Titin in component parameters (dimensionless).
C CONSTS(15) is k_3c in component parameters (per_msec).
C CONSTS(16) is k_3f in component parameters (dimensionless).
C CONSTS(17) is k_3xb in component parameters (dimensionless).
C CONSTS(18) is k_minus3 in component parameters (per_msec).
C CONSTS(19) is k_4 in component parameters (per_msec).
C CONSTS(20) is k_minus4 in component parameters (per_msec).
C CONSTS(21) is k_5 in component parameters (per_msec).
C ALGBRC(1) is F_b in component active_force (mN_per_mm2).
C STATES(1) is RUA_MADPPi in component RU_species (mM).
C STATES(2) is RUA_MADP in component RU_species (mM).
C ALGBRC(6) is Q_CaB in component Ca_binding_rate (mM_per_msec).
C ALGBRC(4) is RUNA in component RU_species (mM).
C STATES(3) is Ca in component RU_species (mM).
C STATES(4) is RUTCa_off in component RU_species (mM).
C ALGBRC(12) is Q_TCaA in component TnI_conf_change_rate (mM_per_msec).
C ALGBRC(10) is k_offI in component TnI_conf_change_rate (per_msec).
C STATES(5) is RUTCa_on in component RU_species (mM).
C ALGBRC(8) is RUA in component RU_species (mM).
C CONSTS(22) is RU_total in component RU_species (mM).
C ALGBRC(3) is Q_TMA in component Tm_conf_change_rate (mM_per_msec).
C ALGBRC(2) is k_tmon in component Tm_conf_change_rate (per_msec).
C STATES(6) is RUTM_on in component RU_species (mM).
C ALGBRC(11) is Q_MB in component myosin_binding_rate (mM_per_msec).
C CONSTS(29) is KTitin in component KTitin (dimensionless).
C CONSTS(25) is L in component length (mmeter).
C CONSTS(23) is L_0 in component length (mmeter).
C CONSTS(24) is MADPPi in component RU_species (mM).
C ALGBRC(9) is RUTM_on_eff in component RU_species (mM).
C STATES(7) is RUAMADPPi in component RU_species (mM).
C CONSTS(26) is f_titinNormal in component KTitin (dimensionless).
C CONSTS(28) is f_titinDamaged in component KTitin (dimensionless).
C ALGBRC(15) is Q_Fgen in component Force_Generating_rate (mM_per_msec).
C ALGBRC(14) is k_3 in component Force_Generating_rate (per_msec).
C ALGBRC(13) is RUTCa_on_eff in component RU_species (mM).
C ALGBRC(5) is Q_PiR in component Pi_release_rate (mM_per_msec).
C ALGBRC(7) is Q_ADPR in component ADP_release_rate (mM_per_msec).
C CONSTS(27) is alpha in component alpha (dimensionless).
C RATES(4) is d/dt RUTCa_off in component RU_species (mM).
C RATES(5) is d/dt RUTCa_on in component RU_species (mM).
C RATES(6) is d/dt RUTM_on in component RU_species (mM).
C RATES(7) is d/dt RUAMADPPi in component RU_species (mM).
C RATES(1) is d/dt RUA_MADPPi in component RU_species (mM).
C RATES(2) is d/dt RUA_MADP in component RU_species (mM).
C RATES(3) is d/dt Ca in component RU_species (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 9000
      CONSTS(2) = 0.2
      CONSTS(3) = 17.3
      CONSTS(4) = 0.2
      CONSTS(5) = 0.075
      CONSTS(6) = -0.32
      CONSTS(7) = -1.37
      CONSTS(8) = 0.067
      CONSTS(9) = 0.014
      CONSTS(10) = 10
      CONSTS(11) = -1.86
      CONSTS(12) = 2
      CONSTS(13) = 0.6
      CONSTS(14) = 0.2
      CONSTS(15) = 0.025
      CONSTS(16) = 50
      CONSTS(17) = -1.23
      CONSTS(18) = 0.008
      CONSTS(19) = 0.077
      CONSTS(20) = 0.001
      CONSTS(21) = 0.03723
      STATES(1) = 0
      STATES(2) = 0
      STATES(3) = 0.01
      STATES(4) = 0
      STATES(5) = 0
      CONSTS(22) = 0.0726
      STATES(6) = 0
      CONSTS(23) = 0.001
      CONSTS(24) = 0.1375
      STATES(7) = 0
      CONSTS(25) = CONSTS(23)
      CONSTS(26) =  CONSTS(14)*(9.96630/(1.00000+EXP((CONSTS(25)/CONSTS(23)+1.06390)/0.0696000)))
      CONSTS(27) = TERNRY(CONSTS(25).LT.1.00000,  1.50000*(CONSTS(25)/CONSTS(23)) - 0.500000, TERNRY(1.00000.LE.CONSTS(25).AND.CONSTS(25).LT.1.10000, 1.00000,  - 1.60000*(CONSTS(25)/CONSTS(23))+2.76000)
      CONSTS(28) =  CONSTS(14)*( - 31.0000*(CONSTS(25)/CONSTS(23))+40.0000)
      CONSTS(29) = TERNRY(CONSTS(25).LE.1.10000, CONSTS(26), CONSTS(28))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) =  CONSTS(19)*STATES(1) -  CONSTS(20)*STATES(2)
      ALGBRC(7) =  CONSTS(21)*STATES(2)
      RATES(2) = ALGBRC(5) - ALGBRC(7)
      ALGBRC(4) = CONSTS(22) - (STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(1)+STATES(2))
      ALGBRC(6) =  CONSTS(3)*ALGBRC(4)*STATES(3) -  CONSTS(2)*STATES(4)
      RATES(3) = - ALGBRC(6)
      ALGBRC(2) =  CONSTS(9)*(STATES(5)/CONSTS(22))*1.00000+ CONSTS(10)*(STATES(6)/CONSTS(22)) ** 2.00000*1.00000+ CONSTS(11)*((STATES(1)+STATES(2))/CONSTS(22)) ** 2.00000
      ALGBRC(3) =  ALGBRC(2)*STATES(5) -  CONSTS(8)*STATES(6)
      ALGBRC(9) =  CONSTS(27)*STATES(6)
      ALGBRC(11) =  CONSTS(12)*CONSTS(29)*CONSTS(24)*ALGBRC(9) -  CONSTS(13)*STATES(7)
      RATES(6) = (ALGBRC(3) - ALGBRC(11))+ALGBRC(7)
      ALGBRC(8) = STATES(5)+STATES(6)+STATES(7)+STATES(1)+STATES(2)
      ALGBRC(10) =  CONSTS(5)*1.00000+ CONSTS(6)*(ALGBRC(8)/CONSTS(22)) ** 2.00000*1.00000+ CONSTS(7)*((STATES(1)+STATES(2))/CONSTS(22)) ** 4.40000
      ALGBRC(12) =  CONSTS(4)*STATES(4) -  ALGBRC(10)*STATES(5)
      RATES(4) = ALGBRC(6) - ALGBRC(12)
      RATES(5) = ALGBRC(12) - ALGBRC(3)
      ALGBRC(13) =  CONSTS(27)*STATES(5)
      ALGBRC(14) =  CONSTS(15)*1.00000+( CONSTS(16)*(ALGBRC(13)+ALGBRC(9)+STATES(7)))/CONSTS(22) ** 2.00000*1.00000+( CONSTS(17)*(STATES(1)+STATES(2)))/CONSTS(22) ** 2.00000
      ALGBRC(15) =  ALGBRC(14)*STATES(7) -  CONSTS(18)*STATES(1)
      RATES(7) = ALGBRC(11) - ALGBRC(15)
      RATES(1) = ALGBRC(15) - ALGBRC(5)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) =  CONSTS(19)*STATES(1) -  CONSTS(20)*STATES(2)
      ALGBRC(7) =  CONSTS(21)*STATES(2)
      ALGBRC(4) = CONSTS(22) - (STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(1)+STATES(2))
      ALGBRC(6) =  CONSTS(3)*ALGBRC(4)*STATES(3) -  CONSTS(2)*STATES(4)
      ALGBRC(2) =  CONSTS(9)*(STATES(5)/CONSTS(22))*1.00000+ CONSTS(10)*(STATES(6)/CONSTS(22)) ** 2.00000*1.00000+ CONSTS(11)*((STATES(1)+STATES(2))/CONSTS(22)) ** 2.00000
      ALGBRC(3) =  ALGBRC(2)*STATES(5) -  CONSTS(8)*STATES(6)
      ALGBRC(9) =  CONSTS(27)*STATES(6)
      ALGBRC(11) =  CONSTS(12)*CONSTS(29)*CONSTS(24)*ALGBRC(9) -  CONSTS(13)*STATES(7)
      ALGBRC(8) = STATES(5)+STATES(6)+STATES(7)+STATES(1)+STATES(2)
      ALGBRC(10) =  CONSTS(5)*1.00000+ CONSTS(6)*(ALGBRC(8)/CONSTS(22)) ** 2.00000*1.00000+ CONSTS(7)*((STATES(1)+STATES(2))/CONSTS(22)) ** 4.40000
      ALGBRC(12) =  CONSTS(4)*STATES(4) -  ALGBRC(10)*STATES(5)
      ALGBRC(13) =  CONSTS(27)*STATES(5)
      ALGBRC(14) =  CONSTS(15)*1.00000+( CONSTS(16)*(ALGBRC(13)+ALGBRC(9)+STATES(7)))/CONSTS(22) ** 2.00000*1.00000+( CONSTS(17)*(STATES(1)+STATES(2)))/CONSTS(22) ** 2.00000
      ALGBRC(15) =  ALGBRC(14)*STATES(7) -  CONSTS(18)*STATES(1)
      ALGBRC(1) =  CONSTS(1)*(STATES(1)+STATES(2))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END