Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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