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
