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 8 entries in each of the rate and state variable arrays. C There are a total of 17 entries in the constant variable array. C C C VOI is time in component environment (second). C STATES(1) is C1 in component C1 (dimensionless). C CONSTS(1) is C1C2 in component reaction_constants (second_order_rate_constant). C CONSTS(2) is C2C1 in component reaction_constants (first_order_rate_constant). C STATES(2) is C2 in component C2 (dimensionless). C CONSTS(3) is Ca in component reaction_constants (micromolar). C CONSTS(4) is C2C3 in component reaction_constants (second_order_rate_constant). C CONSTS(5) is C3C2 in component reaction_constants (first_order_rate_constant). C CONSTS(6) is C2C5 in component reaction_constants (first_order_rate_constant). C CONSTS(7) is C5C2 in component reaction_constants (first_order_rate_constant). C STATES(3) is C3 in component C3 (dimensionless). C STATES(4) is C5 in component C5 (dimensionless). C CONSTS(8) is O1C3 in component reaction_constants (first_order_rate_constant). C CONSTS(9) is C3O1 in component reaction_constants (first_order_rate_constant). C CONSTS(10) is O2C3 in component reaction_constants (first_order_rate_constant). C CONSTS(11) is C3O2 in component reaction_constants (first_order_rate_constant). C CONSTS(12) is O3C3 in component reaction_constants (first_order_rate_constant). C CONSTS(13) is C3O3 in component reaction_constants (first_order_rate_constant). C STATES(5) is O1 in component O1 (dimensionless). C STATES(6) is O2 in component O2 (dimensionless). C STATES(7) is O3 in component O3 (dimensionless). C STATES(8) is C4 in component C4 (dimensionless). C CONSTS(14) is O2C4 in component reaction_constants (first_order_rate_constant). C CONSTS(15) is C4O2 in component reaction_constants (first_order_rate_constant). C CONSTS(16) is O3C4 in component reaction_constants (first_order_rate_constant). C CONSTS(17) is C4O3 in component reaction_constants (first_order_rate_constant). C RATES(1) is d/dt C1 in component C1 (dimensionless). C RATES(2) is d/dt C2 in component C2 (dimensionless). C RATES(3) is d/dt C3 in component C3 (dimensionless). C RATES(8) is d/dt C4 in component C4 (dimensionless). C RATES(4) is d/dt C5 in component C5 (dimensionless). C RATES(5) is d/dt O1 in component O1 (dimensionless). C RATES(6) is d/dt O2 in component O2 (dimensionless). C RATES(7) is d/dt O3 in component O3 (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.125 CONSTS(1) = 2.5 CONSTS(2) = 13.3 STATES(2) = 0.125 CONSTS(3) = 20.0 CONSTS(4) = 68.0 CONSTS(5) = 8000.0 CONSTS(6) = 0.13 CONSTS(7) = 3.6 STATES(3) = 0.125 STATES(4) = 0.125 CONSTS(8) = 3400.0 CONSTS(9) = 1100.0 CONSTS(10) = 92.0 CONSTS(11) = 17.0 CONSTS(12) = 138.0 CONSTS(13) = 14.0 STATES(5) = 0.125 STATES(6) = 0.125 STATES(7) = 0.125 STATES(8) = 0.125 CONSTS(14) = 1900.0 CONSTS(15) = 520.0 CONSTS(16) = 300.0 CONSTS(17) = 46.0 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = CONSTS(2)*STATES(2) - CONSTS(1)*CONSTS(3)*STATES(2) RATES(2) = ( CONSTS(1)*CONSTS(3)*STATES(1)+ CONSTS(5)*STATES(3)+ CONSTS(7)*STATES(4)) - ( CONSTS(2)*STATES(2)+ CONSTS(4)*CONSTS(3)*STATES(2)+ CONSTS(6)*STATES(2)) RATES(3) = ( CONSTS(4)*CONSTS(3)*STATES(2)+ CONSTS(8)*STATES(5)+ CONSTS(12)*STATES(7)+ CONSTS(10)*STATES(6)) - ( CONSTS(5)*STATES(3)+ CONSTS(9)*STATES(3)+ CONSTS(11)*STATES(3)+ CONSTS(13)*STATES(3)) RATES(8) = ( CONSTS(14)*STATES(6)+ CONSTS(16)*STATES(7)) - ( CONSTS(15)*STATES(8)+ CONSTS(17)*STATES(8)) RATES(4) = CONSTS(6)*STATES(2) - CONSTS(7)*STATES(4) RATES(5) = CONSTS(9)*STATES(3) - CONSTS(8)*STATES(5) RATES(6) = ( CONSTS(11)*STATES(3)+ CONSTS(15)*STATES(8)) - ( CONSTS(10)*STATES(6)+ CONSTS(14)*STATES(6)) RATES(7) = ( CONSTS(13)*STATES(3)+ CONSTS(17)*STATES(8)) - ( CONSTS(12)*STATES(7)+ CONSTS(16)*STATES(7)) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RETURN END