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