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 0 entries in the algebraic variable array.
C There are a total of 12 entries in each of the rate and state variable arrays.
C There are a total of 38 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is E in component E (millimolar).
C CONSTS(1) is k23 in component reaction_constants (second_order_rate_constant).
C CONSTS(2) is k24 in component reaction_constants (first_order_rate_constant).
C CONSTS(3) is k17 in component reaction_constants (first_order_rate_constant).
C CONSTS(4) is k18 in component reaction_constants (first_order_rate_constant).
C CONSTS(5) is k1 in component reaction_constants (second_order_rate_constant).
C CONSTS(6) is k2 in component reaction_constants (first_order_rate_constant).
C CONSTS(7) is k3 in component reaction_constants (second_order_rate_constant).
C CONSTS(8) is k4 in component reaction_constants (first_order_rate_constant).
C STATES(2) is ED in component ED (millimolar).
C STATES(3) is E_ in component E_ (millimolar).
C STATES(4) is ENa in component ENa (millimolar).
C STATES(5) is ECl in component ECl (millimolar).
C CONSTS(9) is D in component reaction_constants (millimolar).
C CONSTS(10) is Na in component reaction_constants (millimolar).
C CONSTS(11) is Cl in component reaction_constants (millimolar).
C CONSTS(12) is k29 in component reaction_constants (second_order_rate_constant).
C CONSTS(13) is k30 in component reaction_constants (first_order_rate_constant).
C CONSTS(14) is k11 in component reaction_constants (second_order_rate_constant).
C CONSTS(15) is k12 in component reaction_constants (first_order_rate_constant).
C CONSTS(16) is k9 in component reaction_constants (second_order_rate_constant).
C CONSTS(17) is k10 in component reaction_constants (first_order_rate_constant).
C STATES(6) is ED_ in component ED_ (millimolar).
C STATES(7) is ENa_ in component ENa_ (millimolar).
C STATES(8) is ECl_ in component ECl_ (millimolar).
C CONSTS(18) is D_ in component reaction_constants (millimolar).
C CONSTS(19) is Na_ in component reaction_constants (millimolar).
C CONSTS(20) is Cl_ in component reaction_constants (millimolar).
C CONSTS(21) is k21 in component reaction_constants (second_order_rate_constant).
C CONSTS(22) is k22 in component reaction_constants (first_order_rate_constant).
C STATES(9) is ENaD in component ENaD (millimolar).
C CONSTS(23) is k27 in component reaction_constants (second_order_rate_constant).
C CONSTS(24) is k28 in component reaction_constants (first_order_rate_constant).
C STATES(10) is ENaD_ in component ENaD_ (millimolar).
C CONSTS(25) is k25 in component reaction_constants (second_order_rate_constant).
C CONSTS(26) is k26 in component reaction_constants (first_order_rate_constant).
C CONSTS(27) is k31 in component reaction_constants (second_order_rate_constant).
C CONSTS(28) is k32 in component reaction_constants (first_order_rate_constant).
C CONSTS(29) is k5 in component reaction_constants (second_order_rate_constant).
C CONSTS(30) is k6 in component reaction_constants (first_order_rate_constant).
C STATES(11) is ENaCl in component ENaCl (millimolar).
C CONSTS(31) is k13 in component reaction_constants (second_order_rate_constant).
C CONSTS(32) is k14 in component reaction_constants (first_order_rate_constant).
C STATES(12) is ENaCl_ in component ENaCl_ (millimolar).
C CONSTS(33) is k7 in component reaction_constants (second_order_rate_constant).
C CONSTS(34) is k8 in component reaction_constants (first_order_rate_constant).
C CONSTS(35) is k15 in component reaction_constants (second_order_rate_constant).
C CONSTS(36) is k16 in component reaction_constants (first_order_rate_constant).
C CONSTS(37) is k19 in component reaction_constants (first_order_rate_constant).
C CONSTS(38) is k20 in component reaction_constants (first_order_rate_constant).
C RATES(1) is d/dt E in component E (millimolar).
C RATES(3) is d/dt E_ in component E_ (millimolar).
C RATES(2) is d/dt ED in component ED (millimolar).
C RATES(6) is d/dt ED_ in component ED_ (millimolar).
C RATES(9) is d/dt ENaD in component ENaD (millimolar).
C RATES(10) is d/dt ENaD_ in component ENaD_ (millimolar).
C RATES(4) is d/dt ENa in component ENa (millimolar).
C RATES(7) is d/dt ENa_ in component ENa_ (millimolar).
C RATES(5) is d/dt ECl in component ECl (millimolar).
C RATES(8) is d/dt ECl_ in component ECl_ (millimolar).
C RATES(11) is d/dt ENaCl in component ENaCl (millimolar).
C RATES(12) is d/dt ENaCl_ in component ENaCl_ (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.08333
      CONSTS(1) = 1.0E5
      CONSTS(2) = 3.192E1
      CONSTS(3) = 4.587E5
      CONSTS(4) = 1.0E5
      CONSTS(5) = 1.0E5
      CONSTS(6) = 4.183E5
      CONSTS(7) = 1.0E5
      CONSTS(8) = 4.928E6
      STATES(2) = 0.08333
      STATES(3) = 0.08333
      STATES(4) = 0.08333
      STATES(5) = 0.08333
      CONSTS(9) = 1.0E-6
      CONSTS(10) = 50.0
      CONSTS(11) = 96.0
      CONSTS(12) = 1.0E5
      CONSTS(13) = 3.514E-1
      CONSTS(14) = 1.0E5
      CONSTS(15) = 4.982E6
      CONSTS(16) = 1.0E5
      CONSTS(17) = 4.183E5
      STATES(6) = 0.08333
      STATES(7) = 0.08333
      STATES(8) = 0.08333
      CONSTS(18) = 1.0E-6
      CONSTS(19) = 10.0
      CONSTS(20) = 40.0
      CONSTS(21) = 1.0E5
      CONSTS(22) = 4.183E5
      STATES(9) = 0.08333
      CONSTS(23) = 1.0E5
      CONSTS(24) = 1.389E5
      STATES(10) = 0.08333
      CONSTS(25) = 1.0E5
      CONSTS(26) = 3.192E1
      CONSTS(27) = 1.0E5
      CONSTS(28) = 1.166E-1
      CONSTS(29) = 1.0E5
      CONSTS(30) = 1.065E6
      STATES(11) = 0.08333
      CONSTS(31) = 1.0E5
      CONSTS(32) = 1.065E6
      STATES(12) = 0.08333
      CONSTS(33) = 1.0E5
      CONSTS(34) = 8.940E4
      CONSTS(35) = 1.0E5
      CONSTS(36) = 8.940E4
      CONSTS(37) = 1.0E3
      CONSTS(38) = 2.180E2
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(2)*STATES(2)+ CONSTS(4)*STATES(3)+ CONSTS(6)*STATES(4)+ CONSTS(8)*STATES(5)) - ( CONSTS(1)*CONSTS(9)*STATES(1)+ CONSTS(3)*STATES(1)+ CONSTS(5)*CONSTS(10)*STATES(1)+ CONSTS(7)*CONSTS(11)*STATES(1))
      RATES(3) = ( CONSTS(13)*STATES(6)+ CONSTS(3)*STATES(1)+ CONSTS(17)*STATES(7)+ CONSTS(15)*STATES(8)) - ( CONSTS(12)*CONSTS(18)*STATES(3)+ CONSTS(4)*STATES(3)+ CONSTS(16)*CONSTS(19)*STATES(3)+ CONSTS(14)*CONSTS(20)*STATES(3))
      RATES(2) = ( CONSTS(1)*STATES(1)*CONSTS(9)+ CONSTS(22)*STATES(9)) - ( CONSTS(2)*STATES(2)+ CONSTS(21)*CONSTS(10)*STATES(2))
      RATES(6) = ( CONSTS(12)*STATES(3)*CONSTS(18)+ CONSTS(24)*STATES(10)) - ( CONSTS(13)*STATES(6)+ CONSTS(23)*CONSTS(19)*STATES(6))
      RATES(9) = ( CONSTS(21)*CONSTS(10)*STATES(2)+ CONSTS(25)*CONSTS(9)*STATES(4)) - ( CONSTS(22)*STATES(9)+ CONSTS(26)*STATES(9))
      RATES(10) = ( CONSTS(23)*CONSTS(19)*STATES(6)+ CONSTS(27)*CONSTS(18)*STATES(7)) - ( CONSTS(24)*STATES(10)+ CONSTS(28)*STATES(10))
      RATES(4) = ( CONSTS(5)*CONSTS(10)*STATES(1)+ CONSTS(26)*STATES(9)+ CONSTS(30)*STATES(11)) - ( CONSTS(6)*STATES(4)+ CONSTS(29)*CONSTS(11)*STATES(4)+ CONSTS(25)*CONSTS(9)*STATES(4))
      RATES(7) = ( CONSTS(16)*CONSTS(19)*STATES(3)+ CONSTS(28)*STATES(10)+ CONSTS(32)*STATES(12)) - ( CONSTS(17)*STATES(7)+ CONSTS(31)*CONSTS(20)*STATES(7)+ CONSTS(27)*CONSTS(18)*STATES(7))
      RATES(5) = ( CONSTS(7)*CONSTS(11)*STATES(1)+ CONSTS(34)*STATES(11)) - ( CONSTS(33)*CONSTS(10)*STATES(5)+ CONSTS(8)*STATES(5))
      RATES(8) = ( CONSTS(14)*CONSTS(20)*STATES(3)+ CONSTS(36)*STATES(12)) - ( CONSTS(35)*CONSTS(19)*STATES(8)+ CONSTS(15)*STATES(8))
      RATES(11) = ( CONSTS(29)*CONSTS(11)*STATES(4)+ CONSTS(33)*CONSTS(10)*STATES(5)+ CONSTS(38)*STATES(12)) - ( CONSTS(30)*STATES(11)+ CONSTS(34)*STATES(11)+ CONSTS(37)*STATES(11))
      RATES(12) = ( CONSTS(31)*CONSTS(20)*STATES(7)+ CONSTS(35)*CONSTS(19)*STATES(8)+ CONSTS(37)*STATES(11)) - ( CONSTS(32)*STATES(12)+ CONSTS(36)*STATES(12)+ CONSTS(38)*STATES(12))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END
Source
Derived from workspace Chang, Fujita, 1999 at changeset 43068d996765.
This exposure was expired.
Collaboration
To begin collaborating on this work, please use your git client and issue this command:
License
The terms of use/license for this work is unspecified.