C
C There are a total of 4 entries in the algebraic variable array.
C There are a total of 17 entries in each of the rate and state variable arrays.
C There are a total of 43 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is MPF_n in component MPF_n (dimensionless).
C CONSTS(1) is k_in in component parameters (first_order_rate_constant).
C CONSTS(2) is k_out in component parameters (first_order_rate_constant).
C CONSTS(3) is k_dn_ in component parameters (first_order_rate_constant).
C CONSTS(4) is k_dn__ in component parameters (first_order_rate_constant).
C CONSTS(5) is k_wee_ in component parameters (first_order_rate_constant).
C CONSTS(6) is k_wee__ in component parameters (first_order_rate_constant).
C CONSTS(7) is k_stg_ in component parameters (first_order_rate_constant).
C CONSTS(8) is k_stg__ in component parameters (first_order_rate_constant).
C STATES(2) is MPF_c in component MPF_c (dimensionless).
C STATES(3) is FZY in component FZY (dimensionless).
C STATES(4) is Wee1_n in component Wee1_n (dimensionless).
C STATES(5) is StgP_n in component StgP_n (dimensionless).
C STATES(6) is preMPF_n in component preMPF_n (dimensionless).
C STATES(7) is preMPF_c in component preMPF_c (dimensionless).
C CONSTS(9) is k_sc in component parameters (first_order_rate_constant).
C CONSTS(10) is k_dc_ in component parameters (first_order_rate_constant).
C CONSTS(11) is epsilon in component parameters (dimensionless).
C STATES(8) is N in component N (dimensionless).
C STATES(9) is StgP_c in component StgP_c (dimensionless).
C STATES(10) is Wee1_c in component Wee1_c (dimensionless).
C STATES(11) is IE in component IE (dimensionless).
C CONSTS(12) is j_aie in component parameters (dimensionless).
C CONSTS(13) is j_iie in component parameters (dimensionless).
C CONSTS(14) is k_aie in component parameters (first_order_rate_constant).
C CONSTS(15) is k_iie in component parameters (first_order_rate_constant).
C CONSTS(16) is j_afz in component parameters (dimensionless).
C CONSTS(17) is j_ifz in component parameters (dimensionless).
C CONSTS(18) is k_afz in component parameters (first_order_rate_constant).
C CONSTS(19) is k_ifz in component parameters (first_order_rate_constant).
C STATES(12) is Stg_m in component Stg_m (dimensionless).
C CONSTS(20) is k_dm_ in component parameters (first_order_rate_constant).
C CONSTS(21) is k_dm__ in component parameters (first_order_rate_constant).
C CONSTS(22) is j_m in component parameters (dimensionless).
C STATES(13) is Xp in component Xp (dimensionless).
C STATES(14) is Xm in component Xm (dimensionless).
C CONSTS(23) is k_sxm in component parameters (first_order_rate_constant).
C CONSTS(24) is k_sxp in component parameters (first_order_rate_constant).
C CONSTS(25) is k_ins in component parameters (first_order_rate_constant).
C CONSTS(26) is k_outs in component parameters (first_order_rate_constant).
C CONSTS(27) is k_astg_ in component parameters (first_order_rate_constant).
C CONSTS(28) is k_astg__ in component parameters (first_order_rate_constant).
C CONSTS(29) is k_istg in component parameters (first_order_rate_constant).
C CONSTS(30) is k_dstg in component parameters (first_order_rate_constant).
C CONSTS(31) is j_astg in component parameters (dimensionless).
C CONSTS(32) is j_istg in component parameters (dimensionless).
C STATES(15) is Stg_n in component Stg_n (dimensionless).
C STATES(16) is Stg_c in component Stg_c (dimensionless).
C CONSTS(33) is k_sstg in component parameters (first_order_rate_constant).
C CONSTS(34) is k_inw in component parameters (first_order_rate_constant).
C CONSTS(35) is k_outw in component parameters (first_order_rate_constant).
C CONSTS(36) is k_awee in component parameters (first_order_rate_constant).
C CONSTS(37) is k_iwee_ in component parameters (first_order_rate_constant).
C CONSTS(38) is k_iwee__ in component parameters (first_order_rate_constant).
C CONSTS(39) is j_awee in component parameters (dimensionless).
C CONSTS(40) is j_iwee in component parameters (dimensionless).
C STATES(17) is Wee1P_n in component Wee1P_n (dimensionless).
C ALGBRC(1) is Wee1P_c in component Wee1P_c (dimensionless).
C CONSTS(41) is Wee1_T in component Wee1P_c (dimensionless).
C ALGBRC(2) is CycB_T in component CycB_T (dimensionless).
C ALGBRC(3) is Stg_T in component Stg_T (dimensionless).
C ALGBRC(4) is StgP_T in component StgP_T (dimensionless).
C CONSTS(42) is k_ez in component parameters (first_order_rate_constant).
C RATES(1) is d/dt MPF_n in component MPF_n (dimensionless).
C RATES(6) is d/dt preMPF_n in component preMPF_n (dimensionless).
C RATES(2) is d/dt MPF_c in component MPF_c (dimensionless).
C RATES(7) is d/dt preMPF_c in component preMPF_c (dimensionless).
C RATES(11) is d/dt IE in component IE (dimensionless).
C RATES(3) is d/dt FZY in component FZY (dimensionless).
C RATES(12) is d/dt Stg_m in component Stg_m (dimensionless).
C RATES(14) is d/dt Xm in component Xm (dimensionless).
C RATES(13) is d/dt Xp in component Xp (dimensionless).
C RATES(5) is d/dt StgP_n in component StgP_n (dimensionless).
C RATES(15) is d/dt Stg_n in component Stg_n (dimensionless).
C RATES(9) is d/dt StgP_c in component StgP_c (dimensionless).
C RATES(16) is d/dt Stg_c in component Stg_c (dimensionless).
C RATES(4) is d/dt Wee1_n in component Wee1_n (dimensionless).
C RATES(10) is d/dt Wee1_c in component Wee1_c (dimensionless).
C RATES(17) is d/dt Wee1P_n in component Wee1P_n (dimensionless).
C RATES(8) is d/dt N in component N (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0
      CONSTS(1) = 0.15
      CONSTS(2) = 0
      CONSTS(3) = 0.01
      CONSTS(4) = 1.5
      CONSTS(5) = 0.005
      CONSTS(6) = 1
      CONSTS(7) = 0.2
      CONSTS(8) = 2
      STATES(2) = 1
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      CONSTS(9) = 0.01
      CONSTS(10) = 0.01
      CONSTS(11) = 0.00007
      STATES(8) = 1
      STATES(9) = 0
      STATES(10) = 0
      STATES(11) = 0
      CONSTS(12) = 0.01
      CONSTS(13) = 0.01
      CONSTS(14) = 1
      CONSTS(15) = 0.4
      CONSTS(16) = 0.01
      CONSTS(17) = 0.01
      CONSTS(18) = 1
      CONSTS(19) = 0.2
      STATES(12) = 1
      CONSTS(20) = 0.002
      CONSTS(21) = 0.2
      CONSTS(22) = 0.05
      STATES(13) = 0
      STATES(14) = 0
      CONSTS(23) = 0.0005
      CONSTS(24) = 0.001
      CONSTS(25) = 0.08
      CONSTS(26) = 0.02
      CONSTS(27) = 0
      CONSTS(28) = 1
      CONSTS(29) = 0.3
      CONSTS(30) = 0.015
      CONSTS(31) = 0.05
      CONSTS(32) = 0.05
      STATES(15) = 0
      STATES(16) = 1
      CONSTS(33) = 0.02
      CONSTS(34) = 0.04
      CONSTS(35) = 0.01
      CONSTS(36) = 0.3
      CONSTS(37) = 0.01
      CONSTS(38) = 1
      CONSTS(39) = 0.05
      CONSTS(40) = 0.05
      STATES(17) = 0
      CONSTS(41) = 0.8
      CONSTS(42) = 0.5
      CONSTS(42) = 0.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(8) = CONSTS(42)
      RATES(1) = ((( CONSTS(1)*STATES(2) -  CONSTS(2)*STATES(1)) -  (CONSTS(3)+ CONSTS(4)*STATES(3))*STATES(1)) -  (CONSTS(5)+ CONSTS(6)*STATES(4))*STATES(1))+ (CONSTS(7)+ CONSTS(8)*STATES(5))*STATES(6)
      RATES(6) = ((( CONSTS(1)*STATES(7) -  CONSTS(2)*STATES(6)) -  (CONSTS(3)+ CONSTS(4)*STATES(3))*STATES(6))+ (CONSTS(5)+ CONSTS(6)*STATES(4))*STATES(1)) -  (CONSTS(5)+ CONSTS(8)*STATES(5))*STATES(6)
      RATES(2) = (((CONSTS(9) - ( CONSTS(11)*STATES(8)*( CONSTS(1)*STATES(2) -  CONSTS(2)*STATES(1)))/(1.00000 -  STATES(8)*CONSTS(11))) -  CONSTS(10)*STATES(2))+ (CONSTS(7)+ CONSTS(8)*STATES(9))*STATES(7)) -  (CONSTS(5)+ CONSTS(6)*STATES(10))*STATES(2)
      RATES(7) = ((( - CONSTS(11)*STATES(8)*( CONSTS(1)*STATES(7) -  CONSTS(2)*STATES(6)))/(1.00000 -  STATES(8)*CONSTS(11)) -  CONSTS(10)*STATES(7)) -  (CONSTS(7)+ CONSTS(8)*STATES(9))*STATES(7))+ (CONSTS(5)+ CONSTS(6)*STATES(10))*STATES(2)
      RATES(11) = ( CONSTS(14)*(1.00000 - STATES(11))*STATES(1))/((CONSTS(12)+1.00000) - STATES(11)) - ( CONSTS(15)*STATES(11))/(CONSTS(13)+STATES(11))
      RATES(3) = ( CONSTS(18)*STATES(11)*(1.00000 - STATES(3)))/((CONSTS(16)+1.00000) - STATES(3)) - ( CONSTS(19)*STATES(3))/(CONSTS(17)+STATES(3))
      RATES(12) =  - (CONSTS(20)/(CONSTS(22)+STATES(12))+ CONSTS(21)*STATES(13))*STATES(12)
      RATES(14) =  CONSTS(23)*STATES(8)
      RATES(13) =  CONSTS(24)*STATES(14)
      RATES(5) = ((( CONSTS(25)*STATES(9) -  CONSTS(26)*STATES(5))+( (CONSTS(27)+ CONSTS(28)*STATES(1))*STATES(15))/(CONSTS(31)+STATES(15))) - ( CONSTS(29)*STATES(5))/(CONSTS(32)+STATES(5))) -  CONSTS(30)*STATES(5)
      RATES(15) = ((( CONSTS(25)*STATES(16) -  CONSTS(26)*STATES(15)) - ( (CONSTS(27)+ CONSTS(28)*STATES(1))*STATES(15))/(CONSTS(31)+STATES(15)))+( CONSTS(29)*STATES(5))/(CONSTS(32)+STATES(5))) -  CONSTS(30)*STATES(15)
      RATES(9) = (( - CONSTS(30)*STATES(9) - ( CONSTS(11)*STATES(8)*( CONSTS(25)*STATES(9) -  CONSTS(26)*STATES(5)))/(1.00000 -  STATES(8)*CONSTS(11)))+( (CONSTS(27)+ CONSTS(28)*STATES(2))*STATES(16))/(CONSTS(31)+STATES(16))) - ( CONSTS(29)*STATES(9))/(CONSTS(32)+STATES(9))
      RATES(16) = ((( CONSTS(33)*STATES(12) -  CONSTS(30)*STATES(16)) - ( CONSTS(11)*STATES(8)*( CONSTS(25)*STATES(16) -  CONSTS(26)*STATES(15)))/(1.00000 -  STATES(8)*CONSTS(11))) - ( (CONSTS(27)+ CONSTS(28)*STATES(2))*STATES(16))/(CONSTS(31)+STATES(16)))+( CONSTS(29)*STATES(9))/(CONSTS(32)+STATES(9))
      RATES(4) = (( CONSTS(34)*STATES(10) -  CONSTS(35)*STATES(4))+( CONSTS(36)*STATES(17))/(CONSTS(39)+STATES(17))) - ( (CONSTS(37)+ CONSTS(38)*STATES(1))*STATES(4))/(CONSTS(40)+STATES(4))
      ALGBRC(1) = (CONSTS(41) -  STATES(8)*CONSTS(11)*(STATES(4)+STATES(17)))/(1.00000 -  STATES(8)*CONSTS(11)) - STATES(10)
      RATES(10) = (( - ( CONSTS(34)*STATES(10) -  CONSTS(35)*STATES(4))*STATES(8)*CONSTS(11))/(1.00000 -  STATES(8)*CONSTS(11))+( CONSTS(36)*ALGBRC(1))/(CONSTS(39)+ALGBRC(1))) - ( (CONSTS(37)+ CONSTS(38)*STATES(2))*STATES(10))/(CONSTS(40)+STATES(10))
      RATES(17) = (( CONSTS(34)*ALGBRC(1) -  CONSTS(35)*STATES(17)) - ( CONSTS(36)*STATES(17))/(CONSTS(39)+STATES(17)))+( (CONSTS(37)+ CONSTS(38)*STATES(1))*STATES(4))/(CONSTS(40)+STATES(4))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = (CONSTS(41) -  STATES(8)*CONSTS(11)*(STATES(4)+STATES(17)))/(1.00000 -  STATES(8)*CONSTS(11)) - STATES(10)
      ALGBRC(2) =  (1.00000 -  STATES(8)*CONSTS(11))*(STATES(2)+STATES(7))+ STATES(8)*CONSTS(11)*(STATES(1)+STATES(6))
      ALGBRC(3) =  (1.00000 -  STATES(8)*CONSTS(11))*(STATES(16)+STATES(9))+ STATES(8)*CONSTS(11)*(STATES(15)+STATES(5))
      ALGBRC(4) =  (1.00000 -  STATES(8)*CONSTS(11))*STATES(9)+ STATES(8)*CONSTS(11)*STATES(5)
      RETURN
      END