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 9 entries in each of the rate and state variable arrays.
C There are a total of 40 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is pRB in component pRB (CONC).
C CONSTS(1) is k1 in component model_parameters (dimensionless).
C STATES(2) is E2F1 in component E2F1 (CONC).
C CONSTS(2) is km1 in component model_parameters (dimensionless).
C CONSTS(3) is j11 in component model_parameters (dimensionless).
C CONSTS(4) is j61 in component model_parameters (dimensionless).
C CONSTS(5) is k16 in component model_parameters (dimensionless).
C STATES(3) is CycDa in component CycDa (CONC).
C CONSTS(6) is k61 in component model_parameters (dimensionless).
C STATES(4) is pRBp in component pRBp (CONC).
C CONSTS(7) is phi_pRB in component model_parameters (dimensionless).
C CONSTS(8) is kp in component model_parameters (dimensionless).
C CONSTS(9) is k2 in component model_parameters (dimensionless).
C CONSTS(10) is a in component model_parameters (dimensionless).
C CONSTS(11) is km2 in component model_parameters (dimensionless).
C CONSTS(12) is j12 in component model_parameters (dimensionless).
C CONSTS(13) is j62 in component model_parameters (dimensionless).
C CONSTS(14) is phi_E2F1 in component model_parameters (dimensionless).
C CONSTS(15) is k34 in component model_parameters (dimensionless).
C STATES(5) is CycDi in component CycDi (CONC).
C CONSTS(16) is km4 in component model_parameters (dimensionless).
C CONSTS(17) is k43 in component model_parameters (dimensionless).
C CONSTS(18) is phi_CycDa in component model_parameters (dimensionless).
C CONSTS(19) is k3 in component model_parameters (dimensionless).
C STATES(6) is AP1 in component AP1 (CONC).
C CONSTS(20) is k23 in component model_parameters (dimensionless).
C CONSTS(21) is j13 in component model_parameters (dimensionless).
C CONSTS(22) is j63 in component model_parameters (dimensionless).
C CONSTS(23) is phi_CycDi in component model_parameters (dimensionless).
C CONSTS(24) is Fm in component model_parameters (dimensionless).
C CONSTS(25) is k25 in component model_parameters (dimensionless).
C CONSTS(26) is j15 in component model_parameters (dimensionless).
C CONSTS(27) is j65 in component model_parameters (dimensionless).
C CONSTS(28) is phi_AP1 in component model_parameters (dimensionless).
C CONSTS(29) is k67 in component model_parameters (dimensionless).
C STATES(7) is CycEa in component CycEa (CONC).
C CONSTS(30) is k76 in component model_parameters (dimensionless).
C STATES(8) is pRBpp in component pRBpp (CONC).
C CONSTS(31) is phi_pRBp in component model_parameters (dimensionless).
C CONSTS(32) is phi_pRBpp in component model_parameters (dimensionless).
C STATES(9) is CycEi in component CycEi (CONC).
C CONSTS(33) is k28 in component model_parameters (dimensionless).
C CONSTS(34) is j18 in component model_parameters (dimensionless).
C CONSTS(35) is j68 in component model_parameters (dimensionless).
C CONSTS(36) is k98 in component model_parameters (dimensionless).
C CONSTS(37) is k89 in component model_parameters (dimensionless).
C CONSTS(38) is km9 in component model_parameters (dimensionless).
C CONSTS(39) is phi_CycEi in component model_parameters (dimensionless).
C CONSTS(40) is phi_CycEa in component model_parameters (dimensionless).
C RATES(1) is d/dt pRB in component pRB (CONC).
C RATES(2) is d/dt E2F1 in component E2F1 (CONC).
C RATES(3) is d/dt CycDa in component CycDa (CONC).
C RATES(5) is d/dt CycDi in component CycDi (CONC).
C RATES(6) is d/dt AP1 in component AP1 (CONC).
C RATES(4) is d/dt pRBp in component pRBp (CONC).
C RATES(8) is d/dt pRBpp in component pRBpp (CONC).
C RATES(9) is d/dt CycEi in component CycEi (CONC).
C RATES(7) is d/dt CycEa in component CycEa (CONC).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0
CONSTS(1) = 1
STATES(2) = 0
CONSTS(2) = 0.5
CONSTS(3) = 0.5
CONSTS(4) = 5
CONSTS(5) = 0.4
STATES(3) = 0
CONSTS(6) = 0.3
STATES(4) = 0
CONSTS(7) = 0.005
CONSTS(8) = 0.05
CONSTS(9) = 1.6
CONSTS(10) = 0.04
CONSTS(11) = 4
CONSTS(12) = 5
CONSTS(13) = 8
CONSTS(14) = 0.1
CONSTS(15) = 0.04
STATES(5) = 0
CONSTS(16) = 0.3
CONSTS(17) = 0.01
CONSTS(18) = 0.03
CONSTS(19) = 0.05
STATES(6) = 0
CONSTS(20) = 0.3
CONSTS(21) = 0.002
CONSTS(22) = 2
CONSTS(23) = 0.023
CONSTS(24) = 0
CONSTS(25) = 0.9
CONSTS(26) = 0.001
CONSTS(27) = 6
CONSTS(28) = 0.01
CONSTS(29) = 0.7
STATES(7) = 0
CONSTS(30) = 0.1
STATES(8) = 0
CONSTS(31) = 0.06
CONSTS(32) = 0.04
STATES(9) = 0
CONSTS(33) = 0.06
CONSTS(34) = 0.6
CONSTS(35) = 7
CONSTS(36) = 0.01
CONSTS(37) = 0.07
CONSTS(38) = 0.005
CONSTS(39) = 0.06
CONSTS(40) = 0.05
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = (( CONSTS(1)*(STATES(2)/(CONSTS(2)+STATES(2)))*(CONSTS(3)/( CONSTS(3)*STATES(1)))*(CONSTS(4)/( CONSTS(4)*STATES(4))) - CONSTS(5)*STATES(1)*STATES(3))+ CONSTS(6)*STATES(4)) - CONSTS(7)*STATES(1)
RATES(2) = (CONSTS(8)+ CONSTS(9)*((CONSTS(10) ** 2.00000+STATES(2) ** 2.00000)/(CONSTS(11) ** 2.00000+STATES(2) ** 2.00000))*(CONSTS(12)/(CONSTS(12)+STATES(4)))*(CONSTS(13)/(CONSTS(13)+STATES(4)))) - CONSTS(14)*STATES(2)
RATES(3) = (( CONSTS(15)*STATES(5)*STATES(3))/(CONSTS(16)+STATES(3)) - CONSTS(15)*STATES(3)) - CONSTS(18)*STATES(3)
RATES(5) = (( CONSTS(19)*STATES(6)+ CONSTS(20)*STATES(2)*(CONSTS(21)/(CONSTS(21)+STATES(1)))*(CONSTS(22)/(CONSTS(22)+STATES(4)))+ CONSTS(17)*STATES(3)) - CONSTS(15)*STATES(5)*(STATES(3)/(CONSTS(16)+STATES(3)))) - CONSTS(23)*STATES(5)
RATES(6) = (CONSTS(24)+ CONSTS(25)*STATES(2)*(CONSTS(26)/(CONSTS(26)+STATES(1)))*(CONSTS(27)/(CONSTS(27)+STATES(4)))) - CONSTS(28)*STATES(6)
RATES(4) = ((( CONSTS(5)*STATES(1)*STATES(3) - CONSTS(6)*STATES(4)) - CONSTS(29)*STATES(4)*STATES(7))+ CONSTS(30)*STATES(8)) - CONSTS(31)*STATES(4)
RATES(8) = ( CONSTS(29)*STATES(4)*STATES(7) - CONSTS(30)*STATES(8)) - CONSTS(32)*STATES(8)
RATES(9) = (( CONSTS(33)*STATES(2)*(CONSTS(34)/(CONSTS(34)+STATES(1)))*(CONSTS(35)/(CONSTS(35)+STATES(4)))+ CONSTS(36)*STATES(7)) - CONSTS(37)*STATES(9)*(STATES(7)/(CONSTS(38)+STATES(7)))) - CONSTS(39)*STATES(9)
RATES(7) = ( CONSTS(37)*STATES(9)*(STATES(7)/(CONSTS(38)+STATES(7))) - CONSTS(36)*STATES(7)) - CONSTS(40)*STATES(7)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
