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