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 6 entries in the algebraic variable array. C There are a total of 13 entries in each of the rate and state variable arrays. C There are a total of 23 entries in the constant variable array. C C C VOI is time in component environment (minute). C STATES(1) is Cdk2_CycE in component Cdk2_CycE (dimensionless). C CONSTS(1) is kwee in component kinetic_parameters (first_order_rate_constant). C CONSTS(2) is kon in component kinetic_parameters (first_order_rate_constant). C CONSTS(3) is k25A in component kinetic_parameters (first_order_rate_constant). C CONSTS(4) is koff in component kinetic_parameters (first_order_rate_constant). C CONSTS(5) is kassoc in component kinetic_parameters (first_order_rate_constant). C CONSTS(6) is kdissoc in component kinetic_parameters (first_order_rate_constant). C ALGBRC(5) is phi in component kinetic_parameters (dimensionless). C STATES(2) is Wee1_a in component Wee1_a (dimensionless). C STATES(3) is PCdk2_CycE in component PCdk2_CycE (dimensionless). C STATES(4) is Xic in component Xic (dimensionless). C STATES(5) is Cdk2_CycErem in component Cdk2_CycErem (dimensionless). C STATES(6) is Xic_Cdk2_CycE in component Xic_Cdk2_CycE (dimensionless). C STATES(7) is PCdk2_CycErem in component PCdk2_CycErem (dimensionless). C STATES(8) is Xic_PCdk2_CycE in component Xic_PCdk2_CycE (dimensionless). C CONSTS(7) is Jwact in component Wee1_a (dimensionless). C CONSTS(8) is Jwinact in component Wee1_a (dimensionless). C CONSTS(9) is Wee1_total in component Wee1_a (dimensionless). C CONSTS(10) is kwact in component Wee1_a (first_order_rate_constant). C CONSTS(11) is kwinact in component Wee1_a (first_order_rate_constant). C STATES(9) is Kin_a in component Kin_a (dimensionless). C CONSTS(12) is Jiact in component Kin_a (dimensionless). C CONSTS(13) is Jiinact in component Kin_a (dimensionless). C CONSTS(14) is kiact in component Kin_a (first_order_rate_constant). C CONSTS(15) is kiinact in component Kin_a (first_order_rate_constant). C CONSTS(16) is kedeg in component kinetic_parameters (first_order_rate_constant). C CONSTS(17) is kxdeg in component kinetic_parameters (first_order_rate_constant). C STATES(10) is Deg_a in component Deg_a (dimensionless). C STATES(11) is Xic_Cdk2_CycErem in component Xic_Cdk2_CycErem (dimensionless). C STATES(12) is Xic_PCdk2_CycErem in component Xic_PCdk2_CycErem (dimensionless). C ALGBRC(6) is Heav in component Deg_a (dimensionless). C CONSTS(18) is kdact in component Deg_a (first_order_rate_constant). C CONSTS(19) is theta in component Deg_a (dimensionless). C ALGBRC(2) is x in component Deg_a (dimensionless). C CONSTS(20) is Deg in component kinetic_parameters (dimensionless). C STATES(13) is Xicrem in component Xicrem (dimensionless). C ALGBRC(3) is Cyc_total in component Cyc_total (dimensionless). C ALGBRC(4) is Xic_total in component Xic_total (dimensionless). C CONSTS(21) is epsilon in component kinetic_parameters (dimensionless). C ALGBRC(1) is pool in component kinetic_parameters (dimensionless). C CONSTS(22) is n in component kinetic_parameters (dimensionless). C CONSTS(23) is L in component kinetic_parameters (dimensionless). C RATES(1) is d/dt Cdk2_CycE in component Cdk2_CycE (dimensionless). C RATES(3) is d/dt PCdk2_CycE in component PCdk2_CycE (dimensionless). C RATES(2) is d/dt Wee1_a in component Wee1_a (dimensionless). C RATES(9) is d/dt Kin_a in component Kin_a (dimensionless). C RATES(5) is d/dt Cdk2_CycErem in component Cdk2_CycErem (dimensionless). C RATES(7) is d/dt PCdk2_CycErem in component PCdk2_CycErem (dimensionless). C RATES(10) is d/dt Deg_a in component Deg_a (dimensionless). C RATES(4) is d/dt Xic in component Xic (dimensionless). C RATES(6) is d/dt Xic_Cdk2_CycE in component Xic_Cdk2_CycE (dimensionless). C RATES(8) is d/dt Xic_PCdk2_CycE in component Xic_PCdk2_CycE (dimensionless). C RATES(11) is d/dt Xic_Cdk2_CycErem in component Xic_Cdk2_CycErem (dimensionless). C RATES(12) is d/dt Xic_PCdk2_CycErem in component Xic_PCdk2_CycErem (dimensionless). C RATES(13) is d/dt Xicrem in component Xicrem (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.06 CONSTS(1) = 1.5 CONSTS(2) = 0.02 CONSTS(3) = 0.1 CONSTS(4) = 0.0001 CONSTS(5) = 0.1 CONSTS(6) = 0.001 STATES(2) = 1.02 STATES(3) = 0.94 STATES(4) = 3 STATES(5) = 0 STATES(6) = 0 STATES(7) = 0 STATES(8) = 0 CONSTS(7) = 0.01 CONSTS(8) = 0.01 CONSTS(9) = 8 CONSTS(10) = 0.75 CONSTS(11) = 1.5 STATES(9) = 0.6 CONSTS(12) = 0.01 CONSTS(13) = 0.01 CONSTS(14) = 0.15 CONSTS(15) = 0.6 CONSTS(16) = 0.017 CONSTS(17) = 0.01 STATES(10) = 0 STATES(11) = 0 STATES(12) = 0 CONSTS(18) = 0.023 CONSTS(19) = 0.3 CONSTS(20) = 0 STATES(13) = 0 CONSTS(21) = 0.001 CONSTS(22) = 4 CONSTS(23) = 0.4 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(2) = ( CONSTS(10)*(CONSTS(9) - STATES(2)))/((CONSTS(7)+CONSTS(9)) - STATES(2)) - ( CONSTS(11)*STATES(9)*STATES(2))/(CONSTS(8)+STATES(2)) RATES(9) = ( CONSTS(14)*(1.00000 - STATES(9)))/((CONSTS(12)+1.00000) - STATES(9)) - ( CONSTS(15)*STATES(1)*STATES(9))/(CONSTS(13)+STATES(9)) RATES(4) = CONSTS(6)*(STATES(6)+STATES(8)+STATES(11)+STATES(12)) - CONSTS(5)*STATES(4)*(STATES(1)+STATES(3)+STATES(5)+STATES(7)) RATES(13) = ( CONSTS(16)*STATES(10)*STATES(11)+ CONSTS(16)*STATES(10)*STATES(12)) - CONSTS(17)*STATES(13) ALGBRC(1) = STATES(5)+STATES(7)+STATES(11)+STATES(12) ALGBRC(5) = (CONSTS(21)+ALGBRC(1) ** CONSTS(22))/(CONSTS(23) ** CONSTS(22)+ALGBRC(1) ** CONSTS(22)) RATES(1) = ( CONSTS(3)*STATES(3)+ CONSTS(4)*STATES(5)+ CONSTS(6)*STATES(6)) - ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(2)*ALGBRC(5)*STATES(1)+ CONSTS(5)*STATES(4)*STATES(1)) RATES(3) = ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(4)*STATES(7)+ CONSTS(6)*STATES(8)) - ( CONSTS(3)*STATES(3)+ CONSTS(2)*ALGBRC(5)*STATES(3)+ CONSTS(5)*STATES(4)*STATES(3)) RATES(5) = ( CONSTS(2)*ALGBRC(5)*STATES(1)+ CONSTS(17)*STATES(11)+ CONSTS(6)*STATES(11)) - ( CONSTS(4)*STATES(5)+ CONSTS(16)*STATES(10)*STATES(5)+ CONSTS(5)*STATES(4)*STATES(5)) RATES(7) = ( CONSTS(2)*ALGBRC(5)*STATES(3)+ CONSTS(17)*STATES(12)+ CONSTS(6)*STATES(12)) - ( CONSTS(4)*STATES(7)+ CONSTS(16)*STATES(10)*STATES(7)+ CONSTS(5)*STATES(4)*STATES(7)) ALGBRC(2) = STATES(5) - CONSTS(19) ALGBRC(6) = TERNRY(ALGBRC(2).LT.0.00000, 0.00000, 1.00000) RATES(10) = CONSTS(18)*ALGBRC(6) RATES(6) = ( CONSTS(5)*STATES(4)*STATES(1)+ CONSTS(4)*STATES(11)+ CONSTS(3)*STATES(8)) - ( CONSTS(6)*STATES(6)+ CONSTS(2)*ALGBRC(5)*STATES(6)+ CONSTS(1)*STATES(2)*STATES(6)) RATES(8) = ( CONSTS(5)*STATES(4)*STATES(3)+ CONSTS(4)*STATES(12)+ CONSTS(1)*STATES(2)*STATES(6)) - ( CONSTS(6)*STATES(8)+ CONSTS(2)*ALGBRC(5)*STATES(8)+ CONSTS(3)*STATES(8)) RATES(11) = ( CONSTS(5)*STATES(4)*STATES(5)+ CONSTS(2)*ALGBRC(5)*STATES(6)) - ( CONSTS(6)*STATES(11)+ CONSTS(4)*STATES(11)+ CONSTS(16)*STATES(10)*STATES(11)+ CONSTS(17)*STATES(11)) RATES(12) = ( CONSTS(5)*STATES(4)*STATES(7)+ CONSTS(2)*ALGBRC(5)*STATES(8)) - ( CONSTS(6)*STATES(12)+ CONSTS(4)*STATES(12)+ CONSTS(16)*CONSTS(20)*STATES(12)+ CONSTS(17)*STATES(12)) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = STATES(5)+STATES(7)+STATES(11)+STATES(12) ALGBRC(5) = (CONSTS(21)+ALGBRC(1) ** CONSTS(22))/(CONSTS(23) ** CONSTS(22)+ALGBRC(1) ** CONSTS(22)) ALGBRC(2) = STATES(5) - CONSTS(19) ALGBRC(6) = TERNRY(ALGBRC(2).LT.0.00000, 0.00000, 1.00000) ALGBRC(3) = STATES(8)+STATES(6)+STATES(12)+STATES(11)+STATES(7)+STATES(5)+STATES(1)+STATES(3) ALGBRC(4) = (STATES(4)+STATES(8)+STATES(6)+STATES(12)+STATES(11)+STATES(13))/3.00000 RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END