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