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 17 entries in the algebraic variable array. C There are a total of 5 entries in each of the rate and state variable arrays. C There are a total of 25 entries in the constant variable array. C C C VOI is time in component environment (ms). C ALGBRC(2) is Ca_i in component intracellular_ion_concentrations (mM). C ALGBRC(1) is Ca_b in component intracellular_ion_concentrations (mM). C STATES(1) is TRPN in component intracellular_ion_concentrations (mM). C CONSTS(1) is Ca_TRPN_Max in component troponin (mM). C ALGBRC(17) is J_TRPN in component troponin (mM_per_ms). C STATES(2) is z in component tropomyosin (dimensionless). C ALGBRC(10) is z_max in component tropomyosin (dimensionless). C CONSTS(2) is k_on in component troponin (per_mM_per_ms). C CONSTS(3) is k_Ref_off in component troponin (per_ms). C CONSTS(4) is gamma_trpn in component troponin (dimensionless). C CONSTS(5) is alpha_0 in component tropomyosin (per_ms). C CONSTS(6) is alpha_r1 in component tropomyosin (per_ms). C CONSTS(7) is alpha_r2 in component tropomyosin (per_ms). C CONSTS(8) is n_Rel in component tropomyosin (dimensionless). C CONSTS(9) is K_z in component tropomyosin (dimensionless). C CONSTS(10) is n_Hill in component tropomyosin (dimensionless). C CONSTS(11) is Ca_50ref in component tropomyosin (mM). C CONSTS(12) is z_p in component tropomyosin (dimensionless). C CONSTS(13) is beta_1 in component tropomyosin (dimensionless). C ALGBRC(7) is Ca_50 in component tropomyosin (mM). C ALGBRC(8) is Ca_TRPN_50 in component tropomyosin (mM). C CONSTS(23) is K_2 in component tropomyosin (per_ms). C CONSTS(25) is K_1 in component tropomyosin (per_ms). C ALGBRC(9) is alpha_Tm in component tropomyosin (per_ms). C ALGBRC(3) is beta_Tm in component tropomyosin (per_ms). C CONSTS(14) is beta_0 in component filament_overlap (dimensionless). C ALGBRC(5) is lambda in component Myofilaments (dimensionless). C ALGBRC(16) is k_off in component troponin (per_ms). C ALGBRC(15) is Tension in component Cross_Bridges (N_per_mm2). C CONSTS(15) is T_ref in component length_independent_tension (N_per_mm2). C ALGBRC(4) is ExtensionRatio in component Myofilaments (dimensionless). C CONSTS(24) is dExtensionRatiodt in component Myofilaments (per_ms). C ALGBRC(6) is lambda_prev in component Myofilaments (dimensionless). C ALGBRC(11) is overlap in component filament_overlap (dimensionless). C ALGBRC(12) is T_Base in component length_independent_tension (N_per_mm2). C ALGBRC(13) is T_0 in component isometric_tension (N_per_mm2). C ALGBRC(14) is Q in component Cross_Bridges (dimensionless). C CONSTS(16) is a in component Cross_Bridges (dimensionless). C STATES(3) is Q_1 in component Cross_Bridges (dimensionless). C STATES(4) is Q_2 in component Cross_Bridges (dimensionless). C STATES(5) is Q_3 in component Cross_Bridges (dimensionless). C CONSTS(17) is A_1 in component Cross_Bridges (dimensionless). C CONSTS(18) is A_2 in component Cross_Bridges (dimensionless). C CONSTS(19) is A_3 in component Cross_Bridges (dimensionless). C CONSTS(20) is alpha_1 in component Cross_Bridges (per_ms). C CONSTS(21) is alpha_2 in component Cross_Bridges (per_ms). C CONSTS(22) is alpha_3 in component Cross_Bridges (per_ms). C RATES(1) is d/dt TRPN in component intracellular_ion_concentrations (mM). C RATES(2) is d/dt z in component tropomyosin (dimensionless). C RATES(3) is d/dt Q_1 in component Cross_Bridges (dimensionless). C RATES(4) is d/dt Q_2 in component Cross_Bridges (dimensionless). C RATES(5) is d/dt Q_3 in component Cross_Bridges (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.067593139865 CONSTS(1) = 70e-3 STATES(2) = 0.014417937837 CONSTS(2) = 100 CONSTS(3) = 0.2 CONSTS(4) = 2 CONSTS(5) = 8e-3 CONSTS(6) = 2e-3 CONSTS(7) = 1.75e-3 CONSTS(8) = 3 CONSTS(9) = 0.15 CONSTS(10) = 3 CONSTS(11) = 1.05e-3 CONSTS(12) = 0.85 CONSTS(13) = -4 CONSTS(14) = 4.9 CONSTS(15) = 56.2 CONSTS(16) = 0.35 STATES(3) = 0 STATES(4) = 0 STATES(5) = 0 CONSTS(17) = -29 CONSTS(18) = 138 CONSTS(19) = 129 CONSTS(20) = 0.03 CONSTS(21) = 0.13 CONSTS(22) = 0.625 CONSTS(23) = (( CONSTS(7)*CONSTS(12) ** CONSTS(8))/(CONSTS(12) ** CONSTS(8)+CONSTS(9) ** CONSTS(8)))*(1.00000 - ( CONSTS(8)*CONSTS(9) ** CONSTS(8))/(CONSTS(12) ** CONSTS(8)+CONSTS(9) ** CONSTS(8))) CONSTS(24) = 0.00000 CONSTS(25) = ( CONSTS(7)*CONSTS(12) ** CONSTS(8) - 1.00000*CONSTS(8)*CONSTS(9) ** CONSTS(8))/CONSTS(12) ** CONSTS(8)+CONSTS(9) ** CONSTS(8) ** 2.00000 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(3) = CONSTS(17)*CONSTS(24) - CONSTS(20)*STATES(3) RATES(4) = CONSTS(18)*CONSTS(24) - CONSTS(21)*STATES(4) RATES(5) = CONSTS(19)*CONSTS(24) - CONSTS(22)*STATES(5) ALGBRC(1) = CONSTS(1) - STATES(1) ALGBRC(4) = TERNRY(VOI.GT.1000.00, 1.00000, 1.00000) ALGBRC(5) = TERNRY(ALGBRC(4).GT.0.800000.AND.ALGBRC(4).LE.1.15000, ALGBRC(4), TERNRY(ALGBRC(4).GT.1.15000, 1.15000, 0.800000) ALGBRC(7) = CONSTS(11)*(1.00000+ CONSTS(13)*(ALGBRC(5) - 1.00000)) ALGBRC(8) = ( ALGBRC(7)*CONSTS(1))/(ALGBRC(7)+ (CONSTS(3)/CONSTS(2))*(1.00000 - ( (1.00000+ CONSTS(14)*(ALGBRC(5) - 1.00000))*0.500000)/CONSTS(4))) ALGBRC(9) = CONSTS(5)*ALGBRC(1)/ALGBRC(8) ** CONSTS(10) ALGBRC(3) = CONSTS(6)+( CONSTS(7)*STATES(2) ** CONSTS(8) - 1.00000)/(STATES(2) ** CONSTS(8)+CONSTS(9) ** CONSTS(8)) RATES(2) = ALGBRC(9)*(1.00000 - STATES(2)) - ALGBRC(3)*STATES(2) ALGBRC(2) = TERNRY(VOI.LT.1.00000, 1000.00*1.84330e-07, TERNRY(VOI.GE.10.0000.AND.VOI.LT.15.0000, 1000.00*((( 1.05500*VOI/1000.00 ** 3.00000 - 0.0350700*VOI/1000.00 ** 2.00000)+( 0.000399200*VOI)/1000.00) - 1.35600e-06), TERNRY(VOI.GE.15.0000.AND.VOI.LT.55.0000, 1000.00*((( 0.0140000*VOI/1000.00 ** 3.00000 - 0.00255500*VOI/1000.00 ** 2.00000)+( 0.000149400*VOI)/1000.00) - 1.42800e-06), TERNRY(VOI.GE.55.0000.AND.VOI.LT.250.000, 1000.00*((( 1.73900e-05*VOI/1000.00 ** 3.00000 - 3.20900e-06*VOI/1000.00 ** 2.00000) - ( 5.68900e-06*VOI)/1000.00)+1.71900e-06), TERNRY(VOI.GE.250.000.AND.VOI.LT.490.000, 1000.00*(((( 0.000132100*VOI/1000.00 ** 4.00000 - 0.000219700*VOI/1000.00 ** 3.00000)+ 0.000137400*VOI/1000.00 ** 2.00000) - ( 3.89500e-05*VOI)/1000.00)+4.44100e-06), 1000.00*1.21480e-07) ALGBRC(11) = 1.00000+ CONSTS(14)*(ALGBRC(5) - 1.00000) ALGBRC(10) = (CONSTS(5)/ALGBRC(8)/CONSTS(1) ** CONSTS(10) - CONSTS(23))/(CONSTS(6)+CONSTS(25)+CONSTS(5)/ALGBRC(8)/CONSTS(1) ** CONSTS(10)) ALGBRC(12) = ( CONSTS(15)*STATES(2))/ALGBRC(10) ALGBRC(13) = ALGBRC(12)*ALGBRC(11) ALGBRC(14) = STATES(3)+STATES(4)+STATES(5) ALGBRC(15) = TERNRY(ALGBRC(14).LT.0.00000, ( ALGBRC(13)*( CONSTS(16)*ALGBRC(14)+1.00000))/(1.00000 - ALGBRC(14)), ( ALGBRC(13)*(1.00000+ (CONSTS(16)+2.00000)*ALGBRC(14)))/(1.00000+ALGBRC(14))) ALGBRC(16) = TERNRY(1.00000 - ALGBRC(15)/( CONSTS(4)*CONSTS(15)).GT.0.100000, CONSTS(3)*(1.00000 - ALGBRC(15)/( CONSTS(4)*CONSTS(15))), CONSTS(3)*0.100000) ALGBRC(17) = (CONSTS(1) - STATES(1))*ALGBRC(16) - ALGBRC(2)*STATES(1)*CONSTS(2) RATES(1) = ALGBRC(17) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = CONSTS(1) - STATES(1) ALGBRC(4) = TERNRY(VOI.GT.1000.00, 1.00000, 1.00000) ALGBRC(5) = TERNRY(ALGBRC(4).GT.0.800000.AND.ALGBRC(4).LE.1.15000, ALGBRC(4), TERNRY(ALGBRC(4).GT.1.15000, 1.15000, 0.800000) ALGBRC(7) = CONSTS(11)*(1.00000+ CONSTS(13)*(ALGBRC(5) - 1.00000)) ALGBRC(8) = ( ALGBRC(7)*CONSTS(1))/(ALGBRC(7)+ (CONSTS(3)/CONSTS(2))*(1.00000 - ( (1.00000+ CONSTS(14)*(ALGBRC(5) - 1.00000))*0.500000)/CONSTS(4))) ALGBRC(9) = CONSTS(5)*ALGBRC(1)/ALGBRC(8) ** CONSTS(10) ALGBRC(3) = CONSTS(6)+( CONSTS(7)*STATES(2) ** CONSTS(8) - 1.00000)/(STATES(2) ** CONSTS(8)+CONSTS(9) ** CONSTS(8)) ALGBRC(2) = TERNRY(VOI.LT.1.00000, 1000.00*1.84330e-07, TERNRY(VOI.GE.10.0000.AND.VOI.LT.15.0000, 1000.00*((( 1.05500*VOI/1000.00 ** 3.00000 - 0.0350700*VOI/1000.00 ** 2.00000)+( 0.000399200*VOI)/1000.00) - 1.35600e-06), TERNRY(VOI.GE.15.0000.AND.VOI.LT.55.0000, 1000.00*((( 0.0140000*VOI/1000.00 ** 3.00000 - 0.00255500*VOI/1000.00 ** 2.00000)+( 0.000149400*VOI)/1000.00) - 1.42800e-06), TERNRY(VOI.GE.55.0000.AND.VOI.LT.250.000, 1000.00*((( 1.73900e-05*VOI/1000.00 ** 3.00000 - 3.20900e-06*VOI/1000.00 ** 2.00000) - ( 5.68900e-06*VOI)/1000.00)+1.71900e-06), TERNRY(VOI.GE.250.000.AND.VOI.LT.490.000, 1000.00*(((( 0.000132100*VOI/1000.00 ** 4.00000 - 0.000219700*VOI/1000.00 ** 3.00000)+ 0.000137400*VOI/1000.00 ** 2.00000) - ( 3.89500e-05*VOI)/1000.00)+4.44100e-06), 1000.00*1.21480e-07) ALGBRC(11) = 1.00000+ CONSTS(14)*(ALGBRC(5) - 1.00000) ALGBRC(10) = (CONSTS(5)/ALGBRC(8)/CONSTS(1) ** CONSTS(10) - CONSTS(23))/(CONSTS(6)+CONSTS(25)+CONSTS(5)/ALGBRC(8)/CONSTS(1) ** CONSTS(10)) ALGBRC(12) = ( CONSTS(15)*STATES(2))/ALGBRC(10) ALGBRC(13) = ALGBRC(12)*ALGBRC(11) ALGBRC(14) = STATES(3)+STATES(4)+STATES(5) ALGBRC(15) = TERNRY(ALGBRC(14).LT.0.00000, ( ALGBRC(13)*( CONSTS(16)*ALGBRC(14)+1.00000))/(1.00000 - ALGBRC(14)), ( ALGBRC(13)*(1.00000+ (CONSTS(16)+2.00000)*ALGBRC(14)))/(1.00000+ALGBRC(14))) ALGBRC(16) = TERNRY(1.00000 - ALGBRC(15)/( CONSTS(4)*CONSTS(15)).GT.0.100000, CONSTS(3)*(1.00000 - ALGBRC(15)/( CONSTS(4)*CONSTS(15))), CONSTS(3)*0.100000) ALGBRC(17) = (CONSTS(1) - STATES(1))*ALGBRC(16) - ALGBRC(2)*STATES(1)*CONSTS(2) ALGBRC(6) = ALGBRC(4) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END