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 8 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 35 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Y in component Y (dimensionless).
C CONSTS(1) is k1 in component Y (first_order_rate_constant).
C ALGBRC(7) is k2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is k3 in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(1) is C in component C (dimensionless).
C STATES(2) is M in component M (dimensionless).
C ALGBRC(8) is kwee in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(6) is kcdc25 in component kinetic_parameters (first_order_rate_constant).
C STATES(3) is preMPF in component preMPF (dimensionless).
C STATES(4) is Cdc25P in component Cdc25P (dimensionless).
C CONSTS(3) is k25 in component Cdc25P (first_order_rate_constant).
C CONSTS(4) is Km25 in component Cdc25P (dimensionless).
C CONSTS(5) is k25r in component Cdc25P (first_order_rate_constant).
C CONSTS(34) is k25ro in component Cdc25P (first_order_rate_constant).
C CONSTS(6) is Km25r in component Cdc25P (dimensionless).
C CONSTS(33) is p in component p (dimensionless).
C ALGBRC(3) is Cdc25 in component Cdc25 (dimensionless).
C STATES(5) is Wee1 in component Wee1 (dimensionless).
C CONSTS(7) is kw in component Wee1 (first_order_rate_constant).
C CONSTS(8) is Kmw in component Wee1 (dimensionless).
C CONSTS(9) is kwr in component Wee1 (first_order_rate_constant).
C CONSTS(35) is kwro in component Wee1 (first_order_rate_constant).
C CONSTS(10) is Kmwr in component Wee1 (dimensionless).
C ALGBRC(4) is Wee1P in component Wee1P (dimensionless).
C STATES(6) is IEP in component IEP (dimensionless).
C CONSTS(11) is kie in component IEP (first_order_rate_constant).
C CONSTS(12) is Kmie in component IEP (dimensionless).
C CONSTS(13) is kier in component IEP (first_order_rate_constant).
C CONSTS(14) is Kmier in component IEP (dimensionless).
C ALGBRC(2) is IE in component IE (dimensionless).
C STATES(7) is APC_ in component APC_ (dimensionless).
C CONSTS(15) is kap in component APC_ (first_order_rate_constant).
C CONSTS(16) is Kmap in component APC_ (dimensionless).
C CONSTS(17) is kapr in component APC_ (first_order_rate_constant).
C CONSTS(18) is Kmapr in component APC_ (dimensionless).
C ALGBRC(5) is APC in component APC (dimensionless).
C CONSTS(19) is Cdc2tot in component C (dimensionless).
C CONSTS(20) is Cdc25tot in component Cdc25 (dimensionless).
C CONSTS(21) is Wee1tot in component Wee1P (dimensionless).
C CONSTS(22) is IEtot in component IE (dimensionless).
C CONSTS(23) is APCtot in component APC (dimensionless).
C CONSTS(24) is R in component p (dimensionless).
C CONSTS(25) is s in component p (dimensionless).
C CONSTS(26) is q in component p (dimensionless).
C CONSTS(27) is V2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(28) is V2_ in component kinetic_parameters (first_order_rate_constant).
C CONSTS(29) is V25 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(30) is V25_ in component kinetic_parameters (first_order_rate_constant).
C CONSTS(31) is Vwee in component kinetic_parameters (first_order_rate_constant).
C CONSTS(32) is Vwee_ in component kinetic_parameters (first_order_rate_constant).
C RATES(1) is d/dt Y in component Y (dimensionless).
C RATES(2) is d/dt M in component M (dimensionless).
C RATES(3) is d/dt preMPF in component preMPF (dimensionless).
C RATES(4) is d/dt Cdc25P in component Cdc25P (dimensionless).
C RATES(5) is d/dt Wee1 in component Wee1 (dimensionless).
C RATES(6) is d/dt IEP in component IEP (dimensionless).
C RATES(7) is d/dt APC_ in component APC_ (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.01577552
CONSTS(1) = 0.024
CONSTS(2) = 1.15
STATES(2) = 0.0142077
STATES(3) = 0.02541723
STATES(4) = 0.4844362
CONSTS(3) = 18
CONSTS(4) = 0.1
CONSTS(5) = 0.8
CONSTS(6) = 1
STATES(5) = 0.5155638
CONSTS(7) = 18
CONSTS(8) = 0.1
CONSTS(9) = 0.8
CONSTS(10) = 1
STATES(6) = 0.002287817
CONSTS(11) = 4.5
CONSTS(12) = 0.01
CONSTS(13) = 0.34
CONSTS(14) = 0.01
STATES(7) = 0.5051103
CONSTS(15) = 0.3
CONSTS(16) = 0.01
CONSTS(17) = 0.3
CONSTS(18) = 1
CONSTS(19) = 1
CONSTS(20) = 1
CONSTS(21) = 1
CONSTS(22) = 1
CONSTS(23) = 1
CONSTS(24) = 1
CONSTS(25) = 0.021
CONSTS(26) = 1
CONSTS(27) = 0.01
CONSTS(28) = 0.6
CONSTS(29) = 0.04
CONSTS(30) = 0.4
CONSTS(31) = 0.025
CONSTS(32) = 2.5
CONSTS(33) = CONSTS(25)*CONSTS(24)+CONSTS(26)
CONSTS(34) = CONSTS(5)/CONSTS(33)
CONSTS(35) = CONSTS(9)/CONSTS(33)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = CONSTS(22) - STATES(6)
RATES(6) = ( CONSTS(11)*STATES(2)*ALGBRC(2))/(CONSTS(12)+ALGBRC(2)) - ( CONSTS(13)*STATES(6))/(CONSTS(14)+STATES(6))
ALGBRC(3) = CONSTS(20) - STATES(4)
RATES(4) = ( CONSTS(3)*STATES(2)*ALGBRC(3))/(CONSTS(4)+ALGBRC(3)) - ( CONSTS(5)*STATES(4))/(CONSTS(6)+STATES(4))
ALGBRC(4) = CONSTS(21) - STATES(5)
RATES(5) = ( - CONSTS(7)*STATES(2)*STATES(5))/(CONSTS(8)+STATES(5))+( CONSTS(9)*ALGBRC(4))/(CONSTS(10)+ALGBRC(4))
ALGBRC(5) = CONSTS(23) - STATES(7)
RATES(7) = ( CONSTS(15)*STATES(6)*ALGBRC(5))/(CONSTS(16)+ALGBRC(5)) - ( CONSTS(17)*STATES(7))/(CONSTS(18)+STATES(7))
ALGBRC(7) = CONSTS(27)*ALGBRC(5)+ CONSTS(28)*STATES(7)
ALGBRC(1) = CONSTS(19) - (STATES(2)+STATES(3))
RATES(1) = CONSTS(1) - ( ALGBRC(7)*STATES(1)+ CONSTS(2)*STATES(1)*ALGBRC(1))
ALGBRC(8) = CONSTS(31)*ALGBRC(4)+ CONSTS(32)*STATES(5)
ALGBRC(6) = CONSTS(29)*ALGBRC(3)+ CONSTS(30)*STATES(4)
RATES(2) = ( ALGBRC(6)*STATES(3)+ CONSTS(2)*STATES(1)*ALGBRC(1)) - ( ALGBRC(7)*STATES(2)+ ALGBRC(8)*STATES(2))
RATES(3) = ALGBRC(8)*STATES(2) - ( ALGBRC(7)*STATES(3)+ ALGBRC(6)*STATES(3))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = CONSTS(22) - STATES(6)
ALGBRC(3) = CONSTS(20) - STATES(4)
ALGBRC(4) = CONSTS(21) - STATES(5)
ALGBRC(5) = CONSTS(23) - STATES(7)
ALGBRC(7) = CONSTS(27)*ALGBRC(5)+ CONSTS(28)*STATES(7)
ALGBRC(1) = CONSTS(19) - (STATES(2)+STATES(3))
ALGBRC(8) = CONSTS(31)*ALGBRC(4)+ CONSTS(32)*STATES(5)
ALGBRC(6) = CONSTS(29)*ALGBRC(3)+ CONSTS(30)*STATES(4)
RETURN
END
