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 6 entries in each of the rate and state variable arrays.
C There are a total of 25 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is j_1 in component parameters (flux).
C CONSTS(2) is j_2 in component parameters (flux).
C CONSTS(3) is j_3 in component parameters (flux).
C CONSTS(4) is v_d1 in component parameters (flux).
C CONSTS(5) is v_d2 in component parameters (flux).
C CONSTS(6) is v_d3 in component parameters (flux).
C CONSTS(7) is k_d1 in component parameters (first_order_rate_constant).
C CONSTS(8) is k_d2 in component parameters (first_order_rate_constant).
C CONSTS(9) is k_d3 in component parameters (first_order_rate_constant).
C CONSTS(10) is k_c1 in component parameters (first_order_rate_constant).
C CONSTS(11) is k_c2 in component parameters (first_order_rate_constant).
C CONSTS(12) is k_c3 in component parameters (first_order_rate_constant).
C CONSTS(13) is k_m1 in component parameters (nanomolar).
C CONSTS(14) is k_m2 in component parameters (nanomolar).
C CONSTS(15) is k_m3 in component parameters (nanomolar).
C CONSTS(16) is v_12 in component parameters (flux).
C CONSTS(17) is v_11 in component parameters (flux).
C CONSTS(18) is v_10 in component parameters (flux).
C CONSTS(19) is k_120 in component parameters (nanomolar).
C CONSTS(20) is k_110 in component parameters (nanomolar).
C CONSTS(21) is k_100 in component parameters (nanomolar).
C CONSTS(22) is k_d4 in component parameters (first_order_rate_constant).
C CONSTS(23) is k_d5 in component parameters (first_order_rate_constant).
C CONSTS(24) is k_d6 in component parameters (first_order_rate_constant).
C CONSTS(25) is n in component parameters (dimensionless).
C STATES(1) is C_1 in component C_1 (nanomolar).
C STATES(2) is C_2 in component C_2 (nanomolar).
C STATES(3) is T_1 in component T_1 (nanomolar).
C STATES(4) is C_3 in component C_3 (nanomolar).
C STATES(5) is T_2 in component T_2 (nanomolar).
C STATES(6) is T_3 in component T_3 (nanomolar).
C RATES(1) is d/dt C_1 in component C_1 (nanomolar).
C RATES(2) is d/dt C_2 in component C_2 (nanomolar).
C RATES(4) is d/dt C_3 in component C_3 (nanomolar).
C RATES(3) is d/dt T_1 in component T_1 (nanomolar).
C RATES(5) is d/dt T_2 in component T_2 (nanomolar).
C RATES(6) is d/dt T_3 in component T_3 (nanomolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.9
CONSTS(2) = 0.5
CONSTS(3) = 0.6
CONSTS(4) = 6
CONSTS(5) = 1.052
CONSTS(6) = 3
CONSTS(7) = 0.8
CONSTS(8) = 0.9
CONSTS(9) = 0.8
CONSTS(10) = 0.2
CONSTS(11) = 0.22
CONSTS(12) = 0.6
CONSTS(13) = 5
CONSTS(14) = 5
CONSTS(15) = 5
CONSTS(16) = 15
CONSTS(17) = 15
CONSTS(18) = 15
CONSTS(19) = 10
CONSTS(20) = 10
CONSTS(21) = 10
CONSTS(22) = 0.16
CONSTS(23) = 0.16
CONSTS(24) = 0.16
CONSTS(25) = 2
STATES(1) = 0
STATES(2) = 0
STATES(3) = 6
STATES(4) = 0
STATES(5) = 5
STATES(6) = 1
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = ( CONSTS(16)*STATES(3) ** CONSTS(25))/(CONSTS(19) ** CONSTS(25)+STATES(3) ** CONSTS(25)+STATES(2) ** CONSTS(25)) - CONSTS(22)*STATES(1)
RATES(2) = ( CONSTS(17)*STATES(5) ** CONSTS(25))/(CONSTS(20) ** CONSTS(25)+STATES(5) ** CONSTS(25)+STATES(4) ** CONSTS(25)) - CONSTS(23)*STATES(2)
RATES(4) = ( CONSTS(18)*STATES(6) ** CONSTS(25))/(CONSTS(21) ** CONSTS(25)+STATES(6) ** CONSTS(25)+STATES(1) ** CONSTS(25)) - CONSTS(24)*STATES(4)
RATES(3) = (CONSTS(1)+( CONSTS(4)*STATES(6) ** CONSTS(25))/(CONSTS(13) ** CONSTS(25)+STATES(6) ** CONSTS(25))+ CONSTS(10)*STATES(1)) - CONSTS(7)*STATES(3)
RATES(5) = (CONSTS(2)+( CONSTS(5)*STATES(3) ** CONSTS(25))/(CONSTS(14) ** CONSTS(25)+STATES(3) ** CONSTS(25))+ CONSTS(11)*STATES(2)) - CONSTS(8)*STATES(5)
RATES(6) = (CONSTS(3)+( CONSTS(6)*STATES(5) ** CONSTS(25))/(CONSTS(15) ** CONSTS(25)+STATES(5) ** CONSTS(25))+ CONSTS(12)*STATES(4)) - CONSTS(9)*STATES(6)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
