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