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 18 entries in the algebraic variable array. C There are a total of 16 entries in each of the rate and state variable arrays. C There are a total of 76 entries in the constant variable array. C C C VOI is time in component environment (minute). C STATES(1) is N in component N (nanomolar). C CONSTS(1) is KdN in component N (nanomolar). C CONSTS(2) is vsN in component N (flux). C CONSTS(3) is vdN in component N (flux). C CONSTS(4) is kc in component model_parameters (first_order_rate_constant). C CONSTS(5) is KIF in component model_parameters (nanomolar). C STATES(2) is F in component F (nanomolar). C CONSTS(6) is epsilon in component model_parameters (dimensionless). C CONSTS(7) is j in component model_parameters (dimensionless). C STATES(3) is Na in component Na (nanomolar). C ALGBRC(1) is Vtr in component Na (flux). C CONSTS(8) is KdNa in component Na (nanomolar). C CONSTS(9) is VdNa in component Na (flux). C CONSTS(10) is kt1 in component Na (first_order_rate_constant). C CONSTS(11) is kt2 in component Na (first_order_rate_constant). C STATES(4) is Nan in component Nan (nanomolar). C CONSTS(12) is KdNan in component Nan (nanomolar). C CONSTS(13) is VdNan in component Nan (flux). C STATES(5) is MF in component MF (nanomolar). C CONSTS(14) is KdMF in component MF (nanomolar). C CONSTS(15) is KIG1 in component MF (nanomolar). C ALGBRC(2) is vsFK in component MF (flux). C CONSTS(16) is vsF in component MF (flux). C CONSTS(17) is vmF in component MF (flux). C CONSTS(18) is KA in component MF (nanomolar). C STATES(6) is K in component K (nanomolar). C CONSTS(19) is p in component model_parameters (dimensionless). C CONSTS(20) is KdF in component F (nanomolar). C CONSTS(21) is vdF in component F (flux). C CONSTS(22) is ksF in component F (first_order_rate_constant). C ALGBRC(8) is V1 in component Wnt_parameters (flux). C CONSTS(23) is theta in component model_parameters (dimensionless). C STATES(7) is B in component B (nanomolar). C CONSTS(24) is kd1 in component B (first_order_rate_constant). C CONSTS(25) is vsB in component B (flux). C ALGBRC(14) is VK in component Wnt_parameters (flux). C ALGBRC(18) is VP in component Wnt_parameters (flux). C ALGBRC(9) is V2 in component Wnt_parameters (flux). C CONSTS(26) is Kt in component Wnt_parameters (nanomolar). C ALGBRC(3) is AK in component Wnt_parameters (nanomolar). C STATES(8) is Bp in component Bp (nanomolar). C CONSTS(27) is kd2 in component Bp (first_order_rate_constant). C STATES(9) is BN in component BN (nanomolar). C STATES(10) is MAx in component MAx (nanomolar). C CONSTS(28) is v0 in component MAx (flux). C CONSTS(29) is vMB in component MAx (flux). C CONSTS(30) is vmd in component MAx (flux). C CONSTS(31) is KaB in component MAx (nanomolar). C CONSTS(32) is KaXa in component MAx (nanomolar). C CONSTS(33) is Kmd in component MAx (nanomolar). C CONSTS(34) is n in component MAx (dimensionless). C CONSTS(35) is m in component MAx (dimensionless). C CONSTS(36) is vMXa in component MAx (flux). C STATES(11) is Xa in component Xa (nanomolar). C STATES(12) is A in component A (nanomolar). C CONSTS(37) is ksAx in component A (first_order_rate_constant). C CONSTS(38) is vdAx in component A (flux). C CONSTS(39) is KdAx in component A (nanomolar). C CONSTS(40) is d1 in component Wnt_parameters (first_order_rate_constant). C CONSTS(41) is a1 in component Wnt_parameters (second_order_rate_constant). C CONSTS(42) is K1 in component Wnt_parameters (nanomolar). C CONSTS(43) is K2 in component Wnt_parameters (nanomolar). C CONSTS(44) is D in component Wnt_parameters (nanomolar). C CONSTS(45) is KID in component Wnt_parameters (nanomolar). C CONSTS(46) is kt3 in component Wnt_parameters (first_order_rate_constant). C CONSTS(47) is kt4 in component Wnt_parameters (first_order_rate_constant). C CONSTS(48) is VMK in component Wnt_parameters (flux). C CONSTS(49) is VMP in component Wnt_parameters (flux). C STATES(13) is Rasa in component Rasa (nanomolar). C ALGBRC(10) is VaRas in component FGF_parameters (flux). C ALGBRC(15) is VdRas in component FGF_parameters (flux). C CONSTS(50) is eta in component model_parameters (dimensionless). C STATES(14) is ERKa in component ERKa (nanomolar). C ALGBRC(11) is VaErk in component FGF_parameters (flux). C ALGBRC(16) is VdErk in component FGF_parameters (flux). C ALGBRC(12) is VaX in component FGF_parameters (flux). C ALGBRC(17) is VdX in component FGF_parameters (flux). C STATES(15) is MDusp in component MDusp (nanomolar). C ALGBRC(7) is VsMDusp in component FGF_parameters (flux). C ALGBRC(13) is VdMDusp in component FGF_parameters (flux). C STATES(16) is Dusp in component Dusp (nanomolar). C CONSTS(51) is ksDusp in component Dusp (first_order_rate_constant). C CONSTS(52) is vdDusp in component Dusp (flux). C CONSTS(53) is KdDusp in component Dusp (nanomolar). C ALGBRC(4) is Rasi in component FGF_parameters (nanomolar). C ALGBRC(5) is ERKi in component FGF_parameters (nanomolar). C ALGBRC(6) is Xi in component FGF_parameters (nanomolar). C CONSTS(54) is Rast in component FGF_parameters (nanomolar). C CONSTS(55) is ERKt in component FGF_parameters (nanomolar). C CONSTS(56) is Xt in component FGF_parameters (nanomolar). C CONSTS(57) is kcDusp in component FGF_parameters (first_order_rate_constant). C CONSTS(58) is VMaRas in component FGF_parameters (flux). C CONSTS(59) is VMdRas in component FGF_parameters (flux). C CONSTS(60) is VMaErk in component FGF_parameters (flux). C CONSTS(61) is VMaX in component FGF_parameters (flux). C CONSTS(62) is VMdX in component FGF_parameters (flux). C CONSTS(63) is VMsMDusp in component FGF_parameters (flux). C CONSTS(64) is VMdMDusp in component FGF_parameters (flux). C CONSTS(65) is Fgf in component FGF_parameters (nanomolar). C CONSTS(66) is KaFgf in component FGF_parameters (nanomolar). C CONSTS(67) is KaRas in component FGF_parameters (nanomolar). C CONSTS(68) is KdRas in component FGF_parameters (nanomolar). C CONSTS(69) is KdErk in component FGF_parameters (nanomolar). C CONSTS(70) is KaErk in component FGF_parameters (nanomolar). C CONSTS(71) is KaX in component FGF_parameters (nanomolar). C CONSTS(72) is KdX in component FGF_parameters (nanomolar). C CONSTS(73) is KaMDusp in component FGF_parameters (nanomolar). C CONSTS(74) is KdMDusp in component FGF_parameters (nanomolar). C CONSTS(75) is q in component FGF_parameters (dimensionless). C CONSTS(76) is r in component FGF_parameters (dimensionless). C RATES(1) is d/dt N in component N (nanomolar). C RATES(3) is d/dt Na in component Na (nanomolar). C RATES(4) is d/dt Nan in component Nan (nanomolar). C RATES(5) is d/dt MF in component MF (nanomolar). C RATES(2) is d/dt F in component F (nanomolar). C RATES(6) is d/dt K in component K (nanomolar). C RATES(7) is d/dt B in component B (nanomolar). C RATES(8) is d/dt Bp in component Bp (nanomolar). C RATES(9) is d/dt BN in component BN (nanomolar). C RATES(10) is d/dt MAx in component MAx (nanomolar). C RATES(12) is d/dt A in component A (nanomolar). C RATES(13) is d/dt Rasa in component Rasa (nanomolar). C RATES(14) is d/dt ERKa in component ERKa (nanomolar). C RATES(11) is d/dt Xa in component Xa (nanomolar). C RATES(15) is d/dt MDusp in component MDusp (nanomolar). C RATES(16) is d/dt Dusp in component Dusp (nanomolar). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.5 CONSTS(1) = 1.4 CONSTS(2) = 0.23 CONSTS(3) = 2.82 CONSTS(4) = 3.45 CONSTS(5) = 0.5 STATES(2) = 0.001 CONSTS(6) = 0.3 CONSTS(7) = 2.0 STATES(3) = 0.2 CONSTS(8) = 0.001 CONSTS(9) = 0.01 CONSTS(10) = 0.1 CONSTS(11) = 0.1 STATES(4) = 0.0 CONSTS(12) = 0.001 CONSTS(13) = 0.1 STATES(5) = 0.1 CONSTS(14) = 0.768 CONSTS(15) = 2.5 CONSTS(16) = 3.0 CONSTS(17) = 1.92 CONSTS(18) = 0.05 STATES(6) = 3.0 CONSTS(19) = 2.0 CONSTS(20) = 0.37 CONSTS(21) = 0.39 CONSTS(22) = 0.3 CONSTS(23) = 1.5 STATES(7) = 0.1 CONSTS(24) = 0.0 CONSTS(25) = 0.087 CONSTS(26) = 3.0 STATES(8) = 0.1 CONSTS(27) = 7.062 STATES(9) = 0.001 STATES(10) = 0.1 CONSTS(28) = 0.06 CONSTS(29) = 1.64 CONSTS(30) = 0.8 CONSTS(31) = 0.7 CONSTS(32) = 0.05 CONSTS(33) = 0.48 CONSTS(34) = 2.0 CONSTS(35) = 2.0 CONSTS(36) = 0.5 STATES(11) = 0.1 STATES(12) = 0.1 CONSTS(37) = 0.02 CONSTS(38) = 0.6 CONSTS(39) = 0.63 CONSTS(40) = 0.1 CONSTS(41) = 1.8 CONSTS(42) = 0.28 CONSTS(43) = 0.03 CONSTS(44) = 2.0 CONSTS(45) = 0.5 CONSTS(46) = 0.7 CONSTS(47) = 1.5 CONSTS(48) = 5.08 CONSTS(49) = 1.0 STATES(13) = 0.5 CONSTS(50) = 0.3 STATES(14) = 0.2 STATES(15) = 0.1 STATES(16) = 0.1 CONSTS(51) = 0.5 CONSTS(52) = 2.0 CONSTS(53) = 0.5 CONSTS(54) = 2.0 CONSTS(55) = 2.0 CONSTS(56) = 2.0 CONSTS(57) = 1.35 CONSTS(58) = 4.968 CONSTS(59) = 0.41 CONSTS(60) = 3.30 CONSTS(61) = 1.6 CONSTS(62) = 0.5 CONSTS(63) = 0.9 CONSTS(64) = 0.5 CONSTS(65) = 1.0 CONSTS(66) = 0.5 CONSTS(67) = 0.103 CONSTS(68) = 0.1 CONSTS(69) = 0.05 CONSTS(70) = 0.05 CONSTS(71) = 0.05 CONSTS(72) = 0.05 CONSTS(73) = 0.5 CONSTS(74) = 0.5 CONSTS(75) = 2.0 CONSTS(76) = 2.0 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = CONSTS(6)*(CONSTS(2) - ( CONSTS(3)*(STATES(1)/(CONSTS(1)+STATES(1)))+ CONSTS(4)*STATES(1)*(CONSTS(5) ** CONSTS(7)/(CONSTS(5) ** CONSTS(7)+STATES(2) ** CONSTS(7))))) RATES(2) = CONSTS(6)*( CONSTS(22)*STATES(5) - CONSTS(21)*(STATES(2)/(CONSTS(20)+STATES(2)))) RATES(10) = CONSTS(23)*((CONSTS(28)+ CONSTS(29)*(STATES(9) ** CONSTS(34)/(CONSTS(31) ** CONSTS(34)+STATES(9) ** CONSTS(34)))+ CONSTS(36)*(STATES(11) ** CONSTS(35)/(CONSTS(32) ** CONSTS(35)+STATES(11) ** CONSTS(35)))) - CONSTS(30)*(STATES(10)/(CONSTS(33)+STATES(10)))) RATES(16) = CONSTS(50)*( CONSTS(51)*STATES(15) - CONSTS(52)*(STATES(16)/(CONSTS(53)+STATES(16)))) ALGBRC(1) = CONSTS(10)*STATES(3) - CONSTS(11)*STATES(4) RATES(3) = CONSTS(6)*( CONSTS(4)*STATES(1)*(CONSTS(5) ** CONSTS(7)/(CONSTS(5) ** CONSTS(7)+STATES(2) ** CONSTS(7))) - ( CONSTS(9)*(STATES(3)/(CONSTS(8)+STATES(3)))+ALGBRC(1))) RATES(4) = CONSTS(6)*(ALGBRC(1) - CONSTS(13)*(STATES(4)/(CONSTS(12)+STATES(4)))) ALGBRC(2) = CONSTS(16)*(CONSTS(15)/(CONSTS(15)+STATES(6))) RATES(5) = CONSTS(6)*( ALGBRC(2)*(STATES(4) ** CONSTS(19)/(CONSTS(18) ** CONSTS(19)+STATES(4) ** CONSTS(19))) - CONSTS(17)*(STATES(5)/(CONSTS(14)+STATES(5)))) ALGBRC(3) = CONSTS(26) - STATES(6) ALGBRC(8) = CONSTS(40)*ALGBRC(3) - CONSTS(41)*STATES(12)*STATES(6) RATES(6) = CONSTS(23)*ALGBRC(8) ALGBRC(9) = CONSTS(47)*STATES(9) - CONSTS(46)*STATES(7) RATES(9) = - ( CONSTS(23)*ALGBRC(9)) RATES(12) = CONSTS(23)*(( CONSTS(37)*STATES(10)+ALGBRC(8)) - CONSTS(38)*(STATES(12)/(CONSTS(39)+STATES(12)))) ALGBRC(7) = CONSTS(63)*(STATES(11) ** CONSTS(75)/(CONSTS(73) ** CONSTS(75)+STATES(11) ** CONSTS(75))) ALGBRC(13) = CONSTS(64)*(STATES(15)/(CONSTS(74)+STATES(15))) RATES(15) = CONSTS(50)*(ALGBRC(7) - ALGBRC(13)) ALGBRC(4) = CONSTS(54) - STATES(13) ALGBRC(10) = CONSTS(58)*(CONSTS(65) ** CONSTS(76)/(CONSTS(66) ** CONSTS(76)+CONSTS(65) ** CONSTS(76)))*(ALGBRC(4)/(CONSTS(67)+ALGBRC(4))) ALGBRC(15) = CONSTS(59)*(STATES(13)/(CONSTS(68)+STATES(13))) RATES(13) = CONSTS(50)*(ALGBRC(10) - ALGBRC(15)) ALGBRC(5) = CONSTS(55) - STATES(14) ALGBRC(11) = CONSTS(60)*(STATES(13)/CONSTS(54))*(ALGBRC(5)/(CONSTS(70)+ALGBRC(5))) ALGBRC(16) = CONSTS(57)*STATES(16)*(STATES(14)/(CONSTS(69)+STATES(14))) RATES(14) = CONSTS(50)*(ALGBRC(11) - ALGBRC(16)) ALGBRC(6) = CONSTS(56) - STATES(11) ALGBRC(12) = CONSTS(61)*(STATES(14)/CONSTS(55))*(ALGBRC(6)/(CONSTS(71)+ALGBRC(6))) ALGBRC(17) = CONSTS(62)*(STATES(11)/(CONSTS(72)+STATES(11))) RATES(11) = CONSTS(50)*(ALGBRC(12) - ALGBRC(17)) ALGBRC(14) = CONSTS(48)*(CONSTS(45)/(CONSTS(45)+CONSTS(44)))*(STATES(7)/(CONSTS(42)+STATES(7))) ALGBRC(18) = CONSTS(49)*(STATES(8)/(CONSTS(43)+STATES(8))) RATES(7) = CONSTS(23)*((CONSTS(25)+ALGBRC(18)+ALGBRC(9)) - ( ALGBRC(14)*(ALGBRC(3)/CONSTS(26))+ CONSTS(24)*STATES(7))) RATES(8) = CONSTS(23)*( ALGBRC(14)*(ALGBRC(3)/CONSTS(26)) - (ALGBRC(18)+ CONSTS(27)*STATES(8))) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = CONSTS(10)*STATES(3) - CONSTS(11)*STATES(4) ALGBRC(2) = CONSTS(16)*(CONSTS(15)/(CONSTS(15)+STATES(6))) ALGBRC(3) = CONSTS(26) - STATES(6) ALGBRC(8) = CONSTS(40)*ALGBRC(3) - CONSTS(41)*STATES(12)*STATES(6) ALGBRC(9) = CONSTS(47)*STATES(9) - CONSTS(46)*STATES(7) ALGBRC(7) = CONSTS(63)*(STATES(11) ** CONSTS(75)/(CONSTS(73) ** CONSTS(75)+STATES(11) ** CONSTS(75))) ALGBRC(13) = CONSTS(64)*(STATES(15)/(CONSTS(74)+STATES(15))) ALGBRC(4) = CONSTS(54) - STATES(13) ALGBRC(10) = CONSTS(58)*(CONSTS(65) ** CONSTS(76)/(CONSTS(66) ** CONSTS(76)+CONSTS(65) ** CONSTS(76)))*(ALGBRC(4)/(CONSTS(67)+ALGBRC(4))) ALGBRC(15) = CONSTS(59)*(STATES(13)/(CONSTS(68)+STATES(13))) ALGBRC(5) = CONSTS(55) - STATES(14) ALGBRC(11) = CONSTS(60)*(STATES(13)/CONSTS(54))*(ALGBRC(5)/(CONSTS(70)+ALGBRC(5))) ALGBRC(16) = CONSTS(57)*STATES(16)*(STATES(14)/(CONSTS(69)+STATES(14))) ALGBRC(6) = CONSTS(56) - STATES(11) ALGBRC(12) = CONSTS(61)*(STATES(14)/CONSTS(55))*(ALGBRC(6)/(CONSTS(71)+ALGBRC(6))) ALGBRC(17) = CONSTS(62)*(STATES(11)/(CONSTS(72)+STATES(11))) ALGBRC(14) = CONSTS(48)*(CONSTS(45)/(CONSTS(45)+CONSTS(44)))*(STATES(7)/(CONSTS(42)+STATES(7))) ALGBRC(18) = CONSTS(49)*(STATES(8)/(CONSTS(43)+STATES(8))) RETURN END