Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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