C C There are a total of 2 entries in the algebraic variable array. C There are a total of 34 entries in each of the rate and state variable arrays. C There are a total of 55 entries in the constant variable array. C C C VOI is time in component environment (minute). C STATES(1) is v_c in component v_c (nanomolar). C ALGBRC(1) is dv_c_dt in component dv_c_dt (flux). C CONSTS(1) is k_growth in component rate_constants (first_order_rate_constant). C STATES(2) is v_n in component v_n (nanomolar). C ALGBRC(2) is k_volume in component k_volume (dimensionless). C STATES(3) is mcln2_c in component mcln2_c (nanomolar). C CONSTS(2) is k_50 in component rate_constants (first_order_rate_constant). C CONSTS(3) is k_10 in component rate_constants (first_order_rate_constant). C STATES(4) is mcln2_n in component mcln2_n (nanomolar). C STATES(5) is mclb5_c in component mclb5_c (nanomolar). C CONSTS(4) is k_51 in component rate_constants (first_order_rate_constant). C CONSTS(5) is k_11 in component rate_constants (first_order_rate_constant). C STATES(6) is mclb5_n in component mclb5_n (nanomolar). C CONSTS(6) is k_1 in component rate_constants (first_order_rate_constant). C STATES(7) is sbf_n in component sbf_n (nanomolar). C CONSTS(7) is k_2 in component rate_constants (first_order_rate_constant). C CONSTS(8) is k_39 in component rate_constants (first_order_rate_constant). C CONSTS(9) is k_34 in component rate_constants (second_order_rate_constant). C CONSTS(10) is k_35 in component rate_constants (first_order_rate_constant). C STATES(8) is sbfwhi5p_n in component sbfwhi5p_n (nanomolar). C STATES(9) is whi5_n in component whi5_n (nanomolar). C STATES(10) is sbfwhi5_n in component sbfwhi5_n (nanomolar). C STATES(11) is cln3_c in component cln3_c (nanomolar). C CONSTS(11) is k_6 in component rate_constants (flux). C CONSTS(12) is k_43 in component rate_constants (first_order_rate_constant). C CONSTS(13) is k_15 in component rate_constants (first_order_rate_constant). C STATES(12) is cln2_c in component cln2_c (nanomolar). C CONSTS(14) is k_3 in component rate_constants (first_order_rate_constant). C CONSTS(15) is k_26 in component rate_constants (second_order_rate_constant). C CONSTS(16) is k_27 in component rate_constants (first_order_rate_constant). C CONSTS(17) is k_12 in component rate_constants (first_order_rate_constant). C STATES(13) is cdk1_c in component cdk1_c (nanomolar). C STATES(14) is cdk1cln2_c in component cdk1cln2_c (nanomolar). C STATES(15) is clb5_c in component clb5_c (nanomolar). C CONSTS(18) is k_4 in component rate_constants (first_order_rate_constant). C CONSTS(19) is k_28 in component rate_constants (second_order_rate_constant). C CONSTS(20) is k_29 in component rate_constants (first_order_rate_constant). C CONSTS(21) is k_13 in component rate_constants (first_order_rate_constant). C STATES(16) is cdk1clb5_c in component cdk1clb5_c (nanomolar). C CONSTS(22) is k_7 in component rate_constants (flux). C CONSTS(23) is k_44 in component rate_constants (first_order_rate_constant). C CONSTS(24) is k_49 in component rate_constants (first_order_rate_constant). C CONSTS(25) is k_16 in component rate_constants (first_order_rate_constant). C STATES(17) is cdk1_n in component cdk1_n (nanomolar). C STATES(18) is cln3_n in component cln3_n (nanomolar). C CONSTS(26) is k_24 in component rate_constants (second_order_rate_constant). C CONSTS(27) is k_25 in component rate_constants (first_order_rate_constant). C CONSTS(28) is k_20 in component rate_constants (first_order_rate_constant). C STATES(19) is cdk1cln3_n in component cdk1cln3_n (nanomolar). C CONSTS(29) is k_21 in component rate_constants (first_order_rate_constant). C CONSTS(30) is k_46 in component rate_constants (first_order_rate_constant). C CONSTS(31) is k_53 in component rate_constants (first_order_rate_constant). C STATES(20) is cdk1cln2_n in component cdk1cln2_n (nanomolar). C CONSTS(32) is k_33 in component rate_constants (first_order_rate_constant). C CONSTS(33) is k_32 in component rate_constants (second_order_rate_constant). C CONSTS(34) is k_48 in component rate_constants (first_order_rate_constant). C STATES(21) is cdk1clb5sic1_c in component cdk1clb5sic1_c (nanomolar). C STATES(22) is sic1_c in component sic1_c (nanomolar). C STATES(23) is cdk1clb5_n in component cdk1clb5_n (nanomolar). C CONSTS(35) is k_41 in component rate_constants (first_order_rate_constant). C STATES(24) is cdk1clb5sic1p_n in component cdk1clb5sic1p_n (nanomolar). C CONSTS(36) is k_31 in component rate_constants (first_order_rate_constant). C CONSTS(37) is k_30 in component rate_constants (second_order_rate_constant). C CONSTS(38) is k_40 in component rate_constants (first_order_rate_constant). C STATES(25) is cdk1cln3far1_n in component cdk1cln3far1_n (nanomolar). C STATES(26) is far1_n in component far1_n (nanomolar). C STATES(27) is cdk1cln3far1p_n in component cdk1cln3far1p_n (nanomolar). C CONSTS(39) is k_36 in component rate_constants (second_order_rate_constant). C CONSTS(40) is k_37 in component rate_constants (second_order_rate_constant). C STATES(28) is cdk1clb5sic1_n in component cdk1clb5sic1_n (nanomolar). C CONSTS(41) is k_47 in component rate_constants (first_order_rate_constant). C CONSTS(42) is k_38 in component rate_constants (second_order_rate_constant). C CONSTS(43) is k_9 in component rate_constants (flux). C CONSTS(44) is k_18 in component rate_constants (first_order_rate_constant). C STATES(29) is whi5_c in component whi5_c (nanomolar). C CONSTS(45) is k_8 in component rate_constants (flux). C CONSTS(46) is k_45 in component rate_constants (first_order_rate_constant). C CONSTS(47) is k_17 in component rate_constants (first_order_rate_constant). C STATES(30) is far1_c in component far1_c (nanomolar). C CONSTS(48) is k_5 in component rate_constants (flux). C CONSTS(49) is k_42 in component rate_constants (first_order_rate_constant). C CONSTS(50) is k_14 in component rate_constants (first_order_rate_constant). C STATES(31) is whi5p_c in component whi5p_c (nanomolar). C CONSTS(51) is k_52 in component rate_constants (first_order_rate_constant). C STATES(32) is whi5p_n in component whi5p_n (nanomolar). C STATES(33) is sic1P_n in component sic1P_n (nanomolar). C STATES(34) is far1P_n in component far1P_n (nanomolar). C CONSTS(52) is k_22 in component rate_constants (first_order_rate_constant). C CONSTS(53) is k_19 in component rate_constants (first_order_rate_constant). C CONSTS(54) is k_23 in component rate_constants (first_order_rate_constant). C RATES(1) is d/dt v_c in component v_c (nanomolar). C RATES(2) is d/dt v_n in component v_n (nanomolar). C RATES(3) is d/dt mcln2_c in component mcln2_c (nanomolar). C RATES(5) is d/dt mclb5_c in component mclb5_c (nanomolar). C RATES(4) is d/dt mcln2_n in component mcln2_n (nanomolar). C RATES(6) is d/dt mclb5_n in component mclb5_n (nanomolar). C RATES(7) is d/dt sbf_n in component sbf_n (nanomolar). C RATES(11) is d/dt cln3_c in component cln3_c (nanomolar). C RATES(12) is d/dt cln2_c in component cln2_c (nanomolar). C RATES(15) is d/dt clb5_c in component clb5_c (nanomolar). C RATES(13) is d/dt cdk1_c in component cdk1_c (nanomolar). C RATES(18) is d/dt cln3_n in component cln3_n (nanomolar). C RATES(17) is d/dt cdk1_n in component cdk1_n (nanomolar). C RATES(14) is d/dt cdk1cln2_c in component cdk1cln2_c (nanomolar). C RATES(16) is d/dt cdk1clb5_c in component cdk1clb5_c (nanomolar). C RATES(20) is d/dt cdk1cln2_n in component cdk1cln2_n (nanomolar). C RATES(23) is d/dt cdk1clb5_n in component cdk1clb5_n (nanomolar). C RATES(19) is d/dt cdk1cln3_n in component cdk1cln3_n (nanomolar). C RATES(10) is d/dt sbfwhi5_n in component sbfwhi5_n (nanomolar). C RATES(8) is d/dt sbfwhi5p_n in component sbfwhi5p_n (nanomolar). C RATES(25) is d/dt cdk1cln3far1_n in component cdk1cln3far1_n (nanomolar). C RATES(28) is d/dt cdk1clb5sic1_n in component cdk1clb5sic1_n (nanomolar). C RATES(24) is d/dt cdk1clb5sic1p_n in component cdk1clb5sic1p_n (nanomolar). C RATES(27) is d/dt cdk1cln3far1p_n in component cdk1cln3far1p_n (nanomolar). C RATES(21) is d/dt cdk1clb5sic1_c in component cdk1clb5sic1_c (nanomolar). C RATES(22) is d/dt sic1_c in component sic1_c (nanomolar). C RATES(29) is d/dt whi5_c in component whi5_c (nanomolar). C RATES(30) is d/dt far1_c in component far1_c (nanomolar). C RATES(31) is d/dt whi5p_c in component whi5p_c (nanomolar). C RATES(33) is d/dt sic1P_n in component sic1P_n (nanomolar). C RATES(34) is d/dt far1P_n in component far1P_n (nanomolar). C RATES(9) is d/dt whi5_n in component whi5_n (nanomolar). C RATES(26) is d/dt far1_n in component far1_n (nanomolar). C RATES(32) is d/dt whi5p_n in component whi5p_n (nanomolar). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.5 CONSTS(1) = 0.0051 STATES(2) = 0.5 STATES(3) = 0 CONSTS(2) = 0.6 CONSTS(3) = 0.12 STATES(4) = 0 STATES(5) = 0 CONSTS(4) = 0.6 CONSTS(5) = 0.12 STATES(6) = 0 CONSTS(6) = 0.03523 STATES(7) = 0 CONSTS(7) = 0.03523 CONSTS(8) = 1 CONSTS(9) = 8.46 CONSTS(10) = 0.0005 STATES(8) = 0 STATES(9) = 0 STATES(10) = 0.025544 STATES(11) = 0.000485 CONSTS(11) = 0.00001 CONSTS(12) = 0.005 CONSTS(13) = 0.01 STATES(12) = 0 CONSTS(14) = 0.32 CONSTS(15) = 2.82 CONSTS(16) = 0.55 CONSTS(17) = 0.1 STATES(13) = 0.333333 STATES(14) = 0 STATES(15) = 0 CONSTS(18) = 0.32 CONSTS(19) = 2.82 CONSTS(20) = 0.55 CONSTS(21) = 0.35 STATES(16) = 0 CONSTS(22) = 0.01 CONSTS(23) = 0.005 CONSTS(24) = 0.001 CONSTS(25) = 0.03 STATES(17) = 0.0074127 STATES(18) = 0 CONSTS(26) = 2.82 CONSTS(27) = 0.55 CONSTS(28) = 0.01 STATES(19) = 0 CONSTS(29) = 0 CONSTS(30) = 0.1 CONSTS(31) = 0.001 STATES(20) = 0 CONSTS(32) = 0.55 CONSTS(33) = 84.6 CONSTS(34) = 0.012 STATES(21) = 0 STATES(22) = 0.039234 STATES(23) = 0 CONSTS(35) = 1 STATES(24) = 0 CONSTS(36) = 0.55 CONSTS(37) = 42300 CONSTS(38) = 1 STATES(25) = 0 STATES(26) = 0 STATES(27) = 0 CONSTS(39) = 4363.6 CONSTS(40) = 4363.6 STATES(28) = 0 CONSTS(41) = 1 CONSTS(42) = 4363.6 CONSTS(43) = 0.00005 CONSTS(44) = 0.0008 STATES(29) = 0.073564 CONSTS(45) = 0.00004 CONSTS(46) = 0.005 CONSTS(47) = 0.01 STATES(30) = 0.0037926 CONSTS(48) = 0.000042 CONSTS(49) = 0.005 CONSTS(50) = 0.01 STATES(31) = 0 CONSTS(51) = 0.005 STATES(32) = 0 STATES(33) = 0 STATES(34) = 0 CONSTS(52) = 0.01 CONSTS(53) = 0.01 CONSTS(54) = 0.01 CONSTS(54) = 0.00000 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(2) = CONSTS(54) RATES(4) = CONSTS(6)*STATES(7) - CONSTS(2)*STATES(4) RATES(6) = CONSTS(7)*STATES(7) - CONSTS(4)*STATES(6) RATES(7) = ( CONSTS(8)*STATES(8) - CONSTS(9)*STATES(7)*STATES(9))+ CONSTS(10)*STATES(10) RATES(19) = ((( CONSTS(26)*STATES(18)*STATES(17) - CONSTS(27)*STATES(19))+ CONSTS(36)*STATES(25)) - CONSTS(37)*STATES(26)*STATES(19))+ CONSTS(38)*STATES(27) RATES(10) = ( CONSTS(9)*STATES(7)*STATES(9) - CONSTS(10)*STATES(10)) - CONSTS(39)*STATES(10)*STATES(19) RATES(8) = CONSTS(39)*STATES(10)*STATES(19) - CONSTS(8)*STATES(8) RATES(25) = ( CONSTS(37)*STATES(26)*STATES(19) - CONSTS(36)*STATES(25)) - CONSTS(40)*STATES(25)*STATES(20) RATES(24) = CONSTS(42)*STATES(28)*STATES(20) - CONSTS(35)*STATES(24) RATES(27) = CONSTS(40)*STATES(25)*STATES(20) - CONSTS(38)*STATES(27) RATES(33) = CONSTS(35)*STATES(24) RATES(34) = CONSTS(38)*STATES(27) RATES(32) = ( CONSTS(8)*STATES(8) - CONSTS(51)*STATES(32)) - CONSTS(54)*STATES(32) ALGBRC(1) = CONSTS(1)*STATES(1) RATES(1) = ALGBRC(1) RATES(11) = ((CONSTS(11) - CONSTS(12)*STATES(11)) - CONSTS(13)*STATES(11)) - (ALGBRC(1)/STATES(1))*STATES(11) RATES(12) = ((( CONSTS(14)*STATES(3) - CONSTS(15)*STATES(13)*STATES(12))+ CONSTS(16)*STATES(14)) - CONSTS(17)*STATES(12)) - (ALGBRC(1)/STATES(1))*STATES(12) RATES(15) = ((( CONSTS(18)*STATES(5) - CONSTS(19)*STATES(15)*STATES(13))+ CONSTS(20)*STATES(16)) - CONSTS(21)*STATES(15)) - (ALGBRC(1)/STATES(1))*STATES(15) RATES(13) = ((((((CONSTS(22) - CONSTS(23)*STATES(13))+ CONSTS(24)*STATES(17)+ CONSTS(16)*STATES(14)) - CONSTS(15)*STATES(13)*STATES(12)) - CONSTS(19)*STATES(13)*STATES(15))+ CONSTS(20)*STATES(16)) - CONSTS(25)*STATES(13)) - (ALGBRC(1)/STATES(1))*STATES(13) RATES(16) = (((( CONSTS(19)*STATES(13)*STATES(15) - CONSTS(20)*STATES(16))+ CONSTS(32)*STATES(21)) - CONSTS(33)*STATES(22)*STATES(16)) - CONSTS(34)*STATES(16)) - (ALGBRC(1)/STATES(1))*STATES(16) RATES(21) = (( CONSTS(33)*STATES(22)*STATES(16) - CONSTS(32)*STATES(21)) - CONSTS(41)*STATES(21)) - (ALGBRC(1)/STATES(1))*STATES(21) RATES(22) = (((CONSTS(43) - CONSTS(33)*STATES(22)*STATES(16))+ CONSTS(32)*STATES(21)) - CONSTS(44)*STATES(22)) - (ALGBRC(1)/STATES(1))*STATES(22) RATES(29) = ((CONSTS(45) - CONSTS(46)*STATES(29)) - CONSTS(47)*STATES(29)) - (ALGBRC(1)/STATES(1))*STATES(29) RATES(30) = ((CONSTS(48) - CONSTS(49)*STATES(30)) - CONSTS(50)*STATES(30)) - (ALGBRC(1)/STATES(1))*STATES(30) ALGBRC(2) = STATES(2)/STATES(1) RATES(3) = ( CONSTS(2)*STATES(4)*ALGBRC(2) - CONSTS(3)*STATES(3)) - (ALGBRC(1)/STATES(1))*STATES(3) RATES(5) = ( CONSTS(4)*STATES(6)*ALGBRC(2) - CONSTS(5)*STATES(5)) - (ALGBRC(1)/STATES(1))*STATES(5) RATES(18) = ((( CONSTS(12)*STATES(11))/ALGBRC(2) - CONSTS(26)*STATES(18)*STATES(17))+ CONSTS(27)*STATES(19)) - CONSTS(28)*STATES(18) RATES(17) = (((( CONSTS(23)*STATES(13))/ALGBRC(2) - CONSTS(24)*STATES(17)) - CONSTS(26)*STATES(18)*STATES(17))+ CONSTS(27)*STATES(19)) - CONSTS(29)*STATES(17) RATES(14) = ((( CONSTS(15)*STATES(13)*STATES(12) - CONSTS(16)*STATES(14)) - CONSTS(30)*STATES(14))+ CONSTS(31)*STATES(20)*ALGBRC(2)) - (ALGBRC(1)/STATES(1))*STATES(14) RATES(20) = ( CONSTS(30)*STATES(14))/ALGBRC(2) - CONSTS(31)*STATES(20) RATES(23) = CONSTS(35)*STATES(24)+( CONSTS(34)*STATES(16))/ALGBRC(2) RATES(28) = ( CONSTS(41)*STATES(21))/ALGBRC(2) - CONSTS(42)*STATES(28)*STATES(20) RATES(31) = CONSTS(51)*STATES(32)*ALGBRC(2) RATES(9) = (( CONSTS(46)*STATES(29))/ALGBRC(2) - CONSTS(9)*STATES(9)*STATES(7)) - CONSTS(52)*STATES(9) RATES(26) = ((( CONSTS(49)*STATES(30))/ALGBRC(2) - CONSTS(37)*STATES(26)*STATES(19))+ CONSTS(36)*STATES(25)) - CONSTS(53)*STATES(26) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = CONSTS(1)*STATES(1) ALGBRC(2) = STATES(2)/STATES(1) RETURN END