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