C
C There are a total of 8 entries in the algebraic variable array.
C There are a total of 7 entries in each of the rate and state variable arrays.
C There are a total of 35 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Y in component Y (dimensionless).
C CONSTS(1) is k1 in component Y (first_order_rate_constant).
C ALGBRC(7) is k2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is k3 in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(1) is C in component C (dimensionless).
C STATES(2) is M in component M (dimensionless).
C ALGBRC(8) is kwee in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(6) is kcdc25 in component kinetic_parameters (first_order_rate_constant).
C STATES(3) is preMPF in component preMPF (dimensionless).
C STATES(4) is Cdc25P in component Cdc25P (dimensionless).
C CONSTS(3) is k25 in component Cdc25P (first_order_rate_constant).
C CONSTS(4) is Km25 in component Cdc25P (dimensionless).
C CONSTS(5) is k25r in component Cdc25P (first_order_rate_constant).
C CONSTS(34) is k25ro in component Cdc25P (first_order_rate_constant).
C CONSTS(6) is Km25r in component Cdc25P (dimensionless).
C CONSTS(33) is p in component p (dimensionless).
C ALGBRC(3) is Cdc25 in component Cdc25 (dimensionless).
C STATES(5) is Wee1 in component Wee1 (dimensionless).
C CONSTS(7) is kw in component Wee1 (first_order_rate_constant).
C CONSTS(8) is Kmw in component Wee1 (dimensionless).
C CONSTS(9) is kwr in component Wee1 (first_order_rate_constant).
C CONSTS(35) is kwro in component Wee1 (first_order_rate_constant).
C CONSTS(10) is Kmwr in component Wee1 (dimensionless).
C ALGBRC(4) is Wee1P in component Wee1P (dimensionless).
C STATES(6) is IEP in component IEP (dimensionless).
C CONSTS(11) is kie in component IEP (first_order_rate_constant).
C CONSTS(12) is Kmie in component IEP (dimensionless).
C CONSTS(13) is kier in component IEP (first_order_rate_constant).
C CONSTS(14) is Kmier in component IEP (dimensionless).
C ALGBRC(2) is IE in component IE (dimensionless).
C STATES(7) is APC_ in component APC_ (dimensionless).
C CONSTS(15) is kap in component APC_ (first_order_rate_constant).
C CONSTS(16) is Kmap in component APC_ (dimensionless).
C CONSTS(17) is kapr in component APC_ (first_order_rate_constant).
C CONSTS(18) is Kmapr in component APC_ (dimensionless).
C ALGBRC(5) is APC in component APC (dimensionless).
C CONSTS(19) is Cdc2tot in component C (dimensionless).
C CONSTS(20) is Cdc25tot in component Cdc25 (dimensionless).
C CONSTS(21) is Wee1tot in component Wee1P (dimensionless).
C CONSTS(22) is IEtot in component IE (dimensionless).
C CONSTS(23) is APCtot in component APC (dimensionless).
C CONSTS(24) is R in component p (dimensionless).
C CONSTS(25) is s in component p (dimensionless).
C CONSTS(26) is q in component p (dimensionless).
C CONSTS(27) is V2 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(28) is V2_ in component kinetic_parameters (first_order_rate_constant).
C CONSTS(29) is V25 in component kinetic_parameters (first_order_rate_constant).
C CONSTS(30) is V25_ in component kinetic_parameters (first_order_rate_constant).
C CONSTS(31) is Vwee in component kinetic_parameters (first_order_rate_constant).
C CONSTS(32) is Vwee_ in component kinetic_parameters (first_order_rate_constant).
C RATES(1) is d/dt Y in component Y (dimensionless).
C RATES(2) is d/dt M in component M (dimensionless).
C RATES(3) is d/dt preMPF in component preMPF (dimensionless).
C RATES(4) is d/dt Cdc25P in component Cdc25P (dimensionless).
C RATES(5) is d/dt Wee1 in component Wee1 (dimensionless).
C RATES(6) is d/dt IEP in component IEP (dimensionless).
C RATES(7) is d/dt APC_ in component APC_ (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.01577552
      CONSTS(1) = 0.024
      CONSTS(2) = 1.15
      STATES(2) = 0.0142077
      STATES(3) = 0.02541723
      STATES(4) = 0.4844362
      CONSTS(3) = 18
      CONSTS(4) = 0.1
      CONSTS(5) = 0.8
      CONSTS(6) = 1
      STATES(5) = 0.5155638
      CONSTS(7) = 18
      CONSTS(8) = 0.1
      CONSTS(9) = 0.8
      CONSTS(10) = 1
      STATES(6) = 0.002287817
      CONSTS(11) = 4.5
      CONSTS(12) = 0.01
      CONSTS(13) = 0.34
      CONSTS(14) = 0.01
      STATES(7) = 0.5051103
      CONSTS(15) = 0.3
      CONSTS(16) = 0.01
      CONSTS(17) = 0.3
      CONSTS(18) = 1
      CONSTS(19) = 1
      CONSTS(20) = 1
      CONSTS(21) = 1
      CONSTS(22) = 1
      CONSTS(23) = 1
      CONSTS(24) = 1
      CONSTS(25) = 0.021
      CONSTS(26) = 1
      CONSTS(27) = 0.01
      CONSTS(28) = 0.6
      CONSTS(29) = 0.04
      CONSTS(30) = 0.4
      CONSTS(31) = 0.025
      CONSTS(32) = 2.5
      CONSTS(33) =  CONSTS(25)*CONSTS(24)+CONSTS(26)
      CONSTS(34) = CONSTS(5)/CONSTS(33)
      CONSTS(35) = CONSTS(9)/CONSTS(33)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = CONSTS(22) - STATES(6)
      RATES(6) = ( CONSTS(11)*STATES(2)*ALGBRC(2))/(CONSTS(12)+ALGBRC(2)) - ( CONSTS(13)*STATES(6))/(CONSTS(14)+STATES(6))
      ALGBRC(3) = CONSTS(20) - STATES(4)
      RATES(4) = ( CONSTS(3)*STATES(2)*ALGBRC(3))/(CONSTS(4)+ALGBRC(3)) - ( CONSTS(5)*STATES(4))/(CONSTS(6)+STATES(4))
      ALGBRC(4) = CONSTS(21) - STATES(5)
      RATES(5) = ( - CONSTS(7)*STATES(2)*STATES(5))/(CONSTS(8)+STATES(5))+( CONSTS(9)*ALGBRC(4))/(CONSTS(10)+ALGBRC(4))
      ALGBRC(5) = CONSTS(23) - STATES(7)
      RATES(7) = ( CONSTS(15)*STATES(6)*ALGBRC(5))/(CONSTS(16)+ALGBRC(5)) - ( CONSTS(17)*STATES(7))/(CONSTS(18)+STATES(7))
      ALGBRC(7) =  CONSTS(27)*ALGBRC(5)+ CONSTS(28)*STATES(7)
      ALGBRC(1) = CONSTS(19) - (STATES(2)+STATES(3))
      RATES(1) = CONSTS(1) - ( ALGBRC(7)*STATES(1)+ CONSTS(2)*STATES(1)*ALGBRC(1))
      ALGBRC(8) =  CONSTS(31)*ALGBRC(4)+ CONSTS(32)*STATES(5)
      ALGBRC(6) =  CONSTS(29)*ALGBRC(3)+ CONSTS(30)*STATES(4)
      RATES(2) = ( ALGBRC(6)*STATES(3)+ CONSTS(2)*STATES(1)*ALGBRC(1)) - ( ALGBRC(7)*STATES(2)+ ALGBRC(8)*STATES(2))
      RATES(3) =  ALGBRC(8)*STATES(2) - ( ALGBRC(7)*STATES(3)+ ALGBRC(6)*STATES(3))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = CONSTS(22) - STATES(6)
      ALGBRC(3) = CONSTS(20) - STATES(4)
      ALGBRC(4) = CONSTS(21) - STATES(5)
      ALGBRC(5) = CONSTS(23) - STATES(7)
      ALGBRC(7) =  CONSTS(27)*ALGBRC(5)+ CONSTS(28)*STATES(7)
      ALGBRC(1) = CONSTS(19) - (STATES(2)+STATES(3))
      ALGBRC(8) =  CONSTS(31)*ALGBRC(4)+ CONSTS(32)*STATES(5)
      ALGBRC(6) =  CONSTS(29)*ALGBRC(3)+ CONSTS(30)*STATES(4)
      RETURN
      END