C
C There are a total of 9 entries in the algebraic variable array.
C There are a total of 19 entries in each of the rate and state variable arrays.
C There are a total of 69 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Clb2 in component Clb2 (dimensionless).
C CONSTS(1) is epsilon in component Clb2 (dimensionless).
C CONSTS(2) is Jm in component Clb2 (dimensionless).
C CONSTS(3) is ks_clb in component Clb2 (first_order_rate_constant).
C CONSTS(4) is kd_clb_1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(5) is kd_clb_2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(6) is kd_clb_3 in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(9) is kswe in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(8) is kmih in component kinetic_parameters (first_order_rate_constant).
C CONSTS(7) is kass in component kinetic_parameters (first_order_rate_constant).
C CONSTS(8) is kdiss in component kinetic_parameters (first_order_rate_constant).
C CONSTS(9) is kd_sic in component kinetic_parameters (first_order_rate_constant).
C CONSTS(10) is kd_sic_1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(11) is kd_sic_2 in component kinetic_parameters (first_order_rate_constant).
C STATES(2) is M in component M (dimensionless).
C STATES(3) is Mcm_a in component Mcm_a (dimensionless).
C STATES(4) is Cdh1_a in component Cdh1_a (dimensionless).
C STATES(5) is Cdc20_a in component Cdc20_a (dimensionless).
C STATES(6) is Sic1 in component Sic1 (dimensionless).
C STATES(7) is PClb2 in component PClb2 (dimensionless).
C STATES(8) is Trim in component Trim (dimensionless).
C STATES(9) is Cln in component Cln (dimensionless).
C STATES(10) is PTrim in component PTrim (dimensionless).
C CONSTS(12) is Ji_mcm in component Mcm_a (dimensionless).
C CONSTS(13) is Ja_mcm in component Mcm_a (dimensionless).
C CONSTS(14) is ki_mcm in component Mcm_a (first_order_rate_constant).
C CONSTS(15) is ka_mcm in component Mcm_a (first_order_rate_constant).
C ALGBRC(1) is Mcm in component Mcm (dimensionless).
C CONSTS(16) is ks_sic in component Sic1 (first_order_rate_constant).
C STATES(11) is Mih1_a in component Mih1_a (dimensionless).
C CONSTS(17) is Ji_mih in component Mih1_a (dimensionless).
C CONSTS(18) is Ja_mih in component Mih1_a (dimensionless).
C CONSTS(19) is Vi_mih in component Mih1_a (first_order_rate_constant).
C CONSTS(20) is Va_mih in component Mih1_a (first_order_rate_constant).
C ALGBRC(2) is Mih1 in component Mih1 (dimensionless).
C STATES(12) is IE_a in component IE_a (dimensionless).
C CONSTS(21) is Ji_ie in component IE_a (dimensionless).
C CONSTS(22) is Ja_ie in component IE_a (dimensionless).
C CONSTS(23) is ki_ie in component IE_a (first_order_rate_constant).
C CONSTS(24) is ka_ie in component IE_a (first_order_rate_constant).
C ALGBRC(3) is IE in component IE (dimensionless).
C CONSTS(25) is Ji_cdc20 in component kinetic_parameters (dimensionless).
C CONSTS(26) is Ja_cdc20 in component kinetic_parameters (dimensionless).
C CONSTS(27) is ki_cdc20 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(28) is kd_cdc20 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(29) is ka_cdc20 in component kinetic_parameters (first_order_rate_constant).
C STATES(13) is Cdc20 in component Cdc20 (dimensionless).
C CONSTS(30) is ks_cdc20_1 in component Cdc20 (first_order_rate_constant).
C CONSTS(31) is ks_cdc20_2 in component Cdc20 (first_order_rate_constant).
C CONSTS(32) is Js_cdc20 in component Cdc20 (dimensionless).
C CONSTS(33) is Ji_cdh in component Cdh1_a (dimensionless).
C CONSTS(34) is Ja_cdh in component Cdh1_a (dimensionless).
C CONSTS(35) is ki_cdh in component Cdh1_a (first_order_rate_constant).
C CONSTS(36) is ka_cdh_1 in component Cdh1_a (first_order_rate_constant).
C CONSTS(37) is ki_cdh_1 in component Cdh1_a (first_order_rate_constant).
C CONSTS(38) is ka_cdh_2 in component Cdh1_a (first_order_rate_constant).
C ALGBRC(4) is Cdh1 in component Cdh1 (dimensionless).
C CONSTS(39) is ks_cln in component Cln (first_order_rate_constant).
C CONSTS(40) is kd_cln in component Cln (first_order_rate_constant).
C STATES(14) is SBF_a in component SBF_a (dimensionless).
C CONSTS(41) is Ji_sbf in component SBF_a (dimensionless).
C CONSTS(42) is Ja_sbf in component SBF_a (dimensionless).
C CONSTS(43) is ki_sbf_1 in component SBF_a (first_order_rate_constant).
C CONSTS(44) is ka_sbf_1 in component SBF_a (first_order_rate_constant).
C CONSTS(45) is ki_sbf_2 in component SBF_a (first_order_rate_constant).
C CONSTS(46) is ka_sbf_2 in component SBF_a (first_order_rate_constant).
C ALGBRC(5) is SBF in component SBF (dimensionless).
C STATES(15) is Swe1 in component Swe1 (dimensionless).
C CONSTS(47) is ks_swe in component Swe1 (first_order_rate_constant).
C CONSTS(48) is ks_sweC in component Swe1 (first_order_rate_constant).
C CONSTS(49) is khsl1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(50) is khsl1r in component kinetic_parameters (first_order_rate_constant).
C CONSTS(51) is kd_swe_1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(52) is Ji_wee in component kinetic_parameters (dimensionless).
C CONSTS(53) is Ja_wee in component kinetic_parameters (dimensionless).
C CONSTS(54) is Vi_wee in component kinetic_parameters (first_order_rate_constant).
C CONSTS(55) is Va_wee in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(6) is BUD in component BUD (dimensionless).
C STATES(16) is PSwe1 in component PSwe1 (dimensionless).
C STATES(17) is Swe1M in component Swe1M (dimensionless).
C STATES(18) is PSwe1M in component PSwe1M (dimensionless).
C CONSTS(56) is kd_swe_2 in component PSwe1M (first_order_rate_constant).
C STATES(19) is BE in component BE (dimensionless).
C CONSTS(57) is ks_bud in component BE (first_order_rate_constant).
C CONSTS(58) is kd_bud in component BE (first_order_rate_constant).
C CONSTS(59) is mu in component M (first_order_rate_constant).
C CONSTS(60) is IE_total in component IE (dimensionless).
C CONSTS(61) is Cdh1_total in component Cdh1 (dimensionless).
C CONSTS(62) is Mih1_total in component Mih1 (dimensionless).
C CONSTS(63) is Mcm_total in component Mcm (dimensionless).
C CONSTS(64) is SBF_total in component SBF (dimensionless).
C ALGBRC(7) is Swe1_total in component Swe1_total (dimensionless).
C CONSTS(65) is kswe_1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(66) is kswe_2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(67) is kswe_3 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(68) is kmih_1 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(69) is kmih_2 in component kinetic_parameters (first_order_rate_constant).
C RATES(1) is d/dt Clb2 in component Clb2 (dimensionless).
C RATES(7) is d/dt PClb2 in component PClb2 (dimensionless).
C RATES(8) is d/dt Trim in component Trim (dimensionless).
C RATES(10) is d/dt PTrim in component PTrim (dimensionless).
C RATES(3) is d/dt Mcm_a in component Mcm_a (dimensionless).
C RATES(6) is d/dt Sic1 in component Sic1 (dimensionless).
C RATES(11) is d/dt Mih1_a in component Mih1_a (dimensionless).
C RATES(12) is d/dt IE_a in component IE_a (dimensionless).
C RATES(5) is d/dt Cdc20_a in component Cdc20_a (dimensionless).
C RATES(13) is d/dt Cdc20 in component Cdc20 (dimensionless).
C RATES(4) is d/dt Cdh1_a in component Cdh1_a (dimensionless).
C RATES(9) is d/dt Cln in component Cln (dimensionless).
C RATES(14) is d/dt SBF_a in component SBF_a (dimensionless).
C RATES(15) is d/dt Swe1 in component Swe1 (dimensionless).
C RATES(16) is d/dt PSwe1 in component PSwe1 (dimensionless).
C RATES(17) is d/dt Swe1M in component Swe1M (dimensionless).
C RATES(18) is d/dt PSwe1M in component PSwe1M (dimensionless).
C RATES(19) is d/dt BE in component BE (dimensionless).
C RATES(2) is d/dt M in component M (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.184
      CONSTS(1) = 0.5
      CONSTS(2) = 10
      CONSTS(3) = 0.015
      CONSTS(4) = 0.015
      CONSTS(5) = 1
      CONSTS(6) = 0.1
      CONSTS(7) = 300
      CONSTS(8) = 0.1
      CONSTS(9) = 0.01
      CONSTS(10) = 1
      CONSTS(11) = 3
      STATES(2) = 0.802
      STATES(3) = 0.933
      STATES(4) = 0.993
      STATES(5) = 1.438
      STATES(6) = 0.003
      STATES(7) = 0
      STATES(8) = 0.084
      STATES(9) = 0.054
      STATES(10) = 0
      CONSTS(12) = 0.1
      CONSTS(13) = 0.1
      CONSTS(14) = 0.15
      CONSTS(15) = 1
      CONSTS(16) = 0.1
      STATES(11) = 0.808
      CONSTS(17) = 0.1
      CONSTS(18) = 0.1
      CONSTS(19) = 0.3
      CONSTS(20) = 1
      STATES(12) = 0.522
      CONSTS(21) = 0.01
      CONSTS(22) = 0.01
      CONSTS(23) = 0.04
      CONSTS(24) = 0.1
      CONSTS(25) = 0.001
      CONSTS(26) = 0.001
      CONSTS(27) = 0.25
      CONSTS(28) = 0.1
      CONSTS(29) = 1
      STATES(13) = 1.172
      CONSTS(30) = 0.005
      CONSTS(31) = 0.3
      CONSTS(32) = 0.3
      CONSTS(33) = 0.01
      CONSTS(34) = 0.01
      CONSTS(35) = 35
      CONSTS(36) = 1
      CONSTS(37) = 2
      CONSTS(38) = 10
      CONSTS(39) = 0.1
      CONSTS(40) = 0.1
      STATES(14) = 0.124
      CONSTS(41) = 0.01
      CONSTS(42) = 0.01
      CONSTS(43) = 1
      CONSTS(44) = 1
      CONSTS(45) = 2
      CONSTS(46) = 0
      STATES(15) = 0
      CONSTS(47) = 0.0025
      CONSTS(48) = 0
      CONSTS(49) = 1
      CONSTS(50) = 0.01
      CONSTS(51) = 0.007
      CONSTS(52) = 0.05
      CONSTS(53) = 0.05
      CONSTS(54) = 1
      CONSTS(55) = 0.3
      STATES(16) = 0
      STATES(17) = 0.018
      STATES(18) = 0.013
      CONSTS(56) = 0.05
      STATES(19) = 0
      CONSTS(57) = 0.1
      CONSTS(58) = 0.1
      CONSTS(59) = 0.005
      CONSTS(60) = 1
      CONSTS(61) = 1
      CONSTS(62) = 1
      CONSTS(63) = 1
      CONSTS(64) = 1
      CONSTS(65) = 2
      CONSTS(66) = 0.01
      CONSTS(67) = 0.2
      CONSTS(68) = 5
      CONSTS(69) = 0.5
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(6) = (CONSTS(16)+ CONSTS(8)*(STATES(10)+STATES(8))+ (CONSTS(4)+ CONSTS(5)*STATES(4)+ CONSTS(6)*STATES(5))*(STATES(10)+STATES(8))) - ( (CONSTS(9)+ CONSTS(10)*STATES(9)+ CONSTS(11)*STATES(1))*STATES(6)+ CONSTS(7)*STATES(6)*(STATES(1)+STATES(7)))
      RATES(5) = ( CONSTS(29)*STATES(12)*STATES(13))/(CONSTS(26)+STATES(13)) - (( CONSTS(27)*STATES(5))/(CONSTS(25)+STATES(5))+ CONSTS(28)*STATES(5))
      RATES(13) = (CONSTS(30)+( CONSTS(31)*STATES(1) ** 4.00000)/(CONSTS(32) ** 4.00000+STATES(1) ** 4.00000)+( CONSTS(27)*STATES(5))/(CONSTS(25)+STATES(5))) - (( CONSTS(29)*STATES(12)*STATES(13))/(CONSTS(26)+STATES(13))+ CONSTS(28)*STATES(13))
      RATES(9) =  CONSTS(39)*STATES(14) -  CONSTS(40)*STATES(9)
      RATES(19) =  CONSTS(57)*STATES(9) -  CONSTS(58)*STATES(19)
      RATES(2) =  CONSTS(59)*STATES(2)
      ALGBRC(1) = CONSTS(63) - STATES(3)
      RATES(3) = ( - CONSTS(14)*STATES(3))/(CONSTS(12)+STATES(3))+( CONSTS(15)*ALGBRC(1)*STATES(1))/(CONSTS(13)+ALGBRC(1))
      ALGBRC(2) = CONSTS(62) - STATES(11)
      RATES(11) = ( - CONSTS(19)*STATES(11))/(CONSTS(17)+STATES(11))+( CONSTS(20)*STATES(1)*ALGBRC(2))/(CONSTS(18)+ALGBRC(2))
      ALGBRC(3) = CONSTS(60) - STATES(12)
      RATES(12) = ( - CONSTS(23)*STATES(12))/(CONSTS(21)+STATES(12))+( CONSTS(24)*STATES(1)*ALGBRC(3))/(CONSTS(22)+ALGBRC(3))
      ALGBRC(4) = CONSTS(61) - STATES(4)
      RATES(4) = ( - ( CONSTS(35)*STATES(1)+ CONSTS(37)*STATES(9))*STATES(4))/(CONSTS(33)+STATES(4))+( (CONSTS(36)+ CONSTS(38)*STATES(5))*ALGBRC(4))/(CONSTS(34)+ALGBRC(4))
      ALGBRC(5) = CONSTS(64) - STATES(14)
      RATES(14) = ( - (CONSTS(43)+ CONSTS(45)*STATES(1))*STATES(14))/(CONSTS(41)+STATES(14))+( ( CONSTS(44)*STATES(2)+ CONSTS(46)*STATES(9))*ALGBRC(5))/(CONSTS(42)+ALGBRC(5))
      ALGBRC(6) = TERNRY(STATES(19).LE.0.600000.AND.STATES(1).LT.0.200000, 0.00000, 1.00000)
      RATES(15) = ( CONSTS(47)*STATES(14)+CONSTS(48)+ CONSTS(50)*STATES(17)+( CONSTS(55)*STATES(16))/(CONSTS(53)+STATES(16))) - ( CONSTS(49)*ALGBRC(6)*STATES(15)+( CONSTS(54)*STATES(1)*STATES(15))/(CONSTS(52)+STATES(15))+ CONSTS(51)*STATES(15))
      RATES(16) = ( CONSTS(50)*STATES(18)+( CONSTS(54)*STATES(15)*STATES(1))/(CONSTS(52)+STATES(15))) - ( CONSTS(49)*ALGBRC(6)*STATES(16)+( CONSTS(55)*STATES(16))/(CONSTS(53)+STATES(16))+ CONSTS(51)*STATES(16))
      RATES(17) = ( CONSTS(49)*ALGBRC(6)*STATES(15)+( CONSTS(55)*STATES(18))/(CONSTS(53)+STATES(18))) - ( CONSTS(50)*STATES(17)+( CONSTS(54)*STATES(1)*STATES(17))/(CONSTS(52)+STATES(17))+ CONSTS(51)*STATES(17))
      RATES(18) = ( CONSTS(49)*ALGBRC(6)*STATES(16)+( CONSTS(54)*STATES(1)*STATES(17))/(CONSTS(52)+STATES(17))) - ( CONSTS(50)*STATES(18)+( CONSTS(55)*STATES(18))/(CONSTS(53)+STATES(18))+ CONSTS(56)*STATES(18))
      ALGBRC(9) =  CONSTS(65)*STATES(15)+ CONSTS(66)*STATES(17)+ CONSTS(67)*STATES(16)
      ALGBRC(8) =  CONSTS(68)*STATES(11)+ CONSTS(69)*ALGBRC(2)
      RATES(1) = (( CONSTS(3)*(CONSTS(1)+STATES(3))*STATES(2))/(1.00000+STATES(2)/CONSTS(2))+ ALGBRC(8)*STATES(7)+ CONSTS(8)*STATES(8)+ (CONSTS(9)+ CONSTS(10)*STATES(9)+ CONSTS(11)*STATES(1))*STATES(8)) - ( (CONSTS(4)+ CONSTS(5)*STATES(4)+ CONSTS(6)*STATES(5))*STATES(1)+ ALGBRC(9)*STATES(1)+ CONSTS(7)*STATES(6)*STATES(1))
      RATES(7) = ( ALGBRC(9)*STATES(1)+ CONSTS(8)*STATES(10)+ (CONSTS(9)+ CONSTS(10)*STATES(9)+ CONSTS(11)*STATES(1))*STATES(10)) - ( (CONSTS(4)+ CONSTS(5)*STATES(4)+ CONSTS(6)*STATES(5))*STATES(7)+ ALGBRC(8)*STATES(7)+ CONSTS(7)*STATES(6)*STATES(7))
      RATES(8) = ( CONSTS(7)*STATES(6)*STATES(1)+ ALGBRC(8)*STATES(10)) - ( CONSTS(8)*STATES(8)+ (CONSTS(9)+ CONSTS(10)*STATES(9)+ CONSTS(11)*STATES(1))*STATES(8)+ (CONSTS(4)+ CONSTS(5)*STATES(4)+ CONSTS(6)*STATES(5))*STATES(8)+ ALGBRC(9)*STATES(8))
      RATES(10) = ( ALGBRC(9)*STATES(8)+ CONSTS(7)*STATES(6)*STATES(7)) - ( CONSTS(8)*STATES(10)+ (CONSTS(9)+ CONSTS(10)*STATES(9)+ CONSTS(11)*STATES(1))*STATES(10)+ (CONSTS(4)+ CONSTS(5)*STATES(4)+ CONSTS(6)*STATES(5))*STATES(10)+ ALGBRC(8)*STATES(10))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = CONSTS(63) - STATES(3)
      ALGBRC(2) = CONSTS(62) - STATES(11)
      ALGBRC(3) = CONSTS(60) - STATES(12)
      ALGBRC(4) = CONSTS(61) - STATES(4)
      ALGBRC(5) = CONSTS(64) - STATES(14)
      ALGBRC(6) = TERNRY(STATES(19).LE.0.600000.AND.STATES(1).LT.0.200000, 0.00000, 1.00000)
      ALGBRC(9) =  CONSTS(65)*STATES(15)+ CONSTS(66)*STATES(17)+ CONSTS(67)*STATES(16)
      ALGBRC(8) =  CONSTS(68)*STATES(11)+ CONSTS(69)*ALGBRC(2)
      ALGBRC(7) = STATES(15)+STATES(16)+STATES(17)+STATES(18)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END