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