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 2 entries in the algebraic variable array. C There are a total of 13 entries in each of the rate and state variable arrays. C There are a total of 65 entries in the constant variable array. C C C VOI is time in component environment (ms). C CONSTS(1) is V in component membrane (V). C ALGBRC(2) is i_Na in component sodium_current (dimensionless). C CONSTS(2) is E_Na in component sodium_current (V). C CONSTS(3) is g_Na in component sodium_current (per_V). C ALGBRC(1) is P_open in component Markov_model (dimensionless). C STATES(1) is O1 in component Markov_model (dimensionless). C STATES(2) is O2 in component Markov_model (dimensionless). C STATES(3) is C0 in component Markov_model (dimensionless). C STATES(4) is C1 in component Markov_model (dimensionless). C STATES(5) is C2 in component Markov_model (dimensionless). C STATES(6) is C3 in component Markov_model (dimensionless). C STATES(7) is C4 in component Markov_model (dimensionless). C STATES(8) is C0I in component Markov_model (dimensionless). C STATES(9) is C1I in component Markov_model (dimensionless). C STATES(10) is C2I in component Markov_model (dimensionless). C STATES(11) is C3I in component Markov_model (dimensionless). C STATES(12) is C4I in component Markov_model (dimensionless). C STATES(13) is I in component Markov_model (dimensionless). C CONSTS(4) is a in component Markov_model (dimensionless). C CONSTS(52) is alpha in component rate_constants (first_order_rate_constant). C CONSTS(53) is beta in component rate_constants (first_order_rate_constant). C CONSTS(54) is cf in component rate_constants (first_order_rate_constant). C CONSTS(55) is cn in component rate_constants (first_order_rate_constant). C CONSTS(56) is of in component rate_constants (first_order_rate_constant). C CONSTS(57) is on in component rate_constants (first_order_rate_constant). C CONSTS(58) is eta in component rate_constants (first_order_rate_constant). C CONSTS(59) is gamma in component rate_constants (first_order_rate_constant). C CONSTS(60) is delta in component rate_constants (first_order_rate_constant). C CONSTS(61) is epsilon in component rate_constants (first_order_rate_constant). C CONSTS(62) is omega in component rate_constants (first_order_rate_constant). C CONSTS(63) is v in component rate_constants (first_order_rate_constant). C CONSTS(64) is gamma_gamma in component rate_constants (first_order_rate_constant). C CONSTS(65) is delta_delta in component rate_constants (first_order_rate_constant). C CONSTS(5) is R in component rate_constants (J_mol_K). C CONSTS(6) is T in component rate_constants (K). C CONSTS(7) is F in component rate_constants (C_mol). C CONSTS(8) is k in component rate_constants (J_K). C CONSTS(9) is h in component rate_constants (J_ms). C CONSTS(10) is z_alpha in component rate_constants (dimensionless). C CONSTS(11) is z_beta in component rate_constants (dimensionless). C CONSTS(12) is z_gamma in component rate_constants (dimensionless). C CONSTS(13) is z_delta in component rate_constants (dimensionless). C CONSTS(14) is z_on in component rate_constants (dimensionless). C CONSTS(15) is z_of in component rate_constants (dimensionless). C CONSTS(16) is z_gamma_gamma in component rate_constants (dimensionless). C CONSTS(17) is z_delta_delta in component rate_constants (dimensionless). C CONSTS(18) is z_epsilon in component rate_constants (dimensionless). C CONSTS(19) is z_omega in component rate_constants (dimensionless). C CONSTS(20) is z_eta in component rate_constants (dimensionless). C CONSTS(21) is z_v in component rate_constants (dimensionless). C CONSTS(22) is z_cn in component rate_constants (dimensionless). C CONSTS(23) is z_cf in component rate_constants (dimensionless). C CONSTS(24) is delta_H_alpha in component rate_constants (J_mol). C CONSTS(25) is delta_H_beta in component rate_constants (J_mol). C CONSTS(26) is delta_H_cf in component rate_constants (J_mol). C CONSTS(27) is delta_H_cn in component rate_constants (J_mol). C CONSTS(28) is delta_H_of in component rate_constants (J_mol). C CONSTS(29) is delta_H_on in component rate_constants (J_mol). C CONSTS(30) is delta_H_eta in component rate_constants (J_mol). C CONSTS(31) is delta_H_gamma in component rate_constants (J_mol). C CONSTS(32) is delta_H_delta in component rate_constants (J_mol). C CONSTS(33) is delta_H_epsilon in component rate_constants (J_mol). C CONSTS(34) is delta_H_omega in component rate_constants (J_mol). C CONSTS(35) is delta_H_v in component rate_constants (J_mol). C CONSTS(36) is delta_H_gamma_gamma in component rate_constants (J_mol). C CONSTS(37) is delta_H_delta_delta in component rate_constants (J_mol). C CONSTS(38) is delta_S_alpha in component rate_constants (J_mol_K). C CONSTS(39) is delta_S_beta in component rate_constants (J_mol_K). C CONSTS(40) is delta_S_cf in component rate_constants (J_mol_K). C CONSTS(41) is delta_S_cn in component rate_constants (J_mol_K). C CONSTS(42) is delta_S_of in component rate_constants (J_mol_K). C CONSTS(43) is delta_S_on in component rate_constants (J_mol_K). C CONSTS(44) is delta_S_eta in component rate_constants (J_mol_K). C CONSTS(45) is delta_S_gamma in component rate_constants (J_mol_K). C CONSTS(46) is delta_S_delta in component rate_constants (J_mol_K). C CONSTS(47) is delta_S_epsilon in component rate_constants (J_mol_K). C CONSTS(48) is delta_S_omega in component rate_constants (J_mol_K). C CONSTS(49) is delta_S_v in component rate_constants (J_mol_K). C CONSTS(50) is delta_S_gamma_gamma in component rate_constants (J_mol_K). C CONSTS(51) is delta_S_delta_delta in component rate_constants (J_mol_K). C RATES(3) is d/dt C0 in component Markov_model (dimensionless). C RATES(4) is d/dt C1 in component Markov_model (dimensionless). C RATES(5) is d/dt C2 in component Markov_model (dimensionless). C RATES(6) is d/dt C3 in component Markov_model (dimensionless). C RATES(7) is d/dt C4 in component Markov_model (dimensionless). C RATES(1) is d/dt O1 in component Markov_model (dimensionless). C RATES(2) is d/dt O2 in component Markov_model (dimensionless). C RATES(8) is d/dt C0I in component Markov_model (dimensionless). C RATES(9) is d/dt C1I in component Markov_model (dimensionless). C RATES(10) is d/dt C2I in component Markov_model (dimensionless). C RATES(11) is d/dt C3I in component Markov_model (dimensionless). C RATES(12) is d/dt C4I in component Markov_model (dimensionless). C RATES(13) is d/dt I in component Markov_model (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = -0.12 CONSTS(2) = 0.044675 CONSTS(3) = 0.0131 STATES(1) = 0.0 STATES(2) = 0.0 STATES(3) = 1.0 STATES(4) = 0.0 STATES(5) = 0.0 STATES(6) = 0.0 STATES(7) = 0.0 STATES(8) = 0.0 STATES(9) = 0.0 STATES(10) = 0.0 STATES(11) = 0.0 STATES(12) = 0.0 STATES(13) = 0.0 CONSTS(4) = 2.5218 CONSTS(5) = 8.314472 CONSTS(6) = 286.0 CONSTS(7) = 96500.0 CONSTS(8) = 1.3806504E-23 CONSTS(9) = 6.62607095E-31 CONSTS(10) = 0 CONSTS(11) = -0.9701 CONSTS(12) = 1.5703 CONSTS(13) = -1.3266 CONSTS(14) = 0.6625 CONSTS(15) = 0 CONSTS(16) = 0 CONSTS(17) = -3.5596 CONSTS(18) = 0 CONSTS(19) = 0 CONSTS(20) = 1.5717 CONSTS(21) = -1.3281 CONSTS(22) = 0 CONSTS(23) = 0 CONSTS(24) = 116900 CONSTS(25) = 263870 CONSTS(26) = 57533 CONSTS(27) = 293270 CONSTS(28) = 79035 CONSTS(29) = 62385 CONSTS(30) = 150333 CONSTS(31) = 200240 CONSTS(32) = 127970 CONSTS(33) = 79183 CONSTS(34) = 123020 CONSTS(35) = 121900 CONSTS(36) = -99967 CONSTS(37) = 62555 CONSTS(38) = 224.114 CONSTS(39) = 708.146 CONSTS(40) = 0.00711 CONSTS(41) = 786.217 CONSTS(42) = 1.510 CONSTS(43) = 39.295 CONSTS(44) = 338.915 CONSTS(45) = 529.952 CONSTS(46) = 229.205 CONSTS(47) = 70.078 CONSTS(48) = 225.175 CONSTS(49) = 193.265 CONSTS(50) = -578.317 CONSTS(51) = -130.639 CONSTS(52) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(24)/( CONSTS(5)*CONSTS(6))+CONSTS(38)/CONSTS(5)+( CONSTS(10)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(53) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(25)/( CONSTS(5)*CONSTS(6))+CONSTS(39)/CONSTS(5)+( CONSTS(11)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(54) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(26)/( CONSTS(5)*CONSTS(6))+CONSTS(40)/CONSTS(5)+( CONSTS(23)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(55) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(27)/( CONSTS(5)*CONSTS(6))+CONSTS(41)/CONSTS(5)+( CONSTS(22)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(56) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(28)/( CONSTS(5)*CONSTS(6))+CONSTS(42)/CONSTS(5)+( CONSTS(15)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(57) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(29)/( CONSTS(5)*CONSTS(6))+CONSTS(43)/CONSTS(5)+( CONSTS(14)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(58) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(30)/( CONSTS(5)*CONSTS(6))+CONSTS(44)/CONSTS(5)+( CONSTS(20)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(59) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(31)/( CONSTS(5)*CONSTS(6))+CONSTS(45)/CONSTS(5)+( CONSTS(12)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(60) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(32)/( CONSTS(5)*CONSTS(6))+CONSTS(46)/CONSTS(5)+( CONSTS(13)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(61) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(33)/( CONSTS(5)*CONSTS(6))+CONSTS(47)/CONSTS(5)+( CONSTS(18)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(62) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(34)/( CONSTS(5)*CONSTS(6))+CONSTS(48)/CONSTS(5)+( CONSTS(19)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(63) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(35)/( CONSTS(5)*CONSTS(6))+CONSTS(49)/CONSTS(5)+( CONSTS(21)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(64) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(36)/( CONSTS(5)*CONSTS(6))+CONSTS(50)/CONSTS(5)+( CONSTS(16)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) CONSTS(65) = (( CONSTS(8)*CONSTS(6))/CONSTS(9))*EXP(- CONSTS(37)/( CONSTS(5)*CONSTS(6))+CONSTS(51)/CONSTS(5)+( CONSTS(17)*CONSTS(7)*CONSTS(1))/( CONSTS(5)*CONSTS(6))) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(3) = ( CONSTS(53)*STATES(4)+ CONSTS(54)*STATES(8)) - (CONSTS(55)+ 4.00000*CONSTS(52))*STATES(3) RATES(4) = ( 2.00000*CONSTS(53)*STATES(5)+ 4.00000*CONSTS(52)*STATES(3)+ (CONSTS(54)/CONSTS(4))*STATES(9)) - (CONSTS(53)+ 3.00000*CONSTS(52)+ CONSTS(55)*CONSTS(4))*STATES(4) RATES(5) = ( 3.00000*CONSTS(53)*STATES(6)+ 3.00000*CONSTS(52)*STATES(4)+ (CONSTS(54)/CONSTS(4) ** 2.00000)*STATES(10)) - ( 2.00000*CONSTS(53)+ 2.00000*CONSTS(52)+ CONSTS(55)*CONSTS(4) ** 2.00000)*STATES(5) RATES(6) = ( 4.00000*CONSTS(53)*STATES(7)+ 2.00000*CONSTS(52)*STATES(5)+ (CONSTS(54)/CONSTS(4) ** 3.00000)*STATES(11)) - ( 3.00000*CONSTS(53)+CONSTS(52)+ CONSTS(55)*CONSTS(4) ** 3.00000)*STATES(6) RATES(7) = ( CONSTS(60)*STATES(1)+ CONSTS(63)*STATES(2)+ CONSTS(52)*STATES(6)+ (CONSTS(54)/CONSTS(4) ** 4.00000)*STATES(12)) - ( 4.00000*CONSTS(53)+CONSTS(59)+CONSTS(58)+ CONSTS(55)*CONSTS(4) ** 4.00000)*STATES(7) RATES(1) = ( CONSTS(59)*STATES(7)+ CONSTS(62)*STATES(2)+ CONSTS(56)*STATES(13)) - (CONSTS(60)+CONSTS(61)+CONSTS(57))*STATES(1) RATES(2) = ( CONSTS(61)*STATES(1)+ CONSTS(58)*STATES(7)) - (CONSTS(63)+CONSTS(62))*STATES(2) RATES(8) = ( (CONSTS(53)/CONSTS(4))*STATES(9)+ CONSTS(55)*STATES(3)) - (CONSTS(54)+ 4.00000*CONSTS(52)*CONSTS(4))*STATES(8) RATES(9) = ( (( 2.00000*CONSTS(53))/CONSTS(4))*STATES(10)+ 4.00000*CONSTS(52)*CONSTS(4)*STATES(8)+ CONSTS(55)*CONSTS(4)*STATES(4)) - (CONSTS(53)/CONSTS(4)+ 3.00000*CONSTS(52)*CONSTS(4)+CONSTS(54)/CONSTS(4))*STATES(9) RATES(10) = ( (( 3.00000*CONSTS(53))/CONSTS(4))*STATES(11)+ 3.00000*CONSTS(52)*CONSTS(4)*STATES(9)+ CONSTS(55)*CONSTS(4) ** 2.00000*STATES(5)) - (( 2.00000*CONSTS(53))/CONSTS(4)+ 2.00000*CONSTS(52)*CONSTS(4)+CONSTS(54)/CONSTS(4) ** 2.00000)*STATES(10) RATES(11) = ( (( 4.00000*CONSTS(53))/CONSTS(4))*STATES(12)+ 2.00000*CONSTS(52)*CONSTS(4)*STATES(10)+ CONSTS(55)*CONSTS(4) ** 3.00000*STATES(6)) - (( 3.00000*CONSTS(53))/CONSTS(4)+ CONSTS(52)*CONSTS(4)+CONSTS(54)/CONSTS(4) ** 3.00000)*STATES(11) RATES(12) = ( CONSTS(65)*STATES(13)+ CONSTS(52)*CONSTS(4)*STATES(11)+ CONSTS(55)*CONSTS(4) ** 4.00000*STATES(7)) - (( 4.00000*CONSTS(53))/CONSTS(4)+CONSTS(64)+CONSTS(54)/CONSTS(4) ** 4.00000)*STATES(12) RATES(13) = ( CONSTS(64)*STATES(12)+ CONSTS(57)*STATES(1)) - (CONSTS(65)+CONSTS(56))*STATES(13) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = STATES(1)+STATES(2) ALGBRC(2) = CONSTS(3)*ALGBRC(1)*(CONSTS(1) - CONSTS(2)) RETURN END