C
C There are a total of 10 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 45 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Cdc13_Cdc2 in component Cdc13_Cdc2 (dimensionless).
C CONSTS(1) is k1 in component Cdc13_Cdc2 (first_order_rate_constant).
C ALGBRC(4) is kwee in component kwee (first_order_rate_constant).
C CONSTS(2) is ki in component parameters (first_order_rate_constant).
C ALGBRC(3) is kcdc25 in component kcdc25 (first_order_rate_constant).
C CONSTS(3) is kir in component parameters (first_order_rate_constant).
C ALGBRC(1) is k2 in component k2 (first_order_rate_constant).
C ALGBRC(6) is k4 in component k4 (first_order_rate_constant).
C STATES(2) is mass in component mass (dimensionless).
C STATES(3) is Cdc13_P_Cdc2 in component Cdc13_P_Cdc2 (dimensionless).
C STATES(4) is Rum1_Cdc13_Cdc2 in component Rum1_Cdc13_Cdc2 (dimensionless).
C STATES(5) is Rum1 in component Rum1 (dimensionless).
C ALGBRC(2) is k2c in component k2c (first_order_rate_constant).
C CONSTS(4) is k3 in component Rum1 (first_order_rate_constant).
C STATES(6) is Cdc25P in component Cdc25P (dimensionless).
C CONSTS(5) is k25 in component Cdc25P (first_order_rate_constant).
C ALGBRC(8) is k25r in component Cdc25P (first_order_rate_constant).
C CONSTS(6) is k25r_ in component Cdc25P (first_order_rate_constant).
C CONSTS(7) is J25 in component Cdc25P (dimensionless).
C CONSTS(8) is J25r in component Cdc25P (dimensionless).
C ALGBRC(7) is ks in component parameters (first_order_rate_constant).
C ALGBRC(5) is MPF in component MPF (dimensionless).
C STATES(7) is Wee1 in component Wee1 (dimensionless).
C CONSTS(9) is kw in component Wee1 (first_order_rate_constant).
C ALGBRC(9) is kwr in component Wee1 (first_order_rate_constant).
C CONSTS(10) is kwr_ in component Wee1 (first_order_rate_constant).
C CONSTS(11) is Jw in component Wee1 (dimensionless).
C CONSTS(12) is Jwr in component Wee1 (dimensionless).
C STATES(8) is Mik1 in component Mik1 (dimensionless).
C CONSTS(13) is km in component Mik1 (first_order_rate_constant).
C ALGBRC(10) is kmr in component Mik1 (first_order_rate_constant).
C CONSTS(14) is kmr_ in component Mik1 (first_order_rate_constant).
C CONSTS(15) is Jm in component Mik1 (dimensionless).
C CONSTS(16) is Jmr in component Mik1 (dimensionless).
C STATES(9) is AAE_total in component AAE_total (dimensionless).
C CONSTS(17) is kas in component AAE_total (first_order_rate_constant).
C CONSTS(18) is kad in component parameters (first_order_rate_constant).
C STATES(10) is AAE in component AAE (dimensionless).
C CONSTS(19) is kaa in component AAE (first_order_rate_constant).
C CONSTS(20) is kaa_ in component AAE (first_order_rate_constant).
C CONSTS(45) is kai in component AAE (first_order_rate_constant).
C CONSTS(21) is kai_ in component AAE (first_order_rate_constant).
C CONSTS(22) is kx in component AAE (first_order_rate_constant).
C CONSTS(23) is Jaa in component AAE (dimensionless).
C CONSTS(24) is Jai in component AAE (dimensionless).
C STATES(11) is APC in component APC (dimensionless).
C CONSTS(25) is kapr in component APC (first_order_rate_constant).
C CONSTS(26) is kapr_ in component APC (first_order_rate_constant).
C CONSTS(27) is kap in component APC (first_order_rate_constant).
C CONSTS(28) is Japr in component APC (dimensionless).
C CONSTS(29) is Jap in component APC (dimensionless).
C CONSTS(30) is Puc1 in component parameters (dimensionless).
C CONSTS(31) is mu in component mass (first_order_rate_constant).
C CONSTS(32) is V2 in component k2 (first_order_rate_constant).
C CONSTS(33) is V2_ in component k2 (first_order_rate_constant).
C CONSTS(34) is V2c in component k2c (first_order_rate_constant).
C CONSTS(35) is V2c_ in component k2c (first_order_rate_constant).
C CONSTS(36) is k4_ in component k4 (first_order_rate_constant).
C CONSTS(37) is k4__ in component k4 (first_order_rate_constant).
C CONSTS(38) is V25 in component kcdc25 (first_order_rate_constant).
C CONSTS(39) is V25_ in component kcdc25 (first_order_rate_constant).
C CONSTS(40) is Vwee in component kwee (first_order_rate_constant).
C CONSTS(41) is Vwee_ in component kwee (first_order_rate_constant).
C CONSTS(42) is Vmik in component kwee (first_order_rate_constant).
C CONSTS(43) is Vmik_ in component kwee (first_order_rate_constant).
C CONSTS(44) is alpha in component MPF (dimensionless).
C RATES(1) is d/dt Cdc13_Cdc2 in component Cdc13_Cdc2 (dimensionless).
C RATES(3) is d/dt Cdc13_P_Cdc2 in component Cdc13_P_Cdc2 (dimensionless).
C RATES(4) is d/dt Rum1_Cdc13_Cdc2 in component Rum1_Cdc13_Cdc2 (dimensionless).
C RATES(5) is d/dt Rum1 in component Rum1 (dimensionless).
C RATES(6) is d/dt Cdc25P in component Cdc25P (dimensionless).
C RATES(7) is d/dt Wee1 in component Wee1 (dimensionless).
C RATES(8) is d/dt Mik1 in component Mik1 (dimensionless).
C RATES(9) is d/dt AAE_total in component AAE_total (dimensionless).
C RATES(10) is d/dt AAE in component AAE (dimensionless).
C RATES(11) is d/dt APC in component APC (dimensionless).
C RATES(2) is d/dt mass in component mass (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 1
      CONSTS(1) = 0.03
      CONSTS(2) = 200
      CONSTS(3) = 1
      STATES(2) = 1
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      CONSTS(4) = 0.15
      STATES(6) = 0
      CONSTS(5) = 0.5
      CONSTS(6) = 0.2
      CONSTS(7) = 0.2
      CONSTS(8) = 0.2
      STATES(7) = 0
      CONSTS(9) = 0.5
      CONSTS(10) = 0.2
      CONSTS(11) = 0.2
      CONSTS(12) = 0.2
      STATES(8) = 0
      CONSTS(13) = 0.1
      CONSTS(14) = 0
      CONSTS(15) = 0.2
      CONSTS(16) = 0.2
      STATES(9) = 2
      CONSTS(17) = 0.25
      CONSTS(18) = 0.1
      STATES(10) = 2
      CONSTS(19) = 0.001
      CONSTS(20) = 1
      CONSTS(21) = 0.25
      CONSTS(22) = 0
      CONSTS(23) = 0.1
      CONSTS(24) = 0.1
      STATES(11) = 0
      CONSTS(25) = 0.04
      CONSTS(26) = 3
      CONSTS(27) = 4
      CONSTS(28) = 0.01
      CONSTS(29) = 0.01
      CONSTS(30) = 0.013
      CONSTS(31) = 0.005776
      CONSTS(32) = 0.03
      CONSTS(33) = 1
      CONSTS(34) = 0.03
      CONSTS(35) = 0.16
      CONSTS(36) = 0.15
      CONSTS(37) = 20
      CONSTS(38) = 0.01
      CONSTS(39) = 0.4
      CONSTS(40) = 0.01
      CONSTS(41) = 0.93
      CONSTS(42) = 0.002
      CONSTS(43) = 0.2
      CONSTS(44) = 0.1
      CONSTS(45) = CONSTS(21)+CONSTS(22)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(2) =  CONSTS(31)*STATES(2)
      ALGBRC(4) =  CONSTS(40)*( (1.00000 - STATES(7))*1.00000+ CONSTS(41)*STATES(7)+ CONSTS(42)*(1.00000 - STATES(8)))*1.00000+ CONSTS(43)*STATES(8)
      ALGBRC(3) =  CONSTS(38)*(1.00000 - STATES(6))+ CONSTS(39)*STATES(6)
      ALGBRC(1) =  CONSTS(32)*(1.00000 - STATES(11))+ CONSTS(33)*STATES(11)
      RATES(3) =  ALGBRC(4)*STATES(1) - ( ALGBRC(3)*STATES(3)+ ALGBRC(1)*STATES(3))
      ALGBRC(5) = STATES(1)+ CONSTS(44)*STATES(3)
      RATES(9) =  CONSTS(17)*ALGBRC(5) -  CONSTS(18)*STATES(9)
      RATES(10) = ( (CONSTS(19)+ CONSTS(20)*ALGBRC(5))*(STATES(9) - STATES(10)))/((CONSTS(23)+STATES(9)) - STATES(10)) - (( CONSTS(45)*STATES(10))/(CONSTS(24)+STATES(10))+ CONSTS(18)*STATES(10))
      RATES(11) = ( (CONSTS(25)+ CONSTS(26)*STATES(10))*(1.00000 - STATES(11)))/((CONSTS(28)+1.00000) - STATES(11)) - ( CONSTS(27)*( CONSTS(30)*STATES(2)+ALGBRC(5))*STATES(11))/(CONSTS(29)+STATES(11))
      ALGBRC(6) = CONSTS(36)+ CONSTS(37)*( CONSTS(30)*STATES(2)+ALGBRC(5))
      RATES(1) = ( CONSTS(1)*STATES(2)+ ALGBRC(3)*STATES(3)+ STATES(4)*(CONSTS(3)+ALGBRC(6))) - ( ALGBRC(4)*STATES(1)+ CONSTS(2)*STATES(1)*STATES(5)+ ALGBRC(1)*STATES(1))
      ALGBRC(2) =  CONSTS(34)*(1.00000 - STATES(11))+ CONSTS(35)*STATES(11)
      RATES(4) =  CONSTS(2)*STATES(5)*STATES(1) -  STATES(4)*(ALGBRC(6)+ALGBRC(2)+CONSTS(3))
      RATES(5) = (CONSTS(4)+ STATES(4)*(CONSTS(3)+ALGBRC(2))) - ( ALGBRC(6)*STATES(5)+ CONSTS(2)*STATES(1)*STATES(5))
      ALGBRC(7) = TERNRY(STATES(11).LE.0.200000, 0.500000, 0.00000)
      ALGBRC(8) = CONSTS(6)+ALGBRC(7)
      RATES(6) = ( CONSTS(5)*ALGBRC(5)*(1.00000 - STATES(6)))/((CONSTS(7)+1.00000) - STATES(6)) - ( ALGBRC(8)*STATES(6))/(CONSTS(8)+STATES(6))
      ALGBRC(9) = CONSTS(10)+ALGBRC(7)
      RATES(7) = ( ALGBRC(9)*(1.00000 - STATES(7)))/((CONSTS(12)+1.00000) - STATES(7)) - ( CONSTS(9)*ALGBRC(5)*STATES(7))/(CONSTS(11)+STATES(7))
      ALGBRC(10) = CONSTS(14)+ALGBRC(7)
      RATES(8) = ( ALGBRC(10)*(1.00000 - STATES(8)))/((CONSTS(16)+1.00000) - STATES(8)) - ( CONSTS(13)*STATES(8))/(CONSTS(15)+STATES(8))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) =  CONSTS(40)*( (1.00000 - STATES(7))*1.00000+ CONSTS(41)*STATES(7)+ CONSTS(42)*(1.00000 - STATES(8)))*1.00000+ CONSTS(43)*STATES(8)
      ALGBRC(3) =  CONSTS(38)*(1.00000 - STATES(6))+ CONSTS(39)*STATES(6)
      ALGBRC(1) =  CONSTS(32)*(1.00000 - STATES(11))+ CONSTS(33)*STATES(11)
      ALGBRC(5) = STATES(1)+ CONSTS(44)*STATES(3)
      ALGBRC(6) = CONSTS(36)+ CONSTS(37)*( CONSTS(30)*STATES(2)+ALGBRC(5))
      ALGBRC(2) =  CONSTS(34)*(1.00000 - STATES(11))+ CONSTS(35)*STATES(11)
      ALGBRC(7) = TERNRY(STATES(11).LE.0.200000, 0.500000, 0.00000)
      ALGBRC(8) = CONSTS(6)+ALGBRC(7)
      ALGBRC(9) = CONSTS(10)+ALGBRC(7)
      ALGBRC(10) = CONSTS(14)+ALGBRC(7)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END