C
C There are a total of 5 entries in the algebraic variable array.
C There are a total of 9 entries in each of the rate and state variable arrays.
C There are a total of 31 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is k1AA in component rate_constants (first_order_rate_constant).
C CONSTS(2) is k3 in component rate_constants (first_order_rate_constant).
C CONSTS(3) is kinh in component rate_constants (first_order_rate_constant).
C CONSTS(4) is kcak in component rate_constants (first_order_rate_constant).
C CONSTS(5) is kbPPase in component rate_constants (first_order_rate_constant).
C CONSTS(6) is kfPPase in component rate_constants (first_order_rate_constant).
C CONSTS(7) is khPPase in component rate_constants (first_order_rate_constant).
C CONSTS(8) is kd_antiIE in component rate_constants (first_order_rate_constant).
C CONSTS(9) is ka in component rate_constants (first_order_rate_constant).
C CONSTS(10) is kc in component rate_constants (first_order_rate_constant).
C CONSTS(11) is ke in component rate_constants (first_order_rate_constant).
C CONSTS(12) is kg in component rate_constants (first_order_rate_constant).
C CONSTS(13) is V2_ in component rate_constants (first_order_rate_constant).
C CONSTS(14) is V2__ in component rate_constants (first_order_rate_constant).
C CONSTS(15) is V25_ in component rate_constants (first_order_rate_constant).
C CONSTS(16) is V25__ in component rate_constants (first_order_rate_constant).
C CONSTS(17) is Vwee_ in component rate_constants (first_order_rate_constant).
C CONSTS(18) is Vwee__ in component rate_constants (first_order_rate_constant).
C CONSTS(19) is K_a in component rate_constants (dimensionless).
C CONSTS(20) is K_b in component rate_constants (dimensionless).
C CONSTS(21) is K_c in component rate_constants (dimensionless).
C CONSTS(22) is K_d in component rate_constants (dimensionless).
C CONSTS(23) is K_e in component rate_constants (dimensionless).
C CONSTS(24) is K_f in component rate_constants (dimensionless).
C CONSTS(25) is K_g in component rate_constants (dimensionless).
C CONSTS(26) is K_h in component rate_constants (dimensionless).
C CONSTS(27) is Cdc25_T in component concentration_variables (dimensionless).
C STATES(1) is Cdc25_P in component Cdc25_P (dimensionless).
C STATES(2) is Wee1_P in component Wee1_P (dimensionless).
C CONSTS(28) is Wee1_T in component concentration_variables (dimensionless).
C CONSTS(29) is UbE_T in component concentration_variables (dimensionless).
C STATES(3) is UbE in component UbE (dimensionless).
C ALGBRC(1) is k25 in component rate_constants (first_order_rate_constant).
C ALGBRC(3) is kwee in component rate_constants (first_order_rate_constant).
C ALGBRC(4) is k2 in component rate_constants (first_order_rate_constant).
C CONSTS(30) is Cdc2_T in component concentration_variables (dimensionless).
C CONSTS(31) is IE_T in component concentration_variables (dimensionless).
C STATES(4) is cyclin in component cyclin (dimensionless).
C ALGBRC(5) is Cdc2 in component Cdc2 (dimensionless).
C STATES(5) is cyclin_Cdc2 in component cyclin_Cdc2 (dimensionless).
C STATES(6) is MPF_active in component MPF_active (dimensionless).
C STATES(7) is Tyr15P_dimer in component Tyr15P_dimer (dimensionless).
C STATES(8) is MPF_pre in component MPF_pre (dimensionless).
C ALGBRC(2) is cyclin_T in component cyclin_T (dimensionless).
C STATES(9) is IE_P in component IE_P (dimensionless).
C RATES(4) is d/dt cyclin in component cyclin (dimensionless).
C RATES(5) is d/dt cyclin_Cdc2 in component cyclin_Cdc2 (dimensionless).
C RATES(7) is d/dt Tyr15P_dimer in component Tyr15P_dimer (dimensionless).
C RATES(8) is d/dt MPF_pre in component MPF_pre (dimensionless).
C RATES(6) is d/dt MPF_active in component MPF_active (dimensionless).
C RATES(1) is d/dt Cdc25_P in component Cdc25_P (dimensionless).
C RATES(2) is d/dt Wee1_P in component Wee1_P (dimensionless).
C RATES(9) is d/dt IE_P in component IE_P (dimensionless).
C RATES(3) is d/dt UbE in component UbE (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 0.01
      CONSTS(3) = 0.025
      CONSTS(4) = 0.25
      CONSTS(5) = 0.125
      CONSTS(6) = 0.1
      CONSTS(7) = 0.087
      CONSTS(8) = 0.095
      CONSTS(9) = 0.01
      CONSTS(10) = 0.1
      CONSTS(11) = 0.0133
      CONSTS(12) = 0.0065
      CONSTS(13) = 0.015
      CONSTS(14) = 1
      CONSTS(15) = 0.1
      CONSTS(16) = 2
      CONSTS(17) = 0.1
      CONSTS(18) = 1
      CONSTS(19) = 0.1
      CONSTS(20) = 0.1
      CONSTS(21) = 0.01
      CONSTS(22) = 0.01
      CONSTS(23) = 0.3
      CONSTS(24) = 0.3
      CONSTS(25) = 0.01
      CONSTS(26) = 0.01
      CONSTS(27) = 1
      STATES(1) = 0
      STATES(2) = 0
      CONSTS(28) = 1
      CONSTS(29) = 1
      STATES(3) = 0
      CONSTS(30) = 100
      CONSTS(31) = 1
      STATES(4) = 1
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(9)*STATES(6)*(CONSTS(27) - STATES(1)))/((CONSTS(19)+CONSTS(27)) - STATES(1)) - ( CONSTS(5)*STATES(1))/(CONSTS(20)+STATES(1))
      RATES(2) = ( CONSTS(11)*STATES(6)*(CONSTS(28) - STATES(2)))/((CONSTS(23)+CONSTS(28)) - STATES(2)) - ( CONSTS(6)*STATES(2))/(CONSTS(24)+STATES(2))
      RATES(9) = ( CONSTS(12)*STATES(6)*(CONSTS(31) - STATES(9)))/((CONSTS(25)+CONSTS(31)) - STATES(9)) - ( CONSTS(7)*STATES(9))/(CONSTS(26)+STATES(9))
      RATES(3) = ( CONSTS(10)*STATES(9)*(CONSTS(29) - STATES(3)))/((CONSTS(21)+CONSTS(29)) - STATES(3)) - ( CONSTS(8)*STATES(3))/(CONSTS(22)+STATES(3))
      ALGBRC(1) =  CONSTS(15)*(CONSTS(27) - STATES(1))+ CONSTS(16)*STATES(1)
      ALGBRC(3) =  CONSTS(17)*STATES(2)+ CONSTS(18)*(CONSTS(28) - STATES(2))
      ALGBRC(4) =  CONSTS(13)*(CONSTS(29) - STATES(3))+ CONSTS(14)*STATES(3)
      RATES(7) = ( ALGBRC(3)*STATES(5) -  (ALGBRC(1)+CONSTS(4)+ALGBRC(4))*STATES(7))+ CONSTS(3)*STATES(8)
      RATES(8) = ( ALGBRC(3)*STATES(6) -  (CONSTS(3)+ALGBRC(1)+ALGBRC(4))*STATES(8))+ CONSTS(4)*STATES(7)
      RATES(6) = ( CONSTS(4)*STATES(5) -  (CONSTS(3)+ALGBRC(3)+ALGBRC(4))*STATES(6))+ ALGBRC(1)*STATES(8)
      ALGBRC(5) = CONSTS(30) - (STATES(5)+STATES(6)+STATES(8)+STATES(7))
      RATES(4) = (CONSTS(1) -  ALGBRC(4)*STATES(4)) -  CONSTS(2)*STATES(4)*ALGBRC(5)
      RATES(5) = ( CONSTS(3)*STATES(6) -  (ALGBRC(3)+CONSTS(4)+ALGBRC(4))*STATES(5))+ ALGBRC(1)*STATES(7)+ CONSTS(2)*STATES(4)*ALGBRC(5)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(15)*(CONSTS(27) - STATES(1))+ CONSTS(16)*STATES(1)
      ALGBRC(3) =  CONSTS(17)*STATES(2)+ CONSTS(18)*(CONSTS(28) - STATES(2))
      ALGBRC(4) =  CONSTS(13)*(CONSTS(29) - STATES(3))+ CONSTS(14)*STATES(3)
      ALGBRC(5) = CONSTS(30) - (STATES(5)+STATES(6)+STATES(8)+STATES(7))
      ALGBRC(2) = STATES(4)+STATES(5)+STATES(6)+STATES(8)+STATES(7)
      RETURN
      END