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 233 entries in the algebraic variable array.
C There are a total of 44 entries in each of the rate and state variable arrays.
C There are a total of 186 entries in the constant variable array.
C
C
C CONSTS(1) is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec).
C CONSTS(2) is kappa_R_BU in component BG_parameters (fmol_per_sec).
C CONSTS(3) is kappa_R_UW in component BG_parameters (fmol_per_sec).
C CONSTS(4) is kappa_R_WS in component BG_parameters (fmol_per_sec).
C CONSTS(5) is kappa_R_SU in component BG_parameters (fmol_per_sec).
C CONSTS(6) is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec).
C CONSTS(7) is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec).
C CONSTS(8) is kappa_LCC_K1 in component BG_parameters (fmol_per_sec).
C CONSTS(9) is kappa_LCC_K2 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_PLBph1 in component BG_parameters (fmol_per_sec).
C CONSTS(35) is kappa_PLBph2 in component BG_parameters (fmol_per_sec).
C CONSTS(36) is kappa_PLBd1 in component BG_parameters (fmol_per_sec).
C CONSTS(37) is kappa_PLBd2 in component BG_parameters (fmol_per_sec).
C CONSTS(38) is kappa_Inh in component BG_parameters (fmol_per_sec).
C CONSTS(39) is kappa_RyR in component BG_parameters (fmol_per_sec).
C CONSTS(40) is kappa_OC in component BG_parameters (fmol_per_sec).
C CONSTS(41) is kappa_CCI in component BG_parameters (fmol_per_sec).
C CONSTS(42) is kappa_CII in component BG_parameters (fmol_per_sec).
C CONSTS(43) is kappa_IO in component BG_parameters (fmol_per_sec).
C CONSTS(44) is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec).
C CONSTS(45) is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec).
C CONSTS(46) is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec).
C CONSTS(47) is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec).
C CONSTS(48) is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec).
C CONSTS(49) is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec).
C CONSTS(50) is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec).
C CONSTS(51) is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec).
C CONSTS(52) is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec).
C CONSTS(53) is K_TRPN in component BG_parameters (per_fmol).
C CONSTS(54) is K_Cai in component BG_parameters (per_fmol).
C CONSTS(55) is K_Ca_TRPN in component BG_parameters (per_fmol).
C CONSTS(56) is K_B in component BG_parameters (per_fmol).
C CONSTS(57) is K_U in component BG_parameters (per_fmol).
C CONSTS(58) is K_W in component BG_parameters (per_fmol).
C CONSTS(59) is K_S in component BG_parameters (per_fmol).
C CONSTS(60) is K_Ca_i in component BG_parameters (per_fmol).
C CONSTS(61) is K_Ca_o in component BG_parameters (per_fmol).
C CONSTS(62) is K_K_i in component BG_parameters (per_fmol).
C CONSTS(63) is K_K_o in component BG_parameters (per_fmol).
C CONSTS(64) is K_000_LCC in component BG_parameters (per_fmol).
C CONSTS(65) is K_010_LCC in component BG_parameters (per_fmol).
C CONSTS(66) is K_020_LCC in component BG_parameters (per_fmol).
C CONSTS(67) is K_100_LCC in component BG_parameters (per_fmol).
C CONSTS(68) is K_110_LCC in component BG_parameters (per_fmol).
C CONSTS(69) is K_120_LCC in component BG_parameters (per_fmol).
C CONSTS(70) is K_001_LCC in component BG_parameters (per_fmol).
C CONSTS(71) is K_011_LCC in component BG_parameters (per_fmol).
C CONSTS(72) is K_021_LCC in component BG_parameters (per_fmol).
C CONSTS(73) is K_101_LCC in component BG_parameters (per_fmol).
C CONSTS(74) is K_111_LCC in component BG_parameters (per_fmol).
C CONSTS(75) is K_121_LCC in component BG_parameters (per_fmol).
C CONSTS(76) is K_PLB in component BG_parameters (per_fmol).
C CONSTS(77) is K_PKACI in component BG_parameters (per_fmol).
C CONSTS(78) is K_PLB_PKACI in component BG_parameters (per_fmol).
C CONSTS(79) is K_PP1 in component BG_parameters (per_fmol).
C CONSTS(80) is K_PLBp_PP1 in component BG_parameters (per_fmol).
C CONSTS(81) is K_PLBp in component BG_parameters (per_fmol).
C CONSTS(82) is K_Ip in component BG_parameters (per_fmol).
C CONSTS(83) is K_Ip_PP1 in component BG_parameters (per_fmol).
C CONSTS(84) is K_Ca_SR in component BG_parameters (per_fmol).
C CONSTS(85) is K_C_RyR in component BG_parameters (per_fmol).
C CONSTS(86) is K_CI_RyR in component BG_parameters (per_fmol).
C CONSTS(87) is K_I_RyR in component BG_parameters (per_fmol).
C CONSTS(88) is K_O_RyR in component BG_parameters (per_fmol).
C CONSTS(89) is K_P1_SERCA in component BG_parameters (per_fmol).
C CONSTS(90) is K_P2_SERCA in component BG_parameters (per_fmol).
C CONSTS(91) is K_P2a_SERCA in component BG_parameters (per_fmol).
C CONSTS(92) is K_P4_SERCA in component BG_parameters (per_fmol).
C CONSTS(93) is K_P5_SERCA in component BG_parameters (per_fmol).
C CONSTS(94) is K_P6_SERCA in component BG_parameters (per_fmol).
C CONSTS(95) is K_P8_SERCA in component BG_parameters (per_fmol).
C CONSTS(96) is K_P9_SERCA in component BG_parameters (per_fmol).
C CONSTS(97) is K_P10_SERCA in component BG_parameters (per_fmol).
C CONSTS(98) is K_H in component BG_parameters (per_fmol).
C CONSTS(99) is K_MgATP in component BG_parameters (per_fmol).
C CONSTS(100) is K_MgADP in component BG_parameters (per_fmol).
C CONSTS(101) is K_Pi in component BG_parameters (per_fmol).
C CONSTS(102) is z_rCa in component BG_parameters (dimensionless).
C CONSTS(103) is z_fd in component BG_parameters (dimensionless).
C CONSTS(104) is z_ff1 in component BG_parameters (dimensionless).
C CONSTS(105) is z_ff2 in component BG_parameters (dimensionless).
C CONSTS(106) is z_ff3 in component BG_parameters (dimensionless).
C CONSTS(107) is z_rd in component BG_parameters (dimensionless).
C CONSTS(108) is z_rf1 in component BG_parameters (dimensionless).
C CONSTS(109) is z_rf2 in component BG_parameters (dimensionless).
C CONSTS(110) is z_rf3 in component BG_parameters (dimensionless).
C VOI is time in component environment (second).
C CONSTS(111) is vol_myo in component environment (pL).
C CONSTS(112) is freq in component environment (dimensionless).
C CONSTS(113) is C_m in component environment (fF).
C CONSTS(114) is F in component constants (C_per_mol).
C STATES(1) is q_membrane in component environment (fC).
C ALGBRC(1) is V_m in component environment (volt).
C STATES(2) is q_TRPN in component environment (fmol).
C STATES(3) is q_Ca_TRPN in component environment (fmol).
C STATES(4) is q_Cai in component environment (fmol).
C STATES(5) is q_B in component environment (fmol).
C STATES(6) is q_U in component environment (fmol).
C STATES(7) is q_W in component environment (fmol).
C STATES(8) is q_S in component environment (fmol).
C STATES(9) is q_Cao in component environment (fmol).
C STATES(10) is q_Ki in component environment (fmol).
C STATES(11) is q_Ko in component environment (fmol).
C STATES(12) is q_PLB in component environment (fmol).
C STATES(13) is q_PKACI in component environment (fmol).
C STATES(14) is q_PLB_PKACI in component environment (fmol).
C STATES(15) is q_PP1 in component environment (fmol).
C STATES(16) is q_PLBp_PP1 in component environment (fmol).
C STATES(17) is q_PLBp in component environment (fmol).
C STATES(18) is q_Ip in component environment (fmol).
C STATES(19) is q_Ip_PP1 in component environment (fmol).
C STATES(20) is q_Ca_SR in component environment (fmol).
C STATES(21) is q_H in component environment (fmol).
C STATES(22) is q_P in component environment (fmol).
C STATES(23) is q_MgADP in component environment (fmol).
C STATES(24) is q_MgATP in component environment (fmol).
C STATES(25) is q_000_LCC in component environment (fmol).
C STATES(26) is q_010_LCC in component environment (fmol).
C STATES(27) is q_020_LCC in component environment (fmol).
C STATES(28) is q_100_LCC in component environment (fmol).
C STATES(29) is q_110_LCC in component environment (fmol).
C STATES(30) is q_120_LCC in component environment (fmol).
C STATES(31) is q_001_LCC in component environment (fmol).
C STATES(32) is q_011_LCC in component environment (fmol).
C STATES(33) is q_021_LCC in component environment (fmol).
C STATES(34) is q_101_LCC in component environment (fmol).
C STATES(35) is q_111_LCC in component environment (fmol).
C STATES(36) is q_121_LCC in component environment (fmol).
C STATES(37) is q_C_RyR in component environment (fmol).
C STATES(38) is q_CI_RyR in component environment (fmol).
C STATES(39) is q_I_RyR in component environment (fmol).
C STATES(40) is q_O_RyR in component environment (fmol).
C CONSTS(115) is q_P1_SERCA in component environment (fmol).
C CONSTS(116) is q_P2_SERCA in component environment (fmol).
C CONSTS(117) is q_P2a_SERCA in component environment (fmol).
C CONSTS(118) is q_P4_SERCA in component environment (fmol).
C CONSTS(119) is q_P5_SERCA in component environment (fmol).
C CONSTS(120) is q_P6_SERCA in component environment (fmol).
C CONSTS(121) is q_P8_SERCA in component environment (fmol).
C CONSTS(122) is q_P9_SERCA in component environment (fmol).
C CONSTS(123) is q_P10_SERCA in component environment (fmol).
C ALGBRC(233) is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(44) is v_R_BU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(47) is v_R_UW in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(50) is v_R_WS in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(53) is v_R_SU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(230) is T_total in component crossbridge_TRPN (kPa).
C STATES(41) is SL in component crossbridge_TRPN (metre).
C ALGBRC(164) is v_Ca_i_LCC in component LCC (fmol_per_sec).
C ALGBRC(162) is v_Ca_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(168) is v_K_i_LCC in component LCC (fmol_per_sec).
C ALGBRC(169) is v_K_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(170) is I_mem_LCC in component LCC (fA).
C ALGBRC(45) is v_PLBph1 in component PLB (fmol_per_sec).
C ALGBRC(48) is v_PLBph2 in component PLB (fmol_per_sec).
C ALGBRC(51) is v_PLBd1 in component PLB (fmol_per_sec).
C ALGBRC(54) is v_PLBd2 in component PLB (fmol_per_sec).
C ALGBRC(56) is v_Inh in component PLB (fmol_per_sec).
C ALGBRC(178) is v_RyR in component RyR (fmol_per_sec).
C ALGBRC(183) is v_RyRgate_Ca_i in component RyR (fmol_per_sec).
C ALGBRC(210) is v_Ca_i_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(220) is v_Ca_SR_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(227) is v_H_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(205) is v_MgADP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(202) is v_MgATP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(200) is v_Pi_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(10) is I_pulse in component environment (fA).
C CONSTS(124) is pulse_start in component environment (second).
C CONSTS(125) is pulse_end in component environment (second).
C CONSTS(126) is pulseMag in component environment (fA).
C CONSTS(127) is pulseHolding in component environment (fA).
C ALGBRC(173) is sum_I in component environment (fA).
C CONSTS(128) is zCa in component ion_valences (dimensionless).
C CONSTS(129) is zK in component ion_valences (dimensionless).
C ALGBRC(3) is PLB_T in component environment (fmol).
C ALGBRC(4) is PKACI_T in component environment (fmol).
C ALGBRC(5) is Ip_T in component environment (fmol).
C ALGBRC(6) is K_i_T in component environment (fmol).
C ALGBRC(7) is K_o_T in component environment (fmol).
C CONSTS(130) is R in component constants (J_per_K_per_mol).
C CONSTS(131) is T in component constants (kelvin).
C CONSTS(132) is zNa in component ion_valences (dimensionless).
C CONSTS(133) is zCl in component ion_valences (dimensionless).
C CONSTS(134) is n_Tm in component crossbridge_TRPN (dimensionless).
C ALGBRC(25) is mu_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(29) is mu_Cai in component crossbridge_TRPN (J_per_mol).
C ALGBRC(31) is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(32) is mu_B in component crossbridge_TRPN (J_per_mol).
C ALGBRC(35) is mu_U in component crossbridge_TRPN (J_per_mol).
C ALGBRC(38) is mu_W in component crossbridge_TRPN (J_per_mol).
C ALGBRC(41) is mu_S in component crossbridge_TRPN (J_per_mol).
C ALGBRC(231) is tension in component crossbridge_TRPN (N_per_mm2).
C ALGBRC(232) is mu_tension in component crossbridge_TRPN (J_per_mol).
C CONSTS(135) is kf_coeff in component crossbridge_TRPN (dimensionless).
C CONSTS(161) is n in component crossbridge_TRPN (J_per_mol).
C CONSTS(136) is hh in component crossbridge_TRPN (mm2_per_N).
C CONSTS(137) is SL_0 in component crossbridge_TRPN (metre).
C CONSTS(138) is q_MS in component crossbridge_TRPN (fmol).
C CONSTS(139) is r_s in component crossbridge_TRPN (dimensionless).
C CONSTS(140) is r_w in component crossbridge_TRPN (dimensionless).
C CONSTS(141) is A_eff in component crossbridge_TRPN (dimensionless).
C CONSTS(142) is phi in component crossbridge_TRPN (dimensionless).
C CONSTS(143) is k_uw in component crossbridge_TRPN (per_sec).
C CONSTS(144) is k_ws in component crossbridge_TRPN (per_sec).
C STATES(42) is G_w in component crossbridge_TRPN (metre).
C STATES(43) is G_s in component crossbridge_TRPN (metre).
C CONSTS(158) is c_w in component crossbridge_TRPN (per_sec).
C CONSTS(159) is c_s in component crossbridge_TRPN (per_sec).
C CONSTS(160) is A_w in component crossbridge_TRPN (dimensionless).
C CONSTS(164) is A_s in component crossbridge_TRPN (dimensionless).
C ALGBRC(23) is mu_1 in component crossbridge_TRPN (J_per_m).
C CONSTS(165) is mu_2 in component crossbridge_TRPN (J_per_m).
C ALGBRC(8) is mu_3 in component crossbridge_TRPN (J_per_m).
C ALGBRC(24) is mu_4 in component crossbridge_TRPN (J_per_m).
C CONSTS(168) is mu_5 in component crossbridge_TRPN (J_per_m).
C ALGBRC(9) is mu_6 in component crossbridge_TRPN (J_per_m).
C ALGBRC(27) is v_1 in component crossbridge_TRPN (m_per_s).
C CONSTS(145) is v_2 in component crossbridge_TRPN (m_per_s).
C ALGBRC(28) is v_4 in component crossbridge_TRPN (m_per_s).
C CONSTS(146) is v_to_mu in component crossbridge_TRPN (Js_per_m2).
C ALGBRC(43) is T_active in component crossbridge_TRPN (kPa).
C CONSTS(147) is T_ref in component crossbridge_TRPN (kPa).
C ALGBRC(40) is mu_T_a in component crossbridge_TRPN (kPa).
C ALGBRC(34) is mu_T_S in component crossbridge_TRPN (kPa).
C ALGBRC(37) is mu_T_W in component crossbridge_TRPN (kPa).
C STATES(44) is Cdd in component crossbridge_TRPN (metre).
C ALGBRC(213) is eta in component crossbridge_TRPN (per_sec).
C CONSTS(148) is eta_l in component crossbridge_TRPN (per_sec).
C CONSTS(149) is eta_s in component crossbridge_TRPN (per_sec).
C CONSTS(150) is k in component crossbridge_TRPN (dimensionless).
C CONSTS(151) is alpha in component crossbridge_TRPN (kPa).
C ALGBRC(219) is v_Cdd in component crossbridge_TRPN (m_per_s).
C ALGBRC(214) is mu_d in component crossbridge_TRPN (J_per_m).
C ALGBRC(215) is mu_k in component crossbridge_TRPN (J_per_m).
C ALGBRC(216) is mu_f in component crossbridge_TRPN (J_per_m).
C ALGBRC(223) is mu_T_passive in component crossbridge_TRPN (kPa).
C ALGBRC(226) is T_passive in component crossbridge_TRPN (kPa).
C ALGBRC(49) is mu_Ca_o in component LCC (J_per_mol).
C ALGBRC(52) is mu_Ca_i in component LCC (J_per_mol).
C ALGBRC(55) is mu_K_o in component LCC (J_per_mol).
C ALGBRC(57) is mu_K_i in component LCC (J_per_mol).
C ALGBRC(46) is V_mem in component LCC (volt).
C ALGBRC(66) is Af_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(98) is Ar_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(94) is Am_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(159) is v_LCC_Ca1 in component LCC (fmol_per_sec).
C ALGBRC(67) is Af_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(99) is Ar_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(95) is Am_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(161) is v_LCC_Ca2 in component LCC (fmol_per_sec).
C ALGBRC(68) is Af_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(100) is Ar_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(96) is Am_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(163) is v_LCC_K1 in component LCC (fmol_per_sec).
C ALGBRC(69) is Af_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(101) is Ar_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(97) is Am_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(165) is v_LCC_K2 in component LCC (fmol_per_sec).
C ALGBRC(88) is Af_fCa00 in component LCC (J_per_mol).
C ALGBRC(138) is Ar_fCa00 in component LCC (J_per_mol).
C ALGBRC(139) is v_fCa00 in component LCC (fmol_per_sec).
C ALGBRC(89) is Af_fCa01 in component LCC (J_per_mol).
C ALGBRC(140) is Ar_fCa01 in component LCC (J_per_mol).
C ALGBRC(143) is v_fCa01 in component LCC (fmol_per_sec).
C ALGBRC(90) is Af_fCa02 in component LCC (J_per_mol).
C ALGBRC(144) is Ar_fCa02 in component LCC (J_per_mol).
C ALGBRC(147) is v_fCa02 in component LCC (fmol_per_sec).
C ALGBRC(91) is Af_fCa10 in component LCC (J_per_mol).
C ALGBRC(148) is Ar_fCa10 in component LCC (J_per_mol).
C ALGBRC(151) is v_fCa10 in component LCC (fmol_per_sec).
C ALGBRC(92) is Af_fCa11 in component LCC (J_per_mol).
C ALGBRC(152) is Ar_fCa11 in component LCC (J_per_mol).
C ALGBRC(155) is v_fCa11 in component LCC (fmol_per_sec).
C ALGBRC(93) is Af_fCa12 in component LCC (J_per_mol).
C ALGBRC(156) is Ar_fCa12 in component LCC (J_per_mol).
C ALGBRC(158) is v_fCa12 in component LCC (fmol_per_sec).
C ALGBRC(58) is mu_S000_LCC in component LCC (J_per_mol).
C ALGBRC(141) is v_S000_LCC in component LCC (fmol_per_sec).
C ALGBRC(59) is mu_S010_LCC in component LCC (J_per_mol).
C ALGBRC(145) is v_S010_LCC in component LCC (fmol_per_sec).
C ALGBRC(60) is mu_S020_LCC in component LCC (J_per_mol).
C ALGBRC(149) is v_S020_LCC in component LCC (fmol_per_sec).
C ALGBRC(61) is mu_S100_LCC in component LCC (J_per_mol).
C ALGBRC(153) is v_S100_LCC in component LCC (fmol_per_sec).
C ALGBRC(62) is mu_S110_LCC in component LCC (J_per_mol).
C ALGBRC(157) is v_S110_LCC in component LCC (fmol_per_sec).
C ALGBRC(64) is mu_S120_LCC in component LCC (J_per_mol).
C ALGBRC(160) is v_S120_LCC in component LCC (fmol_per_sec).
C ALGBRC(11) is mu_S001_LCC in component LCC (J_per_mol).
C ALGBRC(142) is v_S001_LCC in component LCC (fmol_per_sec).
C ALGBRC(12) is mu_S011_LCC in component LCC (J_per_mol).
C ALGBRC(146) is v_S011_LCC in component LCC (fmol_per_sec).
C ALGBRC(13) is mu_S021_LCC in component LCC (J_per_mol).
C ALGBRC(150) is v_S021_LCC in component LCC (fmol_per_sec).
C ALGBRC(14) is mu_S101_LCC in component LCC (J_per_mol).
C ALGBRC(154) is v_S101_LCC in component LCC (fmol_per_sec).
C ALGBRC(63) is mu_S111_LCC in component LCC (J_per_mol).
C ALGBRC(167) is v_S111_LCC in component LCC (fmol_per_sec).
C ALGBRC(65) is mu_S121_LCC in component LCC (J_per_mol).
C ALGBRC(172) is v_S121_LCC in component LCC (fmol_per_sec).
C ALGBRC(70) is Af_d000 in component LCC (J_per_mol).
C ALGBRC(102) is Ar_d000 in component LCC (J_per_mol).
C ALGBRC(103) is v_d000 in component LCC (fmol_per_sec).
C ALGBRC(72) is Af_d010 in component LCC (J_per_mol).
C ALGBRC(106) is Ar_d010 in component LCC (J_per_mol).
C ALGBRC(107) is v_d010 in component LCC (fmol_per_sec).
C ALGBRC(74) is Af_d020 in component LCC (J_per_mol).
C ALGBRC(110) is Ar_d020 in component LCC (J_per_mol).
C ALGBRC(111) is v_d020 in component LCC (fmol_per_sec).
C ALGBRC(76) is Af_f1_000 in component LCC (J_per_mol).
C ALGBRC(114) is Ar_f1_000 in component LCC (J_per_mol).
C ALGBRC(115) is v_f1_000 in component LCC (fmol_per_sec).
C ALGBRC(80) is Af_f2_000 in component LCC (J_per_mol).
C ALGBRC(122) is Ar_f2_000 in component LCC (J_per_mol).
C ALGBRC(123) is v_f2_000 in component LCC (fmol_per_sec).
C ALGBRC(84) is Af_f3_010 in component LCC (J_per_mol).
C ALGBRC(130) is Ar_f3_010 in component LCC (J_per_mol).
C ALGBRC(131) is v_f3_010 in component LCC (fmol_per_sec).
C ALGBRC(78) is Af_f1_100 in component LCC (J_per_mol).
C ALGBRC(118) is Ar_f1_100 in component LCC (J_per_mol).
C ALGBRC(119) is v_f1_100 in component LCC (fmol_per_sec).
C ALGBRC(82) is Af_f2_100 in component LCC (J_per_mol).
C ALGBRC(126) is Ar_f2_100 in component LCC (J_per_mol).
C ALGBRC(127) is v_f2_100 in component LCC (fmol_per_sec).
C ALGBRC(86) is Af_f3_110 in component LCC (J_per_mol).
C ALGBRC(134) is Ar_f3_110 in component LCC (J_per_mol).
C ALGBRC(135) is v_f3_110 in component LCC (fmol_per_sec).
C ALGBRC(71) is Af_d001 in component LCC (J_per_mol).
C ALGBRC(104) is Ar_d001 in component LCC (J_per_mol).
C ALGBRC(105) is v_d001 in component LCC (fmol_per_sec).
C ALGBRC(73) is Af_d011 in component LCC (J_per_mol).
C ALGBRC(108) is Ar_d011 in component LCC (J_per_mol).
C ALGBRC(109) is v_d011 in component LCC (fmol_per_sec).
C ALGBRC(75) is Af_d021 in component LCC (J_per_mol).
C ALGBRC(112) is Ar_d021 in component LCC (J_per_mol).
C ALGBRC(113) is v_d021 in component LCC (fmol_per_sec).
C ALGBRC(77) is Af_f1_001 in component LCC (J_per_mol).
C ALGBRC(116) is Ar_f1_001 in component LCC (J_per_mol).
C ALGBRC(117) is v_f1_001 in component LCC (fmol_per_sec).
C ALGBRC(81) is Af_f2_001 in component LCC (J_per_mol).
C ALGBRC(124) is Ar_f2_001 in component LCC (J_per_mol).
C ALGBRC(125) is v_f2_001 in component LCC (fmol_per_sec).
C ALGBRC(85) is Af_f3_011 in component LCC (J_per_mol).
C ALGBRC(132) is Ar_f3_011 in component LCC (J_per_mol).
C ALGBRC(133) is v_f3_011 in component LCC (fmol_per_sec).
C ALGBRC(79) is Af_f1_101 in component LCC (J_per_mol).
C ALGBRC(120) is Ar_f1_101 in component LCC (J_per_mol).
C ALGBRC(121) is v_f1_101 in component LCC (fmol_per_sec).
C ALGBRC(83) is Af_f2_101 in component LCC (J_per_mol).
C ALGBRC(128) is Ar_f2_101 in component LCC (J_per_mol).
C ALGBRC(129) is v_f2_101 in component LCC (fmol_per_sec).
C ALGBRC(87) is Af_f3_111 in component LCC (J_per_mol).
C ALGBRC(136) is Ar_f3_111 in component LCC (J_per_mol).
C ALGBRC(137) is v_f3_111 in component LCC (fmol_per_sec).
C ALGBRC(2) is mu_PLB in component PLB (J_per_mol).
C ALGBRC(22) is mu_PKACI in component PLB (J_per_mol).
C ALGBRC(26) is mu_PLB_PKACI in component PLB (J_per_mol).
C ALGBRC(30) is mu_PP1 in component PLB (J_per_mol).
C ALGBRC(33) is mu_PLBp_PP1 in component PLB (J_per_mol).
C ALGBRC(36) is mu_PLBp in component PLB (J_per_mol).
C ALGBRC(39) is mu_Ip in component PLB (J_per_mol).
C ALGBRC(42) is mu_Ip_PP1 in component PLB (J_per_mol).
C CONSTS(152) is nCa_1 in component RyR (dimensionless).
C CONSTS(153) is nCa_2 in component RyR (dimensionless).
C ALGBRC(166) is mu_Ca_SR in component RyR (J_per_mol).
C ALGBRC(171) is mu_Ca_i in component RyR (J_per_mol).
C ALGBRC(15) is q_Ca_gate_complexes in component RyR (fmol).
C ALGBRC(174) is mu_C_RyR in component RyR (J_per_mol).
C ALGBRC(175) is mu_CI_RyR in component RyR (J_per_mol).
C ALGBRC(176) is mu_I_RyR in component RyR (J_per_mol).
C ALGBRC(177) is mu_O_RyR in component RyR (J_per_mol).
C ALGBRC(179) is v_OC in component RyR (fmol_per_sec).
C ALGBRC(180) is v_CCI in component RyR (fmol_per_sec).
C ALGBRC(181) is v_CII in component RyR (fmol_per_sec).
C ALGBRC(182) is v_IO in component RyR (fmol_per_sec).
C CONSTS(154) is n_Ca_i in component SERCA (dimensionless).
C CONSTS(155) is n_Ca_SR in component SERCA (dimensionless).
C CONSTS(156) is n_H in component SERCA (dimensionless).
C CONSTS(162) is q_SERCA_Ca_complexes in component SERCA (fmol).
C ALGBRC(16) is c_Ca_i in component SERCA (mM).
C ALGBRC(17) is c_Ca_SR in component SERCA (mM).
C ALGBRC(18) is c_H in component SERCA (mM).
C ALGBRC(19) is c_MgADP in component SERCA (mM).
C ALGBRC(20) is c_MgATP in component SERCA (mM).
C ALGBRC(21) is c_Pi in component SERCA (mM).
C CONSTS(157) is vol_i in component SERCA (pL).
C CONSTS(163) is vol_sr in component SERCA (pL).
C CONSTS(167) is vol_isr in component SERCA (pL).
C ALGBRC(184) is mu_Ca_i in component SERCA (J_per_mol).
C ALGBRC(185) is mu_Ca_SR in component SERCA (J_per_mol).
C ALGBRC(186) is mu_H in component SERCA (J_per_mol).
C ALGBRC(187) is mu_MgADP in component SERCA (J_per_mol).
C ALGBRC(188) is mu_MgATP in component SERCA (J_per_mol).
C CONSTS(166) is mu_P1 in component SERCA (J_per_mol).
C ALGBRC(204) is v_P1 in component SERCA (fmol_per_sec).
C ALGBRC(189) is mu_Pi in component SERCA (J_per_mol).
C CONSTS(169) is mu_P2 in component SERCA (J_per_mol).
C ALGBRC(211) is v_P2 in component SERCA (fmol_per_sec).
C CONSTS(170) is mu_P2a in component SERCA (J_per_mol).
C ALGBRC(208) is v_P2a in component SERCA (fmol_per_sec).
C CONSTS(171) is mu_P4 in component SERCA (J_per_mol).
C ALGBRC(217) is v_P4 in component SERCA (fmol_per_sec).
C CONSTS(172) is mu_P5 in component SERCA (J_per_mol).
C ALGBRC(218) is v_P5 in component SERCA (fmol_per_sec).
C CONSTS(174) is mu_P6 in component SERCA (J_per_mol).
C ALGBRC(222) is v_P6 in component SERCA (fmol_per_sec).
C CONSTS(175) is mu_P8 in component SERCA (J_per_mol).
C ALGBRC(225) is v_P8 in component SERCA (fmol_per_sec).
C CONSTS(176) is mu_P9 in component SERCA (J_per_mol).
C ALGBRC(228) is v_P9 in component SERCA (fmol_per_sec).
C CONSTS(173) is mu_P10 in component SERCA (J_per_mol).
C ALGBRC(229) is v_P10 in component SERCA (fmol_per_sec).
C ALGBRC(191) is Af_R1_2 in component SERCA (J_per_mol).
C CONSTS(178) is Ar_R1_2 in component SERCA (J_per_mol).
C ALGBRC(201) is v_SERCA_R1_2 in component SERCA (fmol_per_sec).
C CONSTS(179) is Af_R5_6 in component SERCA (J_per_mol).
C ALGBRC(192) is Ar_R5_6 in component SERCA (J_per_mol).
C ALGBRC(203) is v_SERCA_R5_6 in component SERCA (fmol_per_sec).
C ALGBRC(193) is Af_R2_2a in component SERCA (J_per_mol).
C CONSTS(180) is Ar_R2_2a in component SERCA (J_per_mol).
C ALGBRC(206) is v_SERCA_R2_2a in component SERCA (fmol_per_sec).
C ALGBRC(194) is Af_R2_4 in component SERCA (J_per_mol).
C CONSTS(181) is Ar_R2_4 in component SERCA (J_per_mol).
C ALGBRC(207) is v_SERCA_R2_4 in component SERCA (fmol_per_sec).
C CONSTS(182) is Af_R4_5 in component SERCA (J_per_mol).
C ALGBRC(195) is Ar_R4_5 in component SERCA (J_per_mol).
C ALGBRC(209) is v_SERCA_R4_5 in component SERCA (fmol_per_sec).
C CONSTS(183) is Af_R6_8 in component SERCA (J_per_mol).
C ALGBRC(196) is Ar_R6_8 in component SERCA (J_per_mol).
C ALGBRC(212) is v_SERCA_R6_8 in component SERCA (fmol_per_sec).
C ALGBRC(197) is Af_R8_9 in component SERCA (J_per_mol).
C CONSTS(184) is Ar_R8_9 in component SERCA (J_per_mol).
C ALGBRC(221) is v_SERCA_R8_9 in component SERCA (fmol_per_sec).
C CONSTS(185) is Af_R9_10 in component SERCA (J_per_mol).
C ALGBRC(198) is Ar_R9_10 in component SERCA (J_per_mol).
C ALGBRC(224) is v_SERCA_R9_10 in component SERCA (fmol_per_sec).
C CONSTS(177) is Af_R10_1 in component SERCA (J_per_mol).
C ALGBRC(190) is Ar_R10_1 in component SERCA (J_per_mol).
C ALGBRC(199) is v_SERCA_R10_1 in component SERCA (fmol_per_sec).
C RATES(4) is d/dt q_Cai in component environment (fmol).
C RATES(2) is d/dt q_TRPN in component environment (fmol).
C RATES(3) is d/dt q_Ca_TRPN in component environment (fmol).
C RATES(5) is d/dt q_B in component environment (fmol).
C RATES(6) is d/dt q_U in component environment (fmol).
C RATES(7) is d/dt q_W in component environment (fmol).
C RATES(8) is d/dt q_S in component environment (fmol).
C RATES(9) is d/dt q_Cao in component environment (fmol).
C RATES(10) is d/dt q_Ki in component environment (fmol).
C RATES(11) is d/dt q_Ko in component environment (fmol).
C RATES(12) is d/dt q_PLB in component environment (fmol).
C RATES(13) is d/dt q_PKACI in component environment (fmol).
C RATES(14) is d/dt q_PLB_PKACI in component environment (fmol).
C RATES(15) is d/dt q_PP1 in component environment (fmol).
C RATES(16) is d/dt q_PLBp_PP1 in component environment (fmol).
C RATES(17) is d/dt q_PLBp in component environment (fmol).
C RATES(18) is d/dt q_Ip in component environment (fmol).
C RATES(19) is d/dt q_Ip_PP1 in component environment (fmol).
C RATES(20) is d/dt q_Ca_SR in component environment (fmol).
C RATES(21) is d/dt q_H in component environment (fmol).
C RATES(23) is d/dt q_MgADP in component environment (fmol).
C RATES(24) is d/dt q_MgATP in component environment (fmol).
C RATES(22) is d/dt q_P in component environment (fmol).
C RATES(1) is d/dt q_membrane in component environment (fC).
C RATES(42) is d/dt G_w in component crossbridge_TRPN (metre).
C RATES(41) is d/dt SL in component crossbridge_TRPN (metre).
C RATES(43) is d/dt G_s in component crossbridge_TRPN (metre).
C RATES(44) is d/dt Cdd in component crossbridge_TRPN (metre).
C RATES(25) is d/dt q_000_LCC in component environment (fmol).
C RATES(26) is d/dt q_010_LCC in component environment (fmol).
C RATES(27) is d/dt q_020_LCC in component environment (fmol).
C RATES(28) is d/dt q_100_LCC in component environment (fmol).
C RATES(29) is d/dt q_110_LCC in component environment (fmol).
C RATES(30) is d/dt q_120_LCC in component environment (fmol).
C RATES(31) is d/dt q_001_LCC in component environment (fmol).
C RATES(32) is d/dt q_011_LCC in component environment (fmol).
C RATES(33) is d/dt q_021_LCC in component environment (fmol).
C RATES(34) is d/dt q_101_LCC in component environment (fmol).
C RATES(35) is d/dt q_111_LCC in component environment (fmol).
C RATES(36) is d/dt q_121_LCC in component environment (fmol).
C RATES(40) is d/dt q_O_RyR in component environment (fmol).
C RATES(37) is d/dt q_C_RyR in component environment (fmol).
C RATES(38) is d/dt q_CI_RyR in component environment (fmol).
C RATES(39) is d/dt q_I_RyR in component environment (fmol).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 52.0453
CONSTS(2) = 2.90875
CONSTS(3) = 0.332501
CONSTS(4) = 0.00782356
CONSTS(5) = 140824
CONSTS(6) = 0.756518
CONSTS(7) = 0.996063
CONSTS(8) = 0.0247377
CONSTS(9) = 0.0325707
CONSTS(10) = 329.349
CONSTS(11) = 16.6321
CONSTS(12) = 21.8985
CONSTS(13) = 0.141736
CONSTS(14) = 0.00715766
CONSTS(15) = 0.00942408
CONSTS(16) = 31.7244
CONSTS(17) = 157.535
CONSTS(18) = 0.0136527
CONSTS(19) = 0.0677958
CONSTS(20) = 2.06197
CONSTS(21) = 10.2392
CONSTS(22) = 0.000887375
CONSTS(23) = 0.00440648
CONSTS(24) = 4498.81
CONSTS(25) = 22339.9
CONSTS(26) = 1.93608
CONSTS(27) = 9.61406
CONSTS(28) = 48530.3
CONSTS(29) = 240989
CONSTS(30) = 2450.77
CONSTS(31) = 12169.9
CONSTS(32) = 3226.79
CONSTS(33) = 16023.4
CONSTS(34) = 45.5263
CONSTS(35) = 6.55904
CONSTS(36) = 0.386674
CONSTS(37) = 1.21269
CONSTS(38) = 431.435
CONSTS(39) = 2.85549e+06
CONSTS(40) = 8.83262
CONSTS(41) = 0.12618
CONSTS(42) = 883.262
CONSTS(43) = 73.6052
CONSTS(44) = 1.57289e-05
CONSTS(45) = 175018
CONSTS(46) = 2.44679e+06
CONSTS(47) = 175018
CONSTS(48) = 3487.81
CONSTS(49) = 4.39557e+06
CONSTS(50) = 1.01676e+07
CONSTS(51) = 1.01676e+07
CONSTS(52) = 0.00441447
CONSTS(53) = 1.27424
CONSTS(54) = 1.27424
CONSTS(55) = 0.11171
CONSTS(56) = 6.06164
CONSTS(57) = 2.27312
CONSTS(58) = 14.8627
CONSTS(59) = 3.71567e-06
CONSTS(60) = 0.0343245
CONSTS(61) = 0.227859
CONSTS(62) = 0.000825154
CONSTS(63) = 0.00547767
CONSTS(64) = 1.47795
CONSTS(65) = 29.2665
CONSTS(66) = 22.2281
CONSTS(67) = 0.297629
CONSTS(68) = 5.89367
CONSTS(69) = 4.47629
CONSTS(70) = 3434.28
CONSTS(71) = 68005.7
CONSTS(72) = 51650.9
CONSTS(73) = 691.593
CONSTS(74) = 13695
CONSTS(75) = 10401.4
CONSTS(76) = 0.00235741
CONSTS(77) = 0.395191
CONSTS(78) = 0.638527
CONSTS(79) = 0.361989
CONSTS(80) = 0.203757
CONSTS(81) = 0.326014
CONSTS(82) = 5.41093
CONSTS(83) = 0.0673793
CONSTS(84) = 0.980701
CONSTS(85) = 3962.58
CONSTS(86) = 39.6258
CONSTS(87) = 0.06793
CONSTS(88) = 6.793
CONSTS(89) = 38162.5
CONSTS(90) = 128177
CONSTS(91) = 10585.3
CONSTS(92) = 147985
CONSTS(93) = 0.728251
CONSTS(94) = 29565.2
CONSTS(95) = 4226.28
CONSTS(96) = 2547.31
CONSTS(97) = 4643.82
CONSTS(98) = 212.798
CONSTS(99) = 1254.32
CONSTS(100) = 1.89443e-05
CONSTS(101) = 0.0257106
CONSTS(102) = 2
CONSTS(103) = 2.1404
CONSTS(104) = -1.1495
CONSTS(105) = 0.72162
CONSTS(106) = 4.2933
CONSTS(107) = -2.1404
CONSTS(108) = 1.8993
CONSTS(109) = -0.52288
CONSTS(110) = 0
CONSTS(111) = 34.4
CONSTS(112) = 500
CONSTS(113) = 1.381e5
CONSTS(114) = 96485
STATES(1) = -13039
STATES(2) = 2.57
STATES(3) = 1e-18
STATES(4) = 6.82e-3
STATES(5) = 1
STATES(6) = 0
STATES(7) = 0
STATES(8) = 0
STATES(9) = 6.84
STATES(10) = 5.51E+03
STATES(11) = 2.05E+02
STATES(12) = 4.028E+00
STATES(13) = 2.234E-03
STATES(14) = 1e-18
STATES(15) = 3.382E-02
STATES(16) = 1e-18
STATES(17) = 1e-18
STATES(18) = 1.999E-03
STATES(19) = 1e-18
STATES(20) = 0.641
STATES(21) = 0.004028
STATES(22) = 570
STATES(23) = 1.3794
STATES(24) = 3.8
STATES(25) = 1.3798e-05
STATES(26) = 1.3798e-05
STATES(27) = 1.3798e-05
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) = 4.13941e-05
STATES(38) = 4.13941e-05
STATES(39) = 4.13941e-05
STATES(40) = 4.13941e-05
CONSTS(115) = 2.83036e-14
CONSTS(116) = 2.83036e-14
CONSTS(117) = 2.83036e-14
CONSTS(118) = 2.83036e-14
CONSTS(119) = 2.83036e-14
CONSTS(120) = 2.83036e-14
CONSTS(121) = 2.83036e-14
CONSTS(122) = 2.83036e-14
CONSTS(123) = 2.83036e-14
STATES(41) = 1e-6
CONSTS(124) = 4e-2
CONSTS(125) = 4.1e-2
CONSTS(126) = 1e8
CONSTS(127) = 0
CONSTS(128) = 2
CONSTS(129) = 1
CONSTS(130) = 8.31
CONSTS(131) = 310
CONSTS(132) = 1
CONSTS(133) = -1
CONSTS(134) = 1.1
CONSTS(135) = -0.0118
CONSTS(136) = 1
CONSTS(137) = 2e-6
CONSTS(138) = 1e-6
CONSTS(139) = 0.25
CONSTS(140) = 0.5
CONSTS(141) = 25
CONSTS(142) = 2.23
CONSTS(143) = 26
CONSTS(144) = 4
STATES(42) = 1e-6
STATES(43) = 1e-6
CONSTS(145) = 0
CONSTS(146) = 1
CONSTS(147) = 40.5
STATES(44) = 0
CONSTS(148) = 200e3
CONSTS(149) = 200e3
CONSTS(150) = 7
CONSTS(151) = 2.1
CONSTS(152) = 1
CONSTS(153) = 2
CONSTS(154) = 2
CONSTS(155) = 2
CONSTS(156) = 2
CONSTS(157) = 34.0
CONSTS(158) = ( CONSTS(142)*CONSTS(143)*(1.00000 - CONSTS(140)))/CONSTS(140)
CONSTS(159) = ( CONSTS(142)*CONSTS(144)*(1.00000 - CONSTS(139))*CONSTS(140))/CONSTS(139)
CONSTS(160) = ( CONSTS(141)*CONSTS(139))/( (1.00000 - CONSTS(139))*CONSTS(140)+CONSTS(139))
CONSTS(161) = CONSTS(135)*CONSTS(130)*CONSTS(131)
CONSTS(162) = CONSTS(154)*CONSTS(118)+ CONSTS(154)*CONSTS(119)+ CONSTS(155)*CONSTS(120)
CONSTS(163) = CONSTS(157)*0.0350000
CONSTS(185) = CONSTS(145)
CONSTS(164) = CONSTS(160)
CONSTS(165) = - CONSTS(160)*CONSTS(145)*CONSTS(146)
CONSTS(166) = CONSTS(130)*CONSTS(131)*log( CONSTS(89)*CONSTS(115))
CONSTS(167) = CONSTS(157)+CONSTS(163)
CONSTS(168) = - CONSTS(164)*CONSTS(145)*CONSTS(146)
CONSTS(169) = CONSTS(130)*CONSTS(131)*log( CONSTS(90)*CONSTS(116))
CONSTS(170) = CONSTS(130)*CONSTS(131)*log( CONSTS(91)*CONSTS(117))
CONSTS(171) = CONSTS(130)*CONSTS(131)*log( CONSTS(92)*CONSTS(118))
CONSTS(172) = CONSTS(130)*CONSTS(131)*log( CONSTS(93)*CONSTS(119))
CONSTS(173) = CONSTS(130)*CONSTS(131)*log( CONSTS(97)*CONSTS(123))
CONSTS(174) = CONSTS(130)*CONSTS(131)*log( CONSTS(94)*CONSTS(120))
CONSTS(175) = CONSTS(130)*CONSTS(131)*log( CONSTS(95)*CONSTS(121))
CONSTS(176) = CONSTS(130)*CONSTS(131)*log( CONSTS(96)*CONSTS(122))
CONSTS(177) = CONSTS(173)
CONSTS(178) = CONSTS(169)
CONSTS(179) = CONSTS(172)
CONSTS(180) = CONSTS(170)
CONSTS(181) = CONSTS(171)
CONSTS(182) = CONSTS(171)
CONSTS(183) = CONSTS(174)
CONSTS(184) = CONSTS(176)
CONSTS(185) = CONSTS(176)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(41) = CONSTS(185)
ALGBRC(8) = CONSTS(158)*STATES(42)*CONSTS(146)
ALGBRC(23) = - CONSTS(165) - ALGBRC(8)
ALGBRC(27) = ALGBRC(23)/CONSTS(146)
RATES(42) = ALGBRC(27)
ALGBRC(9) = CONSTS(159)*STATES(43)*CONSTS(146)
ALGBRC(24) = - CONSTS(168) - ALGBRC(9)
ALGBRC(28) = ALGBRC(24)/CONSTS(146)
RATES(43) = ALGBRC(28)
ALGBRC(31) = CONSTS(130)*CONSTS(131)*log( CONSTS(55)*STATES(3))
ALGBRC(32) = CONSTS(130)*CONSTS(131)*log( CONSTS(56)*STATES(5))
ALGBRC(35) = CONSTS(130)*CONSTS(131)*log( CONSTS(57)*STATES(6))
ALGBRC(44) = CONSTS(2)*(EXP((ALGBRC(32)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(35)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))))
RATES(5) = - ALGBRC(44)
ALGBRC(2) = CONSTS(130)*CONSTS(131)*log( CONSTS(76)*STATES(12))
ALGBRC(22) = CONSTS(130)*CONSTS(131)*log( CONSTS(77)*STATES(13))
ALGBRC(26) = CONSTS(130)*CONSTS(131)*log( CONSTS(78)*STATES(14))
ALGBRC(45) = CONSTS(34)*(EXP((ALGBRC(2)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))))
ALGBRC(36) = CONSTS(130)*CONSTS(131)*log( CONSTS(81)*STATES(17))
ALGBRC(48) = CONSTS(35)*(EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(36)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))))
RATES(13) = ALGBRC(48) - ALGBRC(45)
RATES(14) = ALGBRC(45) - ALGBRC(48)
ALGBRC(38) = CONSTS(130)*CONSTS(131)*log( CONSTS(58)*STATES(7))
ALGBRC(47) = CONSTS(3)*(EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))))
ALGBRC(41) = CONSTS(130)*CONSTS(131)*log( CONSTS(59)*STATES(8))
ALGBRC(50) = CONSTS(4)*(EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))))
RATES(7) = ALGBRC(47) - ALGBRC(50)
ALGBRC(30) = CONSTS(130)*CONSTS(131)*log( CONSTS(79)*STATES(15))
ALGBRC(33) = CONSTS(130)*CONSTS(131)*log( CONSTS(80)*STATES(16))
ALGBRC(51) = CONSTS(36)*(EXP((ALGBRC(36)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))))
RATES(17) = ALGBRC(48) - ALGBRC(51)
ALGBRC(53) = CONSTS(5)*(EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))))
RATES(6) = (ALGBRC(44) - ALGBRC(47))+ALGBRC(53)
RATES(8) = ALGBRC(50) - ALGBRC(53)
ALGBRC(54) = CONSTS(37)*(EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(2)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))))
RATES(12) = ALGBRC(54) - ALGBRC(45)
RATES(16) = ALGBRC(51) - ALGBRC(54)
ALGBRC(39) = CONSTS(130)*CONSTS(131)*log( CONSTS(82)*STATES(18))
ALGBRC(42) = CONSTS(130)*CONSTS(131)*log( CONSTS(83)*STATES(19))
ALGBRC(56) = CONSTS(38)*(EXP((ALGBRC(30)+ALGBRC(39))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(42)/( CONSTS(130)*CONSTS(131))))
RATES(15) = (ALGBRC(54) - ALGBRC(51)) - ALGBRC(56)
RATES(18) = - ALGBRC(56)
RATES(19) = ALGBRC(56)
ALGBRC(58) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))
ALGBRC(88) = ALGBRC(58)
ALGBRC(52) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(138) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(102)*ALGBRC(52)
ALGBRC(139) = CONSTS(28)*(EXP(ALGBRC(88)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(138)/( CONSTS(130)*CONSTS(131))))
ALGBRC(46) = STATES(1)/CONSTS(113)
ALGBRC(70) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(102) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(103) = CONSTS(10)*(EXP(ALGBRC(70)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(102)/( CONSTS(130)*CONSTS(131))))
ALGBRC(76) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(114) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(115) = CONSTS(16)*(EXP(ALGBRC(76)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(114)/( CONSTS(130)*CONSTS(131))))
ALGBRC(80) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(122) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(123) = CONSTS(20)*(EXP(ALGBRC(80)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(122)/( CONSTS(130)*CONSTS(131))))
ALGBRC(141) = ((- ALGBRC(103) - ALGBRC(115)) - ALGBRC(123)) - ALGBRC(139)
RATES(25) = ALGBRC(141)
ALGBRC(71) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(104) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(105) = CONSTS(13)*(EXP(ALGBRC(71)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(104)/( CONSTS(130)*CONSTS(131))))
ALGBRC(77) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(116) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(117) = CONSTS(18)*(EXP(ALGBRC(77)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(116)/( CONSTS(130)*CONSTS(131))))
ALGBRC(81) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(124) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(125) = CONSTS(22)*(EXP(ALGBRC(81)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(124)/( CONSTS(130)*CONSTS(131))))
ALGBRC(142) = ((- ALGBRC(105) - ALGBRC(117)) - ALGBRC(125)) - - ALGBRC(139)
RATES(31) = ALGBRC(142)
ALGBRC(59) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))
ALGBRC(89) = ALGBRC(59)
ALGBRC(140) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(102)*ALGBRC(52)
ALGBRC(143) = CONSTS(30)*(EXP(ALGBRC(89)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(140)/( CONSTS(130)*CONSTS(131))))
ALGBRC(72) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(106) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(107) = CONSTS(11)*(EXP(ALGBRC(72)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(106)/( CONSTS(130)*CONSTS(131))))
ALGBRC(84) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(130) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(131) = CONSTS(24)*(EXP(ALGBRC(84)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(130)/( CONSTS(130)*CONSTS(131))))
ALGBRC(145) = ((- ALGBRC(107)+ALGBRC(115)) - ALGBRC(143)) - ALGBRC(131)
RATES(26) = ALGBRC(145)
ALGBRC(73) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(108) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(109) = CONSTS(14)*(EXP(ALGBRC(73)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(108)/( CONSTS(130)*CONSTS(131))))
ALGBRC(85) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(132) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(133) = CONSTS(26)*(EXP(ALGBRC(85)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(132)/( CONSTS(130)*CONSTS(131))))
ALGBRC(146) = ((- ALGBRC(109)+ALGBRC(117)) - - ALGBRC(143)) - ALGBRC(133)
RATES(32) = ALGBRC(146)
ALGBRC(60) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))
ALGBRC(90) = ALGBRC(60)
ALGBRC(144) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(102)*ALGBRC(52)
ALGBRC(147) = CONSTS(32)*(EXP(ALGBRC(90)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(144)/( CONSTS(130)*CONSTS(131))))
ALGBRC(74) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(110) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(111) = CONSTS(12)*(EXP(ALGBRC(74)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(110)/( CONSTS(130)*CONSTS(131))))
ALGBRC(149) = ((- ALGBRC(111)+ALGBRC(123)) - ALGBRC(147))+ALGBRC(131)
RATES(27) = ALGBRC(149)
ALGBRC(75) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(112) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(113) = CONSTS(15)*(EXP(ALGBRC(75)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(112)/( CONSTS(130)*CONSTS(131))))
ALGBRC(150) = ((- ALGBRC(113)+ALGBRC(125)) - - ALGBRC(147))+ALGBRC(133)
RATES(33) = ALGBRC(150)
ALGBRC(61) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))
ALGBRC(91) = ALGBRC(61)
ALGBRC(148) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(102)*ALGBRC(52)
ALGBRC(151) = CONSTS(29)*(EXP(ALGBRC(91)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(148)/( CONSTS(130)*CONSTS(131))))
ALGBRC(78) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(118) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(119) = CONSTS(17)*(EXP(ALGBRC(78)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(118)/( CONSTS(130)*CONSTS(131))))
ALGBRC(82) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(126) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(127) = CONSTS(21)*(EXP(ALGBRC(82)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(126)/( CONSTS(130)*CONSTS(131))))
ALGBRC(153) = ((- ALGBRC(127)+ALGBRC(103)) - ALGBRC(119)) - ALGBRC(151)
RATES(28) = ALGBRC(153)
ALGBRC(79) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(120) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(121) = CONSTS(19)*(EXP(ALGBRC(79)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(120)/( CONSTS(130)*CONSTS(131))))
ALGBRC(83) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(128) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(129) = CONSTS(23)*(EXP(ALGBRC(83)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(128)/( CONSTS(130)*CONSTS(131))))
ALGBRC(154) = ((- ALGBRC(129)+ALGBRC(105)) - ALGBRC(121)) - - ALGBRC(151)
RATES(34) = ALGBRC(154)
ALGBRC(62) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))
ALGBRC(92) = ALGBRC(62)
ALGBRC(152) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(102)*ALGBRC(52)
ALGBRC(155) = CONSTS(31)*(EXP(ALGBRC(92)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(152)/( CONSTS(130)*CONSTS(131))))
ALGBRC(86) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(134) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(135) = CONSTS(25)*(EXP(ALGBRC(86)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(134)/( CONSTS(130)*CONSTS(131))))
ALGBRC(157) = ((ALGBRC(119)+ALGBRC(107)) - ALGBRC(155)) - ALGBRC(135)
RATES(29) = ALGBRC(157)
ALGBRC(64) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))
ALGBRC(93) = ALGBRC(64)
ALGBRC(156) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(102)*ALGBRC(52)
ALGBRC(158) = CONSTS(33)*(EXP(ALGBRC(93)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(156)/( CONSTS(130)*CONSTS(131))))
ALGBRC(160) = ((ALGBRC(127)+ALGBRC(111)) - ALGBRC(158))+ALGBRC(135)
RATES(30) = ALGBRC(160)
ALGBRC(63) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))
ALGBRC(66) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
ALGBRC(49) = CONSTS(130)*CONSTS(131)*log( CONSTS(61)*STATES(9))
ALGBRC(98) = ALGBRC(49)+ALGBRC(63)
ALGBRC(94) = CONSTS(128)*CONSTS(114)*ALGBRC(46)
ALGBRC(159) = TERNRY(ALGBRC(94).EQ.0.00000, CONSTS(6)*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(6)*ALGBRC(94))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(94)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(65) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))
ALGBRC(67) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
ALGBRC(99) = ALGBRC(49)+ALGBRC(65)
ALGBRC(95) = CONSTS(128)*CONSTS(114)*ALGBRC(46)
ALGBRC(161) = TERNRY(ALGBRC(95).EQ.0.00000, CONSTS(7)*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(7)*ALGBRC(95))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(95)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(162) = ALGBRC(161)+ALGBRC(159)
RATES(9) = ALGBRC(162)
ALGBRC(57) = CONSTS(130)*CONSTS(131)*log( CONSTS(62)*STATES(10))
ALGBRC(68) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
ALGBRC(55) = CONSTS(130)*CONSTS(131)*log( CONSTS(63)*STATES(11))
ALGBRC(100) = ALGBRC(55)+ALGBRC(63)
ALGBRC(96) = CONSTS(129)*CONSTS(114)*ALGBRC(46)
ALGBRC(163) = TERNRY(ALGBRC(96).EQ.0.00000, CONSTS(8)*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(8)*ALGBRC(96))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(96)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(87) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(136) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(137) = CONSTS(27)*(EXP(ALGBRC(87)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(136)/( CONSTS(130)*CONSTS(131))))
ALGBRC(167) = ((ALGBRC(121)+ALGBRC(109)) - ((((- ALGBRC(155)+ALGBRC(159)) - ALGBRC(159))+ALGBRC(163)) - ALGBRC(163))) - ALGBRC(137)
RATES(35) = ALGBRC(167)
ALGBRC(69) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
ALGBRC(101) = ALGBRC(55)+ALGBRC(65)
ALGBRC(97) = CONSTS(129)*CONSTS(114)*ALGBRC(46)
ALGBRC(165) = TERNRY(ALGBRC(97).EQ.0.00000, CONSTS(9)*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(9)*ALGBRC(97))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(97)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(168) = - ALGBRC(163) - ALGBRC(165)
RATES(10) = ALGBRC(168)
ALGBRC(169) = ALGBRC(165)+ALGBRC(163)
RATES(11) = ALGBRC(169)
ALGBRC(172) = ((ALGBRC(129)+ALGBRC(113)) - ((((- ALGBRC(158)+ALGBRC(161)) - ALGBRC(161))+ALGBRC(165)) - ALGBRC(165)))+ALGBRC(137)
RATES(36) = ALGBRC(172)
ALGBRC(170) = CONSTS(114)*((((((((((((((((((((((((((((((((((((( CONSTS(109)*ALGBRC(125) - CONSTS(105)*ALGBRC(125))+ CONSTS(110)*ALGBRC(133)) - CONSTS(106)*ALGBRC(133))+ CONSTS(108)*ALGBRC(121)) - CONSTS(104)*ALGBRC(121))+ CONSTS(109)*ALGBRC(129)) - CONSTS(105)*ALGBRC(129))+ CONSTS(110)*ALGBRC(137)) - CONSTS(106)*ALGBRC(137))+ CONSTS(108)*ALGBRC(117)) - CONSTS(104)*ALGBRC(117))+ CONSTS(107)*ALGBRC(113)) - CONSTS(103)*ALGBRC(113))+ CONSTS(107)*ALGBRC(109)) - CONSTS(103)*ALGBRC(109))+ CONSTS(107)*ALGBRC(105)) - CONSTS(103)*ALGBRC(105))+ CONSTS(110)*ALGBRC(135)) - CONSTS(106)*ALGBRC(135))+ CONSTS(109)*ALGBRC(127)) - CONSTS(105)*ALGBRC(127))+ CONSTS(108)*ALGBRC(119)) - CONSTS(104)*ALGBRC(119))+ CONSTS(110)*ALGBRC(131)) - CONSTS(106)*ALGBRC(131))+ CONSTS(109)*ALGBRC(123)) - CONSTS(105)*ALGBRC(123))+ CONSTS(108)*ALGBRC(115)) - CONSTS(104)*ALGBRC(115))+ CONSTS(107)*ALGBRC(111)) - CONSTS(103)*ALGBRC(111))+ CONSTS(107)*ALGBRC(107)) - CONSTS(103)*ALGBRC(107))+ CONSTS(107)*ALGBRC(103)) - CONSTS(103)*ALGBRC(103)) - ( CONSTS(128)*ALGBRC(159)+ CONSTS(128)*ALGBRC(161))) - ( CONSTS(129)*ALGBRC(163)+ CONSTS(129)*ALGBRC(165)))
ALGBRC(10) = TERNRY(VOI.GT.CONSTS(124).AND.VOI.LT.CONSTS(125), CONSTS(126), CONSTS(127))
ALGBRC(173) = - ALGBRC(170)+ALGBRC(10)
RATES(1) = ALGBRC(173)
ALGBRC(171) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(174) = CONSTS(130)*CONSTS(131)*log( CONSTS(85)*STATES(37))
ALGBRC(177) = CONSTS(130)*CONSTS(131)*log( CONSTS(88)*STATES(40))
ALGBRC(179) = CONSTS(40)*(EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(174)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
ALGBRC(175) = CONSTS(130)*CONSTS(131)*log( CONSTS(86)*STATES(38))
ALGBRC(180) = CONSTS(41)*(EXP((ALGBRC(174)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(175)/( CONSTS(130)*CONSTS(131))))
RATES(37) = ALGBRC(179) - ALGBRC(180)
ALGBRC(176) = CONSTS(130)*CONSTS(131)*log( CONSTS(87)*STATES(39))
ALGBRC(181) = CONSTS(42)*(EXP((ALGBRC(175)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))))
RATES(38) = ALGBRC(180) - ALGBRC(181)
ALGBRC(182) = CONSTS(43)*(EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(177)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
RATES(40) = ALGBRC(182) - ALGBRC(179)
RATES(39) = ALGBRC(181) - ALGBRC(182)
ALGBRC(189) = CONSTS(130)*CONSTS(131)*log( CONSTS(101)*STATES(22))
ALGBRC(190) = CONSTS(166)+ALGBRC(189)
ALGBRC(199) = CONSTS(52)*(EXP(CONSTS(177)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(190)/( CONSTS(130)*CONSTS(131))))
ALGBRC(200) = ALGBRC(199)
RATES(22) = ALGBRC(200)
ALGBRC(188) = CONSTS(130)*CONSTS(131)*log( CONSTS(99)*STATES(24))
ALGBRC(191) = CONSTS(166)+ALGBRC(188)
ALGBRC(201) = CONSTS(44)*(EXP(ALGBRC(191)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(178)/( CONSTS(130)*CONSTS(131))))
ALGBRC(202) = - ALGBRC(201)
RATES(24) = ALGBRC(202)
ALGBRC(187) = CONSTS(130)*CONSTS(131)*log( CONSTS(100)*STATES(23))
ALGBRC(192) = ALGBRC(187)+CONSTS(174)
ALGBRC(203) = CONSTS(48)*(EXP(CONSTS(179)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(192)/( CONSTS(130)*CONSTS(131))))
ALGBRC(205) = ALGBRC(203)
RATES(23) = ALGBRC(205)
ALGBRC(166) = CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
ALGBRC(178) = CONSTS(39)*EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131)))*(EXP(ALGBRC(166)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(171)/( CONSTS(130)*CONSTS(131))))
ALGBRC(185) = CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
ALGBRC(196) = CONSTS(175)+ CONSTS(155)*ALGBRC(185)
ALGBRC(212) = CONSTS(49)*(EXP(CONSTS(183)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(196)/( CONSTS(130)*CONSTS(131))))
ALGBRC(220) = CONSTS(155)*ALGBRC(212)
RATES(20) = - ALGBRC(178)+ALGBRC(220)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
ALGBRC(219) = ALGBRC(214)/CONSTS(146)
RATES(44) = ALGBRC(219)
ALGBRC(186) = CONSTS(130)*CONSTS(131)*log( CONSTS(98)*STATES(21))
ALGBRC(193) = CONSTS(169)+ALGBRC(186)
ALGBRC(206) = CONSTS(46)*(EXP(ALGBRC(193)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(180)/( CONSTS(130)*CONSTS(131))))
ALGBRC(195) = CONSTS(172)+ CONSTS(156)*ALGBRC(186)
ALGBRC(209) = CONSTS(47)*(EXP(CONSTS(182)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(195)/( CONSTS(130)*CONSTS(131))))
ALGBRC(197) = CONSTS(175)+ CONSTS(156)*ALGBRC(186)
ALGBRC(221) = CONSTS(50)*(EXP(ALGBRC(197)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(184)/( CONSTS(130)*CONSTS(131))))
ALGBRC(198) = ALGBRC(186)+CONSTS(173)
ALGBRC(224) = CONSTS(51)*(EXP(CONSTS(185)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(198)/( CONSTS(130)*CONSTS(131))))
ALGBRC(227) = (( CONSTS(156)*ALGBRC(209) - ALGBRC(206)) - CONSTS(156)*ALGBRC(221))+ALGBRC(224)
RATES(21) = ALGBRC(227)
ALGBRC(25) = CONSTS(130)*CONSTS(131)*log( CONSTS(53)*STATES(2))
ALGBRC(29) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(34) = (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(8)*(STATES(43)+CONSTS(137))
ALGBRC(37) = (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(7)*STATES(42)
ALGBRC(40) = ALGBRC(34)+ALGBRC(37)
ALGBRC(43) = ALGBRC(40)
ALGBRC(223) = (( CONSTS(151)*ALGBRC(213))/CONSTS(137))*ALGBRC(219)
ALGBRC(226) = ALGBRC(223)
ALGBRC(230) = ALGBRC(43)+ALGBRC(226)
ALGBRC(231) = ALGBRC(230)*1000.00
ALGBRC(232) = CONSTS(136)*CONSTS(161)*ALGBRC(231)
ALGBRC(233) = CONSTS(1)*(EXP((ALGBRC(29)+ALGBRC(25))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(31)+ALGBRC(232))/( CONSTS(130)*CONSTS(131))))
ALGBRC(164) = (( (CONSTS(102))*ALGBRC(151)+ CONSTS(102)*ALGBRC(155)+ CONSTS(102)*ALGBRC(147)+ CONSTS(102)*ALGBRC(143)+ CONSTS(102)*ALGBRC(139)+ CONSTS(102)*ALGBRC(158)) - ALGBRC(159)) - ALGBRC(161)
ALGBRC(183) = (( CONSTS(153)*ALGBRC(179) - CONSTS(152)*ALGBRC(180)) - CONSTS(153)*ALGBRC(181))+ CONSTS(152)*ALGBRC(182)
ALGBRC(184) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(194) = CONSTS(169)+ CONSTS(154)*ALGBRC(184)
ALGBRC(207) = CONSTS(45)*(EXP(ALGBRC(194)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(181)/( CONSTS(130)*CONSTS(131))))
ALGBRC(210) = - CONSTS(154)*ALGBRC(207)
RATES(4) = - ALGBRC(233)+ALGBRC(164)+ALGBRC(178)+ALGBRC(183)+ALGBRC(210)
RATES(2) = - ALGBRC(233)
RATES(3) = ALGBRC(233)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(8) = CONSTS(158)*STATES(42)*CONSTS(146)
ALGBRC(23) = - CONSTS(165) - ALGBRC(8)
ALGBRC(27) = ALGBRC(23)/CONSTS(146)
ALGBRC(9) = CONSTS(159)*STATES(43)*CONSTS(146)
ALGBRC(24) = - CONSTS(168) - ALGBRC(9)
ALGBRC(28) = ALGBRC(24)/CONSTS(146)
ALGBRC(31) = CONSTS(130)*CONSTS(131)*log( CONSTS(55)*STATES(3))
ALGBRC(32) = CONSTS(130)*CONSTS(131)*log( CONSTS(56)*STATES(5))
ALGBRC(35) = CONSTS(130)*CONSTS(131)*log( CONSTS(57)*STATES(6))
ALGBRC(44) = CONSTS(2)*(EXP((ALGBRC(32)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(35)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))))
ALGBRC(2) = CONSTS(130)*CONSTS(131)*log( CONSTS(76)*STATES(12))
ALGBRC(22) = CONSTS(130)*CONSTS(131)*log( CONSTS(77)*STATES(13))
ALGBRC(26) = CONSTS(130)*CONSTS(131)*log( CONSTS(78)*STATES(14))
ALGBRC(45) = CONSTS(34)*(EXP((ALGBRC(2)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))))
ALGBRC(36) = CONSTS(130)*CONSTS(131)*log( CONSTS(81)*STATES(17))
ALGBRC(48) = CONSTS(35)*(EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(36)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))))
ALGBRC(38) = CONSTS(130)*CONSTS(131)*log( CONSTS(58)*STATES(7))
ALGBRC(47) = CONSTS(3)*(EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))))
ALGBRC(41) = CONSTS(130)*CONSTS(131)*log( CONSTS(59)*STATES(8))
ALGBRC(50) = CONSTS(4)*(EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))))
ALGBRC(30) = CONSTS(130)*CONSTS(131)*log( CONSTS(79)*STATES(15))
ALGBRC(33) = CONSTS(130)*CONSTS(131)*log( CONSTS(80)*STATES(16))
ALGBRC(51) = CONSTS(36)*(EXP((ALGBRC(36)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))))
ALGBRC(53) = CONSTS(5)*(EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))))
ALGBRC(54) = CONSTS(37)*(EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(2)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))))
ALGBRC(39) = CONSTS(130)*CONSTS(131)*log( CONSTS(82)*STATES(18))
ALGBRC(42) = CONSTS(130)*CONSTS(131)*log( CONSTS(83)*STATES(19))
ALGBRC(56) = CONSTS(38)*(EXP((ALGBRC(30)+ALGBRC(39))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(42)/( CONSTS(130)*CONSTS(131))))
ALGBRC(58) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))
ALGBRC(88) = ALGBRC(58)
ALGBRC(52) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(138) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(102)*ALGBRC(52)
ALGBRC(139) = CONSTS(28)*(EXP(ALGBRC(88)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(138)/( CONSTS(130)*CONSTS(131))))
ALGBRC(46) = STATES(1)/CONSTS(113)
ALGBRC(70) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(102) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(103) = CONSTS(10)*(EXP(ALGBRC(70)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(102)/( CONSTS(130)*CONSTS(131))))
ALGBRC(76) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(114) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(115) = CONSTS(16)*(EXP(ALGBRC(76)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(114)/( CONSTS(130)*CONSTS(131))))
ALGBRC(80) = CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(122) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(123) = CONSTS(20)*(EXP(ALGBRC(80)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(122)/( CONSTS(130)*CONSTS(131))))
ALGBRC(141) = ((- ALGBRC(103) - ALGBRC(115)) - ALGBRC(123)) - ALGBRC(139)
ALGBRC(71) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(104) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(105) = CONSTS(13)*(EXP(ALGBRC(71)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(104)/( CONSTS(130)*CONSTS(131))))
ALGBRC(77) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(116) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(117) = CONSTS(18)*(EXP(ALGBRC(77)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(116)/( CONSTS(130)*CONSTS(131))))
ALGBRC(81) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(124) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(125) = CONSTS(22)*(EXP(ALGBRC(81)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(124)/( CONSTS(130)*CONSTS(131))))
ALGBRC(142) = ((- ALGBRC(105) - ALGBRC(117)) - ALGBRC(125)) - - ALGBRC(139)
ALGBRC(59) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))
ALGBRC(89) = ALGBRC(59)
ALGBRC(140) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(102)*ALGBRC(52)
ALGBRC(143) = CONSTS(30)*(EXP(ALGBRC(89)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(140)/( CONSTS(130)*CONSTS(131))))
ALGBRC(72) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(106) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(107) = CONSTS(11)*(EXP(ALGBRC(72)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(106)/( CONSTS(130)*CONSTS(131))))
ALGBRC(84) = CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(130) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(131) = CONSTS(24)*(EXP(ALGBRC(84)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(130)/( CONSTS(130)*CONSTS(131))))
ALGBRC(145) = ((- ALGBRC(107)+ALGBRC(115)) - ALGBRC(143)) - ALGBRC(131)
ALGBRC(73) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(108) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(109) = CONSTS(14)*(EXP(ALGBRC(73)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(108)/( CONSTS(130)*CONSTS(131))))
ALGBRC(85) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(132) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(133) = CONSTS(26)*(EXP(ALGBRC(85)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(132)/( CONSTS(130)*CONSTS(131))))
ALGBRC(146) = ((- ALGBRC(109)+ALGBRC(117)) - - ALGBRC(143)) - ALGBRC(133)
ALGBRC(60) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))
ALGBRC(90) = ALGBRC(60)
ALGBRC(144) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(102)*ALGBRC(52)
ALGBRC(147) = CONSTS(32)*(EXP(ALGBRC(90)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(144)/( CONSTS(130)*CONSTS(131))))
ALGBRC(74) = CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(110) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(111) = CONSTS(12)*(EXP(ALGBRC(74)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(110)/( CONSTS(130)*CONSTS(131))))
ALGBRC(149) = ((- ALGBRC(111)+ALGBRC(123)) - ALGBRC(147))+ALGBRC(131)
ALGBRC(75) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
ALGBRC(112) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
ALGBRC(113) = CONSTS(15)*(EXP(ALGBRC(75)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(112)/( CONSTS(130)*CONSTS(131))))
ALGBRC(150) = ((- ALGBRC(113)+ALGBRC(125)) - - ALGBRC(147))+ALGBRC(133)
ALGBRC(61) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))
ALGBRC(91) = ALGBRC(61)
ALGBRC(148) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(102)*ALGBRC(52)
ALGBRC(151) = CONSTS(29)*(EXP(ALGBRC(91)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(148)/( CONSTS(130)*CONSTS(131))))
ALGBRC(78) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(118) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(119) = CONSTS(17)*(EXP(ALGBRC(78)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(118)/( CONSTS(130)*CONSTS(131))))
ALGBRC(82) = CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(126) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(127) = CONSTS(21)*(EXP(ALGBRC(82)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(126)/( CONSTS(130)*CONSTS(131))))
ALGBRC(153) = ((- ALGBRC(127)+ALGBRC(103)) - ALGBRC(119)) - ALGBRC(151)
ALGBRC(79) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
ALGBRC(120) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
ALGBRC(121) = CONSTS(19)*(EXP(ALGBRC(79)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(120)/( CONSTS(130)*CONSTS(131))))
ALGBRC(83) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
ALGBRC(128) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
ALGBRC(129) = CONSTS(23)*(EXP(ALGBRC(83)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(128)/( CONSTS(130)*CONSTS(131))))
ALGBRC(154) = ((- ALGBRC(129)+ALGBRC(105)) - ALGBRC(121)) - - ALGBRC(151)
ALGBRC(62) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))
ALGBRC(92) = ALGBRC(62)
ALGBRC(152) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(102)*ALGBRC(52)
ALGBRC(155) = CONSTS(31)*(EXP(ALGBRC(92)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(152)/( CONSTS(130)*CONSTS(131))))
ALGBRC(86) = CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(134) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(135) = CONSTS(25)*(EXP(ALGBRC(86)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(134)/( CONSTS(130)*CONSTS(131))))
ALGBRC(157) = ((ALGBRC(119)+ALGBRC(107)) - ALGBRC(155)) - ALGBRC(135)
ALGBRC(64) = CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))
ALGBRC(93) = ALGBRC(64)
ALGBRC(156) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(102)*ALGBRC(52)
ALGBRC(158) = CONSTS(33)*(EXP(ALGBRC(93)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(156)/( CONSTS(130)*CONSTS(131))))
ALGBRC(160) = ((ALGBRC(127)+ALGBRC(111)) - ALGBRC(158))+ALGBRC(135)
ALGBRC(63) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))
ALGBRC(66) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
ALGBRC(49) = CONSTS(130)*CONSTS(131)*log( CONSTS(61)*STATES(9))
ALGBRC(98) = ALGBRC(49)+ALGBRC(63)
ALGBRC(94) = CONSTS(128)*CONSTS(114)*ALGBRC(46)
ALGBRC(159) = TERNRY(ALGBRC(94).EQ.0.00000, CONSTS(6)*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(6)*ALGBRC(94))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(94)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(65) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))
ALGBRC(67) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
ALGBRC(99) = ALGBRC(49)+ALGBRC(65)
ALGBRC(95) = CONSTS(128)*CONSTS(114)*ALGBRC(46)
ALGBRC(161) = TERNRY(ALGBRC(95).EQ.0.00000, CONSTS(7)*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(7)*ALGBRC(95))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(95)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(162) = ALGBRC(161)+ALGBRC(159)
ALGBRC(57) = CONSTS(130)*CONSTS(131)*log( CONSTS(62)*STATES(10))
ALGBRC(68) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
ALGBRC(55) = CONSTS(130)*CONSTS(131)*log( CONSTS(63)*STATES(11))
ALGBRC(100) = ALGBRC(55)+ALGBRC(63)
ALGBRC(96) = CONSTS(129)*CONSTS(114)*ALGBRC(46)
ALGBRC(163) = TERNRY(ALGBRC(96).EQ.0.00000, CONSTS(8)*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(8)*ALGBRC(96))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(96)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(87) = CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
ALGBRC(136) = CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
ALGBRC(137) = CONSTS(27)*(EXP(ALGBRC(87)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(136)/( CONSTS(130)*CONSTS(131))))
ALGBRC(167) = ((ALGBRC(121)+ALGBRC(109)) - ((((- ALGBRC(155)+ALGBRC(159)) - ALGBRC(159))+ALGBRC(163)) - ALGBRC(163))) - ALGBRC(137)
ALGBRC(69) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
ALGBRC(101) = ALGBRC(55)+ALGBRC(65)
ALGBRC(97) = CONSTS(129)*CONSTS(114)*ALGBRC(46)
ALGBRC(165) = TERNRY(ALGBRC(97).EQ.0.00000, CONSTS(9)*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))), ((( CONSTS(9)*ALGBRC(97))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(97)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))))
ALGBRC(168) = - ALGBRC(163) - ALGBRC(165)
ALGBRC(169) = ALGBRC(165)+ALGBRC(163)
ALGBRC(172) = ((ALGBRC(129)+ALGBRC(113)) - ((((- ALGBRC(158)+ALGBRC(161)) - ALGBRC(161))+ALGBRC(165)) - ALGBRC(165)))+ALGBRC(137)
ALGBRC(170) = CONSTS(114)*((((((((((((((((((((((((((((((((((((( CONSTS(109)*ALGBRC(125) - CONSTS(105)*ALGBRC(125))+ CONSTS(110)*ALGBRC(133)) - CONSTS(106)*ALGBRC(133))+ CONSTS(108)*ALGBRC(121)) - CONSTS(104)*ALGBRC(121))+ CONSTS(109)*ALGBRC(129)) - CONSTS(105)*ALGBRC(129))+ CONSTS(110)*ALGBRC(137)) - CONSTS(106)*ALGBRC(137))+ CONSTS(108)*ALGBRC(117)) - CONSTS(104)*ALGBRC(117))+ CONSTS(107)*ALGBRC(113)) - CONSTS(103)*ALGBRC(113))+ CONSTS(107)*ALGBRC(109)) - CONSTS(103)*ALGBRC(109))+ CONSTS(107)*ALGBRC(105)) - CONSTS(103)*ALGBRC(105))+ CONSTS(110)*ALGBRC(135)) - CONSTS(106)*ALGBRC(135))+ CONSTS(109)*ALGBRC(127)) - CONSTS(105)*ALGBRC(127))+ CONSTS(108)*ALGBRC(119)) - CONSTS(104)*ALGBRC(119))+ CONSTS(110)*ALGBRC(131)) - CONSTS(106)*ALGBRC(131))+ CONSTS(109)*ALGBRC(123)) - CONSTS(105)*ALGBRC(123))+ CONSTS(108)*ALGBRC(115)) - CONSTS(104)*ALGBRC(115))+ CONSTS(107)*ALGBRC(111)) - CONSTS(103)*ALGBRC(111))+ CONSTS(107)*ALGBRC(107)) - CONSTS(103)*ALGBRC(107))+ CONSTS(107)*ALGBRC(103)) - CONSTS(103)*ALGBRC(103)) - ( CONSTS(128)*ALGBRC(159)+ CONSTS(128)*ALGBRC(161))) - ( CONSTS(129)*ALGBRC(163)+ CONSTS(129)*ALGBRC(165)))
ALGBRC(10) = TERNRY(VOI.GT.CONSTS(124).AND.VOI.LT.CONSTS(125), CONSTS(126), CONSTS(127))
ALGBRC(173) = - ALGBRC(170)+ALGBRC(10)
ALGBRC(171) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(174) = CONSTS(130)*CONSTS(131)*log( CONSTS(85)*STATES(37))
ALGBRC(177) = CONSTS(130)*CONSTS(131)*log( CONSTS(88)*STATES(40))
ALGBRC(179) = CONSTS(40)*(EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(174)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
ALGBRC(175) = CONSTS(130)*CONSTS(131)*log( CONSTS(86)*STATES(38))
ALGBRC(180) = CONSTS(41)*(EXP((ALGBRC(174)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(175)/( CONSTS(130)*CONSTS(131))))
ALGBRC(176) = CONSTS(130)*CONSTS(131)*log( CONSTS(87)*STATES(39))
ALGBRC(181) = CONSTS(42)*(EXP((ALGBRC(175)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))))
ALGBRC(182) = CONSTS(43)*(EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(177)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
ALGBRC(189) = CONSTS(130)*CONSTS(131)*log( CONSTS(101)*STATES(22))
ALGBRC(190) = CONSTS(166)+ALGBRC(189)
ALGBRC(199) = CONSTS(52)*(EXP(CONSTS(177)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(190)/( CONSTS(130)*CONSTS(131))))
ALGBRC(200) = ALGBRC(199)
ALGBRC(188) = CONSTS(130)*CONSTS(131)*log( CONSTS(99)*STATES(24))
ALGBRC(191) = CONSTS(166)+ALGBRC(188)
ALGBRC(201) = CONSTS(44)*(EXP(ALGBRC(191)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(178)/( CONSTS(130)*CONSTS(131))))
ALGBRC(202) = - ALGBRC(201)
ALGBRC(187) = CONSTS(130)*CONSTS(131)*log( CONSTS(100)*STATES(23))
ALGBRC(192) = ALGBRC(187)+CONSTS(174)
ALGBRC(203) = CONSTS(48)*(EXP(CONSTS(179)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(192)/( CONSTS(130)*CONSTS(131))))
ALGBRC(205) = ALGBRC(203)
ALGBRC(166) = CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
ALGBRC(178) = CONSTS(39)*EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131)))*(EXP(ALGBRC(166)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(171)/( CONSTS(130)*CONSTS(131))))
ALGBRC(185) = CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
ALGBRC(196) = CONSTS(175)+ CONSTS(155)*ALGBRC(185)
ALGBRC(212) = CONSTS(49)*(EXP(CONSTS(183)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(196)/( CONSTS(130)*CONSTS(131))))
ALGBRC(220) = CONSTS(155)*ALGBRC(212)
ALGBRC(219) = ALGBRC(214)/CONSTS(146)
ALGBRC(186) = CONSTS(130)*CONSTS(131)*log( CONSTS(98)*STATES(21))
ALGBRC(193) = CONSTS(169)+ALGBRC(186)
ALGBRC(206) = CONSTS(46)*(EXP(ALGBRC(193)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(180)/( CONSTS(130)*CONSTS(131))))
ALGBRC(195) = CONSTS(172)+ CONSTS(156)*ALGBRC(186)
ALGBRC(209) = CONSTS(47)*(EXP(CONSTS(182)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(195)/( CONSTS(130)*CONSTS(131))))
ALGBRC(197) = CONSTS(175)+ CONSTS(156)*ALGBRC(186)
ALGBRC(221) = CONSTS(50)*(EXP(ALGBRC(197)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(184)/( CONSTS(130)*CONSTS(131))))
ALGBRC(198) = ALGBRC(186)+CONSTS(173)
ALGBRC(224) = CONSTS(51)*(EXP(CONSTS(185)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(198)/( CONSTS(130)*CONSTS(131))))
ALGBRC(227) = (( CONSTS(156)*ALGBRC(209) - ALGBRC(206)) - CONSTS(156)*ALGBRC(221))+ALGBRC(224)
ALGBRC(25) = CONSTS(130)*CONSTS(131)*log( CONSTS(53)*STATES(2))
ALGBRC(29) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(34) = (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(8)*(STATES(43)+CONSTS(137))
ALGBRC(37) = (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(7)*STATES(42)
ALGBRC(40) = ALGBRC(34)+ALGBRC(37)
ALGBRC(43) = ALGBRC(40)
ALGBRC(223) = (( CONSTS(151)*ALGBRC(213))/CONSTS(137))*ALGBRC(219)
ALGBRC(226) = ALGBRC(223)
ALGBRC(230) = ALGBRC(43)+ALGBRC(226)
ALGBRC(231) = ALGBRC(230)*1000.00
ALGBRC(232) = CONSTS(136)*CONSTS(161)*ALGBRC(231)
ALGBRC(233) = CONSTS(1)*(EXP((ALGBRC(29)+ALGBRC(25))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(31)+ALGBRC(232))/( CONSTS(130)*CONSTS(131))))
ALGBRC(164) = (( (CONSTS(102))*ALGBRC(151)+ CONSTS(102)*ALGBRC(155)+ CONSTS(102)*ALGBRC(147)+ CONSTS(102)*ALGBRC(143)+ CONSTS(102)*ALGBRC(139)+ CONSTS(102)*ALGBRC(158)) - ALGBRC(159)) - ALGBRC(161)
ALGBRC(183) = (( CONSTS(153)*ALGBRC(179) - CONSTS(152)*ALGBRC(180)) - CONSTS(153)*ALGBRC(181))+ CONSTS(152)*ALGBRC(182)
ALGBRC(184) = CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
ALGBRC(194) = CONSTS(169)+ CONSTS(154)*ALGBRC(184)
ALGBRC(207) = CONSTS(45)*(EXP(ALGBRC(194)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(181)/( CONSTS(130)*CONSTS(131))))
ALGBRC(210) = - CONSTS(154)*ALGBRC(207)
ALGBRC(1) = STATES(1)/CONSTS(113)
ALGBRC(3) = STATES(12)+STATES(14)+STATES(17)+STATES(16)
ALGBRC(4) = STATES(13)+STATES(14)
ALGBRC(5) = STATES(18)+STATES(19)
ALGBRC(6) = STATES(10)
ALGBRC(7) = STATES(11)
ALGBRC(11) = CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))
ALGBRC(12) = CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))
ALGBRC(13) = CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))
ALGBRC(14) = CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))
ALGBRC(15) = CONSTS(153)*STATES(40)+ CONSTS(152)*STATES(38)+ (CONSTS(152)+CONSTS(153))*STATES(39)
ALGBRC(16) = STATES(4)/CONSTS(157)
ALGBRC(17) = STATES(20)/CONSTS(163)
ALGBRC(18) = STATES(21)/CONSTS(167)
ALGBRC(19) = STATES(23)/CONSTS(157)
ALGBRC(20) = STATES(24)/CONSTS(157)
ALGBRC(21) = STATES(22)/CONSTS(157)
ALGBRC(204) = ALGBRC(199) - ALGBRC(201)
ALGBRC(208) = ALGBRC(206)
ALGBRC(211) = ALGBRC(201) - ALGBRC(207)
ALGBRC(217) = ALGBRC(207) - ALGBRC(209)
ALGBRC(218) = ALGBRC(209) - ALGBRC(203)
ALGBRC(222) = ALGBRC(203) - ALGBRC(212)
ALGBRC(225) = ALGBRC(212) - ALGBRC(221)
ALGBRC(228) = ALGBRC(221) - ALGBRC(224)
ALGBRC(229) = ALGBRC(224) - ALGBRC(199)
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0[1] = abs(ALGBRC(213) - (TERNRY(ALGBRC(214).GT.0.00000, CONSTS(148), CONSTS(149))))
minfunc_0[2] = abs(ALGBRC(215) - (CONSTS(150)/ALGBRC(213))*STATES(44)*CONSTS(146))
minfunc_0[3] = abs(ALGBRC(216) - (CONSTS(150)/ALGBRC(213))*(STATES(41) - CONSTS(137))*CONSTS(146))
minfunc_0[4] = abs(ALGBRC(214) - (- ALGBRC(215)+ALGBRC(216)))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END