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 64 entries in the algebraic variable array.
C There are a total of 36 entries in each of the rate and state variable arrays.
C There are a total of 125 entries in the constant variable array.
C
C
C VOI is time in component environment (s).
C STATES(1) is NAn in component NAn (mM).
C ALGBRC(1) is Vn_leak_Na in component Vn_leak_Na (mM_per_s).
C ALGBRC(4) is Vn_pump in component Vn_pump (mM_per_s).
C ALGBRC(40) is Vn_stim in component Vn_stim (mM_per_s).
C STATES(2) is GLCn in component GLCn (mM).
C ALGBRC(5) is V_en_GLC in component V_en_GLC (mM_per_s).
C ALGBRC(6) is Vn_hk in component Vn_hk (mM_per_s).
C STATES(3) is G6Pn in component G6Pn (mM).
C ALGBRC(7) is Vn_pgi in component Vn_pgi (mM_per_s).
C STATES(4) is F6Pn in component F6Pn (mM).
C ALGBRC(8) is Vn_pfk in component Vn_pfk (mM_per_s).
C STATES(5) is GAPn in component GAPn (mM).
C ALGBRC(48) is Vn_pgk in component Vn_pgk (mM_per_s).
C STATES(6) is PEPn in component PEPn (mM).
C ALGBRC(50) is Vn_pk in component Vn_pk (mM_per_s).
C STATES(7) is PYRn in component PYRn (mM).
C ALGBRC(41) is Vn_ldh in component Vn_ldh (mM_per_s).
C ALGBRC(51) is Vn_mito in component Vn_mito (mM_per_s).
C STATES(8) is LACn in component LACn (mM).
C ALGBRC(9) is Vne_LAC in component Vne_LAC (mM_per_s).
C STATES(9) is NADHn in component NADHn (mM).
C STATES(10) is ATPn in component ATPn (mM).
C CONSTS(1) is nOP in component model_parameters (dimensionless).
C ALGBRC(10) is Vn_ATPase in component Vn_ATPase (mM_per_s).
C ALGBRC(52) is Vn_ck in component Vn_ck (mM_per_s).
C ALGBRC(56) is dAMP_dATPn in component dAMP_dATPn (dimensionless).
C STATES(11) is PCrn in component PCrn (mM).
C STATES(12) is O2n in component O2n (mM).
C CONSTS(2) is NAero in component model_parameters (dimensionless).
C ALGBRC(11) is Vcn_O2 in component Vcn_O2 (mM_per_s).
C STATES(13) is GLUn in component GLUn (mM).
C CONSTS(3) is Rng in component model_parameters (dimensionless).
C ALGBRC(29) is Vg_gs in component Vg_gs (mM_per_s).
C ALGBRC(43) is Vn_stim_GLU in component Vn_stim_GLU (mM_per_s).
C STATES(14) is NAg in component NAg (mM).
C ALGBRC(12) is Vg_leak_Na in component Vg_leak_Na (mM_per_s).
C ALGBRC(13) is Vg_pump in component Vg_pump (mM_per_s).
C ALGBRC(30) is Veg_GLU in component Veg_GLU (mM_per_s).
C STATES(15) is GLCg in component GLCg (mM).
C ALGBRC(15) is Vcg_GLC in component Vcg_GLC (mM_per_s).
C ALGBRC(14) is Veg_GLC in component Veg_GLC (mM_per_s).
C ALGBRC(16) is Vg_hk in component Vg_hk (mM_per_s).
C STATES(16) is G6Pg in component G6Pg (mM).
C ALGBRC(17) is Vg_pgi in component Vg_pgi (mM_per_s).
C ALGBRC(19) is Vg_glys in component Vg_glys (mM_per_s).
C ALGBRC(25) is Vg_glyp in component Vg_glyp (mM_per_s).
C STATES(17) is F6Pg in component F6Pg (mM).
C ALGBRC(18) is Vg_pfk in component Vg_pfk (mM_per_s).
C STATES(18) is GAPg in component GAPg (mM).
C ALGBRC(57) is Vg_pgk in component Vg_pgk (mM_per_s).
C STATES(19) is PEPg in component PEPg (mM).
C ALGBRC(59) is Vg_pk in component Vg_pk (mM_per_s).
C STATES(20) is PYRg in component PYRg (mM).
C ALGBRC(44) is Vg_ldh in component Vg_ldh (mM_per_s).
C ALGBRC(60) is Vg_mito in component Vg_mito (mM_per_s).
C STATES(21) is LACg in component LACg (mM).
C ALGBRC(20) is Vge_LAC in component Vge_LAC (mM_per_s).
C ALGBRC(22) is Vgc_LAC in component Vgc_LAC (mM_per_s).
C STATES(22) is NADHg in component NADHg (mM).
C STATES(23) is ATPg in component ATPg (mM).
C ALGBRC(24) is Vg_ATPase in component Vg_ATPase (mM_per_s).
C ALGBRC(61) is Vg_ck in component Vg_ck (mM_per_s).
C ALGBRC(64) is dAMP_dATPg in component dAMP_dATPg (dimensionless).
C STATES(24) is PCrg in component PCrg (mM).
C STATES(25) is O2g in component O2g (mM).
C ALGBRC(26) is Vcg_O2 in component Vcg_O2 (mM_per_s).
C STATES(26) is GLYg in component GLYg (mM).
C STATES(27) is GLUg in component GLUg (mM).
C STATES(28) is GLCe in component GLCe (mM).
C CONSTS(4) is Reg in component model_parameters (dimensionless).
C CONSTS(5) is Ren in component model_parameters (dimensionless).
C ALGBRC(27) is Vce_GLC in component Vce_GLC (mM_per_s).
C STATES(29) is LACe in component LACe (mM).
C ALGBRC(28) is Vec_LAC in component Vec_LAC (mM_per_s).
C STATES(30) is GLUe in component GLUe (mM).
C STATES(31) is O2c in component O2c (mM).
C CONSTS(6) is Rcn in component model_parameters (dimensionless).
C CONSTS(7) is Rcg in component model_parameters (dimensionless).
C ALGBRC(34) is Vc_O2 in component Vc_O2 (mM_per_s).
C STATES(32) is GLCc in component GLCc (mM).
C CONSTS(8) is Rce in component model_parameters (dimensionless).
C ALGBRC(35) is Vc_GLC in component Vc_GLC (mM_per_s).
C STATES(33) is LACc in component LACc (mM).
C ALGBRC(36) is Vc_LAC in component Vc_LAC (mM_per_s).
C STATES(34) is CO2c in component CO2c (mM).
C ALGBRC(53) is Vnc_CO2 in component Vnc_CO2 (mM_per_s).
C ALGBRC(33) is Vc_CO2 in component Vc_CO2 (mM_per_s).
C ALGBRC(62) is Vgc_CO2 in component Vgc_CO2 (mM_per_s).
C STATES(35) is Vv in component Vv (dimensionless).
C ALGBRC(31) is Fin_t in component Fin_t (per_s).
C ALGBRC(37) is Fout_t in component Fout_t (per_s).
C STATES(36) is dHb in component dHb (mM).
C CONSTS(9) is O2a in component model_parameters (mM).
C CONSTS(10) is gn_NA in component Vn_leak_Na (mS_per_cm2).
C CONSTS(11) is Sm_n in component model_parameters (per_cm).
C CONSTS(12) is Vm in component model_parameters (mV).
C CONSTS(13) is Vn in component model_parameters (dimensionless).
C CONSTS(14) is RT in component model_parameters (mV_C_per_mol).
C CONSTS(15) is F in component model_parameters (C_per_mole).
C CONSTS(16) is NAe in component model_parameters (mM).
C CONSTS(17) is kpump in component model_parameters (cm_per_mM_per_s).
C CONSTS(18) is Km_pump in component model_parameters (mM).
C ALGBRC(38) is v_stim in component v_stim (mM_per_s).
C CONSTS(19) is Km_en_GLC in component V_en_GLC (mM).
C CONSTS(20) is Vm_en_GLC in component V_en_GLC (mM_per_s).
C CONSTS(21) is Vmax_n_hk in component Vn_hk (mM_per_s).
C CONSTS(22) is Km_GLC in component model_parameters (mM).
C CONSTS(23) is G6P_inh_hk in component model_parameters (mM).
C CONSTS(24) is aG6P_inh_hk in component model_parameters (dimensionless).
C CONSTS(25) is Vmaxf_n_pgi in component Vn_pgi (mM_per_s).
C CONSTS(26) is Vmaxr_n_pgi in component Vn_pgi (mM_per_s).
C CONSTS(27) is Km_G6P in component model_parameters (mM).
C CONSTS(28) is Km_F6P_pgi in component model_parameters (mM).
C CONSTS(29) is kn_pfk in component Vn_pfk (per_s).
C CONSTS(30) is Km_F6P_pfk in component model_parameters (mM).
C CONSTS(31) is Ki_ATP in component model_parameters (mM).
C CONSTS(32) is nH in component model_parameters (dimensionless).
C CONSTS(33) is kn_pgk in component Vn_pgk (per_mM_per_s).
C ALGBRC(47) is ADPn in component ADPn (mM).
C ALGBRC(39) is NADn in component NADn (mM).
C CONSTS(34) is kn_pk in component Vn_pk (per_mM_per_s).
C CONSTS(35) is kfn_ldh in component Vn_ldh (per_mM_per_s).
C CONSTS(36) is krn_ldh in component Vn_ldh (per_mM_per_s).
C CONSTS(37) is Vmax_n_mito in component Vn_mito (mM_per_s).
C CONSTS(38) is Km_O2 in component model_parameters (mM).
C CONSTS(39) is Km_ADP in component model_parameters (mM).
C CONSTS(40) is Km_PYR in component model_parameters (mM).
C CONSTS(41) is rATP_mito in component model_parameters (dimensionless).
C CONSTS(42) is aATP_mito in component model_parameters (dimensionless).
C CONSTS(43) is Vmax_ne_LAC in component Vne_LAC (mM_per_s).
C CONSTS(44) is Km_ne_LAC in component Vne_LAC (mM).
C CONSTS(45) is Vmax_n_ATPase in component Vn_ATPase (mM_per_s).
C CONSTS(46) is krn_ck in component Vn_ck (per_mM_per_s).
C CONSTS(47) is kfn_ck in component Vn_ck (per_mM_per_s).
C ALGBRC(45) is CRn in component CRn (mM).
C CONSTS(48) is nh_O2 in component Vcn_O2 (dimensionless).
C CONSTS(49) is PScapn in component Vcn_O2 (per_s).
C CONSTS(50) is Ko2 in component model_parameters (mM).
C CONSTS(51) is HbOP in component model_parameters (mM).
C CONSTS(52) is gg_NA in component Vg_leak_Na (mS_per_cm2).
C CONSTS(53) is Sm_g in component model_parameters (per_cm).
C CONSTS(54) is Vg in component model_parameters (dimensionless).
C CONSTS(55) is Km_eg_GLC in component Veg_GLC (mM).
C CONSTS(56) is Vm_eg_GLC in component Veg_GLC (mM_per_s).
C CONSTS(57) is KO1 in component model_parameters (dimensionless).
C CONSTS(58) is Km_cg_GLC in component Vcg_GLC (mM).
C CONSTS(59) is Vm_cg_GLC in component Vcg_GLC (mM_per_s).
C CONSTS(60) is Vmax_g_hk in component Vg_hk (mM_per_s).
C CONSTS(61) is Vmaxf_g_pgi in component Vg_pgi (mM_per_s).
C CONSTS(62) is Vmaxr_g_pgi in component Vg_pgi (mM_per_s).
C CONSTS(63) is kg_pfk in component Vg_pfk (per_s).
C CONSTS(64) is kg_pgk in component Vg_pgk (per_mM_per_s).
C ALGBRC(55) is ADPg in component ADPg (mM).
C ALGBRC(42) is NADg in component NADg (mM).
C CONSTS(65) is kg_pk in component Vg_pk (per_mM_per_s).
C CONSTS(66) is kfg_ldh in component Vg_ldh (per_mM_per_s).
C CONSTS(67) is krg_ldh in component Vg_ldh (per_mM_per_s).
C CONSTS(68) is Vmax_g_mito in component Vg_mito (mM_per_s).
C CONSTS(69) is Vmax_ge_LAC in component Vge_LAC (mM_per_s).
C CONSTS(70) is Km_ge_LAC in component Vge_LAC (mM).
C CONSTS(71) is Vmax_gc_LAC in component Vgc_LAC (mM_per_s).
C CONSTS(72) is Km_gc_LAC in component Vgc_LAC (mM).
C CONSTS(73) is Vmax_g_ATPase in component Vg_ATPase (mM_per_s).
C CONSTS(74) is krg_ck in component Vg_ck (per_mM_per_s).
C CONSTS(75) is kfg_ck in component Vg_ck (per_mM_per_s).
C ALGBRC(46) is CRg in component CRg (mM).
C CONSTS(76) is PScapg in component Vcg_O2 (per_s).
C CONSTS(77) is nh_O2 in component model_parameters (dimensionless).
C CONSTS(78) is Vc in component model_parameters (dimensionless).
C CONSTS(79) is GLCa in component model_parameters (mM).
C CONSTS(80) is Km_ce_GLC in component Vce_GLC (mM).
C CONSTS(81) is Vm_ce_GLC in component Vce_GLC (mM_per_s).
C CONSTS(82) is LACa in component model_parameters (mM).
C CONSTS(83) is Km_ec_LAC in component Vec_LAC (mM).
C CONSTS(84) is Vm_ec_LAC in component Vec_LAC (mM_per_s).
C CONSTS(85) is R_GLU_NA in component model_parameters (dimensionless).
C CONSTS(86) is Km_GLU in component model_parameters (mM).
C CONSTS(87) is KO2 in component model_parameters (dimensionless).
C CONSTS(88) is Vmax_g_gs in component Vg_gs (mM_per_s).
C CONSTS(89) is Km_ATP in component model_parameters (mM).
C CONSTS(90) is Vmax_eg_GLU in component Veg_GLU (mM_per_s).
C CONSTS(91) is CO2a in component model_parameters (mM).
C CONSTS(92) is Vmax_glys in component Vg_glys (mM_per_s).
C CONSTS(93) is Km_G6P_glys in component Vg_glys (mM).
C CONSTS(94) is GLY_inh in component model_parameters (mM).
C CONSTS(95) is aGLY_inh in component model_parameters (dimensionless).
C CONSTS(96) is Vmax_glyp in component Vg_glyp (mM_per_s).
C CONSTS(97) is Km_GLY in component Vg_glyp (mM).
C ALGBRC(23) is deltaVt_GLY in component Vg_glyp (dimensionless).
C ALGBRC(21) is unitstepSB2 in component unitstepSB2 (dimensionless).
C CONSTS(98) is stim in component model_parameters (dimensionless).
C CONSTS(99) is to in component model_parameters (s).
C CONSTS(100) is to_GLY in component model_parameters (s).
C CONSTS(101) is tend_GLY in component model_parameters (s).
C CONSTS(102) is sr_GLY in component model_parameters (dimensionless).
C CONSTS(103) is t1 in component model_parameters (s).
C CONSTS(104) is delta_GLY in component model_parameters (dimensionless).
C CONSTS(105) is KO3 in component model_parameters (dimensionless).
C CONSTS(106) is CBF0 in component Fin_t (per_s).
C CONSTS(107) is tend in component model_parameters (s).
C CONSTS(108) is sr in component model_parameters (dimensionless).
C CONSTS(109) is deltaf in component model_parameters (dimensionless).
C CONSTS(110) is CBF0 in component model_parameters (per_s).
C CONSTS(111) is Vv0 in component model_parameters (dimensionless).
C CONSTS(112) is tv in component model_parameters (s).
C CONSTS(113) is NADH_n_tot in component NADn (mM).
C CONSTS(114) is NADH_g_tot in component NADg (mM).
C CONSTS(115) is PCrn_tot in component CRn (mM).
C CONSTS(116) is PCrg_tot in component CRg (mM).
C CONSTS(117) is ATPtot in component model_parameters (mM).
C CONSTS(118) is qak in component model_parameters (dimensionless).
C ALGBRC(54) is u_n in component u_n (dimensionless).
C ALGBRC(63) is u_g in component u_g (dimensionless).
C ALGBRC(49) is AMPn in component AMPn (mM).
C ALGBRC(58) is AMPg in component AMPg (mM).
C ALGBRC(2) is BOLD in component BOLD (dimensionless).
C CONSTS(119) is k1 in component model_parameters (dimensionless).
C CONSTS(120) is k2 in component model_parameters (dimensionless).
C CONSTS(121) is k3 in component model_parameters (dimensionless).
C CONSTS(122) is dHb0 in component model_parameters (mM).
C ALGBRC(32) is unitpulseSB in component v_stim (dimensionless).
C CONSTS(123) is t_n_stim in component model_parameters (s).
C CONSTS(124) is v1_n in component model_parameters (mM_per_s).
C CONSTS(125) is v2_n in component model_parameters (mM_per_s).
C ALGBRC(3) is unitstepSB in component unitstepSB (dimensionless).
C RATES(1) is d/dt NAn in component NAn (mM).
C RATES(2) is d/dt GLCn in component GLCn (mM).
C RATES(3) is d/dt G6Pn in component G6Pn (mM).
C RATES(4) is d/dt F6Pn in component F6Pn (mM).
C RATES(5) is d/dt GAPn in component GAPn (mM).
C RATES(6) is d/dt PEPn in component PEPn (mM).
C RATES(7) is d/dt PYRn in component PYRn (mM).
C RATES(8) is d/dt LACn in component LACn (mM).
C RATES(9) is d/dt NADHn in component NADHn (mM).
C RATES(10) is d/dt ATPn in component ATPn (mM).
C RATES(11) is d/dt PCrn in component PCrn (mM).
C RATES(12) is d/dt O2n in component O2n (mM).
C RATES(13) is d/dt GLUn in component GLUn (mM).
C RATES(14) is d/dt NAg in component NAg (mM).
C RATES(15) is d/dt GLCg in component GLCg (mM).
C RATES(16) is d/dt G6Pg in component G6Pg (mM).
C RATES(17) is d/dt F6Pg in component F6Pg (mM).
C RATES(18) is d/dt GAPg in component GAPg (mM).
C RATES(19) is d/dt PEPg in component PEPg (mM).
C RATES(20) is d/dt PYRg in component PYRg (mM).
C RATES(21) is d/dt LACg in component LACg (mM).
C RATES(22) is d/dt NADHg in component NADHg (mM).
C RATES(23) is d/dt ATPg in component ATPg (mM).
C RATES(24) is d/dt PCrg in component PCrg (mM).
C RATES(25) is d/dt O2g in component O2g (mM).
C RATES(26) is d/dt GLYg in component GLYg (mM).
C RATES(27) is d/dt GLUg in component GLUg (mM).
C RATES(28) is d/dt GLCe in component GLCe (mM).
C RATES(29) is d/dt LACe in component LACe (mM).
C RATES(30) is d/dt GLUe in component GLUe (mM).
C RATES(31) is d/dt O2c in component O2c (mM).
C RATES(32) is d/dt GLCc in component GLCc (mM).
C RATES(33) is d/dt LACc in component LACc (mM).
C RATES(34) is d/dt CO2c in component CO2c (mM).
C RATES(35) is d/dt Vv in component Vv (dimensionless).
C RATES(36) is d/dt dHb in component dHb (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 15.533
      STATES(2) = 0.2633
      STATES(3) = 0.7275
      STATES(4) = 0.1091
      STATES(5) = 0.0418
      STATES(6) = 0.0037
      STATES(7) = 0.0388
      STATES(8) = 0.3856
      STATES(9) = 0.0319
      STATES(10) = 2.2592
      CONSTS(1) = 15.0
      STATES(11) = 4.2529
      STATES(12) = 0.0975
      CONSTS(2) = 3.0
      STATES(13) = 3.0
      CONSTS(3) = 1.8
      STATES(14) = 13.36
      STATES(15) = 0.1656
      STATES(16) = 0.7326
      STATES(17) = 0.1116
      STATES(18) = 0.0698
      STATES(19) = 0.0254
      STATES(20) = 0.1711
      STATES(21) = 0.4651
      STATES(22) = 0.0445
      STATES(23) = 2.24
      STATES(24) = 4.6817
      STATES(25) = 0.1589
      STATES(26) = 2.5
      STATES(27) = 0.0
      STATES(28) = 0.3339
      CONSTS(4) = 0.8
      CONSTS(5) = 0.4444444444444444
      STATES(29) = 0.3986
      STATES(30) = 0.0
      STATES(31) = 7.4201
      CONSTS(6) = 0.01222
      CONSTS(7) = 0.022
      STATES(32) = 4.6401
      CONSTS(8) = 0.0275
      STATES(33) = 0.3251
      STATES(34) = 2.12
      STATES(35) = 0.0237
      STATES(36) = 0.0218
      CONSTS(9) = 8.34
      CONSTS(10) = 0.0039
      CONSTS(11) = 40500
      CONSTS(12) = -70
      CONSTS(13) = 0.45
      CONSTS(14) = 2577340
      CONSTS(15) = 96500
      CONSTS(16) = 150.0
      CONSTS(17) = 3.17e-7
      CONSTS(18) = 0.4243
      CONSTS(19) = 5.32
      CONSTS(20) = 0.50417
      CONSTS(21) = 0.0513
      CONSTS(22) = 0.105
      CONSTS(23) = 0.6
      CONSTS(24) = 20.0
      CONSTS(25) = 0.5
      CONSTS(26) = 0.45
      CONSTS(27) = 0.5
      CONSTS(28) = 0.06
      CONSTS(29) = 0.55783
      CONSTS(30) = 0.18
      CONSTS(31) = 0.7595
      CONSTS(32) = 4.0
      CONSTS(33) = 0.4287
      CONSTS(34) = 28.6
      CONSTS(35) = 5.30
      CONSTS(36) = 0.1046
      CONSTS(37) = 0.05557
      CONSTS(38) = 0.0029658
      CONSTS(39) = 0.00107
      CONSTS(40) = 0.0632
      CONSTS(41) = 20.0
      CONSTS(42) = 5.0
      CONSTS(43) = 0.1978
      CONSTS(44) = 0.09314
      CONSTS(45) = 0.04889
      CONSTS(46) = 0.015
      CONSTS(47) = 0.0524681
      CONSTS(48) = 2.7
      CONSTS(49) = 0.2202
      CONSTS(50) = 0.089733
      CONSTS(51) = 8.6
      CONSTS(52) = 0.00325
      CONSTS(53) = 10500
      CONSTS(54) = 0.25
      CONSTS(55) = 3.53
      CONSTS(56) = 0.038089
      CONSTS(57) = 1.0
      CONSTS(58) = 9.92
      CONSTS(59) = 0.0098394
      CONSTS(60) = 0.050461
      CONSTS(61) = 0.5
      CONSTS(62) = 0.45
      CONSTS(63) = 0.403
      CONSTS(64) = 0.2514
      CONSTS(65) = 2.73
      CONSTS(66) = 6.2613
      CONSTS(67) = 0.54682
      CONSTS(68) = 0.008454
      CONSTS(69) = 0.086124
      CONSTS(70) = 0.22163
      CONSTS(71) = 0.00021856
      CONSTS(72) = 0.12862
      CONSTS(73) = 0.035657
      CONSTS(74) = 0.02073
      CONSTS(75) = 0.0243
      CONSTS(76) = 0.2457
      CONSTS(77) = 2.7
      CONSTS(78) = 0.0055
      CONSTS(79) = 4.8
      CONSTS(80) = 8.4568
      CONSTS(81) = 0.0489
      CONSTS(82) = 0.313
      CONSTS(83) = 0.764818
      CONSTS(84) = 0.0325
      CONSTS(85) = 0.075
      CONSTS(86) = 0.05
      CONSTS(87) = 1
      CONSTS(88) = 0.3
      CONSTS(89) = 0.01532
      CONSTS(90) = 0.0208
      CONSTS(91) = 1.2
      CONSTS(92) = 0.0001528
      CONSTS(93) = 0.5
      CONSTS(94) = 4.2
      CONSTS(95) = 20.0
      CONSTS(96) = 4.922e-5
      CONSTS(97) = 1.0
      CONSTS(98) = 1
      CONSTS(99) = 200
      CONSTS(100) = 83
      CONSTS(101) = 440
      CONSTS(102) = 4
      CONSTS(103) = 2
      CONSTS(104) = 62
      CONSTS(105) = 1
      CONSTS(106) = 0.012
      CONSTS(107) = 300
      CONSTS(108) = 4.59186
      CONSTS(109) = 0.42
      CONSTS(110) = 0.012
      CONSTS(111) = 0.0236
      CONSTS(112) = 35.0
      CONSTS(113) = 0.22
      CONSTS(114) = 0.22
      CONSTS(115) = 5.0
      CONSTS(116) = 5.0
      CONSTS(117) = 2.379
      CONSTS(118) = 0.92
      CONSTS(119) = 2.22
      CONSTS(120) = 0.46
      CONSTS(121) = 0.43
      CONSTS(122) = 0.064
      CONSTS(123) = 2
      CONSTS(124) = 0.041
      CONSTS(125) = 2.55
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) =  CONSTS(20)*(STATES(28)/(STATES(28)+CONSTS(19)) - STATES(2)/(STATES(2)+CONSTS(19)))
      ALGBRC(6) =  CONSTS(21)*STATES(10)*(STATES(2)/(STATES(2)+CONSTS(22)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(24)*( 1.00000*(STATES(3) - CONSTS(23))))))
      RATES(2) = ALGBRC(5) - ALGBRC(6)
      ALGBRC(7) =  CONSTS(25)*(STATES(3)/(STATES(3)+CONSTS(27))) -  CONSTS(26)*(STATES(4)/(STATES(4)+CONSTS(28)))
      RATES(3) = ALGBRC(6) - ALGBRC(7)
      ALGBRC(8) =  CONSTS(29)*STATES(10)*(STATES(4)/(STATES(4)+CONSTS(30)))*1.00000+STATES(10)/CONSTS(31) ** CONSTS(32) ** -1.00000
      RATES(4) = ALGBRC(7) - ALGBRC(8)
      ALGBRC(15) =  CONSTS(59)*(STATES(32)/(STATES(32)+CONSTS(58)) - STATES(15)/(STATES(15)+CONSTS(58)))
      ALGBRC(14) =  CONSTS(57)*CONSTS(56)*(STATES(28)/(STATES(28)+CONSTS(55)) - STATES(15)/(STATES(15)+CONSTS(55)))
      ALGBRC(16) =  CONSTS(60)*STATES(23)*(STATES(15)/(STATES(15)+CONSTS(22)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(24)*( 1.00000*(STATES(16) - CONSTS(23))))))
      RATES(15) = (ALGBRC(15)+ALGBRC(14)) - ALGBRC(16)
      ALGBRC(17) =  CONSTS(61)*(STATES(16)/(STATES(16)+CONSTS(27))) -  CONSTS(62)*(STATES(17)/(STATES(17)+CONSTS(28)))
      ALGBRC(18) =  CONSTS(63)*STATES(23)*(STATES(17)/(STATES(17)+CONSTS(30)))*1.00000+STATES(23)/CONSTS(31) ** CONSTS(32) ** -1.00000
      RATES(17) = ALGBRC(17) - ALGBRC(18)
      ALGBRC(19) =  CONSTS(92)*(STATES(16)/(STATES(16)+CONSTS(93)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(95)*( 1.00000*(STATES(26) - CONSTS(94))))))
      ALGBRC(21) = TERNRY(VOI - (CONSTS(101)+CONSTS(99)+CONSTS(100)).GE.0.00000, 1.00000, 0.00000)
      ALGBRC(23) = 1.00000+ CONSTS(98)*( CONSTS(104)*CONSTS(105)*(1.00000/(1.00000+EXP( 1.00000*- CONSTS(102)*(VOI - (CONSTS(99)+CONSTS(100))))))*(1.00000 - ALGBRC(21)))
      ALGBRC(25) =  CONSTS(96)*(STATES(26)/(STATES(26)+CONSTS(97)))*ALGBRC(23)
      RATES(16) = (ALGBRC(16)+ALGBRC(25)) - (ALGBRC(17)+ALGBRC(19))
      RATES(26) = ALGBRC(19) - ALGBRC(25)
      ALGBRC(27) =  CONSTS(81)*(STATES(32)/(STATES(32)+CONSTS(80)) - STATES(28)/(STATES(28)+CONSTS(80)))
      RATES(28) = ALGBRC(27) - ( ALGBRC(14)*(1.00000/CONSTS(4))+ ALGBRC(5)*(1.00000/CONSTS(5)))
      ALGBRC(9) =  CONSTS(43)*(STATES(8)/(STATES(8)+CONSTS(44)) - STATES(29)/(STATES(29)+CONSTS(44)))
      ALGBRC(20) =  CONSTS(69)*(STATES(21)/(STATES(21)+CONSTS(70)) - STATES(29)/(STATES(29)+CONSTS(70)))
      ALGBRC(28) =  CONSTS(84)*(STATES(29)/(STATES(29)+CONSTS(83)) - STATES(33)/(STATES(33)+CONSTS(83)))
      RATES(29) = ( ALGBRC(9)*(1.00000/CONSTS(5))+ ALGBRC(20)*(1.00000/CONSTS(4))) - ALGBRC(28)
      ALGBRC(12) =  (CONSTS(53)/CONSTS(54))*(CONSTS(52)/CONSTS(15))*( (CONSTS(14)/CONSTS(15))*log(CONSTS(16)/STATES(14)) - CONSTS(12))
      ALGBRC(13) =  (CONSTS(53)/CONSTS(54))*CONSTS(17)*STATES(23)*STATES(14)*1.00000+STATES(23)/CONSTS(18) ** -1.00000
      ALGBRC(30) =  CONSTS(90)*(STATES(30)/(STATES(30)+CONSTS(86)))
      RATES(14) = (ALGBRC(12)+ 3.00000*ALGBRC(30)) -  3.00000*ALGBRC(13)
      ALGBRC(29) =  CONSTS(88)*( (STATES(27)/(STATES(27)+CONSTS(86)))*(STATES(23)/(STATES(23)+CONSTS(89))))
      RATES(27) = ALGBRC(30) - ALGBRC(29)
      ALGBRC(11) =  (CONSTS(49)/CONSTS(13))*( CONSTS(50)*CONSTS(51)/STATES(31) - 1.00000 ** -1.00000/CONSTS(48) - STATES(12))
      ALGBRC(26) =  (CONSTS(76)/CONSTS(54))*( CONSTS(50)*CONSTS(51)/STATES(31) - 1.00000 ** -1.00000/CONSTS(77) - STATES(25))
      ALGBRC(31) = CONSTS(106)+( CONSTS(98)*CONSTS(106)*CONSTS(109)*(1.00000/(1.00000+EXP( ( 1.00000*- CONSTS(108))*(VOI - ((CONSTS(99)+CONSTS(103)) - 3.00000))))) -  CONSTS(98)*CONSTS(106)*CONSTS(109)*(1.00000/(1.00000+EXP( ( 1.00000*- CONSTS(108))*(VOI - (CONSTS(99)+CONSTS(107)+CONSTS(103)+3.00000))))))
      ALGBRC(34) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(9) - STATES(31))
      RATES(31) = ALGBRC(34) - ( ALGBRC(11)*(1.00000/CONSTS(6))+ ALGBRC(26)*(1.00000/CONSTS(7)))
      ALGBRC(35) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(79) - STATES(32))
      RATES(32) = ALGBRC(35) - ( ALGBRC(27)*(1.00000/CONSTS(8))+ ALGBRC(15)*(1.00000/CONSTS(7)))
      ALGBRC(22) =  CONSTS(71)*(STATES(21)/(STATES(21)+CONSTS(72)) - STATES(33)/(STATES(33)+CONSTS(72)))
      ALGBRC(36) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(82) - STATES(33))
      RATES(33) = ALGBRC(36)+( ALGBRC(28)*(1.00000/CONSTS(8))+ ALGBRC(22)*(1.00000/CONSTS(7)))
      ALGBRC(37) =  CONSTS(110)*((STATES(35)/CONSTS(111) ** 2.00000+ CONSTS(112)*STATES(35)/CONSTS(111) ** -0.500000*(ALGBRC(31)/CONSTS(111)))/(1.00000+ CONSTS(110)*CONSTS(112)*STATES(35)/CONSTS(111) ** -0.500000*(1.00000/CONSTS(111))))
      RATES(35) = ALGBRC(31) - ALGBRC(37)
      RATES(36) =  ALGBRC(31)*(CONSTS(9) - STATES(31)) -  ALGBRC(37)*(STATES(36)/STATES(35))
      ALGBRC(1) =  (CONSTS(11)/CONSTS(13))*(CONSTS(10)/CONSTS(15))*( (CONSTS(14)/CONSTS(15))*log(CONSTS(16)/STATES(1)) - CONSTS(12))
      ALGBRC(4) =  (CONSTS(11)/CONSTS(13))*CONSTS(17)*STATES(10)*STATES(1)*1.00000+STATES(10)/CONSTS(18) ** -1.00000
      ALGBRC(32) = TERNRY(VOI.GE.CONSTS(99).AND.VOI.LE.CONSTS(99)+CONSTS(107), 1.00000, 0.00000)
      ALGBRC(38) =  CONSTS(98)*(CONSTS(124)+ CONSTS(125)*((VOI - CONSTS(99))/CONSTS(123))*EXP(- ( (VOI - CONSTS(99))*(ALGBRC(32)/CONSTS(123)))))*ALGBRC(32)
      ALGBRC(40) = ALGBRC(38)
      RATES(1) = (ALGBRC(1)+ALGBRC(40)) -  3.00000*ALGBRC(4)
      ALGBRC(39) = CONSTS(113) - STATES(9)
      ALGBRC(41) =  CONSTS(35)*STATES(7)*STATES(9) -  CONSTS(36)*STATES(8)*ALGBRC(39)
      RATES(8) = ALGBRC(41) - ALGBRC(9)
      ALGBRC(43) =  ALGBRC(40)*CONSTS(85)*CONSTS(87)*(STATES(13)/(STATES(13)+CONSTS(86)))
      RATES(13) =  ALGBRC(29)*(1.00000/CONSTS(3)) - ALGBRC(43)
      RATES(30) =  ALGBRC(43)*(1.00000/CONSTS(5)) -  ALGBRC(30)*(1.00000/CONSTS(4))
      ALGBRC(42) = CONSTS(114) - STATES(22)
      ALGBRC(44) =  CONSTS(66)*STATES(20)*STATES(22) -  CONSTS(67)*STATES(21)*ALGBRC(42)
      RATES(21) = ALGBRC(44) - (ALGBRC(20)+ALGBRC(22))
      ALGBRC(47) =  (STATES(10)/2.00000)*(- CONSTS(118)+ (CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(10) - 1.00000)) ** (1.0 / 2))
      ALGBRC(48) =  CONSTS(33)*STATES(5)*ALGBRC(47)*(ALGBRC(39)/STATES(9))
      RATES(5) =  2.00000*ALGBRC(8) - ALGBRC(48)
      ALGBRC(50) =  CONSTS(34)*STATES(6)*ALGBRC(47)
      RATES(6) = ALGBRC(48) - ALGBRC(50)
      ALGBRC(51) =  CONSTS(37)*(STATES(12)/(STATES(12)+CONSTS(38)))*(ALGBRC(47)/(ALGBRC(47)+CONSTS(39)))*(STATES(7)/(STATES(7)+CONSTS(40)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(42)*( 1.00000*(STATES(10)/ALGBRC(47) -  1.00000*CONSTS(41))))))
      RATES(7) = ALGBRC(50) - (ALGBRC(41)+ALGBRC(51))
      RATES(9) = ALGBRC(48) - (ALGBRC(41)+ALGBRC(51))
      RATES(12) = ALGBRC(11) -  CONSTS(2)*ALGBRC(51)
      ALGBRC(45) = CONSTS(115) - STATES(11)
      ALGBRC(52) =  CONSTS(47)*STATES(11)*ALGBRC(47) -  CONSTS(46)*ALGBRC(45)*STATES(10)
      RATES(11) = - ALGBRC(52)
      ALGBRC(10) =  CONSTS(45)*(STATES(10)/(STATES(10)+0.00100000))
      ALGBRC(54) = CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(10) - 1.00000)
      ALGBRC(56) = (CONSTS(118)/2.00000+ CONSTS(118)*(CONSTS(117)/( STATES(10)* ALGBRC(54) ** (1.0 / 2)))) - (1.00000+ 0.500000* ALGBRC(54) ** (1.0 / 2))
      RATES(10) =  ((ALGBRC(48)+ALGBRC(50)+ CONSTS(1)*ALGBRC(51)+ALGBRC(52)) - (ALGBRC(6)+ALGBRC(8)+ALGBRC(10)+ALGBRC(4)))*1.00000 - ALGBRC(56) ** -1.00000
      ALGBRC(55) =  (STATES(23)/2.00000)*(- CONSTS(118)+ (CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(23) - 1.00000)) ** (1.0 / 2))
      ALGBRC(57) =  CONSTS(64)*STATES(18)*ALGBRC(55)*(ALGBRC(42)/STATES(22))
      RATES(18) =  2.00000*ALGBRC(18) - ALGBRC(57)
      ALGBRC(59) =  CONSTS(65)*STATES(19)*ALGBRC(55)
      RATES(19) = ALGBRC(57) - ALGBRC(59)
      ALGBRC(60) =  CONSTS(68)*(STATES(25)/(STATES(25)+CONSTS(38)))*(ALGBRC(55)/(ALGBRC(55)+CONSTS(39)))*(STATES(20)/(STATES(20)+CONSTS(40)))*(1.00000 - 1.00000/(1.00000+EXP( 1.00000*- CONSTS(42)*(STATES(23)/ALGBRC(55) -  1.00000*CONSTS(41)))))
      RATES(20) = ALGBRC(59) - (ALGBRC(44)+ALGBRC(60))
      RATES(22) = ALGBRC(57) - (ALGBRC(44)+ALGBRC(60))
      RATES(25) = ALGBRC(26) -  CONSTS(2)*ALGBRC(60)
      ALGBRC(46) = CONSTS(116) - STATES(24)
      ALGBRC(61) =  CONSTS(75)*STATES(24)*ALGBRC(55) -  CONSTS(74)*ALGBRC(46)*STATES(23)
      RATES(24) = - ALGBRC(61)
      ALGBRC(53) =  3.00000*ALGBRC(51)
      ALGBRC(33) =  2.00000*(ALGBRC(31)/CONSTS(78))*(STATES(34) - CONSTS(91))
      ALGBRC(62) =  3.00000*ALGBRC(60)
      RATES(34) = ( ALGBRC(53)*(1.00000/CONSTS(6))+ ALGBRC(62)*(1.00000/CONSTS(7))) - ALGBRC(33)
      ALGBRC(24) =  CONSTS(73)*(STATES(23)/(STATES(23)+0.00100000))
      ALGBRC(63) = CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(23) - 1.00000)
      ALGBRC(64) = (CONSTS(118)/2.00000+ CONSTS(118)*(CONSTS(117)/( STATES(23)* ALGBRC(63) ** (1.0 / 2)))) - (1.00000+ 0.500000* ALGBRC(63) ** (1.0 / 2))
      RATES(23) =  ((ALGBRC(57)+ALGBRC(59)+ CONSTS(1)*ALGBRC(60)+ALGBRC(61)) - (ALGBRC(16)+ALGBRC(18)+ALGBRC(24)+ALGBRC(13)+ALGBRC(29)))*1.00000 - ALGBRC(64) ** -1.00000
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) =  CONSTS(20)*(STATES(28)/(STATES(28)+CONSTS(19)) - STATES(2)/(STATES(2)+CONSTS(19)))
      ALGBRC(6) =  CONSTS(21)*STATES(10)*(STATES(2)/(STATES(2)+CONSTS(22)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(24)*( 1.00000*(STATES(3) - CONSTS(23))))))
      ALGBRC(7) =  CONSTS(25)*(STATES(3)/(STATES(3)+CONSTS(27))) -  CONSTS(26)*(STATES(4)/(STATES(4)+CONSTS(28)))
      ALGBRC(8) =  CONSTS(29)*STATES(10)*(STATES(4)/(STATES(4)+CONSTS(30)))*1.00000+STATES(10)/CONSTS(31) ** CONSTS(32) ** -1.00000
      ALGBRC(15) =  CONSTS(59)*(STATES(32)/(STATES(32)+CONSTS(58)) - STATES(15)/(STATES(15)+CONSTS(58)))
      ALGBRC(14) =  CONSTS(57)*CONSTS(56)*(STATES(28)/(STATES(28)+CONSTS(55)) - STATES(15)/(STATES(15)+CONSTS(55)))
      ALGBRC(16) =  CONSTS(60)*STATES(23)*(STATES(15)/(STATES(15)+CONSTS(22)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(24)*( 1.00000*(STATES(16) - CONSTS(23))))))
      ALGBRC(17) =  CONSTS(61)*(STATES(16)/(STATES(16)+CONSTS(27))) -  CONSTS(62)*(STATES(17)/(STATES(17)+CONSTS(28)))
      ALGBRC(18) =  CONSTS(63)*STATES(23)*(STATES(17)/(STATES(17)+CONSTS(30)))*1.00000+STATES(23)/CONSTS(31) ** CONSTS(32) ** -1.00000
      ALGBRC(19) =  CONSTS(92)*(STATES(16)/(STATES(16)+CONSTS(93)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(95)*( 1.00000*(STATES(26) - CONSTS(94))))))
      ALGBRC(21) = TERNRY(VOI - (CONSTS(101)+CONSTS(99)+CONSTS(100)).GE.0.00000, 1.00000, 0.00000)
      ALGBRC(23) = 1.00000+ CONSTS(98)*( CONSTS(104)*CONSTS(105)*(1.00000/(1.00000+EXP( 1.00000*- CONSTS(102)*(VOI - (CONSTS(99)+CONSTS(100))))))*(1.00000 - ALGBRC(21)))
      ALGBRC(25) =  CONSTS(96)*(STATES(26)/(STATES(26)+CONSTS(97)))*ALGBRC(23)
      ALGBRC(27) =  CONSTS(81)*(STATES(32)/(STATES(32)+CONSTS(80)) - STATES(28)/(STATES(28)+CONSTS(80)))
      ALGBRC(9) =  CONSTS(43)*(STATES(8)/(STATES(8)+CONSTS(44)) - STATES(29)/(STATES(29)+CONSTS(44)))
      ALGBRC(20) =  CONSTS(69)*(STATES(21)/(STATES(21)+CONSTS(70)) - STATES(29)/(STATES(29)+CONSTS(70)))
      ALGBRC(28) =  CONSTS(84)*(STATES(29)/(STATES(29)+CONSTS(83)) - STATES(33)/(STATES(33)+CONSTS(83)))
      ALGBRC(12) =  (CONSTS(53)/CONSTS(54))*(CONSTS(52)/CONSTS(15))*( (CONSTS(14)/CONSTS(15))*log(CONSTS(16)/STATES(14)) - CONSTS(12))
      ALGBRC(13) =  (CONSTS(53)/CONSTS(54))*CONSTS(17)*STATES(23)*STATES(14)*1.00000+STATES(23)/CONSTS(18) ** -1.00000
      ALGBRC(30) =  CONSTS(90)*(STATES(30)/(STATES(30)+CONSTS(86)))
      ALGBRC(29) =  CONSTS(88)*( (STATES(27)/(STATES(27)+CONSTS(86)))*(STATES(23)/(STATES(23)+CONSTS(89))))
      ALGBRC(11) =  (CONSTS(49)/CONSTS(13))*( CONSTS(50)*CONSTS(51)/STATES(31) - 1.00000 ** -1.00000/CONSTS(48) - STATES(12))
      ALGBRC(26) =  (CONSTS(76)/CONSTS(54))*( CONSTS(50)*CONSTS(51)/STATES(31) - 1.00000 ** -1.00000/CONSTS(77) - STATES(25))
      ALGBRC(31) = CONSTS(106)+( CONSTS(98)*CONSTS(106)*CONSTS(109)*(1.00000/(1.00000+EXP( ( 1.00000*- CONSTS(108))*(VOI - ((CONSTS(99)+CONSTS(103)) - 3.00000))))) -  CONSTS(98)*CONSTS(106)*CONSTS(109)*(1.00000/(1.00000+EXP( ( 1.00000*- CONSTS(108))*(VOI - (CONSTS(99)+CONSTS(107)+CONSTS(103)+3.00000))))))
      ALGBRC(34) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(9) - STATES(31))
      ALGBRC(35) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(79) - STATES(32))
      ALGBRC(22) =  CONSTS(71)*(STATES(21)/(STATES(21)+CONSTS(72)) - STATES(33)/(STATES(33)+CONSTS(72)))
      ALGBRC(36) =  2.00000*(ALGBRC(31)/CONSTS(78))*(CONSTS(82) - STATES(33))
      ALGBRC(37) =  CONSTS(110)*((STATES(35)/CONSTS(111) ** 2.00000+ CONSTS(112)*STATES(35)/CONSTS(111) ** -0.500000*(ALGBRC(31)/CONSTS(111)))/(1.00000+ CONSTS(110)*CONSTS(112)*STATES(35)/CONSTS(111) ** -0.500000*(1.00000/CONSTS(111))))
      ALGBRC(1) =  (CONSTS(11)/CONSTS(13))*(CONSTS(10)/CONSTS(15))*( (CONSTS(14)/CONSTS(15))*log(CONSTS(16)/STATES(1)) - CONSTS(12))
      ALGBRC(4) =  (CONSTS(11)/CONSTS(13))*CONSTS(17)*STATES(10)*STATES(1)*1.00000+STATES(10)/CONSTS(18) ** -1.00000
      ALGBRC(32) = TERNRY(VOI.GE.CONSTS(99).AND.VOI.LE.CONSTS(99)+CONSTS(107), 1.00000, 0.00000)
      ALGBRC(38) =  CONSTS(98)*(CONSTS(124)+ CONSTS(125)*((VOI - CONSTS(99))/CONSTS(123))*EXP(- ( (VOI - CONSTS(99))*(ALGBRC(32)/CONSTS(123)))))*ALGBRC(32)
      ALGBRC(40) = ALGBRC(38)
      ALGBRC(39) = CONSTS(113) - STATES(9)
      ALGBRC(41) =  CONSTS(35)*STATES(7)*STATES(9) -  CONSTS(36)*STATES(8)*ALGBRC(39)
      ALGBRC(43) =  ALGBRC(40)*CONSTS(85)*CONSTS(87)*(STATES(13)/(STATES(13)+CONSTS(86)))
      ALGBRC(42) = CONSTS(114) - STATES(22)
      ALGBRC(44) =  CONSTS(66)*STATES(20)*STATES(22) -  CONSTS(67)*STATES(21)*ALGBRC(42)
      ALGBRC(47) =  (STATES(10)/2.00000)*(- CONSTS(118)+ (CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(10) - 1.00000)) ** (1.0 / 2))
      ALGBRC(48) =  CONSTS(33)*STATES(5)*ALGBRC(47)*(ALGBRC(39)/STATES(9))
      ALGBRC(50) =  CONSTS(34)*STATES(6)*ALGBRC(47)
      ALGBRC(51) =  CONSTS(37)*(STATES(12)/(STATES(12)+CONSTS(38)))*(ALGBRC(47)/(ALGBRC(47)+CONSTS(39)))*(STATES(7)/(STATES(7)+CONSTS(40)))*(1.00000 - 1.00000/(1.00000+EXP( - CONSTS(42)*( 1.00000*(STATES(10)/ALGBRC(47) -  1.00000*CONSTS(41))))))
      ALGBRC(45) = CONSTS(115) - STATES(11)
      ALGBRC(52) =  CONSTS(47)*STATES(11)*ALGBRC(47) -  CONSTS(46)*ALGBRC(45)*STATES(10)
      ALGBRC(10) =  CONSTS(45)*(STATES(10)/(STATES(10)+0.00100000))
      ALGBRC(54) = CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(10) - 1.00000)
      ALGBRC(56) = (CONSTS(118)/2.00000+ CONSTS(118)*(CONSTS(117)/( STATES(10)* ALGBRC(54) ** (1.0 / 2)))) - (1.00000+ 0.500000* ALGBRC(54) ** (1.0 / 2))
      ALGBRC(55) =  (STATES(23)/2.00000)*(- CONSTS(118)+ (CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(23) - 1.00000)) ** (1.0 / 2))
      ALGBRC(57) =  CONSTS(64)*STATES(18)*ALGBRC(55)*(ALGBRC(42)/STATES(22))
      ALGBRC(59) =  CONSTS(65)*STATES(19)*ALGBRC(55)
      ALGBRC(60) =  CONSTS(68)*(STATES(25)/(STATES(25)+CONSTS(38)))*(ALGBRC(55)/(ALGBRC(55)+CONSTS(39)))*(STATES(20)/(STATES(20)+CONSTS(40)))*(1.00000 - 1.00000/(1.00000+EXP( 1.00000*- CONSTS(42)*(STATES(23)/ALGBRC(55) -  1.00000*CONSTS(41)))))
      ALGBRC(46) = CONSTS(116) - STATES(24)
      ALGBRC(61) =  CONSTS(75)*STATES(24)*ALGBRC(55) -  CONSTS(74)*ALGBRC(46)*STATES(23)
      ALGBRC(53) =  3.00000*ALGBRC(51)
      ALGBRC(33) =  2.00000*(ALGBRC(31)/CONSTS(78))*(STATES(34) - CONSTS(91))
      ALGBRC(62) =  3.00000*ALGBRC(60)
      ALGBRC(24) =  CONSTS(73)*(STATES(23)/(STATES(23)+0.00100000))
      ALGBRC(63) = CONSTS(118) ** 2.00000+ 4.00000*CONSTS(118)*(CONSTS(117)/STATES(23) - 1.00000)
      ALGBRC(64) = (CONSTS(118)/2.00000+ CONSTS(118)*(CONSTS(117)/( STATES(23)* ALGBRC(63) ** (1.0 / 2)))) - (1.00000+ 0.500000* ALGBRC(63) ** (1.0 / 2))
      ALGBRC(2) =  CONSTS(111)*( (CONSTS(119)+CONSTS(120))*(1.00000 - STATES(36)/CONSTS(122)) -  (CONSTS(120)+CONSTS(121))*(1.00000 - STATES(35)/CONSTS(111)))
      ALGBRC(3) = TERNRY(VOI - (CONSTS(107)+CONSTS(99)).GE.0.00000, 1.00000, 0.00000)
      ALGBRC(49) = CONSTS(117) - (STATES(10)+ALGBRC(47))
      ALGBRC(58) = CONSTS(117) - (STATES(23)+ALGBRC(55))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END