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 4 entries in the algebraic variable array.
C There are a total of 20 entries in each of the rate and state variable arrays.
C There are a total of 26 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is S_0 in component S_0 (dimensionless).
C CONSTS(16) is R_0 in component R (first_order_rate_constant).
C CONSTS(17) is u_s in component kinetic_parameters (first_order_rate_constant).
C ALGBRC(1) is phi in component phi (first_order_rate_constant).
C STATES(2) is S_1 in component S_1 (dimensionless).
C CONSTS(18) is R_1 in component R (first_order_rate_constant).
C CONSTS(1) is epsilon_s in component kinetic_parameters (first_order_rate_constant).
C CONSTS(2) is alpha in component kinetic_parameters (dimensionless).
C CONSTS(3) is u in component kinetic_parameters (first_order_rate_constant).
C STATES(3) is S_2 in component S_2 (dimensionless).
C CONSTS(19) is R_2 in component R (first_order_rate_constant).
C STATES(4) is S_3 in component S_3 (dimensionless).
C CONSTS(20) is R_3 in component R (first_order_rate_constant).
C STATES(5) is S_4 in component S_4 (dimensionless).
C CONSTS(21) is R_4 in component R (first_order_rate_constant).
C STATES(6) is S_5 in component S_5 (dimensionless).
C CONSTS(22) is R_5 in component R (first_order_rate_constant).
C STATES(7) is S_6 in component S_6 (dimensionless).
C CONSTS(23) is R_6 in component R (first_order_rate_constant).
C STATES(8) is S_7 in component S_7 (dimensionless).
C CONSTS(24) is R_7 in component R (first_order_rate_constant).
C STATES(9) is S_8 in component S_8 (dimensionless).
C CONSTS(25) is R_8 in component R (first_order_rate_constant).
C STATES(10) is M_0 in component M_0 (dimensionless).
C CONSTS(4) is epsilon_m in component kinetic_parameters (first_order_rate_constant).
C CONSTS(26) is u_m in component kinetic_parameters (first_order_rate_constant).
C STATES(11) is M_1 in component M_1 (dimensionless).
C STATES(12) is M_1_2 in component M_1 (dimensionless).
C STATES(13) is M_2 in component M_2 (dimensionless).
C STATES(14) is M_3 in component M_3 (dimensionless).
C STATES(15) is M_4 in component M_4 (dimensionless).
C STATES(16) is M_5 in component M_5 (dimensionless).
C STATES(17) is M_6 in component M_6 (dimensionless).
C STATES(18) is M_7 in component M_7 (dimensionless).
C STATES(19) is M_8 in component M_8 (dimensionless).
C CONSTS(5) is r_0 in component R (first_order_rate_constant).
C CONSTS(6) is r_1 in component R (first_order_rate_constant).
C CONSTS(7) is r_2 in component R (first_order_rate_constant).
C CONSTS(8) is r_3 in component R (first_order_rate_constant).
C CONSTS(9) is r_4 in component R (first_order_rate_constant).
C CONSTS(10) is r_5 in component R (first_order_rate_constant).
C CONSTS(11) is r_6 in component R (first_order_rate_constant).
C CONSTS(12) is r_7 in component R (first_order_rate_constant).
C CONSTS(13) is r_8 in component R (first_order_rate_constant).
C CONSTS(14) is a in component kinetic_parameters (dimensionless).
C STATES(20) is w in component w (dimensionless).
C ALGBRC(2) is total_cells in component total_cells (dimensionless).
C ALGBRC(3) is stable_total in component total_cells (dimensionless).
C ALGBRC(4) is mutant_total in component total_cells (dimensionless).
C CONSTS(15) is beta in component kinetic_parameters (dimensionless).
C RATES(1) is d/dt S_0 in component S_0 (dimensionless).
C RATES(2) is d/dt S_1 in component S_1 (dimensionless).
C RATES(3) is d/dt S_2 in component S_2 (dimensionless).
C RATES(4) is d/dt S_3 in component S_3 (dimensionless).
C RATES(5) is d/dt S_4 in component S_4 (dimensionless).
C RATES(6) is d/dt S_5 in component S_5 (dimensionless).
C RATES(7) is d/dt S_6 in component S_6 (dimensionless).
C RATES(8) is d/dt S_7 in component S_7 (dimensionless).
C RATES(9) is d/dt S_8 in component S_8 (dimensionless).
C RATES(10) is d/dt M_0 in component M_0 (dimensionless).
C RATES(11) is d/dt M_1 in component M_1 (dimensionless).
C RATES(12) is d/dt M_1_2 in component M_1 (dimensionless).
C RATES(13) is d/dt M_2 in component M_2 (dimensionless).
C RATES(14) is d/dt M_3 in component M_3 (dimensionless).
C RATES(15) is d/dt M_4 in component M_4 (dimensionless).
C RATES(16) is d/dt M_5 in component M_5 (dimensionless).
C RATES(17) is d/dt M_6 in component M_6 (dimensionless).
C RATES(18) is d/dt M_7 in component M_7 (dimensionless).
C RATES(19) is d/dt M_8 in component M_8 (dimensionless).
C RATES(20) is d/dt w in component w (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.5
STATES(2) = 0
CONSTS(1) = 0.99
CONSTS(2) = 0.6
CONSTS(3) = 0.07
STATES(3) = 0
STATES(4) = 0
STATES(5) = 0
STATES(6) = 0
STATES(7) = 0
STATES(8) = 0
STATES(9) = 0
STATES(10) = 0.5
CONSTS(4) = 0.1
STATES(11) = 0
STATES(12) = 0
STATES(13) = 0
STATES(14) = 0
STATES(15) = 0
STATES(16) = 0
STATES(17) = 0
STATES(18) = 0
STATES(19) = 0
CONSTS(5) = 0.5
CONSTS(6) = 0.6
CONSTS(7) = 0.7
CONSTS(8) = 0.8
CONSTS(9) = 0.9
CONSTS(10) = 1
CONSTS(11) = 1.1
CONSTS(12) = 1.2
CONSTS(13) = 1.3
CONSTS(14) = 0.5
STATES(20) = 0
CONSTS(15) = 0.2
CONSTS(16) = CONSTS(5)
CONSTS(17) = CONSTS(3)*(1.00000 - ( CONSTS(15)*CONSTS(1))/1.00000)
CONSTS(18) = CONSTS(6)*(1.00000 - CONSTS(14))
CONSTS(19) = CONSTS(7)*(1.00000 - CONSTS(14))
CONSTS(20) = CONSTS(8)*(1.00000 - CONSTS(14))
CONSTS(21) = CONSTS(9)*(1.00000 - CONSTS(14))
CONSTS(22) = CONSTS(10)*(1.00000 - CONSTS(14))
CONSTS(23) = CONSTS(11)*(1.00000 - CONSTS(14))
CONSTS(24) = CONSTS(12)*(1.00000 - CONSTS(14))
CONSTS(25) = CONSTS(13)*(1.00000 - CONSTS(14))
CONSTS(26) = CONSTS(3)*(1.00000 - ( CONSTS(15)*CONSTS(4))/1.00000)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(11) = STATES(10)*1.00000
ALGBRC(1) = ((1.00000 - CONSTS(17))/1.00000)*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5))+ ((1.00000 - CONSTS(26))/1.00000)*( CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15))
RATES(1) = (( CONSTS(16)*STATES(1))/1.00000)*(1.00000 - CONSTS(17)) - ALGBRC(1)*STATES(1)
RATES(2) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(16)*STATES(1))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(18)*STATES(2))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(2)
RATES(3) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(18)*STATES(2))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(19)*STATES(3))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(3)
RATES(4) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(19)*STATES(3))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(20)*STATES(4))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(4)
RATES(5) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(20)*STATES(4))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(21)*STATES(5))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(5)
RATES(6) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(21)*STATES(5))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(22)*STATES(6))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(6)
RATES(7) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(22)*STATES(6))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(23)*STATES(7))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(7)
RATES(8) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(23)*STATES(7))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(24)*STATES(8))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(8)
RATES(9) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(24)*STATES(8))/1.00000)*(1.00000 - CONSTS(1))+ (( CONSTS(25)*STATES(9))/1.00000)*(1.00000 - CONSTS(17))) - ALGBRC(1)*STATES(9)
RATES(10) = (( CONSTS(16)*STATES(10))/1.00000)*(1.00000 - CONSTS(26)) - ALGBRC(1)*STATES(10)
RATES(12) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(16)*STATES(10))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(18)*STATES(11))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(11)
RATES(13) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(18)*STATES(11))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(19)*STATES(13))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(13)
RATES(14) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(19)*STATES(13))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(20)*STATES(14))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(14)
RATES(15) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(20)*STATES(15))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(21)*STATES(15))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(15)
RATES(16) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(21)*STATES(15))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(22)*STATES(16))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(16)
RATES(17) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(22)*STATES(16))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(23)*STATES(17))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(17)
RATES(18) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(23)*STATES(17))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(24)*STATES(18))/1.00000)*(1.00000 - CONSTS(26))) - ALGBRC(1)*STATES(18)
RATES(19) = ( (( CONSTS(2)*CONSTS(3)*CONSTS(24)*STATES(18))/1.00000)*(1.00000 - CONSTS(4))+ (( CONSTS(25)*STATES(19))/1.00000)*((1.00000 - CONSTS(26))+ (( CONSTS(2)*CONSTS(3))/1.00000)*(1.00000 - CONSTS(4)))) - ALGBRC(1)*STATES(9)
RATES(20) = (( (1.00000 - CONSTS(2))*CONSTS(3))/1.00000)*( (1.00000 - CONSTS(1))*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5)+ CONSTS(22)*STATES(6)+ CONSTS(23)*STATES(7)+ CONSTS(24)*STATES(8)+ CONSTS(25)*STATES(9))+ (1.00000 - CONSTS(4))*( CONSTS(16)*STATES(10)+ CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15)+ CONSTS(22)*STATES(16)+ CONSTS(23)*STATES(17)+ CONSTS(24)*STATES(18)+ CONSTS(25)*STATES(19))) - ALGBRC(1)*STATES(20)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = ((1.00000 - CONSTS(17))/1.00000)*( CONSTS(16)*STATES(1)+ CONSTS(18)*STATES(2)+ CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(4)+ CONSTS(21)*STATES(5))+ ((1.00000 - CONSTS(26))/1.00000)*( CONSTS(18)*STATES(11)+ CONSTS(19)*STATES(13)+ CONSTS(20)*STATES(14)+ CONSTS(21)*STATES(15))
ALGBRC(2) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)
ALGBRC(3) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)
ALGBRC(4) = STATES(10)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)
RETURN
END