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 11 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 t in component environment (s).
C CONSTS(1) is kf5 in component RAF (per_um_s).
C CONSTS(2) is kb5 in component RAF (per_s).
C CONSTS(3) is km1 in component RAF (um).
C CONSTS(4) is km2 in component RAF (um).
C CONSTS(5) is km3 in component RAF (um).
C CONSTS(6) is km4 in component RAF (um).
C CONSTS(7) is Vmax1 in component RAF (per_s).
C CONSTS(8) is Vmax2 in component RAF (per_s).
C CONSTS(9) is Vmax3 in component RAF (per_s).
C CONSTS(10) is Vmax4 in component RAF (per_s).
C CONSTS(11) is PKC in component undefinedvariables (um).
C CONSTS(12) is PP2A in component RAF (um).
C STATES(1) is MAPK_star in component MAPK (um).
C STATES(2) is Raf in component RAF (um).
C STATES(3) is Raf_star in component RAF (um).
C STATES(4) is Raf_star_star in component RAF (um).
C STATES(5) is GTPRas in component RAF (um).
C STATES(6) is GTPRasRaf_star in component RAF (um).
C CONSTS(13) is km6 in component MAPKK (um).
C CONSTS(14) is km7 in component MAPKK (um).
C CONSTS(15) is km8 in component MAPKK (um).
C CONSTS(16) is km9 in component MAPKK (um).
C CONSTS(17) is Vmax6 in component MAPKK (per_s).
C CONSTS(18) is Vmax7 in component MAPKK (per_s).
C CONSTS(19) is Vmax8 in component MAPKK (per_s).
C CONSTS(20) is Vmax9 in component MAPKK (per_s).
C CONSTS(21) is PP2A in component MAPKK (um).
C STATES(7) is MAPKK in component MAPKK (um).
C STATES(8) is MAPKK_star in component MAPKK (um).
C STATES(9) is MAPKK_star_star in component MAPKK (um).
C CONSTS(22) is km10 in component MAPK (um).
C CONSTS(23) is km11 in component MAPK (um).
C CONSTS(24) is km12 in component MAPK (um).
C CONSTS(25) is km13 in component MAPK (um).
C CONSTS(26) is Vmax10 in component MAPK (per_s).
C CONSTS(27) is Vmax11 in component MAPK (per_s).
C CONSTS(28) is Vmax12 in component MAPK (per_s).
C CONSTS(29) is Vmax13 in component MAPK (per_s).
C STATES(10) is MAPK in component MAPK (um).
C STATES(11) is MAPK_tyr in component MAPK (um).
C CONSTS(30) is MKP1 in component MAPK (um).
C RATES(6) is d/dt GTPRasRaf_star in component RAF (um).
C RATES(5) is d/dt GTPRas in component RAF (um).
C RATES(2) is d/dt Raf in component RAF (um).
C RATES(3) is d/dt Raf_star in component RAF (um).
C RATES(4) is d/dt Raf_star_star in component RAF (um).
C RATES(7) is d/dt MAPKK in component MAPKK (um).
C RATES(8) is d/dt MAPKK_star in component MAPKK (um).
C RATES(9) is d/dt MAPKK_star_star in component MAPKK (um).
C RATES(10) is d/dt MAPK in component MAPK (um).
C RATES(11) is d/dt MAPK_tyr in component MAPK (um).
C RATES(1) is d/dt MAPK_star in component MAPK (um).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.00004
CONSTS(2) = 0.5
CONSTS(3) = 66.666666667
CONSTS(4) = 25.64166667
CONSTS(5) = 15.6565
CONSTS(6) = 15.6565
CONSTS(7) = 4
CONSTS(8) = 10
CONSTS(9) = 6
CONSTS(10) = 6
CONSTS(11) = 0.1
CONSTS(12) = 0.224
STATES(1) = 0
STATES(2) = 0.1
STATES(3) = 0
STATES(4) = 0
STATES(5) = 0.1
STATES(6) = 0
CONSTS(13) = 0.159091667
CONSTS(14) = 0.159091667
CONSTS(15) = 15.6565
CONSTS(16) = 15.6565
CONSTS(17) = 0.105
CONSTS(18) = 0.105
CONSTS(19) = 6
CONSTS(20) = 6
CONSTS(21) = 0.224
STATES(7) = 0.18
STATES(8) = 0
STATES(9) = 0
CONSTS(22) = 0.046296667
CONSTS(23) = 0.046296667
CONSTS(24) = 0.066666667
CONSTS(25) = 0.066666667
CONSTS(26) = 0.15
CONSTS(27) = 0.15
CONSTS(28) = 1
CONSTS(29) = 1
STATES(10) = 0.36
STATES(11) = 0
CONSTS(30) = 0.032
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(6) = STATES(5)*STATES(3)*CONSTS(1) - STATES(6)*CONSTS(2)
RATES(5) = - STATES(5)*STATES(3)*CONSTS(1)+ STATES(6)*CONSTS(2)
RATES(2) = ( - CONSTS(11)*STATES(2)*CONSTS(7))/(CONSTS(3)+STATES(2))+( STATES(3)*CONSTS(12)*CONSTS(9))/(CONSTS(5)+STATES(3))
RATES(3) = ((( CONSTS(11)*STATES(2)*CONSTS(7))/(CONSTS(3)+STATES(2)) - ( STATES(3)*CONSTS(12)*CONSTS(9))/(CONSTS(5)+STATES(3))) - ( STATES(3)*STATES(1)*CONSTS(8))/(CONSTS(4)+STATES(3)))+( STATES(4)*CONSTS(12)*CONSTS(10))/(CONSTS(6)+STATES(4))
RATES(4) = ( STATES(3)*STATES(1)*CONSTS(8))/(CONSTS(4)+STATES(3)) - ( STATES(4)*CONSTS(12)*CONSTS(10))/(CONSTS(6)+STATES(4))
RATES(7) = ( - STATES(7)*CONSTS(17)*STATES(6))/(CONSTS(13)+STATES(7))+( STATES(8)*CONSTS(19)*CONSTS(21))/(CONSTS(15)+STATES(8))
RATES(8) = ((( STATES(7)*CONSTS(17)*STATES(6))/(CONSTS(13)+STATES(7)) - ( STATES(8)*CONSTS(19)*CONSTS(21))/(CONSTS(15)+STATES(8))) - ( STATES(8)*CONSTS(18)*STATES(6))/(CONSTS(14)+STATES(8)))+( STATES(9)*CONSTS(20)*CONSTS(21))/(CONSTS(16)+STATES(9))
RATES(9) = ( STATES(8)*CONSTS(18)*STATES(6))/(CONSTS(14)+STATES(8)) - ( STATES(9)*CONSTS(20)*CONSTS(21))/(CONSTS(16)+STATES(9))
RATES(10) = ( - STATES(10)*CONSTS(26)*STATES(9))/(CONSTS(22)+STATES(10))+( STATES(11)*CONSTS(28)*CONSTS(30))/(CONSTS(24)+STATES(11))
RATES(11) = ((( STATES(10)*CONSTS(26)*STATES(9))/(CONSTS(22)+STATES(10)) - ( STATES(11)*CONSTS(28)*CONSTS(30))/(CONSTS(24)+STATES(11))) - ( STATES(11)*CONSTS(27)*STATES(9))/(CONSTS(23)+STATES(11)))+( STATES(1)*CONSTS(29)*CONSTS(30))/(CONSTS(25)+STATES(1))
RATES(1) = ( STATES(11)*CONSTS(27)*STATES(9))/(CONSTS(23)+STATES(11)) - ( STATES(1)*CONSTS(29)*CONSTS(30))/(CONSTS(25)+STATES(1))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
