C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 55 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is Per_m in component Per_m (nanomolar).
C CONSTS(1) is B1 in component Per_m (dimensionless).
C CONSTS(2) is C1 in component Per_m (flux).
C CONSTS(3) is S1 in component Per_m (flux).
C CONSTS(4) is D1 in component Per_m (flux).
C CONSTS(5) is L1 in component Per_m (nanomolar).
C CONSTS(6) is R1 in component Per_m (nanomolar).
C CONSTS(7) is A1 in component Per_m (nanomolar).
C STATES(2) is PT_n in component PT_n (nanomolar).
C STATES(3) is CC_n in component CC_n (nanomolar).
C CONSTS(8) is D0 in component parameters (first_order_rate_constant).
C CONSTS(9) is a in component parameters (dimensionless).
C CONSTS(10) is r in component parameters (dimensionless).
C STATES(4) is Per_c in component Per_c (nanomolar).
C CONSTS(11) is S2 in component Per_c (first_order_rate_constant).
C CONSTS(12) is D2 in component Per_c (first_order_rate_constant).
C CONSTS(13) is L2 in component Per_c (nanomolar).
C CONSTS(14) is Dbt_c in component Per_c (nanomolar).
C CONSTS(15) is V1 in component parameters (second_order_rate_constant).
C CONSTS(16) is V2 in component parameters (first_order_rate_constant).
C STATES(5) is Tim_c in component Tim_c (nanomolar).
C STATES(6) is PT_c in component PT_c (nanomolar).
C STATES(7) is Tim_m in component Tim_m (nanomolar).
C CONSTS(17) is B2 in component Tim_m (dimensionless).
C CONSTS(18) is C2 in component Tim_m (flux).
C CONSTS(19) is S3 in component Tim_m (flux).
C CONSTS(20) is D3 in component Tim_m (flux).
C CONSTS(21) is L3 in component Tim_m (nanomolar).
C CONSTS(22) is R2 in component Tim_m (nanomolar).
C CONSTS(23) is A2 in component Tim_m (nanomolar).
C CONSTS(24) is S4 in component Tim_c (first_order_rate_constant).
C CONSTS(25) is D4 in component Tim_c (flux).
C CONSTS(26) is L4 in component Tim_c (nanomolar).
C CONSTS(27) is D5 in component PT_c (flux).
C CONSTS(28) is L5 in component PT_c (nanomolar).
C CONSTS(29) is K1 in component parameters (nanomolar).
C CONSTS(30) is K2 in component parameters (nanomolar).
C CONSTS(31) is T1 in component parameters (flux).
C CONSTS(32) is T2 in component parameters (flux).
C CONSTS(33) is D6 in component PT_n (flux).
C CONSTS(34) is L6 in component PT_n (nanomolar).
C STATES(8) is Clk_m in component Clk_m (nanomolar).
C CONSTS(35) is B3 in component Clk_m (dimensionless).
C CONSTS(36) is C3 in component Clk_m (flux).
C CONSTS(37) is S5 in component Clk_m (flux).
C CONSTS(38) is D7 in component Clk_m (flux).
C CONSTS(39) is L7 in component Clk_m (nanomolar).
C CONSTS(40) is R3 in component Clk_m (nanomolar).
C CONSTS(41) is A3 in component Clk_m (nanomolar).
C STATES(9) is Clk_c in component Clk_c (nanomolar).
C CONSTS(42) is S6 in component Clk_c (first_order_rate_constant).
C CONSTS(43) is D8 in component Clk_c (flux).
C CONSTS(44) is L8 in component Clk_c (nanomolar).
C CONSTS(45) is V3 in component parameters (second_order_rate_constant).
C CONSTS(46) is V4 in component parameters (first_order_rate_constant).
C CONSTS(47) is Cyc_c in component Cyc_c (nanomolar).
C STATES(10) is CC_c in component CC_c (nanomolar).
C CONSTS(48) is D9 in component CC_c (flux).
C CONSTS(49) is L9 in component CC_c (nanomolar).
C CONSTS(50) is K3 in component parameters (nanomolar).
C CONSTS(51) is K4 in component parameters (nanomolar).
C CONSTS(52) is T3 in component parameters (flux).
C CONSTS(53) is T4 in component parameters (flux).
C CONSTS(54) is D10 in component CC_n (flux).
C CONSTS(55) is L10 in component CC_n (nanomolar).
C RATES(1) is d/dt Per_m in component Per_m (nanomolar).
C RATES(4) is d/dt Per_c in component Per_c (nanomolar).
C RATES(7) is d/dt Tim_m in component Tim_m (nanomolar).
C RATES(5) is d/dt Tim_c in component Tim_c (nanomolar).
C RATES(6) is d/dt PT_c in component PT_c (nanomolar).
C RATES(2) is d/dt PT_n in component PT_n (nanomolar).
C RATES(8) is d/dt Clk_m in component Clk_m (nanomolar).
C RATES(9) is d/dt Clk_c in component Clk_c (nanomolar).
C RATES(10) is d/dt CC_c in component CC_c (nanomolar).
C RATES(3) is d/dt CC_n in component CC_n (nanomolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.5
      CONSTS(1) = 0.0
      CONSTS(2) = 0.0
      CONSTS(3) = 1.45
      CONSTS(4) = 0.94
      CONSTS(5) = 0.3
      CONSTS(6) = 1.02
      CONSTS(7) = 0.45
      STATES(2) = 1.0
      STATES(3) = 0.4
      CONSTS(8) = 0.012
      CONSTS(9) = 1.0
      CONSTS(10) = 4.0
      STATES(4) = 0.6
      CONSTS(11) = 0.48
      CONSTS(12) = 0.44
      CONSTS(13) = 0.2
      CONSTS(14) = 1.0
      CONSTS(15) = 1.45
      CONSTS(16) = 1.45
      STATES(5) = 0.8
      STATES(6) = 0.9
      STATES(7) = 0.7
      CONSTS(17) = 0.0
      CONSTS(18) = 0.0
      CONSTS(19) = 1.45
      CONSTS(20) = 0.94
      CONSTS(21) = 0.3
      CONSTS(22) = 1.02
      CONSTS(23) = 0.45
      CONSTS(24) = 0.48
      CONSTS(25) = 0.44
      CONSTS(26) = 0.2
      CONSTS(27) = 0.44
      CONSTS(28) = 0.2
      CONSTS(29) = 2.0
      CONSTS(30) = 2.0
      CONSTS(31) = 1.73
      CONSTS(32) = 0.72
      CONSTS(33) = 0.29
      CONSTS(34) = 0.2
      STATES(8) = 0.1
      CONSTS(35) = 0.6
      CONSTS(36) = 0.0
      CONSTS(37) = 1.63
      CONSTS(38) = 0.54
      CONSTS(39) = 0.13
      CONSTS(40) = 0.89
      CONSTS(41) = 0.8
      STATES(9) = 0.2
      CONSTS(42) = 0.47
      CONSTS(43) = 0.6
      CONSTS(44) = 0.2
      CONSTS(45) = 1.63
      CONSTS(46) = 1.63
      CONSTS(47) = 1.0
      STATES(10) = 0.3
      CONSTS(48) = 0.6
      CONSTS(49) = 0.2
      CONSTS(50) = 2.0
      CONSTS(51) = 2.0
      CONSTS(52) = 1.63
      CONSTS(53) = 0.52
      CONSTS(54) = 0.3
      CONSTS(55) = 0.2
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (CONSTS(2)+ CONSTS(3)*((STATES(3)/CONSTS(7) ** CONSTS(9)+CONSTS(1))/(1.00000+STATES(2)/CONSTS(6) ** CONSTS(10)+STATES(3)/CONSTS(7) ** CONSTS(9)+CONSTS(1)))) - ( CONSTS(4)*(STATES(1)/(CONSTS(5)+STATES(1)))+ CONSTS(8)*STATES(1))
      RATES(4) = ( CONSTS(11)*STATES(1)+ CONSTS(16)*STATES(6)) - ( CONSTS(15)*STATES(4)*STATES(5)+ CONSTS(12)*CONSTS(14)*(STATES(4)/(CONSTS(13)+STATES(4)))+ CONSTS(8)*STATES(4))
      RATES(7) = (CONSTS(18)+ CONSTS(19)*((STATES(3)/CONSTS(23) ** CONSTS(9)+CONSTS(17))/(1.00000+STATES(2)/CONSTS(22) ** CONSTS(10)+STATES(3)/CONSTS(23) ** CONSTS(9)+CONSTS(17)))) - ( CONSTS(20)*(STATES(7)/(CONSTS(21)+STATES(7)))+ CONSTS(8)*STATES(7))
      RATES(5) = ( CONSTS(24)*STATES(7)+ CONSTS(16)*STATES(6)) - ( CONSTS(15)*STATES(4)*STATES(5)+ CONSTS(25)*(STATES(5)/(CONSTS(26)+STATES(5)))+ CONSTS(8)*STATES(5))
      RATES(6) = ( CONSTS(15)*STATES(4)*STATES(5)+ CONSTS(32)*(STATES(2)/(CONSTS(30)+STATES(2)))) - ( CONSTS(16)*STATES(6)+ CONSTS(31)*(STATES(6)/(CONSTS(29)+STATES(6)))+ CONSTS(27)*(STATES(6)/(CONSTS(28)+STATES(6)))+ CONSTS(8)*STATES(6))
      RATES(2) =  CONSTS(31)*(STATES(6)/(CONSTS(29)+STATES(6))) - ( CONSTS(32)*(STATES(2)/(CONSTS(30)+STATES(2)))+ CONSTS(33)*(STATES(2)/(CONSTS(34)+STATES(2)))+ CONSTS(8)*STATES(2))
      RATES(8) = (CONSTS(36)+ CONSTS(37)*((STATES(2)/CONSTS(41) ** CONSTS(9)+CONSTS(35))/(1.00000+STATES(3)/CONSTS(40) ** CONSTS(10)+STATES(2)/CONSTS(41) ** CONSTS(9)+CONSTS(35)))) - ( CONSTS(38)*(STATES(8)/(CONSTS(39)+STATES(8)))+ CONSTS(8)*STATES(8))
      RATES(9) = ( CONSTS(42)*STATES(8)+ CONSTS(46)*STATES(10)) - ( CONSTS(45)*STATES(9)*CONSTS(47)+ CONSTS(43)*(STATES(9)/(CONSTS(44)+STATES(9)))+ CONSTS(8)*STATES(9))
      RATES(10) = ( CONSTS(45)*STATES(9)*CONSTS(47)+ CONSTS(53)*(STATES(3)/(CONSTS(51)+STATES(3)))) - ( CONSTS(46)*STATES(10)+ CONSTS(52)*(STATES(10)/(CONSTS(50)+STATES(10)))+ CONSTS(48)*(STATES(10)/(CONSTS(49)+STATES(10)))+ CONSTS(8)*STATES(10))
      RATES(3) =  CONSTS(52)*(STATES(10)/(CONSTS(50)+STATES(10))) - ( CONSTS(53)*(STATES(3)/(CONSTS(51)+STATES(3)))+ CONSTS(54)*(STATES(3)/(CONSTS(55)+STATES(3)))+ CONSTS(8)*STATES(3))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END