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 16 entries in each of the rate and state variable arrays. C There are a total of 58 entries in the constant variable array. C C C VOI is time in component environment (minute). C STATES(1) is MPF in component MPF (dimensionless). C CONSTS(1) is k1 in component rate_constants (first_order_rate_constant). C ALGBRC(1) is k2 in component rate_constants (first_order_rate_constant). C ALGBRC(3) is kwee in component rate_constants (first_order_rate_constant). C ALGBRC(5) is kc25 in component rate_constants (first_order_rate_constant). C CONSTS(2) is kj in component rate_constants (first_order_rate_constant). C CONSTS(3) is kjr in component rate_constants (first_order_rate_constant). C CONSTS(4) is k6_ in component rate_constants (first_order_rate_constant). C CONSTS(5) is k6 in component rate_constants (first_order_rate_constant). C STATES(2) is mass in component mass (dimensionless). C STATES(3) is preMPF in component preMPF (dimensionless). C STATES(4) is Rum1 in component Rum1 (dimensionless). C STATES(5) is Rum1P in component Rum1P (dimensionless). C STATES(6) is CR in component CR (dimensionless). C STATES(7) is CRP in component CRP (dimensionless). C STATES(8) is Ste9 in component Ste9 (dimensionless). C CONSTS(6) is kste9r_ in component rate_constants (first_order_rate_constant). C CONSTS(7) is kste9r in component rate_constants (first_order_rate_constant). C CONSTS(8) is kste9 in component rate_constants (first_order_rate_constant). C CONSTS(9) is Jste9r in component Ste9 (dimensionless). C CONSTS(10) is Jste9 in component Ste9 (dimensionless). C ALGBRC(2) is MPF_a in component MPF_a (dimensionless). C ALGBRC(4) is PP in component PP (dimensionless). C CONSTS(11) is SK in component dimensionless_constants (dimensionless). C STATES(9) is Mik1 in component Mik1 (dimensionless). C CONSTS(12) is ks in component rate_constants (first_order_rate_constant). C CONSTS(13) is kmr_ in component rate_constants (first_order_rate_constant). C CONSTS(14) is kmr in component rate_constants (first_order_rate_constant). C CONSTS(15) is km in component rate_constants (first_order_rate_constant). C CONSTS(16) is Jmikr in component Mik1 (dimensionless). C CONSTS(17) is Jmik in component Mik1 (dimensionless). C STATES(10) is Wee1 in component Wee1 (dimensionless). C CONSTS(18) is kwr_ in component rate_constants (first_order_rate_constant). C CONSTS(19) is kwr in component rate_constants (first_order_rate_constant). C CONSTS(20) is kw in component rate_constants (first_order_rate_constant). C CONSTS(21) is Jweer in component Wee1 (dimensionless). C CONSTS(22) is Jwee in component Wee1 (dimensionless). C STATES(11) is Cdc25 in component Cdc25 (dimensionless). C CONSTS(23) is k25 in component rate_constants (first_order_rate_constant). C CONSTS(24) is k5 in component rate_constants (first_order_rate_constant). C CONSTS(25) is k25r_ in component rate_constants (first_order_rate_constant). C CONSTS(26) is k25r in component rate_constants (first_order_rate_constant). C CONSTS(27) is J25 in component Cdc25 (dimensionless). C CONSTS(28) is J25r in component Cdc25 (dimensionless). C STATES(12) is Slp1 in component Slp1 (dimensionless). C CONSTS(29) is kas in component rate_constants (first_order_rate_constant). C CONSTS(30) is kad in component rate_constants (first_order_rate_constant). C STATES(13) is Slp1_a in component Slp1_a (dimensionless). C CONSTS(31) is kaa_ in component rate_constants (first_order_rate_constant). C CONSTS(32) is kaa in component rate_constants (first_order_rate_constant). C CONSTS(33) is kai in component rate_constants (first_order_rate_constant). C STATES(14) is Inh in component Inh (dimensionless). C CONSTS(34) is k3 in component rate_constants (first_order_rate_constant). C CONSTS(35) is ki in component rate_constants (first_order_rate_constant). C CONSTS(36) is kir in component rate_constants (first_order_rate_constant). C ALGBRC(6) is k4 in component rate_constants (first_order_rate_constant). C STATES(15) is PI in component PI (dimensionless). C CONSTS(37) is kp in component rate_constants (first_order_rate_constant). C CONSTS(38) is kpp_ in component rate_constants (first_order_rate_constant). C CONSTS(39) is kpp in component rate_constants (first_order_rate_constant). C ALGBRC(7) is k2c in component rate_constants (first_order_rate_constant). C CONSTS(40) is epsilon_p in component dimensionless_constants (dimensionless). C CONSTS(41) is mu in component mass (first_order_rate_constant). C STATES(16) is R_dna in component R_dna (dimensionless). C CONSTS(42) is K in component R_dna (dimensionless). C CONSTS(43) is Y in component R_dna (dimensionless). C CONSTS(44) is epsilon in component dimensionless_constants (dimensionless). C ALGBRC(8) is ratio in component ratio (dimensionless). C CONSTS(45) is V2_ in component rate_constants (first_order_rate_constant). C CONSTS(46) is V2 in component rate_constants (first_order_rate_constant). C CONSTS(47) is V2c in component rate_constants (first_order_rate_constant). C CONSTS(48) is V2c_ in component rate_constants (first_order_rate_constant). C CONSTS(49) is Vwee in component rate_constants (first_order_rate_constant). C CONSTS(50) is Vwee_ in component rate_constants (first_order_rate_constant). C CONSTS(51) is Vmik in component rate_constants (first_order_rate_constant). C CONSTS(52) is Vmik_ in component rate_constants (first_order_rate_constant). C CONSTS(53) is V25 in component rate_constants (first_order_rate_constant). C CONSTS(54) is V25_ in component rate_constants (first_order_rate_constant). C CONSTS(55) is Vpyp in component rate_constants (first_order_rate_constant). C CONSTS(56) is V4 in component rate_constants (first_order_rate_constant). C CONSTS(57) is V4_ in component rate_constants (first_order_rate_constant). C CONSTS(58) is Pyp3 in component rate_constants (dimensionless). C RATES(1) is d/dt MPF in component MPF (dimensionless). C RATES(3) is d/dt preMPF in component preMPF (dimensionless). C RATES(8) is d/dt Ste9 in component Ste9 (dimensionless). C RATES(9) is d/dt Mik1 in component Mik1 (dimensionless). C RATES(10) is d/dt Wee1 in component Wee1 (dimensionless). C RATES(11) is d/dt Cdc25 in component Cdc25 (dimensionless). C RATES(12) is d/dt Slp1 in component Slp1 (dimensionless). C RATES(13) is d/dt Slp1_a in component Slp1_a (dimensionless). C RATES(14) is d/dt Inh in component Inh (dimensionless). C RATES(15) is d/dt PI in component PI (dimensionless). C RATES(7) is d/dt CRP in component CRP (dimensionless). C RATES(6) is d/dt CR in component CR (dimensionless). C RATES(4) is d/dt Rum1 in component Rum1 (dimensionless). C RATES(5) is d/dt Rum1P in component Rum1P (dimensionless). C RATES(2) is d/dt mass in component mass (dimensionless). C RATES(16) is d/dt R_dna in component R_dna (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0 CONSTS(1) = 0.02 CONSTS(2) = 400 CONSTS(3) = 1 CONSTS(4) = 0.1 CONSTS(5) = 5 STATES(2) = 1 STATES(3) = 1 STATES(4) = 0 STATES(5) = 0 STATES(6) = 0 STATES(7) = 0 STATES(8) = 0 CONSTS(6) = 0.03 CONSTS(7) = 8 CONSTS(8) = 5 CONSTS(9) = 0.01 CONSTS(10) = 0.01 CONSTS(11) = 0.018 STATES(9) = 0 CONSTS(12) = 0.1 CONSTS(13) = 0.01 CONSTS(14) = 5 CONSTS(15) = 1 CONSTS(16) = 0.15 CONSTS(17) = 0.15 STATES(10) = 0 CONSTS(18) = 0.4 CONSTS(19) = 1 CONSTS(20) = 2 CONSTS(21) = 0.2 CONSTS(22) = 0.2 STATES(11) = 0 CONSTS(23) = 1 CONSTS(24) = 0.1 CONSTS(25) = 0.4 CONSTS(26) = 2 CONSTS(27) = 0.05 CONSTS(28) = 0.05 STATES(12) = 0 CONSTS(29) = 0.1 CONSTS(30) = 0.1 STATES(13) = 0 CONSTS(31) = 0.01 CONSTS(32) = 0.1 CONSTS(33) = 0.1 STATES(14) = 0 CONSTS(34) = 0.1 CONSTS(35) = 50 CONSTS(36) = 0.5 STATES(15) = 0.2 CONSTS(37) = 100 CONSTS(38) = 1 CONSTS(39) = 100 CONSTS(40) = 0.025 CONSTS(41) = 0.00462 STATES(16) = 1 CONSTS(42) = 0.06 CONSTS(43) = 0 CONSTS(44) = 0.05 CONSTS(45) = 0.02 CONSTS(46) = 1 CONSTS(47) = 0.5 CONSTS(48) = 0.02 CONSTS(49) = 10 CONSTS(50) = 0.08 CONSTS(51) = 2 CONSTS(52) = 0.04 CONSTS(53) = 10 CONSTS(54) = 0.05 CONSTS(55) = 0.07 CONSTS(56) = 1 CONSTS(57) = 0.01 CONSTS(58) = 1 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(2) = CONSTS(41)*STATES(2) ALGBRC(2) = STATES(1)+ CONSTS(44)*STATES(3) RATES(12) = CONSTS(29)*ALGBRC(2) - CONSTS(30)*STATES(12) RATES(13) = (CONSTS(31)+ CONSTS(32)*ALGBRC(2))*(STATES(12) - STATES(13)) - (CONSTS(33)+CONSTS(30))*STATES(13) RATES(16) = ( CONSTS(42)*1.00000)/(1.00000+ CONSTS(43)*ALGBRC(2)) ALGBRC(4) = 1.00000 - STATES(15) RATES(8) = ( (CONSTS(6)+ CONSTS(7)*ALGBRC(4))*(1.00000 - STATES(8)))/((CONSTS(9)+1.00000) - STATES(8)) - ( CONSTS(8)*(ALGBRC(2)+ CONSTS(11)*STATES(2))*STATES(8))/(CONSTS(10)+STATES(8)) RATES(9) = ( (CONSTS(12)+CONSTS(13)+ CONSTS(14)*ALGBRC(4))*(1.00000 - STATES(9)))/((CONSTS(16)+1.00000) - STATES(9)) - ( CONSTS(15)*ALGBRC(2)*STATES(9))/(CONSTS(17)+STATES(9)) RATES(10) = ( (CONSTS(18)+ CONSTS(19)*ALGBRC(4))*(1.00000 - STATES(10)))/((CONSTS(21)+1.00000) - STATES(10)) - ( CONSTS(20)*ALGBRC(2)*STATES(10))/(CONSTS(22)+STATES(10)) RATES(11) = ( CONSTS(23)*ALGBRC(2)*(1.00000 - STATES(11)))/((CONSTS(27)+1.00000) - STATES(11)) - ( (CONSTS(24)+CONSTS(25)+ CONSTS(26)*ALGBRC(4))*STATES(11))/(CONSTS(28)+STATES(11)) ALGBRC(1) = CONSTS(45)+ CONSTS(46)*STATES(8) ALGBRC(3) = CONSTS(49)*STATES(10)+ CONSTS(50)*(1.00000 - STATES(10))+ CONSTS(51)*STATES(9)+ CONSTS(52)*(1.00000 - STATES(9)) ALGBRC(5) = CONSTS(53)*STATES(11)+ CONSTS(54)*(1.00000 - STATES(11))+ CONSTS(55)*CONSTS(58) RATES(1) = (((( CONSTS(1)*STATES(2) - ALGBRC(1)*STATES(1)) - ALGBRC(3)*STATES(1))+ ALGBRC(5)*STATES(3)) - CONSTS(2)*STATES(1)*(STATES(4)+STATES(5)))+ CONSTS(3)*(STATES(6)+STATES(7))+ CONSTS(4)*STATES(6)+ (CONSTS(4)+CONSTS(5))*STATES(7) RATES(3) = ( ALGBRC(3)*STATES(1) - ALGBRC(5)*STATES(3)) - ALGBRC(1)*STATES(3) ALGBRC(6) = CONSTS(57)+ CONSTS(56)*STATES(13) RATES(14) = ((CONSTS(34) - CONSTS(35)*STATES(14)*ALGBRC(4))+ CONSTS(36)*STATES(15)) - ALGBRC(6)*STATES(14) RATES(15) = ( CONSTS(35)*STATES(14)*ALGBRC(4) - CONSTS(36)*STATES(15)) - ALGBRC(6)*STATES(15) ALGBRC(7) = CONSTS(48)+ CONSTS(47)*STATES(8) RATES(7) = (((( CONSTS(37)*(ALGBRC(2)+ CONSTS(40)*CONSTS(11)*STATES(2))*STATES(6) - (CONSTS(38)+ CONSTS(39)*ALGBRC(4))*STATES(7))+ CONSTS(2)*STATES(1)*STATES(5)) - CONSTS(3)*STATES(7)) - ALGBRC(7)*STATES(7)) - (CONSTS(4)+CONSTS(5))*STATES(7) RATES(6) = (((( CONSTS(2)*STATES(1)*STATES(4) - CONSTS(3)*STATES(6)) - ALGBRC(7)*STATES(6)) - CONSTS(4)*STATES(6)) - CONSTS(37)*(ALGBRC(2)+ CONSTS(40)*CONSTS(11)*STATES(2))*STATES(6))+ (CONSTS(38)+ CONSTS(39)*ALGBRC(4))*STATES(7) RATES(4) = ((((CONSTS(24) - CONSTS(4)*STATES(4)) - CONSTS(37)*(ALGBRC(2)+ CONSTS(40)*CONSTS(11)*STATES(2))*STATES(4))+ (CONSTS(38)+ CONSTS(39)*ALGBRC(4))*STATES(5)) - CONSTS(2)*STATES(1)*STATES(4))+ CONSTS(3)*STATES(6)+ ALGBRC(7)*STATES(6) RATES(5) = ((( CONSTS(37)*(ALGBRC(2)+ CONSTS(40)*CONSTS(11)*STATES(2))*STATES(4) - (CONSTS(38)+ CONSTS(39)*ALGBRC(4))*STATES(5)) - (CONSTS(4)+CONSTS(5))*STATES(5)) - CONSTS(2)*STATES(1)*STATES(5))+ CONSTS(3)*STATES(7)+ ALGBRC(7)*STATES(7) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(2) = STATES(1)+ CONSTS(44)*STATES(3) ALGBRC(4) = 1.00000 - STATES(15) ALGBRC(1) = CONSTS(45)+ CONSTS(46)*STATES(8) ALGBRC(3) = CONSTS(49)*STATES(10)+ CONSTS(50)*(1.00000 - STATES(10))+ CONSTS(51)*STATES(9)+ CONSTS(52)*(1.00000 - STATES(9)) ALGBRC(5) = CONSTS(53)*STATES(11)+ CONSTS(54)*(1.00000 - STATES(11))+ CONSTS(55)*CONSTS(58) ALGBRC(6) = CONSTS(57)+ CONSTS(56)*STATES(13) ALGBRC(7) = CONSTS(48)+ CONSTS(47)*STATES(8) ALGBRC(8) = ALGBRC(2)/ALGBRC(4) RETURN END