C
C There are a total of 6 entries in the algebraic variable array.
C There are a total of 13 entries in each of the rate and state variable arrays.
C There are a total of 23 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Cdk2_CycE in component Cdk2_CycE (dimensionless).
C CONSTS(1) is kwee in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is kon in component kinetic_parameters (first_order_rate_constant).
C CONSTS(3) is k25A in component kinetic_parameters (first_order_rate_constant).
C CONSTS(4) is koff in component kinetic_parameters (first_order_rate_constant).
C CONSTS(5) is kassoc in component kinetic_parameters (first_order_rate_constant).
C CONSTS(6) is kdissoc in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(5) is phi in component kinetic_parameters (dimensionless).
C STATES(2) is Wee1_a in component Wee1_a (dimensionless).
C STATES(3) is PCdk2_CycE in component PCdk2_CycE (dimensionless).
C STATES(4) is Xic in component Xic (dimensionless).
C STATES(5) is Cdk2_CycErem in component Cdk2_CycErem (dimensionless).
C STATES(6) is Xic_Cdk2_CycE in component Xic_Cdk2_CycE (dimensionless).
C STATES(7) is PCdk2_CycErem in component PCdk2_CycErem (dimensionless).
C STATES(8) is Xic_PCdk2_CycE in component Xic_PCdk2_CycE (dimensionless).
C CONSTS(7) is Jwact in component Wee1_a (dimensionless).
C CONSTS(8) is Jwinact in component Wee1_a (dimensionless).
C CONSTS(9) is Wee1_total in component Wee1_a (dimensionless).
C CONSTS(10) is kwact in component Wee1_a (first_order_rate_constant).
C CONSTS(11) is kwinact in component Wee1_a (first_order_rate_constant).
C STATES(9) is Kin_a in component Kin_a (dimensionless).
C CONSTS(12) is Jiact in component Kin_a (dimensionless).
C CONSTS(13) is Jiinact in component Kin_a (dimensionless).
C CONSTS(14) is kiact in component Kin_a (first_order_rate_constant).
C CONSTS(15) is kiinact in component Kin_a (first_order_rate_constant).
C CONSTS(16) is kedeg in component kinetic_parameters (first_order_rate_constant).
C CONSTS(17) is kxdeg in component kinetic_parameters (first_order_rate_constant).
C STATES(10) is Deg_a in component Deg_a (dimensionless).
C STATES(11) is Xic_Cdk2_CycErem in component Xic_Cdk2_CycErem (dimensionless).
C STATES(12) is Xic_PCdk2_CycErem in component Xic_PCdk2_CycErem (dimensionless).
C ALGBRC(6) is Heav in component Deg_a (dimensionless).
C CONSTS(18) is kdact in component Deg_a (first_order_rate_constant).
C CONSTS(19) is theta in component Deg_a (dimensionless).
C ALGBRC(2) is x in component Deg_a (dimensionless).
C CONSTS(20) is Deg in component kinetic_parameters (dimensionless).
C STATES(13) is Xicrem in component Xicrem (dimensionless).
C ALGBRC(3) is Cyc_total in component Cyc_total (dimensionless).
C ALGBRC(4) is Xic_total in component Xic_total (dimensionless).
C CONSTS(21) is epsilon in component kinetic_parameters (dimensionless).
C ALGBRC(1) is pool in component kinetic_parameters (dimensionless).
C CONSTS(22) is n in component kinetic_parameters (dimensionless).
C CONSTS(23) is L in component kinetic_parameters (dimensionless).
C RATES(1) is d/dt Cdk2_CycE in component Cdk2_CycE (dimensionless).
C RATES(3) is d/dt PCdk2_CycE in component PCdk2_CycE (dimensionless).
C RATES(2) is d/dt Wee1_a in component Wee1_a (dimensionless).
C RATES(9) is d/dt Kin_a in component Kin_a (dimensionless).
C RATES(5) is d/dt Cdk2_CycErem in component Cdk2_CycErem (dimensionless).
C RATES(7) is d/dt PCdk2_CycErem in component PCdk2_CycErem (dimensionless).
C RATES(10) is d/dt Deg_a in component Deg_a (dimensionless).
C RATES(4) is d/dt Xic in component Xic (dimensionless).
C RATES(6) is d/dt Xic_Cdk2_CycE in component Xic_Cdk2_CycE (dimensionless).
C RATES(8) is d/dt Xic_PCdk2_CycE in component Xic_PCdk2_CycE (dimensionless).
C RATES(11) is d/dt Xic_Cdk2_CycErem in component Xic_Cdk2_CycErem (dimensionless).
C RATES(12) is d/dt Xic_PCdk2_CycErem in component Xic_PCdk2_CycErem (dimensionless).
C RATES(13) is d/dt Xicrem in component Xicrem (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.06
      CONSTS(1) = 1.5
      CONSTS(2) = 0.02
      CONSTS(3) = 0.1
      CONSTS(4) = 0.0001
      CONSTS(5) = 0.1
      CONSTS(6) = 0.001
      STATES(2) = 1.02
      STATES(3) = 0.94
      STATES(4) = 3
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      CONSTS(7) = 0.01
      CONSTS(8) = 0.01
      CONSTS(9) = 8
      CONSTS(10) = 0.75
      CONSTS(11) = 1.5
      STATES(9) = 0.6
      CONSTS(12) = 0.01
      CONSTS(13) = 0.01
      CONSTS(14) = 0.15
      CONSTS(15) = 0.6
      CONSTS(16) = 0.017
      CONSTS(17) = 0.01
      STATES(10) = 0
      STATES(11) = 0
      STATES(12) = 0
      CONSTS(18) = 0.023
      CONSTS(19) = 0.3
      CONSTS(20) = 0
      STATES(13) = 0
      CONSTS(21) = 0.001
      CONSTS(22) = 4
      CONSTS(23) = 0.4
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(2) = ( CONSTS(10)*(CONSTS(9) - STATES(2)))/((CONSTS(7)+CONSTS(9)) - STATES(2)) - ( CONSTS(11)*STATES(9)*STATES(2))/(CONSTS(8)+STATES(2))
      RATES(9) = ( CONSTS(14)*(1.00000 - STATES(9)))/((CONSTS(12)+1.00000) - STATES(9)) - ( CONSTS(15)*STATES(1)*STATES(9))/(CONSTS(13)+STATES(9))
      RATES(4) =  CONSTS(6)*(STATES(6)+STATES(8)+STATES(11)+STATES(12)) -  CONSTS(5)*STATES(4)*(STATES(1)+STATES(3)+STATES(5)+STATES(7))
      RATES(13) = ( CONSTS(16)*STATES(10)*STATES(11)+ CONSTS(16)*STATES(10)*STATES(12)) -  CONSTS(17)*STATES(13)
      ALGBRC(1) = STATES(5)+STATES(7)+STATES(11)+STATES(12)
      ALGBRC(5) = (CONSTS(21)+ALGBRC(1) ** CONSTS(22))/(CONSTS(23) ** CONSTS(22)+ALGBRC(1) ** CONSTS(22))
      RATES(1) = ( CONSTS(3)*STATES(3)+ CONSTS(4)*STATES(5)+ CONSTS(6)*STATES(6)) - ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(2)*ALGBRC(5)*STATES(1)+ CONSTS(5)*STATES(4)*STATES(1))
      RATES(3) = ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(4)*STATES(7)+ CONSTS(6)*STATES(8)) - ( CONSTS(3)*STATES(3)+ CONSTS(2)*ALGBRC(5)*STATES(3)+ CONSTS(5)*STATES(4)*STATES(3))
      RATES(5) = ( CONSTS(2)*ALGBRC(5)*STATES(1)+ CONSTS(17)*STATES(11)+ CONSTS(6)*STATES(11)) - ( CONSTS(4)*STATES(5)+ CONSTS(16)*STATES(10)*STATES(5)+ CONSTS(5)*STATES(4)*STATES(5))
      RATES(7) = ( CONSTS(2)*ALGBRC(5)*STATES(3)+ CONSTS(17)*STATES(12)+ CONSTS(6)*STATES(12)) - ( CONSTS(4)*STATES(7)+ CONSTS(16)*STATES(10)*STATES(7)+ CONSTS(5)*STATES(4)*STATES(7))
      ALGBRC(2) = STATES(5) - CONSTS(19)
      ALGBRC(6) = TERNRY(ALGBRC(2).LT.0.00000, 0.00000, 1.00000)
      RATES(10) =  CONSTS(18)*ALGBRC(6)
      RATES(6) = ( CONSTS(5)*STATES(4)*STATES(1)+ CONSTS(4)*STATES(11)+ CONSTS(3)*STATES(8)) - ( CONSTS(6)*STATES(6)+ CONSTS(2)*ALGBRC(5)*STATES(6)+ CONSTS(1)*STATES(2)*STATES(6))
      RATES(8) = ( CONSTS(5)*STATES(4)*STATES(3)+ CONSTS(4)*STATES(12)+ CONSTS(1)*STATES(2)*STATES(6)) - ( CONSTS(6)*STATES(8)+ CONSTS(2)*ALGBRC(5)*STATES(8)+ CONSTS(3)*STATES(8))
      RATES(11) = ( CONSTS(5)*STATES(4)*STATES(5)+ CONSTS(2)*ALGBRC(5)*STATES(6)) - ( CONSTS(6)*STATES(11)+ CONSTS(4)*STATES(11)+ CONSTS(16)*STATES(10)*STATES(11)+ CONSTS(17)*STATES(11))
      RATES(12) = ( CONSTS(5)*STATES(4)*STATES(7)+ CONSTS(2)*ALGBRC(5)*STATES(8)) - ( CONSTS(6)*STATES(12)+ CONSTS(4)*STATES(12)+ CONSTS(16)*CONSTS(20)*STATES(12)+ CONSTS(17)*STATES(12))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = STATES(5)+STATES(7)+STATES(11)+STATES(12)
      ALGBRC(5) = (CONSTS(21)+ALGBRC(1) ** CONSTS(22))/(CONSTS(23) ** CONSTS(22)+ALGBRC(1) ** CONSTS(22))
      ALGBRC(2) = STATES(5) - CONSTS(19)
      ALGBRC(6) = TERNRY(ALGBRC(2).LT.0.00000, 0.00000, 1.00000)
      ALGBRC(3) = STATES(8)+STATES(6)+STATES(12)+STATES(11)+STATES(7)+STATES(5)+STATES(1)+STATES(3)
      ALGBRC(4) = (STATES(4)+STATES(8)+STATES(6)+STATES(12)+STATES(11)+STATES(13))/3.00000
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END