C
C There are a total of 5 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 55 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is V_max1 in component V1 (micromolar_per_minute).
C CONSTS(2) is GEFt in component V1 (micromolar).
C CONSTS(44) is V_1 in component V1 (per_minute).
C CONSTS(3) is Str in component V2 (dimensionless).
C CONSTS(4) is V_max2 in component V2 (micromolar_per_minute).
C CONSTS(45) is V_2 in component V2 (per_minute).
C CONSTS(5) is k_c3 in component V3 (per_minute).
C CONSTS(6) is PKAt in component V3 (micromolar).
C CONSTS(7) is GAPt in component V3 (micromolar).
C CONSTS(46) is V_3 in component V3 (per_minute).
C CONSTS(8) is V_max4 in component V4 (micromolar_per_minute).
C CONSTS(47) is V_4 in component V4 (per_minute).
C CONSTS(9) is k_gef in component V5 (per_minute).
C CONSTS(10) is RASt in component V5 (micromolar).
C CONSTS(48) is V_5 in component V5 (per_minute).
C CONSTS(11) is k_gap in component V6 (per_minute).
C CONSTS(49) is V_6 in component V6 (per_minute).
C CONSTS(12) is k_c7 in component V7 (per_minute).
C CONSTS(13) is PDEt in component V7 (micromolar).
C CONSTS(50) is V_7 in component V7 (per_minute).
C CONSTS(14) is V_max8 in component V8 (micromolar_per_minute).
C CONSTS(51) is V_8 in component V8 (per_minute).
C CONSTS(15) is a in component VPKAact (per_micromolar_squared_minute).
C CONSTS(16) is r in component VPKAact (per_minute).
C STATES(1) is R2C2 in component holoenzyme_R_C (dimensionless).
C STATES(2) is cAMP in component cyclic_AMP (micromolar).
C ALGBRC(4) is C in component C_subunit (dimensionless).
C ALGBRC(1) is R2cAMP2 in component holoenzyme_R_cAMP (dimensionless).
C ALGBRC(5) is V_PKAact in component VPKAact (per_minute).
C CONSTS(17) is K_1 in component active_GEF (dimensionless).
C CONSTS(18) is K_2 in component active_GEF (dimensionless).
C STATES(3) is GEFa in component active_GEF (dimensionless).
C CONSTS(19) is K_3 in component active_GAP (dimensionless).
C CONSTS(20) is K_4 in component active_GAP (dimensionless).
C STATES(4) is GAPa in component active_GAP (dimensionless).
C CONSTS(21) is K_5 in component RAS_to_GTP (dimensionless).
C CONSTS(22) is K_6 in component RAS_to_GTP (dimensionless).
C STATES(5) is RGTP in component RAS_to_GTP (dimensionless).
C CONSTS(23) is k_a in component adenylate_cyclase (per_micromolar_minute).
C CONSTS(24) is k_i in component adenylate_cyclase (per_minute).
C STATES(6) is CYCLa in component adenylate_cyclase (dimensionless).
C CONSTS(25) is K_7 in component active_PDE (dimensionless).
C CONSTS(26) is K_8 in component active_PDE (dimensionless).
C STATES(7) is PDEa in component active_PDE (dimensionless).
C CONSTS(27) is k_s in component cyclic_AMP (per_minute).
C CONSTS(28) is k_d in component cyclic_AMP (per_minute).
C CONSTS(29) is CYCLt in component cyclic_AMP (micromolar).
C CONSTS(30) is K_md in component cyclic_AMP (micromolar).
C CONSTS(31) is k_c9 in component V9 (per_minute).
C CONSTS(32) is MSNt in component V9 (micromolar).
C CONSTS(52) is V_9 in component V9 (per_minute).
C CONSTS(33) is V_max10 in component V10 (micromolar_per_minute).
C CONSTS(53) is V_10 in component V10 (per_minute).
C CONSTS(34) is k_c11 in component V11 (per_minute).
C CONSTS(54) is V_11 in component V11 (per_minute).
C CONSTS(35) is V_max12 in component V12 (micromolar_per_minute).
C CONSTS(55) is V_12 in component V12 (per_minute).
C CONSTS(36) is k_t1 in component cytosol (per_minute).
C CONSTS(37) is k_t2 in component cytosol (per_minute).
C CONSTS(38) is K_11 in component cytosol (dimensionless).
C CONSTS(39) is K_12 in component cytosol (dimensionless).
C STATES(8) is MN in component nucleus (dimensionless).
C STATES(9) is MCP in component cytosol_phos (dimensionless).
C STATES(10) is MC in component cytosol (dimensionless).
C CONSTS(40) is K_9 in component nucleus (dimensionless).
C CONSTS(41) is K_10 in component nucleus (dimensionless).
C STATES(11) is MNP in component nucleus_phos (dimensionless).
C CONSTS(42) is k_t3 in component nucleus_phos (per_minute).
C CONSTS(43) is k_t4 in component nucleus_phos (per_minute).
C ALGBRC(2) is M_cyto in component Mcyto (dimensionless).
C ALGBRC(3) is M_nucl in component Mnucl (dimensionless).
C RATES(3) is d/dt GEFa in component active_GEF (dimensionless).
C RATES(4) is d/dt GAPa in component active_GAP (dimensionless).
C RATES(5) is d/dt RGTP in component RAS_to_GTP (dimensionless).
C RATES(6) is d/dt CYCLa in component adenylate_cyclase (dimensionless).
C RATES(7) is d/dt PDEa in component active_PDE (dimensionless).
C RATES(2) is d/dt cAMP in component cyclic_AMP (micromolar).
C RATES(1) is d/dt R2C2 in component holoenzyme_R_C (dimensionless).
C RATES(10) is d/dt MC in component cytosol (dimensionless).
C RATES(8) is d/dt MN in component nucleus (dimensionless).
C RATES(11) is d/dt MNP in component nucleus_phos (dimensionless).
C RATES(9) is d/dt MCP in component cytosol_phos (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 4
      CONSTS(3) = 1
      CONSTS(4) = 1
      CONSTS(5) = 3.5
      CONSTS(6) = 0.3
      CONSTS(7) = 1.5
      CONSTS(8) = 1.3
      CONSTS(9) = 240
      CONSTS(10) = 250
      CONSTS(11) = 600
      CONSTS(12) = 3.333
      CONSTS(13) = 0.5
      CONSTS(14) = 1.5
      CONSTS(15) = 1
      CONSTS(16) = 1
      STATES(1) = 0.5
      STATES(2) = 1
      CONSTS(17) = 0.05
      CONSTS(18) = 0.05
      STATES(3) = 0.36
      CONSTS(19) = 0.01
      CONSTS(20) = 0.01
      STATES(4) = 0.5
      CONSTS(21) = 0.001
      CONSTS(22) = 0.001
      STATES(5) = 0.1
      CONSTS(23) = 0.01
      CONSTS(24) = 1
      STATES(6) = 0.1
      CONSTS(25) = 0.01
      CONSTS(26) = 0.01
      STATES(7) = 0.5
      CONSTS(27) = 4
      CONSTS(28) = 100
      CONSTS(29) = 0.7
      CONSTS(30) = 20
      CONSTS(31) = 3.333
      CONSTS(32) = 1
      CONSTS(33) = 0.6
      CONSTS(34) = 3.333
      CONSTS(35) = 2
      CONSTS(36) = 10
      CONSTS(37) = 0.001
      CONSTS(38) = 0.05
      CONSTS(39) = 0.05
      STATES(8) = 0.25
      STATES(9) = 0.25
      STATES(10) = 0.25
      CONSTS(40) = 0.05
      CONSTS(41) = 0.05
      STATES(11) = 0.25
      CONSTS(42) = 0.001
      CONSTS(43) = 10
      CONSTS(44) = CONSTS(1)/CONSTS(2)
      CONSTS(45) = ( CONSTS(3)*CONSTS(4))/CONSTS(2)
      CONSTS(46) = ( CONSTS(5)*CONSTS(6))/CONSTS(7)
      CONSTS(47) = CONSTS(8)/CONSTS(7)
      CONSTS(48) = ( CONSTS(9)*CONSTS(2))/CONSTS(10)
      CONSTS(49) = ( CONSTS(11)*CONSTS(7))/CONSTS(10)
      CONSTS(50) = ( CONSTS(12)*CONSTS(6))/CONSTS(13)
      CONSTS(51) = CONSTS(14)/CONSTS(13)
      CONSTS(52) = ( CONSTS(31)*CONSTS(6))/CONSTS(32)
      CONSTS(53) = ( CONSTS(3)*CONSTS(33))/CONSTS(32)
      CONSTS(54) = ( CONSTS(34)*CONSTS(6))/CONSTS(32)
      CONSTS(55) = ( CONSTS(3)*CONSTS(35))/CONSTS(32)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(3) = ( CONSTS(44)*(1.00000 - STATES(3)))/(CONSTS(17)+(1.00000 - STATES(3))) - ( CONSTS(45)*STATES(3))/(CONSTS(18)+STATES(3))
      RATES(5) = ( CONSTS(48)*STATES(3)*(1.00000 - STATES(5)))/(CONSTS(21)+(1.00000 - STATES(5))) - ( CONSTS(49)*STATES(4)*STATES(5))/(CONSTS(22)+STATES(5))
      RATES(6) =  CONSTS(23)*STATES(5)*CONSTS(10)*(1.00000 - STATES(6)) -  CONSTS(24)*STATES(6)
      ALGBRC(4) =  2.00000*(1.00000 - STATES(1))
      RATES(4) = ( CONSTS(46)*ALGBRC(4)*(1.00000 - STATES(4)))/(CONSTS(19)+(1.00000 - STATES(4))) - ( CONSTS(47)*STATES(4))/(CONSTS(20)+STATES(4))
      RATES(7) = ( CONSTS(50)*ALGBRC(4)*(1.00000 - STATES(7)))/(CONSTS(25)+(1.00000 - STATES(7))) - ( CONSTS(51)*STATES(7))/(CONSTS(26)+STATES(7))
      ALGBRC(1) = 1.00000 - STATES(1)
      RATES(1) =  - CONSTS(15)*STATES(1)*STATES(2) ** 2.00000+ CONSTS(16)*ALGBRC(4) ** 2.00000*ALGBRC(1)*CONSTS(6) ** 2.00000*1.00000
      RATES(10) = (( - CONSTS(36)*STATES(10)+ CONSTS(37)*STATES(8)) - ( CONSTS(54)*ALGBRC(4)*STATES(10))/(CONSTS(38)+STATES(10)))+( CONSTS(55)*STATES(9))/(CONSTS(39)+STATES(9))
      RATES(8) = (( CONSTS(36)*STATES(10) -  CONSTS(37)*STATES(8)) - ( CONSTS(52)*ALGBRC(4)*STATES(8))/(CONSTS(40)+STATES(8)))+( CONSTS(53)*STATES(11))/(CONSTS(41)+STATES(11))
      RATES(11) = ((( CONSTS(52)*ALGBRC(4)*STATES(8))/(CONSTS(40)+STATES(8)) - ( CONSTS(53)*STATES(11))/(CONSTS(41)+STATES(11)))+ CONSTS(42)*STATES(9)) -  CONSTS(43)*STATES(11)
      RATES(9) = ( - CONSTS(42)*STATES(9)+ CONSTS(43)*STATES(11)+( CONSTS(54)*ALGBRC(4)*STATES(10))/(CONSTS(38)+STATES(10))) - ( CONSTS(55)*STATES(9))/(CONSTS(39)+STATES(9))
      ALGBRC(5) =  CONSTS(15)*STATES(1)*STATES(2) ** 2.00000 -  CONSTS(16)*ALGBRC(4)*ALGBRC(1)*CONSTS(6) ** 2.00000*1.00000
      RATES(2) = ( CONSTS(27)*STATES(6)*CONSTS(29) - ( CONSTS(28)*CONSTS(13)*STATES(7)*STATES(2))/(CONSTS(30)+STATES(2))) -  2.00000*ALGBRC(5)*CONSTS(6)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) =  2.00000*(1.00000 - STATES(1))
      ALGBRC(1) = 1.00000 - STATES(1)
      ALGBRC(5) =  CONSTS(15)*STATES(1)*STATES(2) ** 2.00000 -  CONSTS(16)*ALGBRC(4)*ALGBRC(1)*CONSTS(6) ** 2.00000*1.00000
      ALGBRC(2) = STATES(10)+STATES(9)
      ALGBRC(3) = STATES(8)+STATES(11)
      RETURN
      END