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 303 entries in the algebraic variable array.
C There are a total of 62 entries in each of the rate and state variable arrays.
C There are a total of 177 entries in the constant variable array.
C
C
C CONSTS(1) is kappa_leak in component BG_parameters (fmol_per_sec).
C CONSTS(2) is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec).
C CONSTS(3) is kappa_R_BU in component BG_parameters (fmol_per_sec).
C CONSTS(4) is kappa_R_UW in component BG_parameters (fmol_per_sec).
C CONSTS(5) is kappa_R_WS in component BG_parameters (fmol_per_sec).
C CONSTS(6) is kappa_R_SU in component BG_parameters (fmol_per_sec).
C CONSTS(7) is kappa_diadDf in component BG_parameters (fmol_per_sec).
C CONSTS(8) is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec).
C CONSTS(9) is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec).
C CONSTS(10) is kappa_d000 in component BG_parameters (fmol_per_sec).
C CONSTS(11) is kappa_d010 in component BG_parameters (fmol_per_sec).
C CONSTS(12) is kappa_d020 in component BG_parameters (fmol_per_sec).
C CONSTS(13) is kappa_d001 in component BG_parameters (fmol_per_sec).
C CONSTS(14) is kappa_d011 in component BG_parameters (fmol_per_sec).
C CONSTS(15) is kappa_d021 in component BG_parameters (fmol_per_sec).
C CONSTS(16) is kappa_f1_000 in component BG_parameters (fmol_per_sec).
C CONSTS(17) is kappa_f1_100 in component BG_parameters (fmol_per_sec).
C CONSTS(18) is kappa_f1_001 in component BG_parameters (fmol_per_sec).
C CONSTS(19) is kappa_f1_101 in component BG_parameters (fmol_per_sec).
C CONSTS(20) is kappa_f2_000 in component BG_parameters (fmol_per_sec).
C CONSTS(21) is kappa_f2_100 in component BG_parameters (fmol_per_sec).
C CONSTS(22) is kappa_f2_001 in component BG_parameters (fmol_per_sec).
C CONSTS(23) is kappa_f2_101 in component BG_parameters (fmol_per_sec).
C CONSTS(24) is kappa_f3_010 in component BG_parameters (fmol_per_sec).
C CONSTS(25) is kappa_f3_110 in component BG_parameters (fmol_per_sec).
C CONSTS(26) is kappa_f3_011 in component BG_parameters (fmol_per_sec).
C CONSTS(27) is kappa_f3_111 in component BG_parameters (fmol_per_sec).
C CONSTS(28) is kappa_fCa000 in component BG_parameters (fmol_per_sec).
C CONSTS(29) is kappa_fCa100 in component BG_parameters (fmol_per_sec).
C CONSTS(30) is kappa_fCa010 in component BG_parameters (fmol_per_sec).
C CONSTS(31) is kappa_fCa110 in component BG_parameters (fmol_per_sec).
C CONSTS(32) is kappa_fCa020 in component BG_parameters (fmol_per_sec).
C CONSTS(33) is kappa_fCa120 in component BG_parameters (fmol_per_sec).
C CONSTS(34) is kappa_1_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(35) is kappa_2_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(36) is kappa_3_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(37) is kappa_4_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(38) is kappa_5_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(39) is kappa_6_NCX in component BG_parameters (fmol_per_sec).
C CONSTS(40) is kappa_pCa_R1 in component BG_parameters (fmol_per_sec).
C CONSTS(41) is kappa_pCa_R2 in component BG_parameters (fmol_per_sec).
C CONSTS(42) is kappa_PLBph1 in component BG_parameters (fmol_per_sec).
C CONSTS(43) is kappa_PLBph2 in component BG_parameters (fmol_per_sec).
C CONSTS(44) is kappa_PLBd1 in component BG_parameters (fmol_per_sec).
C CONSTS(45) is kappa_PLBd2 in component BG_parameters (fmol_per_sec).
C CONSTS(46) is kappa_Inh in component BG_parameters (fmol_per_sec).
C CONSTS(47) is kappa_RyR in component BG_parameters (fmol_per_sec).
C CONSTS(48) is kappa_OC in component BG_parameters (fmol_per_sec).
C CONSTS(49) is kappa_CCI in component BG_parameters (fmol_per_sec).
C CONSTS(50) is kappa_CII in component BG_parameters (fmol_per_sec).
C CONSTS(51) is kappa_IO in component BG_parameters (fmol_per_sec).
C CONSTS(52) is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec).
C CONSTS(53) is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec).
C CONSTS(54) is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec).
C CONSTS(55) is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec).
C CONSTS(56) is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec).
C CONSTS(57) is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec).
C CONSTS(58) is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec).
C CONSTS(59) is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec).
C CONSTS(60) is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec).
C CONSTS(61) is K_Ca_SR in component BG_parameters (per_fmol).
C CONSTS(62) is K_Ca_i in component BG_parameters (per_fmol).
C CONSTS(63) is K_TRPN in component BG_parameters (per_fmol).
C CONSTS(64) is K_Ca_TRPN in component BG_parameters (per_fmol).
C CONSTS(65) is K_B_XB in component BG_parameters (per_fmol).
C CONSTS(66) is K_U_XB in component BG_parameters (per_fmol).
C CONSTS(67) is K_W_XB in component BG_parameters (per_fmol).
C CONSTS(68) is K_S_XB in component BG_parameters (per_fmol).
C CONSTS(69) is K_Ca_D in component BG_parameters (per_fmol).
C CONSTS(70) is K_Ca_o in component BG_parameters (per_fmol).
C CONSTS(71) is K_000_LCC in component BG_parameters (per_fmol).
C CONSTS(72) is K_010_LCC in component BG_parameters (per_fmol).
C CONSTS(73) is K_020_LCC in component BG_parameters (per_fmol).
C CONSTS(74) is K_100_LCC in component BG_parameters (per_fmol).
C CONSTS(75) is K_110_LCC in component BG_parameters (per_fmol).
C CONSTS(76) is K_120_LCC in component BG_parameters (per_fmol).
C CONSTS(77) is K_001_LCC in component BG_parameters (per_fmol).
C CONSTS(78) is K_011_LCC in component BG_parameters (per_fmol).
C CONSTS(79) is K_021_LCC in component BG_parameters (per_fmol).
C CONSTS(80) is K_101_LCC in component BG_parameters (per_fmol).
C CONSTS(81) is K_111_LCC in component BG_parameters (per_fmol).
C CONSTS(82) is K_121_LCC in component BG_parameters (per_fmol).
C CONSTS(83) is K_1_NCX in component BG_parameters (per_fmol).
C CONSTS(84) is K_2_NCX in component BG_parameters (per_fmol).
C CONSTS(85) is K_3_NCX in component BG_parameters (per_fmol).
C CONSTS(86) is K_4_NCX in component BG_parameters (per_fmol).
C CONSTS(87) is K_5_NCX in component BG_parameters (per_fmol).
C CONSTS(88) is K_6_NCX in component BG_parameters (per_fmol).
C CONSTS(89) is K_Na_i in component BG_parameters (per_fmol).
C CONSTS(90) is K_Na_o in component BG_parameters (per_fmol).
C CONSTS(91) is K_pCa in component BG_parameters (per_fmol).
C CONSTS(92) is K_pCa_Ca in component BG_parameters (per_fmol).
C CONSTS(93) is K_PLB in component BG_parameters (per_fmol).
C CONSTS(94) is K_PKACI in component BG_parameters (per_fmol).
C CONSTS(95) is K_PLB_PKACI in component BG_parameters (per_fmol).
C CONSTS(96) is K_PP1 in component BG_parameters (per_fmol).
C CONSTS(97) is K_PLBp_PP1 in component BG_parameters (per_fmol).
C CONSTS(98) is K_PLBp in component BG_parameters (per_fmol).
C CONSTS(99) is K_Ip in component BG_parameters (per_fmol).
C CONSTS(100) is K_Ip_PP1 in component BG_parameters (per_fmol).
C CONSTS(101) is K_C_RyR in component BG_parameters (per_fmol).
C CONSTS(102) is K_CI_RyR in component BG_parameters (per_fmol).
C CONSTS(103) is K_I_RyR in component BG_parameters (per_fmol).
C CONSTS(104) is K_O_RyR in component BG_parameters (per_fmol).
C CONSTS(105) is K_P1_SERCA in component BG_parameters (per_fmol).
C CONSTS(106) is K_P2_SERCA in component BG_parameters (per_fmol).
C CONSTS(107) is K_P2a_SERCA in component BG_parameters (per_fmol).
C CONSTS(108) is K_P4_SERCA in component BG_parameters (per_fmol).
C CONSTS(109) is K_P5_SERCA in component BG_parameters (per_fmol).
C CONSTS(110) is K_P6_SERCA in component BG_parameters (per_fmol).
C CONSTS(111) is K_P8_SERCA in component BG_parameters (per_fmol).
C CONSTS(112) is K_P9_SERCA in component BG_parameters (per_fmol).
C CONSTS(113) is K_P10_SERCA in component BG_parameters (per_fmol).
C CONSTS(114) is K_H in component BG_parameters (per_fmol).
C CONSTS(115) is K_MgATP in component BG_parameters (per_fmol).
C CONSTS(116) is K_MgADP in component BG_parameters (per_fmol).
C CONSTS(117) is K_Pi in component BG_parameters (per_fmol).
C CONSTS(118) is z_rCa in component BG_parameters (dimensionless).
C CONSTS(119) is z_fd in component BG_parameters (dimensionless).
C CONSTS(120) is z_ff1 in component BG_parameters (dimensionless).
C CONSTS(121) is z_ff2 in component BG_parameters (dimensionless).
C CONSTS(122) is z_ff3 in component BG_parameters (dimensionless).
C CONSTS(123) is z_rd in component BG_parameters (dimensionless).
C CONSTS(124) is z_rf1 in component BG_parameters (dimensionless).
C CONSTS(125) is z_rf2 in component BG_parameters (dimensionless).
C CONSTS(126) is z_rf3 in component BG_parameters (dimensionless).
C CONSTS(127) is nNa_o_NCX in component BG_parameters (dimensionless).
C CONSTS(128) is nNa_i_NCX in component BG_parameters (dimensionless).
C CONSTS(129) is zf_NCX in component BG_parameters (dimensionless).
C CONSTS(130) is zr_NCX in component BG_parameters (dimensionless).
C CONSTS(131) is n_Tm in component BG_parameters (dimensionless).
C CONSTS(132) is n_CaTRPN in component BG_parameters (dimensionless).
C VOI is time in component environment (second).
C CONSTS(133) is vol_myo in component environment (pL).
C CONSTS(134) is C_m in component environment (fF).
C STATES(1) is q_membrane in component environment (fC).
C STATES(2) is q_Ca_D in component environment (fmol).
C STATES(3) is q_TRPN in component environment (fmol).
C STATES(4) is q_Ca_TRPN in component environment (fmol).
C STATES(5) is q_Ca_i in component environment (fmol).
C STATES(6) is q_Na_i in component environment (fmol).
C STATES(7) is q_Na_o in component environment (fmol).
C STATES(8) is q_B_XB in component environment (fmol).
C STATES(9) is q_U_XB in component environment (fmol).
C STATES(10) is q_W_XB in component environment (fmol).
C STATES(11) is q_S_XB in component environment (fmol).
C STATES(12) is q_Ca_o in component environment (fmol).
C STATES(13) is q_PLB in component environment (fmol).
C STATES(14) is q_PKACI in component environment (fmol).
C STATES(15) is q_PLB_PKACI in component environment (fmol).
C STATES(16) is q_PP1 in component environment (fmol).
C STATES(17) is q_PLBp_PP1 in component environment (fmol).
C STATES(18) is q_PLBp in component environment (fmol).
C CONSTS(135) is q_PLBp_baseline in component environment (fmol).
C STATES(19) is q_Ip in component environment (fmol).
C STATES(20) is q_Ip_PP1 in component environment (fmol).
C STATES(21) is q_Ca_SR in component environment (fmol).
C STATES(22) is q_pCa in component environment (fmol).
C STATES(23) is q_pCa_Ca in component environment (fmol).
C STATES(24) is q_H in component environment (fmol).
C STATES(25) is q_Pi in component environment (fmol).
C STATES(26) is q_MgADP in component environment (fmol).
C STATES(27) is q_MgATP in component environment (fmol).
C STATES(28) is q_000_LCC in component environment (fmol).
C STATES(29) is q_010_LCC in component environment (fmol).
C STATES(30) is q_020_LCC in component environment (fmol).
C STATES(31) is q_100_LCC in component environment (fmol).
C STATES(32) is q_110_LCC in component environment (fmol).
C STATES(33) is q_120_LCC in component environment (fmol).
C STATES(34) is q_001_LCC in component environment (fmol).
C STATES(35) is q_011_LCC in component environment (fmol).
C STATES(36) is q_021_LCC in component environment (fmol).
C STATES(37) is q_101_LCC in component environment (fmol).
C STATES(38) is q_111_LCC in component environment (fmol).
C STATES(39) is q_121_LCC in component environment (fmol).
C STATES(40) is q_P1_NCX in component environment (fmol).
C STATES(41) is q_P2_NCX in component environment (fmol).
C STATES(42) is q_P3_NCX in component environment (fmol).
C STATES(43) is q_P4_NCX in component environment (fmol).
C STATES(44) is q_P5_NCX in component environment (fmol).
C STATES(45) is q_P6_NCX in component environment (fmol).
C STATES(46) is q_C_RyR in component environment (fmol).
C STATES(47) is q_CI_RyR in component environment (fmol).
C STATES(48) is q_I_RyR in component environment (fmol).
C STATES(49) is q_O_RyR in component environment (fmol).
C STATES(50) is q_P1_SERCA in component environment (fmol).
C STATES(51) is q_P2_SERCA in component environment (fmol).
C STATES(52) is q_P2a_SERCA in component environment (fmol).
C STATES(53) is q_P4_SERCA in component environment (fmol).
C STATES(54) is q_P5_SERCA in component environment (fmol).
C STATES(55) is q_P6_SERCA in component environment (fmol).
C STATES(56) is q_P8_SERCA in component environment (fmol).
C STATES(57) is q_P9_SERCA in component environment (fmol).
C STATES(58) is q_P10_SERCA in component environment (fmol).
C ALGBRC(55) is v_leak in component Ca_leak (fmol_per_sec).
C ALGBRC(67) is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(74) is v_R_BU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(76) is v_R_UW in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(78) is v_R_WS in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(80) is v_R_SU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(60) is T_total in component crossbridge_TRPN (kPa).
C STATES(59) is SL in component crossbridge_TRPN (metre).
C ALGBRC(73) is v_diadDf in component diadCaDiffusion (fmol_per_sec).
C ALGBRC(183) is v_Ca_D_LCC in component LCC (fmol_per_sec).
C ALGBRC(182) is v_Ca_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(181) is I_mem_LCC in component LCC (fA).
C ALGBRC(222) is v_Ca_o_NCX in component NCX (fmol_per_sec).
C ALGBRC(223) is v_Ca_i_NCX in component NCX (fmol_per_sec).
C ALGBRC(215) is v_Na_o_NCX in component NCX (fmol_per_sec).
C ALGBRC(204) is v_Na_i_NCX in component NCX (fmol_per_sec).
C ALGBRC(221) is I_mem_NCX in component NCX (fA).
C ALGBRC(228) is v_pCa_R1 in component pCa (fmol_per_sec).
C ALGBRC(229) is v_pCa_R2 in component pCa (fmol_per_sec).
C ALGBRC(231) is I_mem_pCa in component pCa (fA).
C ALGBRC(56) is v_PLBph1 in component PLB (fmol_per_sec).
C ALGBRC(61) is v_PLBph2 in component PLB (fmol_per_sec).
C ALGBRC(63) is v_PLBd1 in component PLB (fmol_per_sec).
C ALGBRC(66) is v_PLBd2 in component PLB (fmol_per_sec).
C ALGBRC(57) is v_Inh in component PLB (fmol_per_sec).
C ALGBRC(238) is v_RyR in component RyR (fmol_per_sec).
C ALGBRC(247) is v_RyRgate_Ca_D in component RyR (fmol_per_sec).
C ALGBRC(290) is v_Ca_i_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(296) is v_Ca_SR_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(301) is v_H_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(287) is v_MgADP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(284) is v_MgATP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(282) is v_Pi_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(2) is V_m in component environment (volt).
C ALGBRC(29) is I_stim in component environment (fA).
C CONSTS(136) is stimPeriod in component environment (second).
C CONSTS(137) is stimDuration in component environment (second).
C ALGBRC(12) is tPeriod in component environment (second).
C ALGBRC(34) is v_Ca_stim in component environment (fmol_per_sec).
C ALGBRC(233) is sum_I in component environment (fA).
C ALGBRC(3) is dSL in component environment (m_per_s).
C CONSTS(138) is V_max in component environment (m_per_s).
C CONSTS(139) is vol_cell in component environment (pL).
C CONSTS(171) is vol_SR in component environment (pL).
C CONSTS(140) is vol_o in component environment (pL).
C CONSTS(141) is vol_D in component environment (pL).
C ALGBRC(4) is cCai in component environment (mM).
C ALGBRC(5) is cCao in component environment (mM).
C ALGBRC(6) is cCaSR in component environment (mM).
C ALGBRC(7) is cCaD in component environment (mM).
C ALGBRC(8) is cNai in component environment (mM).
C ALGBRC(9) is pHi in component environment (dimensionless).
C ALGBRC(10) is Ca_tot in component environment (fmol).
C CONSTS(142) is R in component constants (J_per_K_per_mol).
C CONSTS(143) is T in component constants (kelvin).
C CONSTS(144) is F in component constants (C_per_mol).
C CONSTS(145) is zNa in component ion_valences (dimensionless).
C CONSTS(146) is zCl in component ion_valences (dimensionless).
C CONSTS(147) is zK in component ion_valences (dimensionless).
C CONSTS(148) is zCa in component ion_valences (dimensionless).
C ALGBRC(39) is mu_Ca_SR in component Ca_leak (J_per_mol).
C ALGBRC(49) is mu_Ca_i in component Ca_leak (J_per_mol).
C ALGBRC(59) is mu_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(62) is mu_Ca_i in component crossbridge_TRPN (J_per_mol).
C ALGBRC(64) is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(65) is mu_B_XB in component crossbridge_TRPN (J_per_mol).
C ALGBRC(68) is mu_U_XB in component crossbridge_TRPN (J_per_mol).
C ALGBRC(70) is mu_W_XB in component crossbridge_TRPN (J_per_mol).
C ALGBRC(72) is mu_S_XB in component crossbridge_TRPN (J_per_mol).
C CONSTS(149) is SL_0 in component crossbridge_TRPN (metre).
C CONSTS(150) is q_MS in component crossbridge_TRPN (fmol).
C CONSTS(151) is r_s in component crossbridge_TRPN (dimensionless).
C CONSTS(152) is r_w in component crossbridge_TRPN (dimensionless).
C CONSTS(153) is A_eff in component crossbridge_TRPN (dimensionless).
C CONSTS(154) is phi in component crossbridge_TRPN (dimensionless).
C CONSTS(155) is k_uw in component crossbridge_TRPN (per_sec).
C CONSTS(156) is k_ws in component crossbridge_TRPN (per_sec).
C STATES(60) is G_w in component crossbridge_TRPN (metre).
C STATES(61) is G_s in component crossbridge_TRPN (metre).
C CONSTS(172) is c_w in component crossbridge_TRPN (per_sec).
C CONSTS(173) is c_s in component crossbridge_TRPN (per_sec).
C CONSTS(170) is A_w in component crossbridge_TRPN (dimensionless).
C CONSTS(176) is A_s in component crossbridge_TRPN (dimensionless).
C ALGBRC(42) is nu_1 in component crossbridge_TRPN (m_per_s).
C ALGBRC(31) is nu_2 in component crossbridge_TRPN (m_per_s).
C ALGBRC(36) is nu_3 in component crossbridge_TRPN (m_per_s).
C ALGBRC(43) is nu_4 in component crossbridge_TRPN (m_per_s).
C ALGBRC(32) is nu_5 in component crossbridge_TRPN (m_per_s).
C ALGBRC(37) is nu_6 in component crossbridge_TRPN (m_per_s).
C ALGBRC(52) is v_1 in component crossbridge_TRPN (m_per_s).
C ALGBRC(27) is v_2 in component crossbridge_TRPN (m_per_s).
C ALGBRC(53) is v_4 in component crossbridge_TRPN (m_per_s).
C CONSTS(157) is v_to_mu in component crossbridge_TRPN (Js_per_m2).
C ALGBRC(33) is T_active in component crossbridge_TRPN (kPa).
C CONSTS(158) is T_ref in component crossbridge_TRPN (kPa).
C ALGBRC(11) is F_T_S in component crossbridge_TRPN (kPa).
C ALGBRC(28) is F_T_W in component crossbridge_TRPN (kPa).
C STATES(62) is Cdd in component crossbridge_TRPN (metre).
C ALGBRC(44) is eta in component crossbridge_TRPN (second).
C CONSTS(159) is eta_l in component crossbridge_TRPN (second).
C CONSTS(160) is eta_s in component crossbridge_TRPN (second).
C CONSTS(161) is k in component crossbridge_TRPN (dimensionless).
C CONSTS(162) is alpha in component crossbridge_TRPN (kPa).
C ALGBRC(45) is v_Cdd in component crossbridge_TRPN (m_per_s).
C ALGBRC(46) is nu_d in component crossbridge_TRPN (m_per_s).
C ALGBRC(47) is nu_k in component crossbridge_TRPN (m_per_s).
C ALGBRC(48) is nu_f in component crossbridge_TRPN (m_per_s).
C ALGBRC(54) is F_d in component crossbridge_TRPN (kPa).
C ALGBRC(38) is F_1 in component crossbridge_TRPN (kPa).
C ALGBRC(58) is T_passive in component crossbridge_TRPN (kPa).
C CONSTS(163) is b in component crossbridge_TRPN (dimensionless).
C ALGBRC(69) is mu_Ca_D in component diadCaDiffusion (J_per_mol).
C ALGBRC(71) is mu_Ca_i in component diadCaDiffusion (J_per_mol).
C ALGBRC(77) is mu_Ca_o in component LCC (J_per_mol).
C ALGBRC(79) is mu_Ca_D in component LCC (J_per_mol).
C ALGBRC(75) is V_mem in component LCC (volt).
C ALGBRC(89) is Af_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(117) is Ar_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(115) is Am_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(176) is v_LCC_Ca1 in component LCC (fmol_per_sec).
C ALGBRC(90) is Af_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(118) is Ar_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(116) is Am_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(178) is v_LCC_Ca2 in component LCC (fmol_per_sec).
C ALGBRC(109) is Af_fCa00 in component LCC (J_per_mol).
C ALGBRC(155) is Ar_fCa00 in component LCC (J_per_mol).
C ALGBRC(156) is v_fCa00 in component LCC (fmol_per_sec).
C ALGBRC(110) is Af_fCa01 in component LCC (J_per_mol).
C ALGBRC(157) is Ar_fCa01 in component LCC (J_per_mol).
C ALGBRC(160) is v_fCa01 in component LCC (fmol_per_sec).
C ALGBRC(111) is Af_fCa02 in component LCC (J_per_mol).
C ALGBRC(161) is Ar_fCa02 in component LCC (J_per_mol).
C ALGBRC(164) is v_fCa02 in component LCC (fmol_per_sec).
C ALGBRC(112) is Af_fCa10 in component LCC (J_per_mol).
C ALGBRC(165) is Ar_fCa10 in component LCC (J_per_mol).
C ALGBRC(168) is v_fCa10 in component LCC (fmol_per_sec).
C ALGBRC(113) is Af_fCa11 in component LCC (J_per_mol).
C ALGBRC(169) is Ar_fCa11 in component LCC (J_per_mol).
C ALGBRC(172) is v_fCa11 in component LCC (fmol_per_sec).
C ALGBRC(114) is Af_fCa12 in component LCC (J_per_mol).
C ALGBRC(173) is Ar_fCa12 in component LCC (J_per_mol).
C ALGBRC(175) is v_fCa12 in component LCC (fmol_per_sec).
C ALGBRC(81) is mu_S000_LCC in component LCC (J_per_mol).
C ALGBRC(158) is v_S000_LCC in component LCC (fmol_per_sec).
C ALGBRC(82) is mu_S010_LCC in component LCC (J_per_mol).
C ALGBRC(162) is v_S010_LCC in component LCC (fmol_per_sec).
C ALGBRC(83) is mu_S020_LCC in component LCC (J_per_mol).
C ALGBRC(166) is v_S020_LCC in component LCC (fmol_per_sec).
C ALGBRC(84) is mu_S100_LCC in component LCC (J_per_mol).
C ALGBRC(170) is v_S100_LCC in component LCC (fmol_per_sec).
C ALGBRC(85) is mu_S110_LCC in component LCC (J_per_mol).
C ALGBRC(174) is v_S110_LCC in component LCC (fmol_per_sec).
C ALGBRC(87) is mu_S120_LCC in component LCC (J_per_mol).
C ALGBRC(177) is v_S120_LCC in component LCC (fmol_per_sec).
C ALGBRC(13) is mu_S001_LCC in component LCC (J_per_mol).
C ALGBRC(159) is v_S001_LCC in component LCC (fmol_per_sec).
C ALGBRC(14) is mu_S011_LCC in component LCC (J_per_mol).
C ALGBRC(163) is v_S011_LCC in component LCC (fmol_per_sec).
C ALGBRC(15) is mu_S021_LCC in component LCC (J_per_mol).
C ALGBRC(167) is v_S021_LCC in component LCC (fmol_per_sec).
C ALGBRC(16) is mu_S101_LCC in component LCC (J_per_mol).
C ALGBRC(171) is v_S101_LCC in component LCC (fmol_per_sec).
C ALGBRC(86) is mu_S111_LCC in component LCC (J_per_mol).
C ALGBRC(179) is v_S111_LCC in component LCC (fmol_per_sec).
C ALGBRC(88) is mu_S121_LCC in component LCC (J_per_mol).
C ALGBRC(180) is v_S121_LCC in component LCC (fmol_per_sec).
C ALGBRC(91) is Af_d000 in component LCC (J_per_mol).
C ALGBRC(119) is Ar_d000 in component LCC (J_per_mol).
C ALGBRC(120) is v_d000 in component LCC (fmol_per_sec).
C ALGBRC(93) is Af_d010 in component LCC (J_per_mol).
C ALGBRC(123) is Ar_d010 in component LCC (J_per_mol).
C ALGBRC(124) is v_d010 in component LCC (fmol_per_sec).
C ALGBRC(95) is Af_d020 in component LCC (J_per_mol).
C ALGBRC(127) is Ar_d020 in component LCC (J_per_mol).
C ALGBRC(128) is v_d020 in component LCC (fmol_per_sec).
C ALGBRC(97) is Af_f1_000 in component LCC (J_per_mol).
C ALGBRC(131) is Ar_f1_000 in component LCC (J_per_mol).
C ALGBRC(132) is v_f1_000 in component LCC (fmol_per_sec).
C ALGBRC(101) is Af_f2_000 in component LCC (J_per_mol).
C ALGBRC(139) is Ar_f2_000 in component LCC (J_per_mol).
C ALGBRC(140) is v_f2_000 in component LCC (fmol_per_sec).
C ALGBRC(105) is Af_f3_010 in component LCC (J_per_mol).
C ALGBRC(147) is Ar_f3_010 in component LCC (J_per_mol).
C ALGBRC(148) is v_f3_010 in component LCC (fmol_per_sec).
C ALGBRC(99) is Af_f1_100 in component LCC (J_per_mol).
C ALGBRC(135) is Ar_f1_100 in component LCC (J_per_mol).
C ALGBRC(136) is v_f1_100 in component LCC (fmol_per_sec).
C ALGBRC(103) is Af_f2_100 in component LCC (J_per_mol).
C ALGBRC(143) is Ar_f2_100 in component LCC (J_per_mol).
C ALGBRC(144) is v_f2_100 in component LCC (fmol_per_sec).
C ALGBRC(107) is Af_f3_110 in component LCC (J_per_mol).
C ALGBRC(151) is Ar_f3_110 in component LCC (J_per_mol).
C ALGBRC(152) is v_f3_110 in component LCC (fmol_per_sec).
C ALGBRC(92) is Af_d001 in component LCC (J_per_mol).
C ALGBRC(121) is Ar_d001 in component LCC (J_per_mol).
C ALGBRC(122) is v_d001 in component LCC (fmol_per_sec).
C ALGBRC(94) is Af_d011 in component LCC (J_per_mol).
C ALGBRC(125) is Ar_d011 in component LCC (J_per_mol).
C ALGBRC(126) is v_d011 in component LCC (fmol_per_sec).
C ALGBRC(96) is Af_d021 in component LCC (J_per_mol).
C ALGBRC(129) is Ar_d021 in component LCC (J_per_mol).
C ALGBRC(130) is v_d021 in component LCC (fmol_per_sec).
C ALGBRC(98) is Af_f1_001 in component LCC (J_per_mol).
C ALGBRC(133) is Ar_f1_001 in component LCC (J_per_mol).
C ALGBRC(134) is v_f1_001 in component LCC (fmol_per_sec).
C ALGBRC(102) is Af_f2_001 in component LCC (J_per_mol).
C ALGBRC(141) is Ar_f2_001 in component LCC (J_per_mol).
C ALGBRC(142) is v_f2_001 in component LCC (fmol_per_sec).
C ALGBRC(106) is Af_f3_011 in component LCC (J_per_mol).
C ALGBRC(149) is Ar_f3_011 in component LCC (J_per_mol).
C ALGBRC(150) is v_f3_011 in component LCC (fmol_per_sec).
C ALGBRC(100) is Af_f1_101 in component LCC (J_per_mol).
C ALGBRC(137) is Ar_f1_101 in component LCC (J_per_mol).
C ALGBRC(138) is v_f1_101 in component LCC (fmol_per_sec).
C ALGBRC(104) is Af_f2_101 in component LCC (J_per_mol).
C ALGBRC(145) is Ar_f2_101 in component LCC (J_per_mol).
C ALGBRC(146) is v_f2_101 in component LCC (fmol_per_sec).
C ALGBRC(108) is Af_f3_111 in component LCC (J_per_mol).
C ALGBRC(153) is Ar_f3_111 in component LCC (J_per_mol).
C ALGBRC(154) is v_f3_111 in component LCC (fmol_per_sec).
C ALGBRC(17) is Ca_tot in component LCC (fmol).
C ALGBRC(185) is mu_Ca_o in component NCX (J_per_mol).
C ALGBRC(186) is mu_Ca_i in component NCX (J_per_mol).
C ALGBRC(187) is mu_Na_o in component NCX (J_per_mol).
C ALGBRC(188) is mu_Na_i in component NCX (J_per_mol).
C ALGBRC(189) is mu_P1_NCX in component NCX (J_per_mol).
C ALGBRC(190) is mu_P2_NCX in component NCX (J_per_mol).
C ALGBRC(191) is mu_P3_NCX in component NCX (J_per_mol).
C ALGBRC(192) is mu_P4_NCX in component NCX (J_per_mol).
C ALGBRC(193) is mu_P5_NCX in component NCX (J_per_mol).
C ALGBRC(194) is mu_P6_NCX in component NCX (J_per_mol).
C ALGBRC(219) is v_P1_NCX in component NCX (fmol_per_sec).
C ALGBRC(207) is v_P2_NCX in component NCX (fmol_per_sec).
C ALGBRC(210) is v_P3_NCX in component NCX (fmol_per_sec).
C ALGBRC(213) is v_P4_NCX in component NCX (fmol_per_sec).
C ALGBRC(216) is v_P5_NCX in component NCX (fmol_per_sec).
C ALGBRC(220) is v_P6_NCX in component NCX (fmol_per_sec).
C ALGBRC(195) is Af_r1 in component NCX (J_per_mol).
C ALGBRC(201) is Ar_r1 in component NCX (J_per_mol).
C ALGBRC(202) is v_r1 in component NCX (fmol_per_sec).
C ALGBRC(196) is Af_r2 in component NCX (J_per_mol).
C ALGBRC(203) is Ar_r2 in component NCX (J_per_mol).
C ALGBRC(205) is v_r2 in component NCX (fmol_per_sec).
C ALGBRC(197) is Af_r3 in component NCX (J_per_mol).
C ALGBRC(206) is Ar_r3 in component NCX (J_per_mol).
C ALGBRC(208) is v_r3 in component NCX (fmol_per_sec).
C ALGBRC(198) is Af_r4 in component NCX (J_per_mol).
C ALGBRC(209) is Ar_r4 in component NCX (J_per_mol).
C ALGBRC(211) is v_r4 in component NCX (fmol_per_sec).
C ALGBRC(199) is Af_r5 in component NCX (J_per_mol).
C ALGBRC(212) is Ar_r5 in component NCX (J_per_mol).
C ALGBRC(214) is v_r5 in component NCX (fmol_per_sec).
C ALGBRC(200) is Af_r6 in component NCX (J_per_mol).
C ALGBRC(217) is Ar_r6 in component NCX (J_per_mol).
C ALGBRC(218) is v_r6 in component NCX (fmol_per_sec).
C ALGBRC(184) is V_mem in component NCX (volt).
C ALGBRC(224) is mu_Ca_i in component pCa (J_per_mol).
C ALGBRC(225) is mu_Ca_o in component pCa (J_per_mol).
C ALGBRC(226) is mu_pCa in component pCa (J_per_mol).
C ALGBRC(227) is mu_pCa_Ca in component pCa (J_per_mol).
C ALGBRC(1) is mu_PLB in component PLB (J_per_mol).
C ALGBRC(26) is mu_PKACI in component PLB (J_per_mol).
C ALGBRC(30) is mu_PLB_PKACI in component PLB (J_per_mol).
C ALGBRC(35) is mu_PP1 in component PLB (J_per_mol).
C ALGBRC(40) is mu_PLBp_PP1 in component PLB (J_per_mol).
C ALGBRC(50) is mu_PLBp in component PLB (J_per_mol).
C ALGBRC(41) is mu_Ip in component PLB (J_per_mol).
C ALGBRC(51) is mu_Ip_PP1 in component PLB (J_per_mol).
C CONSTS(164) is nCa_1 in component RyR (dimensionless).
C CONSTS(165) is nCa_2 in component RyR (dimensionless).
C ALGBRC(230) is mu_Ca_SR in component RyR (J_per_mol).
C ALGBRC(232) is mu_Ca_D in component RyR (J_per_mol).
C ALGBRC(18) is q_Ca_gate_complexes in component RyR (fmol).
C ALGBRC(234) is mu_C_RyR in component RyR (J_per_mol).
C ALGBRC(235) is mu_CI_RyR in component RyR (J_per_mol).
C ALGBRC(236) is mu_I_RyR in component RyR (J_per_mol).
C ALGBRC(237) is mu_O_RyR in component RyR (J_per_mol).
C ALGBRC(239) is v_OC in component RyR (fmol_per_sec).
C ALGBRC(241) is v_CCI in component RyR (fmol_per_sec).
C ALGBRC(243) is v_CII in component RyR (fmol_per_sec).
C ALGBRC(245) is v_IO in component RyR (fmol_per_sec).
C CONSTS(166) is n_Ca_i in component SERCA (dimensionless).
C CONSTS(167) is n_Ca_SR in component SERCA (dimensionless).
C CONSTS(168) is n_H in component SERCA (dimensionless).
C ALGBRC(19) is q_SERCA_Ca_complexes in component SERCA (fmol).
C ALGBRC(20) is c_Ca_i in component SERCA (mM).
C ALGBRC(21) is c_Ca_SR in component SERCA (mM).
C ALGBRC(22) is c_H in component SERCA (mM).
C ALGBRC(23) is c_MgADP in component SERCA (mM).
C ALGBRC(24) is c_MgATP in component SERCA (mM).
C ALGBRC(25) is c_Pi in component SERCA (mM).
C CONSTS(169) is vol_i in component SERCA (pL).
C CONSTS(175) is vol_sr in component SERCA (pL).
C CONSTS(177) is vol_isr in component SERCA (pL).
C ALGBRC(240) is mu_Ca_i in component SERCA (J_per_mol).
C ALGBRC(242) is mu_Ca_SR in component SERCA (J_per_mol).
C ALGBRC(244) is mu_H in component SERCA (J_per_mol).
C ALGBRC(246) is mu_MgADP in component SERCA (J_per_mol).
C ALGBRC(248) is mu_MgATP in component SERCA (J_per_mol).
C ALGBRC(250) is mu_P1 in component SERCA (J_per_mol).
C ALGBRC(285) is v_P1 in component SERCA (fmol_per_sec).
C ALGBRC(249) is mu_Pi in component SERCA (J_per_mol).
C ALGBRC(251) is mu_P2 in component SERCA (J_per_mol).
C ALGBRC(292) is v_P2 in component SERCA (fmol_per_sec).
C ALGBRC(252) is mu_P2a in component SERCA (J_per_mol).
C ALGBRC(289) is v_P2a in component SERCA (fmol_per_sec).
C ALGBRC(253) is mu_P4 in component SERCA (J_per_mol).
C ALGBRC(294) is v_P4 in component SERCA (fmol_per_sec).
C ALGBRC(254) is mu_P5 in component SERCA (J_per_mol).
C ALGBRC(295) is v_P5 in component SERCA (fmol_per_sec).
C ALGBRC(256) is mu_P6 in component SERCA (J_per_mol).
C ALGBRC(298) is v_P6 in component SERCA (fmol_per_sec).
C ALGBRC(258) is mu_P8 in component SERCA (J_per_mol).
C ALGBRC(300) is v_P8 in component SERCA (fmol_per_sec).
C ALGBRC(259) is mu_P9 in component SERCA (J_per_mol).
C ALGBRC(302) is v_P9 in component SERCA (fmol_per_sec).
C ALGBRC(255) is mu_P10 in component SERCA (J_per_mol).
C ALGBRC(303) is v_P10 in component SERCA (fmol_per_sec).
C ALGBRC(260) is mu_PLBp_baseline in component SERCA (J_per_mol).
C CONSTS(174) is mu_PLBp in component SERCA (J_per_mol).
C ALGBRC(261) is mu_PLBp_eff in component SERCA (J_per_mol).
C ALGBRC(264) is Af_R1_2 in component SERCA (J_per_mol).
C ALGBRC(265) is Ar_R1_2 in component SERCA (J_per_mol).
C ALGBRC(281) is v_SERCA_R1_2 in component SERCA (fmol_per_sec).
C ALGBRC(266) is Af_R5_6 in component SERCA (J_per_mol).
C ALGBRC(267) is Ar_R5_6 in component SERCA (J_per_mol).
C ALGBRC(283) is v_SERCA_R5_6 in component SERCA (fmol_per_sec).
C ALGBRC(268) is Af_R2_2a in component SERCA (J_per_mol).
C ALGBRC(269) is Ar_R2_2a in component SERCA (J_per_mol).
C ALGBRC(286) is v_SERCA_R2_2a in component SERCA (fmol_per_sec).
C ALGBRC(270) is Af_R2_4 in component SERCA (J_per_mol).
C ALGBRC(271) is Ar_R2_4 in component SERCA (J_per_mol).
C ALGBRC(288) is v_SERCA_R2_4 in component SERCA (fmol_per_sec).
C ALGBRC(272) is Af_R4_5 in component SERCA (J_per_mol).
C ALGBRC(273) is Ar_R4_5 in component SERCA (J_per_mol).
C ALGBRC(291) is v_SERCA_R4_5 in component SERCA (fmol_per_sec).
C ALGBRC(274) is Af_R6_8 in component SERCA (J_per_mol).
C ALGBRC(275) is Ar_R6_8 in component SERCA (J_per_mol).
C ALGBRC(293) is v_SERCA_R6_8 in component SERCA (fmol_per_sec).
C ALGBRC(276) is Af_R8_9 in component SERCA (J_per_mol).
C ALGBRC(277) is Ar_R8_9 in component SERCA (J_per_mol).
C ALGBRC(297) is v_SERCA_R8_9 in component SERCA (fmol_per_sec).
C ALGBRC(278) is Af_R9_10 in component SERCA (J_per_mol).
C ALGBRC(279) is Ar_R9_10 in component SERCA (J_per_mol).
C ALGBRC(299) is v_SERCA_R9_10 in component SERCA (fmol_per_sec).
C ALGBRC(262) is Af_R10_1 in component SERCA (J_per_mol).
C ALGBRC(263) is Ar_R10_1 in component SERCA (J_per_mol).
C ALGBRC(280) is v_SERCA_R10_1 in component SERCA (fmol_per_sec).
C ALGBRC(257) is v_SERCA_R10_1_noPLBp in component SERCA (fmol_per_sec).
C RATES(21) is d/dt q_Ca_SR in component environment (fmol).
C RATES(5) is d/dt q_Ca_i in component environment (fmol).
C RATES(3) is d/dt q_TRPN in component environment (fmol).
C RATES(4) is d/dt q_Ca_TRPN in component environment (fmol).
C RATES(8) is d/dt q_B_XB in component environment (fmol).
C RATES(9) is d/dt q_U_XB in component environment (fmol).
C RATES(10) is d/dt q_W_XB in component environment (fmol).
C RATES(11) is d/dt q_S_XB in component environment (fmol).
C RATES(2) is d/dt q_Ca_D in component environment (fmol).
C RATES(12) is d/dt q_Ca_o in component environment (fmol).
C RATES(6) is d/dt q_Na_i in component environment (fmol).
C RATES(7) is d/dt q_Na_o in component environment (fmol).
C RATES(22) is d/dt q_pCa in component environment (fmol).
C RATES(23) is d/dt q_pCa_Ca in component environment (fmol).
C RATES(13) is d/dt q_PLB in component environment (fmol).
C RATES(14) is d/dt q_PKACI in component environment (fmol).
C RATES(15) is d/dt q_PLB_PKACI in component environment (fmol).
C RATES(16) is d/dt q_PP1 in component environment (fmol).
C RATES(17) is d/dt q_PLBp_PP1 in component environment (fmol).
C RATES(18) is d/dt q_PLBp in component environment (fmol).
C RATES(19) is d/dt q_Ip in component environment (fmol).
C RATES(20) is d/dt q_Ip_PP1 in component environment (fmol).
C RATES(24) is d/dt q_H in component environment (fmol).
C RATES(26) is d/dt q_MgADP in component environment (fmol).
C RATES(27) is d/dt q_MgATP in component environment (fmol).
C RATES(25) is d/dt q_Pi in component environment (fmol).
C RATES(1) is d/dt q_membrane in component environment (fC).
C RATES(60) is d/dt G_w in component crossbridge_TRPN (metre).
C RATES(59) is d/dt SL in component crossbridge_TRPN (metre).
C RATES(61) is d/dt G_s in component crossbridge_TRPN (metre).
C RATES(62) is d/dt Cdd in component crossbridge_TRPN (metre).
C RATES(28) is d/dt q_000_LCC in component environment (fmol).
C RATES(29) is d/dt q_010_LCC in component environment (fmol).
C RATES(30) is d/dt q_020_LCC in component environment (fmol).
C RATES(31) is d/dt q_100_LCC in component environment (fmol).
C RATES(32) is d/dt q_110_LCC in component environment (fmol).
C RATES(33) is d/dt q_120_LCC in component environment (fmol).
C RATES(34) is d/dt q_001_LCC in component environment (fmol).
C RATES(35) is d/dt q_011_LCC in component environment (fmol).
C RATES(36) is d/dt q_021_LCC in component environment (fmol).
C RATES(37) is d/dt q_101_LCC in component environment (fmol).
C RATES(38) is d/dt q_111_LCC in component environment (fmol).
C RATES(39) is d/dt q_121_LCC in component environment (fmol).
C RATES(40) is d/dt q_P1_NCX in component environment (fmol).
C RATES(41) is d/dt q_P2_NCX in component environment (fmol).
C RATES(42) is d/dt q_P3_NCX in component environment (fmol).
C RATES(43) is d/dt q_P4_NCX in component environment (fmol).
C RATES(44) is d/dt q_P5_NCX in component environment (fmol).
C RATES(45) is d/dt q_P6_NCX in component environment (fmol).
C RATES(49) is d/dt q_O_RyR in component environment (fmol).
C RATES(46) is d/dt q_C_RyR in component environment (fmol).
C RATES(47) is d/dt q_CI_RyR in component environment (fmol).
C RATES(48) is d/dt q_I_RyR in component environment (fmol).
C RATES(50) is d/dt q_P1_SERCA in component environment (fmol).
C RATES(51) is d/dt q_P2_SERCA in component environment (fmol).
C RATES(52) is d/dt q_P2a_SERCA in component environment (fmol).
C RATES(53) is d/dt q_P4_SERCA in component environment (fmol).
C RATES(54) is d/dt q_P5_SERCA in component environment (fmol).
C RATES(58) is d/dt q_P10_SERCA in component environment (fmol).
C RATES(55) is d/dt q_P6_SERCA in component environment (fmol).
C RATES(56) is d/dt q_P8_SERCA in component environment (fmol).
C RATES(57) is d/dt q_P9_SERCA in component environment (fmol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1.19278e-07
      CONSTS(2) = 61.975
      CONSTS(3) = 6.70563
      CONSTS(4) = 0.295103
      CONSTS(5) = 0.00694361
      CONSTS(6) = 124985
      CONSTS(7) = 0.0319966
      CONSTS(8) = 1.45746
      CONSTS(9) = 1.91895
      CONSTS(10) = 306.532
      CONSTS(11) = 15.4798
      CONSTS(12) = 20.3814
      CONSTS(13) = 0.0998404
      CONSTS(14) = 0.00504193
      CONSTS(15) = 0.00663841
      CONSTS(16) = 29.5266
      CONSTS(17) = 146.621
      CONSTS(18) = 0.00961708
      CONSTS(19) = 0.047756
      CONSTS(20) = 1.91912
      CONSTS(21) = 9.52987
      CONSTS(22) = 0.000625075
      CONSTS(23) = 0.00310397
      CONSTS(24) = 4187.14
      CONSTS(25) = 20792.3
      CONSTS(26) = 1.36379
      CONSTS(27) = 6.77224
      CONSTS(28) = 34185.2
      CONSTS(29) = 169755
      CONSTS(30) = 1726.35
      CONSTS(31) = 8572.61
      CONSTS(32) = 2272.98
      CONSTS(33) = 11287.1
      CONSTS(34) = 429216
      CONSTS(35) = 1.06923e+08
      CONSTS(36) = 3622.96
      CONSTS(37) = 1.06923e+08
      CONSTS(38) = 429216
      CONSTS(39) = 0.380974
      CONSTS(40) = 10726.4
      CONSTS(41) = 1.40741e-05
      CONSTS(42) = 12.9276
      CONSTS(43) = 1.86249
      CONSTS(44) = 0.0640148
      CONSTS(45) = 0.200763
      CONSTS(46) = 1269.27
      CONSTS(47) = 1.13157e+06
      CONSTS(48) = 8.50173
      CONSTS(49) = 0.164843
      CONSTS(50) = 626.393
      CONSTS(51) = 52.1994
      CONSTS(52) = 3.59662e-07
      CONSTS(53) = 6.73317
      CONSTS(54) = 32063.4
      CONSTS(55) = 6.73317
      CONSTS(56) = 0.408566
      CONSTS(57) = 313.54
      CONSTS(58) = 6675
      CONSTS(59) = 6675
      CONSTS(60) = 5.05704e-06
      CONSTS(61) = 4.69975
      CONSTS(62) = 1.11404
      CONSTS(63) = 12.0566
      CONSTS(64) = 0.062444
      CONSTS(65) = 9.09233
      CONSTS(66) = 3.40963
      CONSTS(67) = 22.2937
      CONSTS(68) = 5.57343e-06
      CONSTS(69) = 66.2656
      CONSTS(70) = 0.0265011
      CONSTS(71) = 1.58796
      CONSTS(72) = 31.4449
      CONSTS(73) = 23.8827
      CONSTS(74) = 0.319783
      CONSTS(75) = 6.33236
      CONSTS(76) = 4.80948
      CONSTS(77) = 4875.4
      CONSTS(78) = 96542.9
      CONSTS(79) = 73325.1
      CONSTS(80) = 981.805
      CONSTS(81) = 19441.8
      CONSTS(82) = 14766.2
      CONSTS(83) = 2.85254e+08
      CONSTS(84) = 16.9249
      CONSTS(85) = 34379.9
      CONSTS(86) = 313010
      CONSTS(87) = 6.01862e+07
      CONSTS(88) = 929.686
      CONSTS(89) = 0.012842
      CONSTS(90) = 0.0107799
      CONSTS(91) = 0.0126813
      CONSTS(92) = 0.0178289
      CONSTS(93) = 0.482706
      CONSTS(94) = 0.0120458
      CONSTS(95) = 2.99358
      CONSTS(96) = 0.0189254
      CONSTS(97) = 1.63848
      CONSTS(98) = 66.755
      CONSTS(99) = 62.347
      CONSTS(100) = 0.0304897
      CONSTS(101) = 2234.81
      CONSTS(102) = 30.332
      CONSTS(103) = 0.0957866
      CONSTS(104) = 7.05739
      CONSTS(105) = 62.115
      CONSTS(106) = 1.88391
      CONSTS(107) = 0.271481
      CONSTS(108) = 1292.79
      CONSTS(109) = 0.00208939
      CONSTS(110) = 139.3
      CONSTS(111) = 0.710561
      CONSTS(112) = 1.30406
      CONSTS(113) = 1.3624
      CONSTS(114) = 486.59
      CONSTS(115) = 15.0787
      CONSTS(116) = 1.53572e-05
      CONSTS(117) = 0.00616948
      CONSTS(118) = 2
      CONSTS(119) = 2.1404
      CONSTS(120) = -1.1495
      CONSTS(121) = 0.72162
      CONSTS(122) = 4.2933
      CONSTS(123) = -2.1404
      CONSTS(124) = 1.8993
      CONSTS(125) = -0.52288
      CONSTS(126) = 0
      CONSTS(127) = 3
      CONSTS(128) = 3
      CONSTS(129) = -0.55679
      CONSTS(130) = 0.44321
      CONSTS(131) = 1.1
      CONSTS(132) = 2
      CONSTS(133) = 25.84
      CONSTS(134) = 153400
      STATES(1) = -13039
      STATES(2) = 3.68876E-07
      STATES(3) = 2.57
      STATES(4) = 1e-18
      STATES(5) = 0.005663912
      STATES(6) = 481.46
      STATES(7) = 725.48
      STATES(8) = 1
      STATES(9) = 0
      STATES(10) = 0
      STATES(11) = 0
      STATES(12) = 6.84
      STATES(13) = 4.028E+00
      STATES(14) = 2.234E-03
      STATES(15) = 1e-18
      STATES(16) = 3.382E-02
      STATES(17) = 1e-18
      STATES(18) = 1e-18
      CONSTS(135) = 1e-18
      STATES(19) = 1.999E-03
      STATES(20) = 1e-18
      STATES(21) = 0.641
      STATES(22) = 0.0032
      STATES(23) = 1e-9
      STATES(24) = 2.584e-6
      STATES(25) = 570
      STATES(26) = 1.3794
      STATES(27) = 3.8
      STATES(28) = 1.3798e-05
      STATES(29) = 1.3798e-05
      STATES(30) = 1.3798e-05
      STATES(31) = 1.3798e-05
      STATES(32) = 1.3798e-05
      STATES(33) = 1.3798e-05
      STATES(34) = 1.3798e-05
      STATES(35) = 1.3798e-05
      STATES(36) = 1.3798e-05
      STATES(37) = 1.3798e-05
      STATES(38) = 1.3798e-05
      STATES(39) = 1.3798e-05
      STATES(40) = 4.24554e-07
      STATES(41) = 4.24554e-07
      STATES(42) = 4.24554e-07
      STATES(43) = 4.24554e-07
      STATES(44) = 4.24554e-07
      STATES(45) = 4.24554e-07
      STATES(46) = 4.13941e-05
      STATES(47) = 4.13941e-05
      STATES(48) = 4.13941e-05
      STATES(49) = 4.13941e-05
      STATES(50) = 2.83036e-08
      STATES(51) = 2.83036e-08
      STATES(52) = 2.83036e-08
      STATES(53) = 2.83036e-08
      STATES(54) = 2.83036e-08
      STATES(55) = 2.83036e-08
      STATES(56) = 2.83036e-08
      STATES(57) = 2.83036e-08
      STATES(58) = 2.83036e-08
      STATES(59) = 2e-6
      CONSTS(136) = 1
      CONSTS(137) = 0.001
      CONSTS(138) = 19e-6
      CONSTS(139) = 38
      CONSTS(140) = 5.182
      CONSTS(141) = 0.0025375
      CONSTS(142) = 8.31
      CONSTS(143) = 310
      CONSTS(144) = 96485
      CONSTS(145) = 1
      CONSTS(146) = -1
      CONSTS(147) = 1
      CONSTS(148) = 2
      CONSTS(149) = 2e-6
      CONSTS(150) = 1e-6
      CONSTS(151) = 0.25
      CONSTS(152) = 0.5
      CONSTS(153) = 25
      CONSTS(154) = 2.23
      CONSTS(155) = 26
      CONSTS(156) = 4
      STATES(60) = 1e-6
      STATES(61) = 1e-6
      CONSTS(157) = 1
      CONSTS(158) = 40.5
      STATES(62) = 0
      CONSTS(159) = 200e3
      CONSTS(160) = 20e3
      CONSTS(161) = 7
      CONSTS(162) = 2.1
      CONSTS(163) = 9.1
      CONSTS(164) = 1
      CONSTS(165) = 2
      CONSTS(166) = 2
      CONSTS(167) = 2
      CONSTS(168) = 2
      CONSTS(169) = 34.0
      CONSTS(170) = ( CONSTS(153)*CONSTS(151))/( (1.00000 - CONSTS(151))*CONSTS(152)+CONSTS(151))
      CONSTS(171) =  CONSTS(139)*0.0350000
      CONSTS(172) = ( CONSTS(154)*CONSTS(155)*(1.00000 - CONSTS(152)))/CONSTS(152)
      CONSTS(173) = ( CONSTS(154)*CONSTS(156)*(1.00000 - CONSTS(151))*CONSTS(152))/CONSTS(151)
      CONSTS(174) =  CONSTS(142)*CONSTS(143)*log( CONSTS(98)*CONSTS(135))
      CONSTS(175) =  CONSTS(169)*0.0350000
      CONSTS(176) = CONSTS(170)
      CONSTS(177) = CONSTS(169)+CONSTS(175)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = TERNRY(VOI.LT.0.0600000,  1.00000*CONSTS(138), TERNRY(VOI.GT.0.500000.AND.VOI.LT.1.10000,  - 0.100000*CONSTS(138),  0.00000*CONSTS(138))
      ALGBRC(27) = ALGBRC(3)
      RATES(59) = ALGBRC(27)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      RATES(62) = ALGBRC(45)
      ALGBRC(31) =  - CONSTS(170)*ALGBRC(27)
      ALGBRC(36) =  CONSTS(172)*STATES(60)
      ALGBRC(42) = - ALGBRC(31) - ALGBRC(36)
      ALGBRC(52) = ALGBRC(42)
      RATES(60) = ALGBRC(52)
      ALGBRC(32) =  - CONSTS(176)*ALGBRC(27)
      ALGBRC(37) =  CONSTS(173)*STATES(61)
      ALGBRC(43) = - ALGBRC(32) - ALGBRC(37)
      ALGBRC(53) = ALGBRC(43)
      RATES(61) = ALGBRC(53)
      ALGBRC(35) =  CONSTS(142)*CONSTS(143)*log( CONSTS(96)*STATES(16))
      ALGBRC(41) =  CONSTS(142)*CONSTS(143)*log( CONSTS(99)*STATES(19))
      ALGBRC(51) =  CONSTS(142)*CONSTS(143)*log( CONSTS(100)*STATES(20))
      ALGBRC(57) =  CONSTS(46)*(EXP((ALGBRC(35)+ALGBRC(41))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(51)/( CONSTS(142)*CONSTS(143))))
      RATES(19) = - ALGBRC(57)
      RATES(20) = ALGBRC(57)
      ALGBRC(1) =  CONSTS(142)*CONSTS(143)*log( CONSTS(93)*STATES(13))
      ALGBRC(26) =  CONSTS(142)*CONSTS(143)*log( CONSTS(94)*STATES(14))
      ALGBRC(30) =  CONSTS(142)*CONSTS(143)*log( CONSTS(95)*STATES(15))
      ALGBRC(56) =  CONSTS(42)*(EXP((ALGBRC(1)+ALGBRC(26))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(30)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(50) =  CONSTS(142)*CONSTS(143)*log( CONSTS(98)*STATES(18))
      ALGBRC(61) =  CONSTS(43)*(EXP(ALGBRC(30)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(50)+ALGBRC(26))/( CONSTS(142)*CONSTS(143))))
      RATES(14) = ALGBRC(61) - ALGBRC(56)
      RATES(15) = ALGBRC(56) - ALGBRC(61)
      RATES(16) = ALGBRC(56) - ALGBRC(61)
      ALGBRC(40) =  CONSTS(142)*CONSTS(143)*log( CONSTS(97)*STATES(17))
      ALGBRC(63) =  CONSTS(44)*(EXP((ALGBRC(50)+ALGBRC(35))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(40)/( CONSTS(142)*CONSTS(143))))
      RATES(18) = ALGBRC(61) - ALGBRC(63)
      ALGBRC(59) =  CONSTS(142)*CONSTS(143)*log( CONSTS(63)*STATES(3))
      ALGBRC(62) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(64) =  CONSTS(142)*CONSTS(143)*log( CONSTS(64)*STATES(4))
      ALGBRC(67) =  CONSTS(2)*(EXP(( CONSTS(132)*ALGBRC(62)+ALGBRC(59))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(64)/( CONSTS(142)*CONSTS(143))))
      RATES(3) = - ALGBRC(67)
      RATES(4) = ALGBRC(67)
      ALGBRC(66) =  CONSTS(45)*(EXP(ALGBRC(40)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(1)+ALGBRC(35))/( CONSTS(142)*CONSTS(143))))
      RATES(13) = ALGBRC(66) - ALGBRC(56)
      RATES(17) = ALGBRC(63) - ALGBRC(66)
      ALGBRC(65) =  CONSTS(142)*CONSTS(143)*log( CONSTS(65)*STATES(8))
      ALGBRC(68) =  CONSTS(142)*CONSTS(143)*log( CONSTS(66)*STATES(9))
      ALGBRC(74) =  CONSTS(3)*(EXP((ALGBRC(65)+ CONSTS(131)*ALGBRC(64))/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(68)+ CONSTS(131)*ALGBRC(64))/( CONSTS(142)*CONSTS(143))))
      RATES(8) = - ALGBRC(74)
      ALGBRC(70) =  CONSTS(142)*CONSTS(143)*log( CONSTS(67)*STATES(10))
      ALGBRC(76) =  CONSTS(4)*(EXP(ALGBRC(68)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(70)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(72) =  CONSTS(142)*CONSTS(143)*log( CONSTS(68)*STATES(11))
      ALGBRC(78) =  CONSTS(5)*(EXP(ALGBRC(70)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(72)/( CONSTS(142)*CONSTS(143))))
      RATES(10) = ALGBRC(76) - ALGBRC(78)
      ALGBRC(80) =  CONSTS(6)*(EXP(ALGBRC(72)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(68)/( CONSTS(142)*CONSTS(143))))
      RATES(9) = (ALGBRC(74) - ALGBRC(76))+ALGBRC(80)
      RATES(11) = ALGBRC(78) - ALGBRC(80)
      ALGBRC(81) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))
      ALGBRC(109) = ALGBRC(81)
      ALGBRC(79) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(155) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(156) =  CONSTS(28)*(EXP(ALGBRC(109)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(155)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(75) = STATES(1)/CONSTS(134)
      ALGBRC(91) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(119) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(120) =  CONSTS(10)*(EXP(ALGBRC(91)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(119)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(97) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(131) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(132) =  CONSTS(16)*(EXP(ALGBRC(97)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(131)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(101) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(139) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(140) =  CONSTS(20)*(EXP(ALGBRC(101)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(139)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(158) = ((- ALGBRC(120) - ALGBRC(132)) - ALGBRC(140)) - ALGBRC(156)
      RATES(28) = ALGBRC(158)
      ALGBRC(92) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(121) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(122) =  CONSTS(13)*(EXP(ALGBRC(92)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(121)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(98) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(133) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(134) =  CONSTS(18)*(EXP(ALGBRC(98)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(133)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(102) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(141) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(142) =  CONSTS(22)*(EXP(ALGBRC(102)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(141)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(159) = ((- ALGBRC(122) - ALGBRC(134)) - ALGBRC(142)) - - ALGBRC(156)
      RATES(34) = ALGBRC(159)
      ALGBRC(82) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))
      ALGBRC(110) = ALGBRC(82)
      ALGBRC(157) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(160) =  CONSTS(30)*(EXP(ALGBRC(110)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(157)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(93) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(123) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(124) =  CONSTS(11)*(EXP(ALGBRC(93)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(123)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(105) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(147) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(148) =  CONSTS(24)*(EXP(ALGBRC(105)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(147)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(162) = ((- ALGBRC(124)+ALGBRC(132)) - ALGBRC(160)) - ALGBRC(148)
      RATES(29) = ALGBRC(162)
      ALGBRC(94) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(125) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(126) =  CONSTS(14)*(EXP(ALGBRC(94)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(125)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(106) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(149) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(150) =  CONSTS(26)*(EXP(ALGBRC(106)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(149)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(163) = ((- ALGBRC(126)+ALGBRC(134)) - - ALGBRC(160)) - ALGBRC(150)
      RATES(35) = ALGBRC(163)
      ALGBRC(83) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))
      ALGBRC(111) = ALGBRC(83)
      ALGBRC(161) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(164) =  CONSTS(32)*(EXP(ALGBRC(111)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(161)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(95) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(127) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(128) =  CONSTS(12)*(EXP(ALGBRC(95)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(127)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(166) = ((- ALGBRC(128)+ALGBRC(140)) - ALGBRC(164))+ALGBRC(148)
      RATES(30) = ALGBRC(166)
      ALGBRC(96) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(129) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(130) =  CONSTS(15)*(EXP(ALGBRC(96)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(129)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(167) = ((- ALGBRC(130)+ALGBRC(142)) - - ALGBRC(164))+ALGBRC(150)
      RATES(36) = ALGBRC(167)
      ALGBRC(84) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))
      ALGBRC(112) = ALGBRC(84)
      ALGBRC(165) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(168) =  CONSTS(29)*(EXP(ALGBRC(112)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(165)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(99) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(135) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(136) =  CONSTS(17)*(EXP(ALGBRC(99)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(135)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(103) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(143) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(144) =  CONSTS(21)*(EXP(ALGBRC(103)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(143)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(170) = ((- ALGBRC(144)+ALGBRC(120)) - ALGBRC(136)) - ALGBRC(168)
      RATES(31) = ALGBRC(170)
      ALGBRC(100) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(137) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(138) =  CONSTS(19)*(EXP(ALGBRC(100)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(137)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(104) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(145) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(146) =  CONSTS(23)*(EXP(ALGBRC(104)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(145)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(171) = ((- ALGBRC(146)+ALGBRC(122)) - ALGBRC(138)) - - ALGBRC(168)
      RATES(37) = ALGBRC(171)
      ALGBRC(85) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))
      ALGBRC(113) = ALGBRC(85)
      ALGBRC(169) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(172) =  CONSTS(31)*(EXP(ALGBRC(113)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(169)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(107) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(151) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(152) =  CONSTS(25)*(EXP(ALGBRC(107)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(151)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(174) = ((ALGBRC(136)+ALGBRC(124)) - ALGBRC(172)) - ALGBRC(152)
      RATES(32) = ALGBRC(174)
      ALGBRC(87) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))
      ALGBRC(114) = ALGBRC(87)
      ALGBRC(173) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(175) =  CONSTS(33)*(EXP(ALGBRC(114)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(173)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(177) = ((ALGBRC(144)+ALGBRC(128)) - ALGBRC(175))+ALGBRC(152)
      RATES(33) = ALGBRC(177)
      ALGBRC(86) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))
      ALGBRC(89) = ALGBRC(79)+ CONSTS(148)*CONSTS(144)*ALGBRC(75)+ALGBRC(86)
      ALGBRC(77) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(117) = ALGBRC(77)+ALGBRC(86)
      ALGBRC(115) =  CONSTS(148)*CONSTS(144)*ALGBRC(75)
      ALGBRC(176) = TERNRY(ALGBRC(115).EQ.0.00000,  CONSTS(8)*(EXP(ALGBRC(89)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(117)/( CONSTS(142)*CONSTS(143)))),  ((( CONSTS(8)*ALGBRC(115))/( CONSTS(142)*CONSTS(143)))/(EXP(ALGBRC(115)/( CONSTS(142)*CONSTS(143))) - 1.00000))*(EXP(ALGBRC(89)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(117)/( CONSTS(142)*CONSTS(143)))))
      ALGBRC(108) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(153) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(154) =  CONSTS(27)*(EXP(ALGBRC(108)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(153)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(179) = ((ALGBRC(138)+ALGBRC(126)) - ((- ALGBRC(172)+ALGBRC(176)) - ALGBRC(176))) - ALGBRC(154)
      RATES(38) = ALGBRC(179)
      ALGBRC(88) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))
      ALGBRC(90) = ALGBRC(79)+ CONSTS(148)*CONSTS(144)*ALGBRC(75)+ALGBRC(88)
      ALGBRC(118) = ALGBRC(77)+ALGBRC(88)
      ALGBRC(116) =  CONSTS(148)*CONSTS(144)*ALGBRC(75)
      ALGBRC(178) = TERNRY(ALGBRC(116).EQ.0.00000,  CONSTS(9)*(EXP(ALGBRC(90)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(118)/( CONSTS(142)*CONSTS(143)))),  ((( CONSTS(9)*ALGBRC(116))/( CONSTS(142)*CONSTS(143)))/(EXP(ALGBRC(116)/( CONSTS(142)*CONSTS(143))) - 1.00000))*(EXP(ALGBRC(90)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(118)/( CONSTS(142)*CONSTS(143)))))
      ALGBRC(180) = ((ALGBRC(146)+ALGBRC(130)) - ((- ALGBRC(175)+ALGBRC(178)) - ALGBRC(178)))+ALGBRC(154)
      RATES(39) = ALGBRC(180)
      ALGBRC(189) =  CONSTS(142)*CONSTS(143)*log( CONSTS(83)*STATES(40))
      ALGBRC(195) = ALGBRC(189)
      ALGBRC(188) =  CONSTS(142)*CONSTS(143)*log( CONSTS(89)*STATES(6))
      ALGBRC(190) =  CONSTS(142)*CONSTS(143)*log( CONSTS(84)*STATES(41))
      ALGBRC(201) =  CONSTS(128)*ALGBRC(188)+ALGBRC(190)
      ALGBRC(202) =  CONSTS(34)*(EXP(ALGBRC(195)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(201)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(204) =  CONSTS(128)*ALGBRC(202)
      RATES(6) = ALGBRC(204)
      ALGBRC(186) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(196) = ALGBRC(190)+ALGBRC(186)
      ALGBRC(191) =  CONSTS(142)*CONSTS(143)*log( CONSTS(85)*STATES(42))
      ALGBRC(203) = ALGBRC(191)
      ALGBRC(205) =  CONSTS(35)*(EXP(ALGBRC(196)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(203)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(207) = ALGBRC(202) - ALGBRC(205)
      RATES(41) = ALGBRC(207)
      ALGBRC(197) = ALGBRC(191)
      ALGBRC(192) =  CONSTS(142)*CONSTS(143)*log( CONSTS(86)*STATES(43))
      ALGBRC(206) = ALGBRC(192)
      ALGBRC(208) =  CONSTS(36)*(EXP(ALGBRC(197)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(206)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(210) = ALGBRC(205) - ALGBRC(208)
      RATES(42) = ALGBRC(210)
      ALGBRC(198) = ALGBRC(192)
      ALGBRC(185) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(193) =  CONSTS(142)*CONSTS(143)*log( CONSTS(87)*STATES(44))
      ALGBRC(209) = ALGBRC(193)+ALGBRC(185)
      ALGBRC(211) =  CONSTS(37)*(EXP(ALGBRC(198)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(209)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(213) = - ALGBRC(211)+ALGBRC(208)
      RATES(43) = ALGBRC(213)
      ALGBRC(187) =  CONSTS(142)*CONSTS(143)*log( CONSTS(90)*STATES(7))
      ALGBRC(199) = ALGBRC(193)+ CONSTS(127)*ALGBRC(187)
      ALGBRC(194) =  CONSTS(142)*CONSTS(143)*log( CONSTS(88)*STATES(45))
      ALGBRC(212) = ALGBRC(194)
      ALGBRC(214) =  CONSTS(38)*(EXP(ALGBRC(199)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(212)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(215) =  - CONSTS(127)*ALGBRC(214)
      RATES(7) = ALGBRC(215)
      ALGBRC(216) = - ALGBRC(214)+ALGBRC(211)
      RATES(44) = ALGBRC(216)
      ALGBRC(184) = STATES(1)/CONSTS(134)
      ALGBRC(200) =  CONSTS(142)*CONSTS(143)*log( CONSTS(88)*STATES(45))+ CONSTS(129)*CONSTS(144)*ALGBRC(184)
      ALGBRC(217) =  CONSTS(142)*CONSTS(143)*log( CONSTS(83)*STATES(40))+ CONSTS(130)*CONSTS(144)*ALGBRC(184)
      ALGBRC(218) =  CONSTS(39)*(EXP(ALGBRC(200)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(217)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(219) = ALGBRC(218) - ALGBRC(202)
      RATES(40) = ALGBRC(219)
      ALGBRC(220) = - ALGBRC(218)+ALGBRC(214)
      RATES(45) = ALGBRC(220)
      ALGBRC(182) = ALGBRC(178)+ALGBRC(176)
      ALGBRC(222) = ALGBRC(211)
      ALGBRC(225) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(226) =  CONSTS(142)*CONSTS(143)*log( CONSTS(91)*STATES(22))
      ALGBRC(227) =  CONSTS(142)*CONSTS(143)*log( CONSTS(92)*STATES(23))
      ALGBRC(229) =  CONSTS(41)*(EXP(ALGBRC(227)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(225)+ALGBRC(226))/( CONSTS(142)*CONSTS(143))))
      RATES(12) = ALGBRC(182)+ALGBRC(222)+ALGBRC(229)
      ALGBRC(224) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(228) =  CONSTS(40)*(EXP((ALGBRC(224)+ALGBRC(226))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(227)/( CONSTS(142)*CONSTS(143))))
      RATES(22) = - ALGBRC(228)+ALGBRC(229)
      RATES(23) = ALGBRC(228) - ALGBRC(229)
      ALGBRC(181) =  CONSTS(144)*(((((((((((((((((((((((((((((((((((( CONSTS(125)*ALGBRC(142) -  CONSTS(121)*ALGBRC(142))+ CONSTS(126)*ALGBRC(150)) -  CONSTS(122)*ALGBRC(150))+ CONSTS(124)*ALGBRC(138)) -  CONSTS(120)*ALGBRC(138))+ CONSTS(125)*ALGBRC(146)) -  CONSTS(121)*ALGBRC(146))+ CONSTS(126)*ALGBRC(154)) -  CONSTS(122)*ALGBRC(154))+ CONSTS(124)*ALGBRC(134)) -  CONSTS(120)*ALGBRC(134))+ CONSTS(123)*ALGBRC(130)) -  CONSTS(119)*ALGBRC(130))+ CONSTS(123)*ALGBRC(126)) -  CONSTS(119)*ALGBRC(126))+ CONSTS(123)*ALGBRC(122)) -  CONSTS(119)*ALGBRC(122))+ CONSTS(126)*ALGBRC(152)) -  CONSTS(122)*ALGBRC(152))+ CONSTS(125)*ALGBRC(144)) -  CONSTS(121)*ALGBRC(144))+ CONSTS(124)*ALGBRC(136)) -  CONSTS(120)*ALGBRC(136))+ CONSTS(126)*ALGBRC(148)) -  CONSTS(122)*ALGBRC(148))+ CONSTS(125)*ALGBRC(140)) -  CONSTS(121)*ALGBRC(140))+ CONSTS(124)*ALGBRC(132)) -  CONSTS(120)*ALGBRC(132))+ CONSTS(123)*ALGBRC(128)) -  CONSTS(119)*ALGBRC(128))+ CONSTS(123)*ALGBRC(124)) -  CONSTS(119)*ALGBRC(124))+ CONSTS(123)*ALGBRC(120)) -  CONSTS(119)*ALGBRC(120)) - ( CONSTS(148)*ALGBRC(176)+ CONSTS(148)*ALGBRC(178)))
      ALGBRC(221) =  CONSTS(144)*( - CONSTS(129)*ALGBRC(218)+ CONSTS(130)*ALGBRC(218))
      ALGBRC(231) =  - CONSTS(148)*CONSTS(144)*ALGBRC(229)
      ALGBRC(12) = VOI -  INT(VOI/CONSTS(136))*CONSTS(136)
      ALGBRC(29) = TERNRY(ALGBRC(12).GE.0.300000.AND.ALGBRC(12).LE.0.300000+CONSTS(137), ( 0.00000*0.0300000*CONSTS(134))/CONSTS(137), 0.00000)
      ALGBRC(233) = ALGBRC(181)+ALGBRC(221)+ALGBRC(231)+ALGBRC(29)
      RATES(1) = ALGBRC(233)
      ALGBRC(232) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(234) =  CONSTS(142)*CONSTS(143)*log( CONSTS(101)*STATES(46))
      ALGBRC(237) =  CONSTS(142)*CONSTS(143)*log( CONSTS(104)*STATES(49))
      ALGBRC(239) =  CONSTS(48)*(EXP(ALGBRC(237)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(234)+ CONSTS(165)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(235) =  CONSTS(142)*CONSTS(143)*log( CONSTS(102)*STATES(47))
      ALGBRC(241) =  CONSTS(49)*(EXP((ALGBRC(234)+ CONSTS(164)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(235)/( CONSTS(142)*CONSTS(143))))
      RATES(46) = ALGBRC(239) - ALGBRC(241)
      ALGBRC(236) =  CONSTS(142)*CONSTS(143)*log( CONSTS(103)*STATES(48))
      ALGBRC(243) =  CONSTS(50)*(EXP((ALGBRC(235)+ CONSTS(165)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(236)/( CONSTS(142)*CONSTS(143))))
      RATES(47) = ALGBRC(241) - ALGBRC(243)
      ALGBRC(245) =  CONSTS(51)*(EXP(ALGBRC(236)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(237)+ CONSTS(164)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))))
      RATES(49) = ALGBRC(245) - ALGBRC(239)
      RATES(48) = ALGBRC(243) - ALGBRC(245)
      ALGBRC(69) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(71) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(73) =  CONSTS(7)*(EXP(ALGBRC(69)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(71)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(183) = (( (CONSTS(118))*ALGBRC(168)+ CONSTS(118)*ALGBRC(172)+ CONSTS(118)*ALGBRC(164)+ CONSTS(118)*ALGBRC(160)+ CONSTS(118)*ALGBRC(156)+ CONSTS(118)*ALGBRC(175)) - ALGBRC(176)) - ALGBRC(178)
      ALGBRC(230) =  CONSTS(142)*CONSTS(143)*log( CONSTS(61)*STATES(21))
      ALGBRC(238) =  CONSTS(47)*EXP(ALGBRC(237)/( CONSTS(142)*CONSTS(143)))*(EXP(ALGBRC(230)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(232)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(247) = (( CONSTS(165)*ALGBRC(239) -  CONSTS(164)*ALGBRC(241)) -  CONSTS(165)*ALGBRC(243))+ CONSTS(164)*ALGBRC(245)
      RATES(2) = - ALGBRC(73)+ALGBRC(183)+ALGBRC(238)+ALGBRC(247)
      ALGBRC(255) =  CONSTS(142)*CONSTS(143)*log( CONSTS(113)*STATES(58))
      ALGBRC(260) =  CONSTS(142)*CONSTS(143)*log( CONSTS(98)*STATES(18))
      ALGBRC(261) = ALGBRC(260) - CONSTS(174)
      ALGBRC(262) = ALGBRC(255)+ALGBRC(261)
      ALGBRC(250) =  CONSTS(142)*CONSTS(143)*log( CONSTS(105)*STATES(50))
      ALGBRC(249) =  CONSTS(142)*CONSTS(143)*log( CONSTS(117)*STATES(25))
      ALGBRC(263) = ALGBRC(250)+ALGBRC(249)+ALGBRC(261)
      ALGBRC(280) =  CONSTS(60)*(EXP(ALGBRC(262)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(263)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(282) = ALGBRC(280)
      RATES(25) = ALGBRC(282)
      ALGBRC(248) =  CONSTS(142)*CONSTS(143)*log( CONSTS(115)*STATES(27))
      ALGBRC(264) = ALGBRC(250)+ALGBRC(248)+ALGBRC(261)
      ALGBRC(251) =  CONSTS(142)*CONSTS(143)*log( CONSTS(106)*STATES(51))
      ALGBRC(265) = ALGBRC(251)+ALGBRC(261)
      ALGBRC(281) =  CONSTS(52)*(EXP(ALGBRC(264)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(265)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(284) = - ALGBRC(281)
      RATES(27) = ALGBRC(284)
      ALGBRC(285) = ALGBRC(280) - ALGBRC(281)
      RATES(50) = ALGBRC(285)
      ALGBRC(254) =  CONSTS(142)*CONSTS(143)*log( CONSTS(109)*STATES(54))
      ALGBRC(266) = ALGBRC(254)+ALGBRC(261)
      ALGBRC(246) =  CONSTS(142)*CONSTS(143)*log( CONSTS(116)*STATES(26))
      ALGBRC(256) =  CONSTS(142)*CONSTS(143)*log( CONSTS(110)*STATES(55))
      ALGBRC(267) = ALGBRC(246)+ALGBRC(256)+ALGBRC(261)
      ALGBRC(283) =  CONSTS(56)*(EXP(ALGBRC(266)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(267)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(287) = ALGBRC(283)
      RATES(26) = ALGBRC(287)
      ALGBRC(244) =  CONSTS(142)*CONSTS(143)*log( CONSTS(114)*STATES(24))
      ALGBRC(268) = ALGBRC(251)+ALGBRC(244)+ALGBRC(261)
      ALGBRC(252) =  CONSTS(142)*CONSTS(143)*log( CONSTS(107)*STATES(52))
      ALGBRC(269) = ALGBRC(252)+ALGBRC(261)
      ALGBRC(286) =  CONSTS(54)*(EXP(ALGBRC(268)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(269)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(289) = ALGBRC(286)
      RATES(52) = ALGBRC(289)
      ALGBRC(39) =  ( CONSTS(142)*CONSTS(143))*log( CONSTS(61)*STATES(21))
      ALGBRC(49) =  ( CONSTS(142)*CONSTS(143))*log( CONSTS(62)*STATES(5))
      ALGBRC(55) =  CONSTS(1)*(EXP(ALGBRC(39)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(49)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(223) = - ALGBRC(205)
      ALGBRC(240) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(270) = ALGBRC(251)+ CONSTS(166)*ALGBRC(240)+ALGBRC(261)
      ALGBRC(253) =  CONSTS(142)*CONSTS(143)*log( CONSTS(108)*STATES(53))
      ALGBRC(271) = ALGBRC(253)+ALGBRC(261)
      ALGBRC(288) =  CONSTS(53)*(EXP(ALGBRC(270)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(271)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(290) =  - CONSTS(166)*ALGBRC(288)
      ALGBRC(34) = TERNRY(ALGBRC(12).GE.0.300000.AND.ALGBRC(12).LE.0.300000+CONSTS(137),  0.00000*8000.00, 0.00000)
      RATES(5) = (( - ALGBRC(67)*CONSTS(132)+ALGBRC(290)+ALGBRC(55)+ALGBRC(73)+ALGBRC(223)) - ALGBRC(228))+ALGBRC(34)
      ALGBRC(292) = ALGBRC(281) - ALGBRC(288)
      RATES(51) = ALGBRC(292)
      ALGBRC(272) = ALGBRC(253)+ALGBRC(261)
      ALGBRC(273) = ALGBRC(254)+ CONSTS(168)*ALGBRC(244)+ALGBRC(261)
      ALGBRC(291) =  CONSTS(55)*(EXP(ALGBRC(272)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(273)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(294) = ALGBRC(288) - ALGBRC(291)
      RATES(53) = ALGBRC(294)
      ALGBRC(295) = ALGBRC(291) - ALGBRC(283)
      RATES(54) = ALGBRC(295)
      ALGBRC(274) = ALGBRC(256)+ALGBRC(261)
      ALGBRC(242) =  CONSTS(142)*CONSTS(143)*log( CONSTS(61)*STATES(21))
      ALGBRC(258) =  CONSTS(142)*CONSTS(143)*log( CONSTS(111)*STATES(56))
      ALGBRC(275) = ALGBRC(258)+ CONSTS(167)*ALGBRC(242)+ALGBRC(261)
      ALGBRC(293) =  CONSTS(57)*(EXP(ALGBRC(274)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(275)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(296) =  CONSTS(167)*ALGBRC(293)
      RATES(21) = (- ALGBRC(55) - ALGBRC(238))+ALGBRC(296)
      ALGBRC(298) = ALGBRC(283) - ALGBRC(293)
      RATES(55) = ALGBRC(298)
      ALGBRC(276) = ALGBRC(258)+ CONSTS(168)*ALGBRC(244)+ALGBRC(261)
      ALGBRC(259) =  CONSTS(142)*CONSTS(143)*log( CONSTS(112)*STATES(57))
      ALGBRC(277) = ALGBRC(259)+ALGBRC(261)
      ALGBRC(297) =  CONSTS(58)*(EXP(ALGBRC(276)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(277)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(300) = ALGBRC(293) - ALGBRC(297)
      RATES(56) = ALGBRC(300)
      ALGBRC(278) = ALGBRC(259)+ALGBRC(261)
      ALGBRC(279) = ALGBRC(244)+ALGBRC(255)+ALGBRC(261)
      ALGBRC(299) =  CONSTS(59)*(EXP(ALGBRC(278)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(279)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(301) = (( CONSTS(168)*ALGBRC(291) - ALGBRC(286)) -  CONSTS(168)*ALGBRC(297))+ALGBRC(299)
      RATES(24) = ALGBRC(301)
      ALGBRC(303) = ALGBRC(299) - ALGBRC(280)
      RATES(58) = ALGBRC(303)
      ALGBRC(302) = ALGBRC(297) - ALGBRC(299)
      RATES(57) = ALGBRC(302)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = TERNRY(VOI.LT.0.0600000,  1.00000*CONSTS(138), TERNRY(VOI.GT.0.500000.AND.VOI.LT.1.10000,  - 0.100000*CONSTS(138),  0.00000*CONSTS(138))
      ALGBRC(27) = ALGBRC(3)
      ALGBRC(31) =  - CONSTS(170)*ALGBRC(27)
      ALGBRC(36) =  CONSTS(172)*STATES(60)
      ALGBRC(42) = - ALGBRC(31) - ALGBRC(36)
      ALGBRC(52) = ALGBRC(42)
      ALGBRC(32) =  - CONSTS(176)*ALGBRC(27)
      ALGBRC(37) =  CONSTS(173)*STATES(61)
      ALGBRC(43) = - ALGBRC(32) - ALGBRC(37)
      ALGBRC(53) = ALGBRC(43)
      ALGBRC(35) =  CONSTS(142)*CONSTS(143)*log( CONSTS(96)*STATES(16))
      ALGBRC(41) =  CONSTS(142)*CONSTS(143)*log( CONSTS(99)*STATES(19))
      ALGBRC(51) =  CONSTS(142)*CONSTS(143)*log( CONSTS(100)*STATES(20))
      ALGBRC(57) =  CONSTS(46)*(EXP((ALGBRC(35)+ALGBRC(41))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(51)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(1) =  CONSTS(142)*CONSTS(143)*log( CONSTS(93)*STATES(13))
      ALGBRC(26) =  CONSTS(142)*CONSTS(143)*log( CONSTS(94)*STATES(14))
      ALGBRC(30) =  CONSTS(142)*CONSTS(143)*log( CONSTS(95)*STATES(15))
      ALGBRC(56) =  CONSTS(42)*(EXP((ALGBRC(1)+ALGBRC(26))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(30)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(50) =  CONSTS(142)*CONSTS(143)*log( CONSTS(98)*STATES(18))
      ALGBRC(61) =  CONSTS(43)*(EXP(ALGBRC(30)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(50)+ALGBRC(26))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(40) =  CONSTS(142)*CONSTS(143)*log( CONSTS(97)*STATES(17))
      ALGBRC(63) =  CONSTS(44)*(EXP((ALGBRC(50)+ALGBRC(35))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(40)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(59) =  CONSTS(142)*CONSTS(143)*log( CONSTS(63)*STATES(3))
      ALGBRC(62) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(64) =  CONSTS(142)*CONSTS(143)*log( CONSTS(64)*STATES(4))
      ALGBRC(67) =  CONSTS(2)*(EXP(( CONSTS(132)*ALGBRC(62)+ALGBRC(59))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(64)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(66) =  CONSTS(45)*(EXP(ALGBRC(40)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(1)+ALGBRC(35))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(65) =  CONSTS(142)*CONSTS(143)*log( CONSTS(65)*STATES(8))
      ALGBRC(68) =  CONSTS(142)*CONSTS(143)*log( CONSTS(66)*STATES(9))
      ALGBRC(74) =  CONSTS(3)*(EXP((ALGBRC(65)+ CONSTS(131)*ALGBRC(64))/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(68)+ CONSTS(131)*ALGBRC(64))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(70) =  CONSTS(142)*CONSTS(143)*log( CONSTS(67)*STATES(10))
      ALGBRC(76) =  CONSTS(4)*(EXP(ALGBRC(68)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(70)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(72) =  CONSTS(142)*CONSTS(143)*log( CONSTS(68)*STATES(11))
      ALGBRC(78) =  CONSTS(5)*(EXP(ALGBRC(70)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(72)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(80) =  CONSTS(6)*(EXP(ALGBRC(72)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(68)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(81) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))
      ALGBRC(109) = ALGBRC(81)
      ALGBRC(79) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(155) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(156) =  CONSTS(28)*(EXP(ALGBRC(109)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(155)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(75) = STATES(1)/CONSTS(134)
      ALGBRC(91) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(119) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(120) =  CONSTS(10)*(EXP(ALGBRC(91)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(119)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(97) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(131) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(132) =  CONSTS(16)*(EXP(ALGBRC(97)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(131)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(101) =  CONSTS(142)*CONSTS(143)*log( CONSTS(71)*STATES(28))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(139) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(140) =  CONSTS(20)*(EXP(ALGBRC(101)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(139)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(158) = ((- ALGBRC(120) - ALGBRC(132)) - ALGBRC(140)) - ALGBRC(156)
      ALGBRC(92) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(121) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(122) =  CONSTS(13)*(EXP(ALGBRC(92)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(121)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(98) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(133) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(134) =  CONSTS(18)*(EXP(ALGBRC(98)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(133)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(102) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(141) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(142) =  CONSTS(22)*(EXP(ALGBRC(102)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(141)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(159) = ((- ALGBRC(122) - ALGBRC(134)) - ALGBRC(142)) - - ALGBRC(156)
      ALGBRC(82) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))
      ALGBRC(110) = ALGBRC(82)
      ALGBRC(157) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(160) =  CONSTS(30)*(EXP(ALGBRC(110)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(157)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(93) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(123) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(124) =  CONSTS(11)*(EXP(ALGBRC(93)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(123)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(105) =  CONSTS(142)*CONSTS(143)*log( CONSTS(72)*STATES(29))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(147) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(148) =  CONSTS(24)*(EXP(ALGBRC(105)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(147)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(162) = ((- ALGBRC(124)+ALGBRC(132)) - ALGBRC(160)) - ALGBRC(148)
      ALGBRC(94) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(125) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(126) =  CONSTS(14)*(EXP(ALGBRC(94)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(125)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(106) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(149) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(150) =  CONSTS(26)*(EXP(ALGBRC(106)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(149)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(163) = ((- ALGBRC(126)+ALGBRC(134)) - - ALGBRC(160)) - ALGBRC(150)
      ALGBRC(83) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))
      ALGBRC(111) = ALGBRC(83)
      ALGBRC(161) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(164) =  CONSTS(32)*(EXP(ALGBRC(111)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(161)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(95) =  CONSTS(142)*CONSTS(143)*log( CONSTS(73)*STATES(30))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(127) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(128) =  CONSTS(12)*(EXP(ALGBRC(95)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(127)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(166) = ((- ALGBRC(128)+ALGBRC(140)) - ALGBRC(164))+ALGBRC(148)
      ALGBRC(96) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))+ CONSTS(119)*CONSTS(144)*ALGBRC(75)
      ALGBRC(129) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(123)*CONSTS(144)*ALGBRC(75)
      ALGBRC(130) =  CONSTS(15)*(EXP(ALGBRC(96)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(129)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(167) = ((- ALGBRC(130)+ALGBRC(142)) - - ALGBRC(164))+ALGBRC(150)
      ALGBRC(84) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))
      ALGBRC(112) = ALGBRC(84)
      ALGBRC(165) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(168) =  CONSTS(29)*(EXP(ALGBRC(112)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(165)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(99) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(135) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(136) =  CONSTS(17)*(EXP(ALGBRC(99)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(135)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(103) =  CONSTS(142)*CONSTS(143)*log( CONSTS(74)*STATES(31))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(143) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(144) =  CONSTS(21)*(EXP(ALGBRC(103)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(143)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(170) = ((- ALGBRC(144)+ALGBRC(120)) - ALGBRC(136)) - ALGBRC(168)
      ALGBRC(100) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(120)*CONSTS(144)*ALGBRC(75)
      ALGBRC(137) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(124)*CONSTS(144)*ALGBRC(75)
      ALGBRC(138) =  CONSTS(19)*(EXP(ALGBRC(100)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(137)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(104) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))+ CONSTS(121)*CONSTS(144)*ALGBRC(75)
      ALGBRC(145) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(125)*CONSTS(144)*ALGBRC(75)
      ALGBRC(146) =  CONSTS(23)*(EXP(ALGBRC(104)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(145)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(171) = ((- ALGBRC(146)+ALGBRC(122)) - ALGBRC(138)) - - ALGBRC(168)
      ALGBRC(85) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))
      ALGBRC(113) = ALGBRC(85)
      ALGBRC(169) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(172) =  CONSTS(31)*(EXP(ALGBRC(113)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(169)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(107) =  CONSTS(142)*CONSTS(143)*log( CONSTS(75)*STATES(32))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(151) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(152) =  CONSTS(25)*(EXP(ALGBRC(107)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(151)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(174) = ((ALGBRC(136)+ALGBRC(124)) - ALGBRC(172)) - ALGBRC(152)
      ALGBRC(87) =  CONSTS(142)*CONSTS(143)*log( CONSTS(76)*STATES(33))
      ALGBRC(114) = ALGBRC(87)
      ALGBRC(173) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(118)*ALGBRC(79)
      ALGBRC(175) =  CONSTS(33)*(EXP(ALGBRC(114)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(173)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(177) = ((ALGBRC(144)+ALGBRC(128)) - ALGBRC(175))+ALGBRC(152)
      ALGBRC(86) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))
      ALGBRC(89) = ALGBRC(79)+ CONSTS(148)*CONSTS(144)*ALGBRC(75)+ALGBRC(86)
      ALGBRC(77) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(117) = ALGBRC(77)+ALGBRC(86)
      ALGBRC(115) =  CONSTS(148)*CONSTS(144)*ALGBRC(75)
      ALGBRC(176) = TERNRY(ALGBRC(115).EQ.0.00000,  CONSTS(8)*(EXP(ALGBRC(89)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(117)/( CONSTS(142)*CONSTS(143)))),  ((( CONSTS(8)*ALGBRC(115))/( CONSTS(142)*CONSTS(143)))/(EXP(ALGBRC(115)/( CONSTS(142)*CONSTS(143))) - 1.00000))*(EXP(ALGBRC(89)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(117)/( CONSTS(142)*CONSTS(143)))))
      ALGBRC(108) =  CONSTS(142)*CONSTS(143)*log( CONSTS(81)*STATES(38))+ CONSTS(122)*CONSTS(144)*ALGBRC(75)
      ALGBRC(153) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))+ CONSTS(126)*CONSTS(144)*ALGBRC(75)
      ALGBRC(154) =  CONSTS(27)*(EXP(ALGBRC(108)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(153)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(179) = ((ALGBRC(138)+ALGBRC(126)) - ((- ALGBRC(172)+ALGBRC(176)) - ALGBRC(176))) - ALGBRC(154)
      ALGBRC(88) =  CONSTS(142)*CONSTS(143)*log( CONSTS(82)*STATES(39))
      ALGBRC(90) = ALGBRC(79)+ CONSTS(148)*CONSTS(144)*ALGBRC(75)+ALGBRC(88)
      ALGBRC(118) = ALGBRC(77)+ALGBRC(88)
      ALGBRC(116) =  CONSTS(148)*CONSTS(144)*ALGBRC(75)
      ALGBRC(178) = TERNRY(ALGBRC(116).EQ.0.00000,  CONSTS(9)*(EXP(ALGBRC(90)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(118)/( CONSTS(142)*CONSTS(143)))),  ((( CONSTS(9)*ALGBRC(116))/( CONSTS(142)*CONSTS(143)))/(EXP(ALGBRC(116)/( CONSTS(142)*CONSTS(143))) - 1.00000))*(EXP(ALGBRC(90)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(118)/( CONSTS(142)*CONSTS(143)))))
      ALGBRC(180) = ((ALGBRC(146)+ALGBRC(130)) - ((- ALGBRC(175)+ALGBRC(178)) - ALGBRC(178)))+ALGBRC(154)
      ALGBRC(189) =  CONSTS(142)*CONSTS(143)*log( CONSTS(83)*STATES(40))
      ALGBRC(195) = ALGBRC(189)
      ALGBRC(188) =  CONSTS(142)*CONSTS(143)*log( CONSTS(89)*STATES(6))
      ALGBRC(190) =  CONSTS(142)*CONSTS(143)*log( CONSTS(84)*STATES(41))
      ALGBRC(201) =  CONSTS(128)*ALGBRC(188)+ALGBRC(190)
      ALGBRC(202) =  CONSTS(34)*(EXP(ALGBRC(195)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(201)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(204) =  CONSTS(128)*ALGBRC(202)
      ALGBRC(186) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(196) = ALGBRC(190)+ALGBRC(186)
      ALGBRC(191) =  CONSTS(142)*CONSTS(143)*log( CONSTS(85)*STATES(42))
      ALGBRC(203) = ALGBRC(191)
      ALGBRC(205) =  CONSTS(35)*(EXP(ALGBRC(196)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(203)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(207) = ALGBRC(202) - ALGBRC(205)
      ALGBRC(197) = ALGBRC(191)
      ALGBRC(192) =  CONSTS(142)*CONSTS(143)*log( CONSTS(86)*STATES(43))
      ALGBRC(206) = ALGBRC(192)
      ALGBRC(208) =  CONSTS(36)*(EXP(ALGBRC(197)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(206)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(210) = ALGBRC(205) - ALGBRC(208)
      ALGBRC(198) = ALGBRC(192)
      ALGBRC(185) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(193) =  CONSTS(142)*CONSTS(143)*log( CONSTS(87)*STATES(44))
      ALGBRC(209) = ALGBRC(193)+ALGBRC(185)
      ALGBRC(211) =  CONSTS(37)*(EXP(ALGBRC(198)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(209)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(213) = - ALGBRC(211)+ALGBRC(208)
      ALGBRC(187) =  CONSTS(142)*CONSTS(143)*log( CONSTS(90)*STATES(7))
      ALGBRC(199) = ALGBRC(193)+ CONSTS(127)*ALGBRC(187)
      ALGBRC(194) =  CONSTS(142)*CONSTS(143)*log( CONSTS(88)*STATES(45))
      ALGBRC(212) = ALGBRC(194)
      ALGBRC(214) =  CONSTS(38)*(EXP(ALGBRC(199)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(212)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(215) =  - CONSTS(127)*ALGBRC(214)
      ALGBRC(216) = - ALGBRC(214)+ALGBRC(211)
      ALGBRC(184) = STATES(1)/CONSTS(134)
      ALGBRC(200) =  CONSTS(142)*CONSTS(143)*log( CONSTS(88)*STATES(45))+ CONSTS(129)*CONSTS(144)*ALGBRC(184)
      ALGBRC(217) =  CONSTS(142)*CONSTS(143)*log( CONSTS(83)*STATES(40))+ CONSTS(130)*CONSTS(144)*ALGBRC(184)
      ALGBRC(218) =  CONSTS(39)*(EXP(ALGBRC(200)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(217)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(219) = ALGBRC(218) - ALGBRC(202)
      ALGBRC(220) = - ALGBRC(218)+ALGBRC(214)
      ALGBRC(182) = ALGBRC(178)+ALGBRC(176)
      ALGBRC(222) = ALGBRC(211)
      ALGBRC(225) =  CONSTS(142)*CONSTS(143)*log( CONSTS(70)*STATES(12))
      ALGBRC(226) =  CONSTS(142)*CONSTS(143)*log( CONSTS(91)*STATES(22))
      ALGBRC(227) =  CONSTS(142)*CONSTS(143)*log( CONSTS(92)*STATES(23))
      ALGBRC(229) =  CONSTS(41)*(EXP(ALGBRC(227)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(225)+ALGBRC(226))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(224) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(228) =  CONSTS(40)*(EXP((ALGBRC(224)+ALGBRC(226))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(227)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(181) =  CONSTS(144)*(((((((((((((((((((((((((((((((((((( CONSTS(125)*ALGBRC(142) -  CONSTS(121)*ALGBRC(142))+ CONSTS(126)*ALGBRC(150)) -  CONSTS(122)*ALGBRC(150))+ CONSTS(124)*ALGBRC(138)) -  CONSTS(120)*ALGBRC(138))+ CONSTS(125)*ALGBRC(146)) -  CONSTS(121)*ALGBRC(146))+ CONSTS(126)*ALGBRC(154)) -  CONSTS(122)*ALGBRC(154))+ CONSTS(124)*ALGBRC(134)) -  CONSTS(120)*ALGBRC(134))+ CONSTS(123)*ALGBRC(130)) -  CONSTS(119)*ALGBRC(130))+ CONSTS(123)*ALGBRC(126)) -  CONSTS(119)*ALGBRC(126))+ CONSTS(123)*ALGBRC(122)) -  CONSTS(119)*ALGBRC(122))+ CONSTS(126)*ALGBRC(152)) -  CONSTS(122)*ALGBRC(152))+ CONSTS(125)*ALGBRC(144)) -  CONSTS(121)*ALGBRC(144))+ CONSTS(124)*ALGBRC(136)) -  CONSTS(120)*ALGBRC(136))+ CONSTS(126)*ALGBRC(148)) -  CONSTS(122)*ALGBRC(148))+ CONSTS(125)*ALGBRC(140)) -  CONSTS(121)*ALGBRC(140))+ CONSTS(124)*ALGBRC(132)) -  CONSTS(120)*ALGBRC(132))+ CONSTS(123)*ALGBRC(128)) -  CONSTS(119)*ALGBRC(128))+ CONSTS(123)*ALGBRC(124)) -  CONSTS(119)*ALGBRC(124))+ CONSTS(123)*ALGBRC(120)) -  CONSTS(119)*ALGBRC(120)) - ( CONSTS(148)*ALGBRC(176)+ CONSTS(148)*ALGBRC(178)))
      ALGBRC(221) =  CONSTS(144)*( - CONSTS(129)*ALGBRC(218)+ CONSTS(130)*ALGBRC(218))
      ALGBRC(231) =  - CONSTS(148)*CONSTS(144)*ALGBRC(229)
      ALGBRC(12) = VOI -  INT(VOI/CONSTS(136))*CONSTS(136)
      ALGBRC(29) = TERNRY(ALGBRC(12).GE.0.300000.AND.ALGBRC(12).LE.0.300000+CONSTS(137), ( 0.00000*0.0300000*CONSTS(134))/CONSTS(137), 0.00000)
      ALGBRC(233) = ALGBRC(181)+ALGBRC(221)+ALGBRC(231)+ALGBRC(29)
      ALGBRC(232) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(234) =  CONSTS(142)*CONSTS(143)*log( CONSTS(101)*STATES(46))
      ALGBRC(237) =  CONSTS(142)*CONSTS(143)*log( CONSTS(104)*STATES(49))
      ALGBRC(239) =  CONSTS(48)*(EXP(ALGBRC(237)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(234)+ CONSTS(165)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(235) =  CONSTS(142)*CONSTS(143)*log( CONSTS(102)*STATES(47))
      ALGBRC(241) =  CONSTS(49)*(EXP((ALGBRC(234)+ CONSTS(164)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(235)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(236) =  CONSTS(142)*CONSTS(143)*log( CONSTS(103)*STATES(48))
      ALGBRC(243) =  CONSTS(50)*(EXP((ALGBRC(235)+ CONSTS(165)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(236)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(245) =  CONSTS(51)*(EXP(ALGBRC(236)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(237)+ CONSTS(164)*ALGBRC(232))/( CONSTS(142)*CONSTS(143))))
      ALGBRC(69) =  CONSTS(142)*CONSTS(143)*log( CONSTS(69)*STATES(2))
      ALGBRC(71) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(73) =  CONSTS(7)*(EXP(ALGBRC(69)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(71)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(183) = (( (CONSTS(118))*ALGBRC(168)+ CONSTS(118)*ALGBRC(172)+ CONSTS(118)*ALGBRC(164)+ CONSTS(118)*ALGBRC(160)+ CONSTS(118)*ALGBRC(156)+ CONSTS(118)*ALGBRC(175)) - ALGBRC(176)) - ALGBRC(178)
      ALGBRC(230) =  CONSTS(142)*CONSTS(143)*log( CONSTS(61)*STATES(21))
      ALGBRC(238) =  CONSTS(47)*EXP(ALGBRC(237)/( CONSTS(142)*CONSTS(143)))*(EXP(ALGBRC(230)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(232)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(247) = (( CONSTS(165)*ALGBRC(239) -  CONSTS(164)*ALGBRC(241)) -  CONSTS(165)*ALGBRC(243))+ CONSTS(164)*ALGBRC(245)
      ALGBRC(255) =  CONSTS(142)*CONSTS(143)*log( CONSTS(113)*STATES(58))
      ALGBRC(260) =  CONSTS(142)*CONSTS(143)*log( CONSTS(98)*STATES(18))
      ALGBRC(261) = ALGBRC(260) - CONSTS(174)
      ALGBRC(262) = ALGBRC(255)+ALGBRC(261)
      ALGBRC(250) =  CONSTS(142)*CONSTS(143)*log( CONSTS(105)*STATES(50))
      ALGBRC(249) =  CONSTS(142)*CONSTS(143)*log( CONSTS(117)*STATES(25))
      ALGBRC(263) = ALGBRC(250)+ALGBRC(249)+ALGBRC(261)
      ALGBRC(280) =  CONSTS(60)*(EXP(ALGBRC(262)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(263)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(282) = ALGBRC(280)
      ALGBRC(248) =  CONSTS(142)*CONSTS(143)*log( CONSTS(115)*STATES(27))
      ALGBRC(264) = ALGBRC(250)+ALGBRC(248)+ALGBRC(261)
      ALGBRC(251) =  CONSTS(142)*CONSTS(143)*log( CONSTS(106)*STATES(51))
      ALGBRC(265) = ALGBRC(251)+ALGBRC(261)
      ALGBRC(281) =  CONSTS(52)*(EXP(ALGBRC(264)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(265)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(284) = - ALGBRC(281)
      ALGBRC(285) = ALGBRC(280) - ALGBRC(281)
      ALGBRC(254) =  CONSTS(142)*CONSTS(143)*log( CONSTS(109)*STATES(54))
      ALGBRC(266) = ALGBRC(254)+ALGBRC(261)
      ALGBRC(246) =  CONSTS(142)*CONSTS(143)*log( CONSTS(116)*STATES(26))
      ALGBRC(256) =  CONSTS(142)*CONSTS(143)*log( CONSTS(110)*STATES(55))
      ALGBRC(267) = ALGBRC(246)+ALGBRC(256)+ALGBRC(261)
      ALGBRC(283) =  CONSTS(56)*(EXP(ALGBRC(266)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(267)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(287) = ALGBRC(283)
      ALGBRC(244) =  CONSTS(142)*CONSTS(143)*log( CONSTS(114)*STATES(24))
      ALGBRC(268) = ALGBRC(251)+ALGBRC(244)+ALGBRC(261)
      ALGBRC(252) =  CONSTS(142)*CONSTS(143)*log( CONSTS(107)*STATES(52))
      ALGBRC(269) = ALGBRC(252)+ALGBRC(261)
      ALGBRC(286) =  CONSTS(54)*(EXP(ALGBRC(268)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(269)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(289) = ALGBRC(286)
      ALGBRC(39) =  ( CONSTS(142)*CONSTS(143))*log( CONSTS(61)*STATES(21))
      ALGBRC(49) =  ( CONSTS(142)*CONSTS(143))*log( CONSTS(62)*STATES(5))
      ALGBRC(55) =  CONSTS(1)*(EXP(ALGBRC(39)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(49)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(223) = - ALGBRC(205)
      ALGBRC(240) =  CONSTS(142)*CONSTS(143)*log( CONSTS(62)*STATES(5))
      ALGBRC(270) = ALGBRC(251)+ CONSTS(166)*ALGBRC(240)+ALGBRC(261)
      ALGBRC(253) =  CONSTS(142)*CONSTS(143)*log( CONSTS(108)*STATES(53))
      ALGBRC(271) = ALGBRC(253)+ALGBRC(261)
      ALGBRC(288) =  CONSTS(53)*(EXP(ALGBRC(270)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(271)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(290) =  - CONSTS(166)*ALGBRC(288)
      ALGBRC(34) = TERNRY(ALGBRC(12).GE.0.300000.AND.ALGBRC(12).LE.0.300000+CONSTS(137),  0.00000*8000.00, 0.00000)
      ALGBRC(292) = ALGBRC(281) - ALGBRC(288)
      ALGBRC(272) = ALGBRC(253)+ALGBRC(261)
      ALGBRC(273) = ALGBRC(254)+ CONSTS(168)*ALGBRC(244)+ALGBRC(261)
      ALGBRC(291) =  CONSTS(55)*(EXP(ALGBRC(272)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(273)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(294) = ALGBRC(288) - ALGBRC(291)
      ALGBRC(295) = ALGBRC(291) - ALGBRC(283)
      ALGBRC(274) = ALGBRC(256)+ALGBRC(261)
      ALGBRC(242) =  CONSTS(142)*CONSTS(143)*log( CONSTS(61)*STATES(21))
      ALGBRC(258) =  CONSTS(142)*CONSTS(143)*log( CONSTS(111)*STATES(56))
      ALGBRC(275) = ALGBRC(258)+ CONSTS(167)*ALGBRC(242)+ALGBRC(261)
      ALGBRC(293) =  CONSTS(57)*(EXP(ALGBRC(274)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(275)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(296) =  CONSTS(167)*ALGBRC(293)
      ALGBRC(298) = ALGBRC(283) - ALGBRC(293)
      ALGBRC(276) = ALGBRC(258)+ CONSTS(168)*ALGBRC(244)+ALGBRC(261)
      ALGBRC(259) =  CONSTS(142)*CONSTS(143)*log( CONSTS(112)*STATES(57))
      ALGBRC(277) = ALGBRC(259)+ALGBRC(261)
      ALGBRC(297) =  CONSTS(58)*(EXP(ALGBRC(276)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(277)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(300) = ALGBRC(293) - ALGBRC(297)
      ALGBRC(278) = ALGBRC(259)+ALGBRC(261)
      ALGBRC(279) = ALGBRC(244)+ALGBRC(255)+ALGBRC(261)
      ALGBRC(299) =  CONSTS(59)*(EXP(ALGBRC(278)/( CONSTS(142)*CONSTS(143))) - EXP(ALGBRC(279)/( CONSTS(142)*CONSTS(143))))
      ALGBRC(301) = (( CONSTS(168)*ALGBRC(291) - ALGBRC(286)) -  CONSTS(168)*ALGBRC(297))+ALGBRC(299)
      ALGBRC(303) = ALGBRC(299) - ALGBRC(280)
      ALGBRC(302) = ALGBRC(297) - ALGBRC(299)
      ALGBRC(2) = STATES(1)/CONSTS(134)
      ALGBRC(4) = STATES(5)/CONSTS(133)
      ALGBRC(5) = STATES(12)/CONSTS(140)
      ALGBRC(6) = STATES(21)/CONSTS(171)
      ALGBRC(7) = STATES(2)/CONSTS(141)
      ALGBRC(8) = STATES(6)/CONSTS(133)
      ALGBRC(9) = - arbitrary_log(STATES(24)/CONSTS(133), 10)
      ALGBRC(10) = STATES(5)+STATES(12)+STATES(4)+STATES(21)+STATES(2)+STATES(23)+ 2.00000*(STATES(53)+STATES(54)+STATES(55))+STATES(42)+STATES(43)
      ALGBRC(11) =  (CONSTS(158)/( CONSTS(149)*CONSTS(150)*CONSTS(151)))*STATES(11)*(STATES(61)+CONSTS(149))
      ALGBRC(13) =  CONSTS(142)*CONSTS(143)*log( CONSTS(77)*STATES(34))
      ALGBRC(14) =  CONSTS(142)*CONSTS(143)*log( CONSTS(78)*STATES(35))
      ALGBRC(15) =  CONSTS(142)*CONSTS(143)*log( CONSTS(79)*STATES(36))
      ALGBRC(16) =  CONSTS(142)*CONSTS(143)*log( CONSTS(80)*STATES(37))
      ALGBRC(17) = STATES(12)+STATES(2)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(38)+STATES(39)
      ALGBRC(18) =  CONSTS(165)*STATES(49)+ CONSTS(164)*STATES(47)+ (CONSTS(164)+CONSTS(165))*STATES(48)
      ALGBRC(19) =  CONSTS(166)*STATES(53)+ CONSTS(166)*STATES(54)+ CONSTS(167)*STATES(55)
      ALGBRC(20) = STATES(5)/CONSTS(169)
      ALGBRC(21) = STATES(21)/CONSTS(175)
      ALGBRC(22) = STATES(24)/CONSTS(177)
      ALGBRC(23) = STATES(26)/CONSTS(169)
      ALGBRC(24) = STATES(27)/CONSTS(169)
      ALGBRC(25) = STATES(25)/CONSTS(169)
      ALGBRC(28) =  (CONSTS(158)/( CONSTS(149)*CONSTS(150)*CONSTS(151)))*STATES(10)*STATES(60)
      ALGBRC(33) = ALGBRC(11)+ALGBRC(28)
      ALGBRC(38) =  CONSTS(162)*(EXP(( CONSTS(163)*STATES(59))/CONSTS(149)) - 1.00000)
      ALGBRC(54) =  (( CONSTS(162)*ALGBRC(44))/CONSTS(149))*ALGBRC(45)
      ALGBRC(58) = ALGBRC(38)+ALGBRC(54)
      ALGBRC(60) = ALGBRC(33)+ALGBRC(58)
      ALGBRC(257) =  CONSTS(60)*(EXP(ALGBRC(255)/( CONSTS(142)*CONSTS(143))) - EXP((ALGBRC(250)+ALGBRC(249))/( CONSTS(142)*CONSTS(143))))
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(44) - (TERNRY(ALGBRC(45).GT.0.00000, CONSTS(159), CONSTS(160))))
      minfunc_0[2] = abs(ALGBRC(45) - ALGBRC(46))
      minfunc_0[3] = abs(ALGBRC(47) -  (CONSTS(161)/ALGBRC(44))*STATES(62))
      minfunc_0[4] = abs(ALGBRC(48) -  (CONSTS(161)/ALGBRC(44))*(STATES(59) - CONSTS(149)))
      minfunc_0[5] = abs(ALGBRC(46) - (- ALGBRC(47)+ALGBRC(48)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END