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 19 entries in the algebraic variable array.
C There are a total of 9 entries in each of the rate and state variable arrays.
C There are a total of 30 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is NO in component NO (micromolar).
C STATES(2) is O_2m in component O_2m (micromolar).
C CONSTS(1) is O_2 in component NO (micromolar).
C STATES(3) is NO_2 in component NO (micromolar).
C STATES(4) is N2O3 in component N2O3 (micromolar).
C STATES(5) is GSNO in component GSNO (micromolar).
C STATES(6) is CcOX in component NO (micromolar).
C STATES(7) is FeLn in component FeLn (micromolar).
C CONSTS(29) is r_1NO in component NO (flux).
C ALGBRC(1) is r_4NO in component NO (flux).
C ALGBRC(2) is r_12aNO in component NO (flux).
C ALGBRC(3) is r_12bNOp in component NO (flux).
C ALGBRC(4) is r_12bNOm in component NO (flux).
C ALGBRC(5) is r_14NO in component NO (flux).
C ALGBRC(6) is r_15NO in component NO (flux).
C ALGBRC(7) is r_16NO in component NO (flux).
C CONSTS(2) is k_1NO in component model_constant (flux).
C CONSTS(3) is k_4NO in component model_constant (second_order_rate_constant).
C CONSTS(4) is k_12aNO in component model_constant (rate2).
C CONSTS(5) is k_12bNOp in component model_constant (second_order_rate_constant).
C CONSTS(6) is k_12bNOm in component model_constant (first_order_rate_constant).
C CONSTS(7) is k_14NO in component model_constant (first_order_rate_constant).
C CONSTS(8) is k_15NO in component model_constant (second_order_rate_constant).
C CONSTS(9) is k_16NO in component model_constant (second_order_rate_constant).
C CONSTS(10) is SOD in component O_2m (micromolar).
C CONSTS(30) is r_2NO in component O_2m (flux).
C ALGBRC(8) is r_5NO in component O_2m (flux).
C ALGBRC(9) is r_10NO in component O_2m (flux).
C CONSTS(11) is k_2NO in component model_constant (flux).
C CONSTS(12) is k_5NO in component model_constant (second_order_rate_constant).
C CONSTS(13) is k_10NO in component model_constant (rate2).
C STATES(8) is ONOO_m in component ONOO_m (micromolar).
C STATES(9) is GSH in component GSH (micromolar).
C CONSTS(14) is GPX in component ONOO_m (micromolar).
C CONSTS(15) is CO_2 in component ONOO_m (micromolar).
C CONSTS(16) is Cyt_c in component ONOO_m (micromolar).
C ALGBRC(10) is r_6NO in component ONOO_m (flux).
C ALGBRC(11) is r_7NO in component ONOO_m (flux).
C ALGBRC(13) is r_8NO in component ONOO_m (flux).
C ALGBRC(15) is r_9NO in component ONOO_m (flux).
C CONSTS(17) is k_6NO in component model_constant (second_order_rate_constant).
C CONSTS(18) is k_7NO in component model_constant (second_order_rate_constant).
C CONSTS(19) is k_8NO in component model_constant (second_order_rate_constant).
C CONSTS(20) is k_9NO in component model_constant (second_order_rate_constant).
C ALGBRC(14) is GSSG in component GSH (micromolar).
C ALGBRC(12) is FeLnNO in component GSH (micromolar).
C CONSTS(21) is FeLn_0 in component model_constant (micromolar).
C CONSTS(22) is GSH_0 in component model_constant (micromolar).
C ALGBRC(16) is r_11NO in component GSH (flux).
C ALGBRC(17) is r_m in component GSH (flux).
C ALGBRC(19) is r_17NO in component GSH (flux).
C CONSTS(23) is k_11NO in component model_constant (second_order_rate_constant).
C CONSTS(24) is v_m in component model_constant (flux).
C CONSTS(25) is k_m in component model_constant (micromolar).
C CONSTS(26) is k_17NO in component model_constant (second_order_rate_constant).
C ALGBRC(18) is r_13NO in component N2O3 (flux).
C CONSTS(27) is k_13NO in component model_constant (first_order_rate_constant).
C CONSTS(28) is k_17bNO in component model_constant (second_order_rate_constant).
C RATES(1) is d/dt NO in component NO (micromolar).
C RATES(6) is d/dt CcOX in component NO (micromolar).
C RATES(3) is d/dt NO_2 in component NO (micromolar).
C RATES(2) is d/dt O_2m in component O_2m (micromolar).
C RATES(8) is d/dt ONOO_m in component ONOO_m (micromolar).
C RATES(9) is d/dt GSH in component GSH (micromolar).
C RATES(5) is d/dt GSNO in component GSNO (micromolar).
C RATES(4) is d/dt N2O3 in component N2O3 (micromolar).
C RATES(7) is d/dt FeLn in component FeLn (micromolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0
STATES(2) = 0
CONSTS(1) = 35
STATES(3) = 0
STATES(4) = 0
STATES(5) = 0
STATES(6) = 0.1
STATES(7) = 0.05
CONSTS(2) = 1
CONSTS(3) = 6700
CONSTS(4) = 0.000006
CONSTS(5) = 1100
CONSTS(6) = 81000
CONSTS(7) = 0.0002
CONSTS(8) = 100
CONSTS(9) = 1.21
CONSTS(10) = 10
CONSTS(11) = 0.1
CONSTS(12) = 2400
CONSTS(13) = 0.0006
STATES(8) = 0
STATES(9) = 10000
CONSTS(14) = 5.8
CONSTS(15) = 1000
CONSTS(16) = 400
CONSTS(17) = 0.00135
CONSTS(18) = 2
CONSTS(19) = 0.058
CONSTS(20) = 0.025
CONSTS(21) = 0.05
CONSTS(22) = 10000
CONSTS(23) = 66
CONSTS(24) = 320
CONSTS(25) = 50
CONSTS(26) = 66
CONSTS(27) = 1600
CONSTS(28) = 0.0002
CONSTS(29) = CONSTS(2)
CONSTS(30) = CONSTS(11)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = CONSTS(4)*STATES(1) ** 2.00000*CONSTS(1)
ALGBRC(3) = CONSTS(5)*STATES(3)*STATES(1)
ALGBRC(4) = CONSTS(6)*STATES(4)
RATES(3) = ( 2.00000*ALGBRC(2) - ALGBRC(3))+ALGBRC(4)
ALGBRC(6) = CONSTS(8)*STATES(6)*STATES(1)
RATES(6) = - ALGBRC(6)
ALGBRC(1) = CONSTS(3)*STATES(1)*STATES(2)
ALGBRC(5) = CONSTS(7)*STATES(5)
ALGBRC(7) = CONSTS(9)*STATES(7)*STATES(1)
RATES(1) = (((((CONSTS(29) - ALGBRC(1)) - 2.00000*ALGBRC(2)) - ALGBRC(3))+ALGBRC(4)+ALGBRC(5)) - ALGBRC(6)) - ALGBRC(7)
ALGBRC(8) = CONSTS(12)*CONSTS(10)*STATES(2)
ALGBRC(9) = CONSTS(13)*STATES(5) ** 2.00000*STATES(2)
RATES(2) = ((CONSTS(30) - ALGBRC(1)) - ALGBRC(8)) - ALGBRC(9)
ALGBRC(10) = CONSTS(17)*STATES(8)*STATES(9)
ALGBRC(11) = CONSTS(18)*STATES(8)*CONSTS(14)
ALGBRC(13) = CONSTS(19)*STATES(8)*CONSTS(15)
ALGBRC(15) = CONSTS(20)*STATES(8)*CONSTS(16)
RATES(8) = (((ALGBRC(1) - ALGBRC(10)) - ALGBRC(11)) - ALGBRC(13)) - ALGBRC(15)
ALGBRC(16) = CONSTS(23)*STATES(4)*STATES(9)
ALGBRC(18) = CONSTS(27)*STATES(4)
RATES(4) = ((- ALGBRC(16)+ALGBRC(3)) - ALGBRC(4)) - ALGBRC(18)
ALGBRC(14) = ((CONSTS(22) - STATES(9)) - STATES(5))/2.00000
ALGBRC(17) = ( CONSTS(24)*ALGBRC(14))/(CONSTS(25)+ALGBRC(14))
ALGBRC(12) = CONSTS(21) - STATES(7)
ALGBRC(19) = CONSTS(26)*ALGBRC(12)*STATES(9)
RATES(9) = ((- ALGBRC(10) - ALGBRC(16))+ 2.00000*ALGBRC(17)) - ALGBRC(19)
RATES(5) = (((ALGBRC(10) - 2.00000*ALGBRC(9))+ALGBRC(16)) - ALGBRC(5))+ALGBRC(19)
RATES(7) = - ALGBRC(7)+ALGBRC(19)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = CONSTS(4)*STATES(1) ** 2.00000*CONSTS(1)
ALGBRC(3) = CONSTS(5)*STATES(3)*STATES(1)
ALGBRC(4) = CONSTS(6)*STATES(4)
ALGBRC(6) = CONSTS(8)*STATES(6)*STATES(1)
ALGBRC(1) = CONSTS(3)*STATES(1)*STATES(2)
ALGBRC(5) = CONSTS(7)*STATES(5)
ALGBRC(7) = CONSTS(9)*STATES(7)*STATES(1)
ALGBRC(8) = CONSTS(12)*CONSTS(10)*STATES(2)
ALGBRC(9) = CONSTS(13)*STATES(5) ** 2.00000*STATES(2)
ALGBRC(10) = CONSTS(17)*STATES(8)*STATES(9)
ALGBRC(11) = CONSTS(18)*STATES(8)*CONSTS(14)
ALGBRC(13) = CONSTS(19)*STATES(8)*CONSTS(15)
ALGBRC(15) = CONSTS(20)*STATES(8)*CONSTS(16)
ALGBRC(16) = CONSTS(23)*STATES(4)*STATES(9)
ALGBRC(18) = CONSTS(27)*STATES(4)
ALGBRC(14) = ((CONSTS(22) - STATES(9)) - STATES(5))/2.00000
ALGBRC(17) = ( CONSTS(24)*ALGBRC(14))/(CONSTS(25)+ALGBRC(14))
ALGBRC(12) = CONSTS(21) - STATES(7)
ALGBRC(19) = CONSTS(26)*ALGBRC(12)*STATES(9)
RETURN
END