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