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 7 entries in each of the rate and state variable arrays. C There are a total of 11 entries in the constant variable array. C C C VOI is time in component environment (hour). C STATES(1) is T in component T (cell_ml). C CONSTS(1) is k1 in component model_parameters (cell_per_ml_molecule2_hour). C STATES(2) is Lm in component Lm (molecule_cell). C STATES(3) is RT in component RT (molecule_cell). C STATES(4) is m in component m (cell_ml). C STATES(5) is LT in component LT (molecule_cell). C CONSTS(2) is k2 in component model_parameters (molecule_ml2_per_cell3_hour). C CONSTS(3) is k3 in component model_parameters (per_hour). C CONSTS(4) is k4 in component model_parameters (per_hour). C CONSTS(5) is k5 in component model_parameters (ml2_per_cell_molecule_hour). C STATES(6) is Rm in component Rm (molecule_cell). C CONSTS(6) is k6 in component model_parameters (molecule_per_cell_hour). C CONSTS(7) is k7 in component model_parameters (per_hour). C CONSTS(8) is k8 in component model_parameters (per_molar_per_minute). C STATES(7) is SL in component SL (molecule_ml). C CONSTS(9) is k10 in component model_parameters (molecule_per_cell_hour). C CONSTS(10) is k11 in component model_parameters (per_hour). C CONSTS(11) is k9 in component model_parameters (per_molar_per_minute). C RATES(1) is d/dt T in component T (cell_ml). C RATES(5) is d/dt LT in component LT (molecule_cell). C RATES(3) is d/dt RT in component RT (molecule_cell). C RATES(4) is d/dt m in component m (cell_ml). C RATES(2) is d/dt Lm in component Lm (molecule_cell). C RATES(6) is d/dt Rm in component Rm (molecule_cell). C RATES(7) is d/dt SL in component SL (molecule_ml). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 500.0 CONSTS(1) = 8.38E-10 STATES(2) = 1E3 STATES(3) = 1E4 STATES(4) = 500.0 STATES(5) = 0.0 CONSTS(2) = 6E-3 CONSTS(3) = 5.9413 CONSTS(4) = 0.35 CONSTS(5) = 2.52E-9 STATES(6) = 1E3 CONSTS(6) = 2.244E3 CONSTS(7) = 0.35 CONSTS(8) = 1.92E10 STATES(7) = 0.0 CONSTS(9) = 3.11E3 CONSTS(10) = 13.9 CONSTS(11) = 87.3E8 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = - ( 1.00000*CONSTS(1)*STATES(4)*STATES(1)*STATES(2)*STATES(3)) RATES(5) = CONSTS(2)*STATES(1)*STATES(4) - ( CONSTS(3)*STATES(5)+ CONSTS(4)*STATES(5)+ CONSTS(5)*STATES(4)*STATES(1)*STATES(5)*STATES(6)) RATES(3) = CONSTS(6) - ( CONSTS(7)*STATES(3)+ 1.00000*CONSTS(8)*STATES(3)*STATES(7)) RATES(4) = - ( 1.00000*CONSTS(1)*STATES(4)*STATES(1)*STATES(5)*STATES(6)) RATES(2) = CONSTS(9) - ( CONSTS(3)*STATES(2)+ CONSTS(4)*STATES(2)+ CONSTS(5)*STATES(4)*STATES(1)*STATES(2)*STATES(3)) RATES(6) = CONSTS(6) - ( CONSTS(7)*STATES(6)+ 1.00000*CONSTS(8)*STATES(6)*STATES(7)) RATES(7) = CONSTS(3)*( STATES(2)*STATES(4)+ STATES(5)*STATES(1)) - ( CONSTS(10)*STATES(7)+ 1.00000*CONSTS(11)*( STATES(6)*STATES(4)+ STATES(3)*STATES(1))*STATES(7)) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RETURN END