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 32 entries in the algebraic variable array.
C There are a total of 23 entries in each of the rate and state variable arrays.
C There are a total of 50 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is EGF in component EGF (nanomolar).
C ALGBRC(1) is v1 in component v1 (flux).
C STATES(2) is R in component R (nanomolar).
C STATES(3) is Ra in component Ra (nanomolar).
C ALGBRC(9) is v2 in component v2 (flux).
C STATES(4) is R2 in component R2 (nanomolar).
C ALGBRC(10) is v3 in component v3 (flux).
C ALGBRC(11) is v4 in component v4 (flux).
C STATES(5) is RP in component RP (nanomolar).
C ALGBRC(14) is v7 in component v7 (flux).
C ALGBRC(19) is v11 in component v11 (flux).
C ALGBRC(23) is v15 in component v15 (flux).
C ALGBRC(26) is v18 in component v18 (flux).
C ALGBRC(28) is v20 in component v20 (flux).
C ALGBRC(12) is v5 in component v5 (flux).
C ALGBRC(15) is v9 in component v9 (flux).
C ALGBRC(21) is v13 in component v13 (flux).
C STATES(6) is R_PL in component R_PL (nanomolar).
C ALGBRC(13) is v6 in component v6 (flux).
C STATES(7) is R_PLP in component R_PLP (nanomolar).
C STATES(8) is R_G in component R_G (nanomolar).
C ALGBRC(17) is v10 in component v10 (flux).
C STATES(9) is R_G_S in component R_G_S (nanomolar).
C STATES(10) is R_Sh in component R_Sh (nanomolar).
C ALGBRC(22) is v14 in component v14 (flux).
C STATES(11) is R_ShP in component R_ShP (nanomolar).
C ALGBRC(32) is v24 in component v24 (flux).
C ALGBRC(25) is v17 in component v17 (flux).
C STATES(12) is R_Sh_G in component R_Sh_G (nanomolar).
C ALGBRC(27) is v19 in component v19 (flux).
C STATES(13) is R_Sh_G_S in component R_Sh_G_S (nanomolar).
C STATES(14) is G_S in component G_S (nanomolar).
C ALGBRC(31) is v23 in component v23 (flux).
C ALGBRC(20) is v12 in component v12 (flux).
C STATES(15) is ShP in component ShP (nanomolar).
C ALGBRC(29) is v21 in component v21 (flux).
C ALGBRC(24) is v16 in component v16 (flux).
C STATES(16) is Sh_G in component Sh_G (nanomolar).
C ALGBRC(30) is v22 in component v22 (flux).
C STATES(17) is Sh_G_S in component Sh_G_S (nanomolar).
C STATES(18) is PLC_gamma in component PLC_gamma (nanomolar).
C ALGBRC(16) is v8 in component v8 (flux).
C STATES(19) is PLC_gamma_P in component PLC_gamma_P (nanomolar).
C ALGBRC(18) is v25 in component v25 (flux).
C STATES(20) is PLC_gamma_P_I in component PLC_gamma_P_I (nanomolar).
C STATES(21) is Grb in component Grb (nanomolar).
C STATES(22) is Shc in component Shc (nanomolar).
C STATES(23) is SOS in component SOS (nanomolar).
C CONSTS(1) is k1 in component v1 (second_order_rate_constant).
C CONSTS(2) is k1_ in component v1 (first_order_rate_constant).
C CONSTS(3) is k2 in component v2 (second_order_rate_constant).
C CONSTS(4) is k2_ in component v2 (first_order_rate_constant).
C CONSTS(5) is k3 in component v3 (first_order_rate_constant).
C CONSTS(6) is k3_ in component v3 (first_order_rate_constant).
C CONSTS(7) is K4 in component v4 (nanomolar).
C CONSTS(8) is V4 in component v4 (flux).
C CONSTS(9) is k5 in component v5 (second_order_rate_constant).
C CONSTS(10) is k5_ in component v5 (first_order_rate_constant).
C CONSTS(11) is k6 in component v6 (first_order_rate_constant).
C CONSTS(12) is k6_ in component v6 (first_order_rate_constant).
C CONSTS(13) is k7 in component v7 (first_order_rate_constant).
C CONSTS(14) is k7_ in component v7 (second_order_rate_constant).
C CONSTS(15) is K8 in component v8 (nanomolar).
C CONSTS(16) is V8 in component v8 (flux).
C CONSTS(17) is k9 in component v9 (second_order_rate_constant).
C CONSTS(18) is k9_ in component v9 (first_order_rate_constant).
C CONSTS(19) is k10 in component v10 (second_order_rate_constant).
C CONSTS(20) is k10_ in component v10 (first_order_rate_constant).
C CONSTS(21) is k11 in component v11 (first_order_rate_constant).
C CONSTS(22) is k11_ in component v11 (second_order_rate_constant).
C CONSTS(23) is k12 in component v12 (first_order_rate_constant).
C CONSTS(24) is k12_ in component v12 (second_order_rate_constant).
C CONSTS(25) is k13 in component v13 (second_order_rate_constant).
C CONSTS(26) is k13_ in component v13 (first_order_rate_constant).
C CONSTS(27) is k14 in component v14 (first_order_rate_constant).
C CONSTS(28) is k14_ in component v14 (first_order_rate_constant).
C CONSTS(29) is k15 in component v15 (first_order_rate_constant).
C CONSTS(30) is k15_ in component v15 (second_order_rate_constant).
C CONSTS(31) is K16 in component v16 (nanomolar).
C CONSTS(32) is V16 in component v16 (flux).
C CONSTS(33) is k17 in component v17 (second_order_rate_constant).
C CONSTS(34) is k17_ in component v17 (first_order_rate_constant).
C CONSTS(35) is k18 in component v18 (first_order_rate_constant).
C CONSTS(36) is k18_ in component v18 (second_order_rate_constant).
C CONSTS(37) is k19 in component v19 (second_order_rate_constant).
C CONSTS(38) is k19_ in component v19 (first_order_rate_constant).
C CONSTS(39) is k20 in component v20 (first_order_rate_constant).
C CONSTS(40) is k20_ in component v20 (second_order_rate_constant).
C CONSTS(41) is k21 in component v21 (second_order_rate_constant).
C CONSTS(42) is k21_ in component v21 (first_order_rate_constant).
C CONSTS(43) is k22 in component v22 (second_order_rate_constant).
C CONSTS(44) is k22_ in component v22 (first_order_rate_constant).
C CONSTS(45) is k23 in component v23 (first_order_rate_constant).
C CONSTS(46) is k23_ in component v23 (second_order_rate_constant).
C CONSTS(47) is k24 in component v24 (second_order_rate_constant).
C CONSTS(48) is k24_ in component v24 (first_order_rate_constant).
C CONSTS(49) is k25 in component v25 (first_order_rate_constant).
C CONSTS(50) is k25_ in component v25 (first_order_rate_constant).
C ALGBRC(2) is totEGFRphos in component combined_concs (nanomolar).
C ALGBRC(3) is totPLCgammaphos in component combined_concs (nanomolar).
C ALGBRC(4) is totGrb_EGFR in component combined_concs (nanomolar).
C ALGBRC(5) is totGrb_Shc in component combined_concs (nanomolar).
C ALGBRC(6) is totShcphos in component combined_concs (nanomolar).
C ALGBRC(7) is totShc_EGFR in component combined_concs (nanomolar).
C ALGBRC(8) is totSOS_EGFR in component combined_concs (nanomolar).
C RATES(1) is d/dt EGF in component EGF (nanomolar).
C RATES(2) is d/dt R in component R (nanomolar).
C RATES(3) is d/dt Ra in component Ra (nanomolar).
C RATES(4) is d/dt R2 in component R2 (nanomolar).
C RATES(5) is d/dt RP in component RP (nanomolar).
C RATES(6) is d/dt R_PL in component R_PL (nanomolar).
C RATES(7) is d/dt R_PLP in component R_PLP (nanomolar).
C RATES(8) is d/dt R_G in component R_G (nanomolar).
C RATES(9) is d/dt R_G_S in component R_G_S (nanomolar).
C RATES(10) is d/dt R_Sh in component R_Sh (nanomolar).
C RATES(11) is d/dt R_ShP in component R_ShP (nanomolar).
C RATES(12) is d/dt R_Sh_G in component R_Sh_G (nanomolar).
C RATES(13) is d/dt R_Sh_G_S in component R_Sh_G_S (nanomolar).
C RATES(14) is d/dt G_S in component G_S (nanomolar).
C RATES(15) is d/dt ShP in component ShP (nanomolar).
C RATES(16) is d/dt Sh_G in component Sh_G (nanomolar).
C RATES(17) is d/dt Sh_G_S in component Sh_G_S (nanomolar).
C RATES(18) is d/dt PLC_gamma in component PLC_gamma (nanomolar).
C RATES(19) is d/dt PLC_gamma_P in component PLC_gamma_P (nanomolar).
C RATES(20) is d/dt PLC_gamma_P_I in component PLC_gamma_P_I (nanomolar).
C RATES(21) is d/dt Grb in component Grb (nanomolar).
C RATES(22) is d/dt Shc in component Shc (nanomolar).
C RATES(23) is d/dt SOS in component SOS (nanomolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 680
STATES(2) = 100
STATES(3) = 0
STATES(4) = 0
STATES(5) = 0
STATES(6) = 0
STATES(7) = 0
STATES(8) = 0
STATES(9) = 0
STATES(10) = 0
STATES(11) = 0
STATES(12) = 0
STATES(13) = 0
STATES(14) = 0
STATES(15) = 0
STATES(16) = 0
STATES(17) = 0
STATES(18) = 105
STATES(19) = 0
STATES(20) = 0
STATES(21) = 85
STATES(22) = 150
STATES(23) = 34
CONSTS(1) = 0.003
CONSTS(2) = 0.06
CONSTS(3) = 0.01
CONSTS(4) = 0.1
CONSTS(5) = 1
CONSTS(6) = 0.01
CONSTS(7) = 50
CONSTS(8) = 450
CONSTS(9) = 0.06
CONSTS(10) = 0.2
CONSTS(11) = 1
CONSTS(12) = 0.05
CONSTS(13) = 0.3
CONSTS(14) = 0.006
CONSTS(15) = 100
CONSTS(16) = 1
CONSTS(17) = 0.003
CONSTS(18) = 0.05
CONSTS(19) = 0.01
CONSTS(20) = 0.06
CONSTS(21) = 0.03
CONSTS(22) = 4.5e-3
CONSTS(23) = 1.5e-3
CONSTS(24) = 1e-4
CONSTS(25) = 0.09
CONSTS(26) = 0.6
CONSTS(27) = 6
CONSTS(28) = 0.06
CONSTS(29) = 0.3
CONSTS(30) = 9e-4
CONSTS(31) = 340
CONSTS(32) = 1.7
CONSTS(33) = 0.003
CONSTS(34) = 0.1
CONSTS(35) = 0.3
CONSTS(36) = 9e-4
CONSTS(37) = 0.01
CONSTS(38) = 2.14e-2
CONSTS(39) = 0.12
CONSTS(40) = 2.4e-4
CONSTS(41) = 0.003
CONSTS(42) = 0.1
CONSTS(43) = 0.03
CONSTS(44) = 0.064
CONSTS(45) = 0.1
CONSTS(46) = 0.021
CONSTS(47) = 0.009
CONSTS(48) = 4.29e-2
CONSTS(49) = 1
CONSTS(50) = 0.03
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = CONSTS(1)*STATES(2)*STATES(1) - CONSTS(2)*STATES(3)
RATES(1) = - ALGBRC(1)
RATES(2) = - ALGBRC(1)
ALGBRC(9) = CONSTS(3)*STATES(3)*STATES(3) - CONSTS(4)*STATES(4)
RATES(3) = ALGBRC(1) - 2.00000*ALGBRC(9)
ALGBRC(10) = CONSTS(5)*STATES(4) - CONSTS(6)*STATES(5)
ALGBRC(11) = ( CONSTS(8)*STATES(5))/(CONSTS(7)+STATES(5))
RATES(4) = (ALGBRC(9)+ALGBRC(11)) - ALGBRC(10)
ALGBRC(12) = CONSTS(9)*STATES(5)*STATES(18) - CONSTS(10)*STATES(6)
ALGBRC(13) = CONSTS(11)*STATES(6) - CONSTS(12)*STATES(7)
RATES(6) = ALGBRC(12) - ALGBRC(13)
ALGBRC(14) = CONSTS(13)*STATES(7) - CONSTS(14)*STATES(5)*STATES(19)
RATES(7) = ALGBRC(13) - ALGBRC(14)
ALGBRC(16) = ( CONSTS(16)*STATES(19))/(CONSTS(15)+STATES(19))
RATES(18) = ALGBRC(16) - ALGBRC(12)
ALGBRC(15) = CONSTS(17)*STATES(5)*STATES(21) - CONSTS(18)*STATES(8)
ALGBRC(17) = CONSTS(19)*STATES(8)*STATES(23) - CONSTS(20)*STATES(9)
RATES(8) = ALGBRC(15) - ALGBRC(17)
ALGBRC(18) = CONSTS(49)*STATES(19) - CONSTS(50)*STATES(20)
RATES(19) = ALGBRC(14) - (ALGBRC(16)+ALGBRC(18))
RATES(20) = ALGBRC(18)
ALGBRC(19) = CONSTS(21)*STATES(9) - CONSTS(22)*STATES(5)*STATES(14)
RATES(9) = ALGBRC(17) - ALGBRC(19)
ALGBRC(21) = CONSTS(25)*STATES(5)*STATES(22) - CONSTS(26)*STATES(10)
ALGBRC(22) = CONSTS(27)*STATES(10) - CONSTS(28)*STATES(11)
RATES(10) = ALGBRC(21) - ALGBRC(22)
ALGBRC(24) = ( CONSTS(32)*STATES(15))/(CONSTS(31)+STATES(15))
RATES(22) = ALGBRC(24) - ALGBRC(21)
ALGBRC(26) = CONSTS(35)*STATES(12) - CONSTS(36)*STATES(5)*STATES(16)
ALGBRC(25) = CONSTS(33)*STATES(11)*STATES(21) - CONSTS(34)*STATES(12)
ALGBRC(27) = CONSTS(37)*STATES(12)*STATES(23) - CONSTS(38)*STATES(13)
RATES(12) = ALGBRC(25) - (ALGBRC(26)+ALGBRC(27))
ALGBRC(23) = CONSTS(29)*STATES(11) - CONSTS(30)*STATES(15)*STATES(5)
ALGBRC(28) = CONSTS(39)*STATES(13) - CONSTS(40)*STATES(17)*STATES(5)
RATES(5) = (ALGBRC(10)+ALGBRC(14)+ALGBRC(19)+ALGBRC(23)+ALGBRC(26)+ALGBRC(28)) - (ALGBRC(11)+ALGBRC(12)+ALGBRC(15)+ALGBRC(21))
ALGBRC(20) = CONSTS(23)*STATES(14) - CONSTS(24)*STATES(21)*STATES(23)
ALGBRC(29) = CONSTS(41)*STATES(15)*STATES(21) - CONSTS(42)*STATES(16)
RATES(21) = ALGBRC(20) - (ALGBRC(15)+ALGBRC(25)+ALGBRC(29))
ALGBRC(30) = CONSTS(43)*STATES(16)*STATES(23) - CONSTS(44)*STATES(17)
RATES(16) = (ALGBRC(26)+ALGBRC(29)) - ALGBRC(30)
RATES(23) = ALGBRC(20) - (ALGBRC(17)+ALGBRC(27)+ALGBRC(30))
ALGBRC(31) = CONSTS(45)*STATES(17) - CONSTS(46)*STATES(15)*STATES(14)
RATES(15) = (ALGBRC(23)+ALGBRC(31)) - (ALGBRC(29)+ALGBRC(24))
RATES(17) = (ALGBRC(28)+ALGBRC(30)) - ALGBRC(31)
ALGBRC(32) = CONSTS(47)*STATES(11)*STATES(14) - CONSTS(48)*STATES(13)
RATES(11) = ALGBRC(22) - (ALGBRC(32)+ALGBRC(23)+ALGBRC(25))
RATES(13) = (ALGBRC(27)+ALGBRC(32)) - ALGBRC(28)
RATES(14) = (ALGBRC(19)+ALGBRC(31)) - (ALGBRC(20)+ALGBRC(32))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = CONSTS(1)*STATES(2)*STATES(1) - CONSTS(2)*STATES(3)
ALGBRC(9) = CONSTS(3)*STATES(3)*STATES(3) - CONSTS(4)*STATES(4)
ALGBRC(10) = CONSTS(5)*STATES(4) - CONSTS(6)*STATES(5)
ALGBRC(11) = ( CONSTS(8)*STATES(5))/(CONSTS(7)+STATES(5))
ALGBRC(12) = CONSTS(9)*STATES(5)*STATES(18) - CONSTS(10)*STATES(6)
ALGBRC(13) = CONSTS(11)*STATES(6) - CONSTS(12)*STATES(7)
ALGBRC(14) = CONSTS(13)*STATES(7) - CONSTS(14)*STATES(5)*STATES(19)
ALGBRC(16) = ( CONSTS(16)*STATES(19))/(CONSTS(15)+STATES(19))
ALGBRC(15) = CONSTS(17)*STATES(5)*STATES(21) - CONSTS(18)*STATES(8)
ALGBRC(17) = CONSTS(19)*STATES(8)*STATES(23) - CONSTS(20)*STATES(9)
ALGBRC(18) = CONSTS(49)*STATES(19) - CONSTS(50)*STATES(20)
ALGBRC(19) = CONSTS(21)*STATES(9) - CONSTS(22)*STATES(5)*STATES(14)
ALGBRC(21) = CONSTS(25)*STATES(5)*STATES(22) - CONSTS(26)*STATES(10)
ALGBRC(22) = CONSTS(27)*STATES(10) - CONSTS(28)*STATES(11)
ALGBRC(24) = ( CONSTS(32)*STATES(15))/(CONSTS(31)+STATES(15))
ALGBRC(26) = CONSTS(35)*STATES(12) - CONSTS(36)*STATES(5)*STATES(16)
ALGBRC(25) = CONSTS(33)*STATES(11)*STATES(21) - CONSTS(34)*STATES(12)
ALGBRC(27) = CONSTS(37)*STATES(12)*STATES(23) - CONSTS(38)*STATES(13)
ALGBRC(23) = CONSTS(29)*STATES(11) - CONSTS(30)*STATES(15)*STATES(5)
ALGBRC(28) = CONSTS(39)*STATES(13) - CONSTS(40)*STATES(17)*STATES(5)
ALGBRC(20) = CONSTS(23)*STATES(14) - CONSTS(24)*STATES(21)*STATES(23)
ALGBRC(29) = CONSTS(41)*STATES(15)*STATES(21) - CONSTS(42)*STATES(16)
ALGBRC(30) = CONSTS(43)*STATES(16)*STATES(23) - CONSTS(44)*STATES(17)
ALGBRC(31) = CONSTS(45)*STATES(17) - CONSTS(46)*STATES(15)*STATES(14)
ALGBRC(32) = CONSTS(47)*STATES(11)*STATES(14) - CONSTS(48)*STATES(13)
ALGBRC(2) = 2.00000*(STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13))
ALGBRC(3) = STATES(7)+STATES(19)
ALGBRC(4) = STATES(8)+STATES(9)+STATES(12)+STATES(13)
ALGBRC(5) = STATES(12)+STATES(16)+STATES(13)+STATES(17)
ALGBRC(6) = STATES(11)+STATES(12)+STATES(13)+STATES(15)+STATES(16)+STATES(17)
ALGBRC(7) = STATES(11)+STATES(12)+STATES(13)
ALGBRC(8) = STATES(9)+STATES(13)
RETURN
END