C
C There are a total of 152 entries in the algebraic variable array.
C There are a total of 56 entries in each of the rate and state variable arrays.
C There are a total of 150 entries in the constant variable array.
C
C
C CONSTS(1) is kappa_1a in component BG_parameters (fmol_per_sec).
C CONSTS(2) is kappa_1b in component BG_parameters (fmol_per_sec).
C CONSTS(3) is kappa_2a in component BG_parameters (fmol_per_sec).
C CONSTS(4) is kappa_2b in component BG_parameters (fmol_per_sec).
C CONSTS(5) is kappa_3a in component BG_parameters (fmol_per_sec).
C CONSTS(6) is kappa_3b in component BG_parameters (fmol_per_sec).
C CONSTS(7) is kappa_4a in component BG_parameters (fmol_per_sec).
C CONSTS(8) is kappa_4b in component BG_parameters (fmol_per_sec).
C CONSTS(9) is kappa_5 in component BG_parameters (fmol_per_sec).
C CONSTS(10) is kappa_6 in component BG_parameters (fmol_per_sec).
C CONSTS(11) is kappa_7 in component BG_parameters (fmol_per_sec).
C CONSTS(12) is kappa_GiAC in component BG_parameters (fmol_per_sec).
C CONSTS(13) is kappa_Rswitch_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(14) is kappa_LRswitch_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(15) is kappa_C_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(16) is kappa_R_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(17) is kappa_L_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(18) is kappa_Act1_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(19) is kappa_Act2_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(20) is kappa_Hyd_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(21) is kappa_Reassoc_Gs in component BG_parameters (fmol_per_sec).
C CONSTS(22) is kappa_InternR_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(23) is kappa_InternLR_B1 in component BG_parameters (fmol_per_sec).
C CONSTS(24) is kappa_Rswitch_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(25) is kappa_LRswitch_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(26) is kappa_C_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(27) is kappa_R_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(28) is kappa_L_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(29) is kappa_Act1_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(30) is kappa_Act2_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(31) is kappa_Hyd_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(32) is kappa_Reassoc_Gi in component BG_parameters (fmol_per_sec).
C CONSTS(33) is kappa_InternR_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(34) is kappa_InternLR_M2 in component BG_parameters (fmol_per_sec).
C CONSTS(35) is K_ATP in component BG_parameters (per_fmol).
C CONSTS(36) is K_cAMP in component BG_parameters (per_fmol).
C CONSTS(37) is K_AC in component BG_parameters (per_fmol).
C CONSTS(38) is K_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(39) is K_Gsa_GTP_AC in component BG_parameters (per_fmol).
C CONSTS(40) is K_Gsa_GTP_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(41) is K_FSK_AC in component BG_parameters (per_fmol).
C CONSTS(42) is K_FSK_AC_ATP in component BG_parameters (per_fmol).
C CONSTS(43) is K_PDE in component BG_parameters (per_fmol).
C CONSTS(44) is K_PDE_cAMP in component BG_parameters (per_fmol).
C CONSTS(45) is K_five_AMP in component BG_parameters (per_fmol).
C CONSTS(46) is K_IBMX in component BG_parameters (per_fmol).
C CONSTS(47) is K_PDEinh in component BG_parameters (per_fmol).
C CONSTS(48) is K_Gsa_GTP in component BG_parameters (per_fmol).
C CONSTS(49) is K_FSK in component BG_parameters (per_fmol).
C CONSTS(50) is K_Gia_GTP in component BG_parameters (per_fmol).
C CONSTS(51) is K_ACinh in component BG_parameters (per_fmol).
C CONSTS(52) is K_PPi in component BG_parameters (per_fmol).
C CONSTS(53) is K_RB1_inactive in component BG_parameters (per_fmol).
C CONSTS(54) is K_L_RB1_inactive in component BG_parameters (per_fmol).
C CONSTS(55) is K_LB1 in component BG_parameters (per_fmol).
C CONSTS(56) is K_RB1 in component BG_parameters (per_fmol).
C CONSTS(57) is K_Gs in component BG_parameters (per_fmol).
C CONSTS(58) is K_RB1_Gs in component BG_parameters (per_fmol).
C CONSTS(59) is K_L_RB1 in component BG_parameters (per_fmol).
C CONSTS(60) is K_L_RB1_Gs in component BG_parameters (per_fmol).
C CONSTS(61) is K_Gsbetagamma in component BG_parameters (per_fmol).
C CONSTS(62) is K_Gsa_GDP in component BG_parameters (per_fmol).
C CONSTS(63) is K_GTP in component BG_parameters (per_fmol).
C CONSTS(64) is K_GDP in component BG_parameters (per_fmol).
C CONSTS(65) is K_Pi in component BG_parameters (per_fmol).
C CONSTS(66) is K_RB1_tag in component BG_parameters (per_fmol).
C CONSTS(67) is K_L_RB1_tag in component BG_parameters (per_fmol).
C CONSTS(68) is K_RB1_GRKArr in component BG_parameters (per_fmol).
C CONSTS(69) is K_L_RB1_GRKArr in component BG_parameters (per_fmol).
C CONSTS(70) is K_GRKArr in component BG_parameters (per_fmol).
C CONSTS(71) is K_RM2_inactive in component BG_parameters (per_fmol).
C CONSTS(72) is K_L_RM2_inactive in component BG_parameters (per_fmol).
C CONSTS(73) is K_LM2 in component BG_parameters (per_fmol).
C CONSTS(74) is K_RM2 in component BG_parameters (per_fmol).
C CONSTS(75) is K_Gi in component BG_parameters (per_fmol).
C CONSTS(76) is K_RM2_Gi in component BG_parameters (per_fmol).
C CONSTS(77) is K_L_RM2 in component BG_parameters (per_fmol).
C CONSTS(78) is K_L_RM2_Gi in component BG_parameters (per_fmol).
C CONSTS(79) is K_Gibetagamma in component BG_parameters (per_fmol).
C CONSTS(80) is K_Gia_GDP in component BG_parameters (per_fmol).
C CONSTS(81) is K_RM2_tag in component BG_parameters (per_fmol).
C CONSTS(82) is K_L_RM2_tag in component BG_parameters (per_fmol).
C CONSTS(83) is K_RM2_GRKArr in component BG_parameters (per_fmol).
C CONSTS(84) is K_L_RM2_GRKArr in component BG_parameters (per_fmol).
C VOI is time in component environment (second).
C CONSTS(85) is vol_myo in component environment (pL).
C CONSTS(86) is freq in component environment (dimensionless).
C CONSTS(87) is stimSt1 in component environment (second).
C CONSTS(88) is stimDur1 in component environment (second).
C CONSTS(89) is tRamp1 in component environment (second).
C CONSTS(90) is stimMag1 in component environment (fmol).
C CONSTS(91) is stimHolding1 in component environment (fmol).
C CONSTS(149) is m1 in component environment (fmol_per_sec).
C CONSTS(92) is stimSt2 in component environment (second).
C CONSTS(93) is stimDur2 in component environment (second).
C CONSTS(94) is tRamp2 in component environment (second).
C CONSTS(95) is stimMag2 in component environment (fmol).
C CONSTS(96) is stimHolding2 in component environment (fmol).
C CONSTS(150) is m2 in component environment (fmol_per_sec).
C CONSTS(97) is q_ATP_init in component environment (fmol).
C CONSTS(98) is q_AC_init in component environment (fmol).
C CONSTS(99) is q_cAMP_init in component environment (fmol).
C CONSTS(100) is q_AC_ATP_init in component environment (fmol).
C CONSTS(101) is q_FSK_init in component environment (fmol).
C CONSTS(102) is q_FSK_AC_init in component environment (fmol).
C CONSTS(103) is q_FSK_AC_ATP_init in component environment (fmol).
C CONSTS(104) is q_Gsa_GTP_init in component environment (fmol).
C CONSTS(105) is q_Gsa_GTP_AC_init in component environment (fmol).
C CONSTS(106) is q_Gsa_GTP_AC_ATP_init in component environment (fmol).
C CONSTS(107) is q_PDE_init in component environment (fmol).
C CONSTS(108) is q_PDEinh_init in component environment (fmol).
C CONSTS(109) is q_PDE_cAMP_init in component environment (fmol).
C CONSTS(110) is q_IBMX_init in component environment (fmol).
C CONSTS(111) is q_five_AMP_init in component environment (fmol).
C CONSTS(112) is q_Gia_GTP_init in component environment (fmol).
C CONSTS(113) is q_ACinh_init in component environment (fmol).
C CONSTS(114) is q_PPi_init in component environment (fmol).
C CONSTS(115) is q_RB1_inactive_init in component environment (fmol).
C CONSTS(116) is q_L_RB1_inactive_init in component environment (fmol).
C ALGBRC(1) is q_LB1_init in component environment (fmol).
C CONSTS(117) is q_RB1_init in component environment (fmol).
C CONSTS(118) is q_Gs_init in component environment (fmol).
C CONSTS(119) is q_RB1_Gs_init in component environment (fmol).
C CONSTS(120) is q_L_RB1_init in component environment (fmol).
C CONSTS(121) is q_L_RB1_Gs_init in component environment (fmol).
C CONSTS(122) is q_Gsbetagamma_init in component environment (fmol).
C CONSTS(123) is q_Gsa_GDP_init in component environment (fmol).
C CONSTS(124) is q_GTP_init in component environment (fmol).
C CONSTS(125) is q_GDP_init in component environment (fmol).
C CONSTS(126) is q_Pi_init in component environment (fmol).
C CONSTS(127) is q_RB1_tag_init in component environment (fmol).
C CONSTS(128) is q_L_RB1_tag_init in component environment (fmol).
C CONSTS(129) is q_RB1_GRKArr_init in component environment (fmol).
C CONSTS(130) is q_L_RB1_GRKArr_init in component environment (fmol).
C CONSTS(131) is q_GRKArr_init in component environment (fmol).
C CONSTS(132) is q_RM2_inactive_init in component environment (fmol).
C CONSTS(133) is q_L_RM2_inactive_init in component environment (fmol).
C ALGBRC(2) is q_LM2_init in component environment (fmol).
C CONSTS(134) is q_RM2_init in component environment (fmol).
C CONSTS(135) is q_Gi_init in component environment (fmol).
C CONSTS(136) is q_RM2_Gi_init in component environment (fmol).
C CONSTS(137) is q_L_RM2_init in component environment (fmol).
C CONSTS(138) is q_L_RM2_Gi_init in component environment (fmol).
C CONSTS(139) is q_Gibetagamma_init in component environment (fmol).
C CONSTS(140) is q_Gia_GDP_init in component environment (fmol).
C CONSTS(141) is q_RM2_tag_init in component environment (fmol).
C CONSTS(142) is q_L_RM2_tag_init in component environment (fmol).
C CONSTS(143) is q_RM2_GRKArr_init in component environment (fmol).
C CONSTS(144) is q_L_RM2_GRKArr_init in component environment (fmol).
C ALGBRC(59) is LB1_T in component environment (fmol).
C ALGBRC(60) is RB1_T in component environment (fmol).
C ALGBRC(44) is Gs_T in component environment (fmol).
C ALGBRC(17) is adenosine_T in component environment (fmol).
C ALGBRC(3) is q_ATP in component environment (fmol).
C ALGBRC(4) is q_cAMP in component environment (fmol).
C ALGBRC(5) is q_AC in component environment (fmol).
C ALGBRC(6) is q_AC_ATP in component environment (fmol).
C ALGBRC(7) is q_Gsa_GTP_AC in component environment (fmol).
C ALGBRC(8) is q_Gsa_GTP_AC_ATP in component environment (fmol).
C ALGBRC(9) is q_FSK_AC in component environment (fmol).
C ALGBRC(10) is q_FSK_AC_ATP in component environment (fmol).
C ALGBRC(11) is q_PDE in component environment (fmol).
C ALGBRC(14) is q_PDE_cAMP in component environment (fmol).
C ALGBRC(15) is q_five_AMP in component environment (fmol).
C ALGBRC(18) is q_IBMX in component environment (fmol).
C ALGBRC(19) is q_PDEinh in component environment (fmol).
C ALGBRC(20) is q_Gsa_GTP in component environment (fmol).
C ALGBRC(21) is q_FSK in component environment (fmol).
C ALGBRC(23) is q_Gia_GTP in component environment (fmol).
C ALGBRC(24) is q_ACinh in component environment (fmol).
C ALGBRC(26) is q_PPi in component environment (fmol).
C ALGBRC(25) is q_RB1_inactive in component environment (fmol).
C ALGBRC(27) is q_L_RB1_inactive in component environment (fmol).
C ALGBRC(28) is q_LB1 in component environment (fmol).
C ALGBRC(30) is q_RB1 in component environment (fmol).
C ALGBRC(32) is q_Gs in component environment (fmol).
C ALGBRC(34) is q_RB1_Gs in component environment (fmol).
C ALGBRC(36) is q_L_RB1 in component environment (fmol).
C ALGBRC(38) is q_L_RB1_Gs in component environment (fmol).
C ALGBRC(40) is q_Gsbetagamma in component environment (fmol).
C ALGBRC(42) is q_Gsa_GDP in component environment (fmol).
C ALGBRC(45) is q_GTP in component environment (fmol).
C ALGBRC(47) is q_GDP in component environment (fmol).
C ALGBRC(49) is q_Pi in component environment (fmol).
C ALGBRC(51) is q_RB1_tag in component environment (fmol).
C ALGBRC(53) is q_L_RB1_tag in component environment (fmol).
C ALGBRC(55) is q_RB1_GRKArr in component environment (fmol).
C ALGBRC(57) is q_L_RB1_GRKArr in component environment (fmol).
C ALGBRC(61) is q_GRKArr in component environment (fmol).
C ALGBRC(63) is q_RM2_inactive in component environment (fmol).
C ALGBRC(66) is q_L_RM2_inactive in component environment (fmol).
C ALGBRC(69) is q_LM2 in component environment (fmol).
C ALGBRC(72) is q_RM2 in component environment (fmol).
C ALGBRC(75) is q_Gi in component environment (fmol).
C ALGBRC(78) is q_RM2_Gi in component environment (fmol).
C ALGBRC(81) is q_L_RM2 in component environment (fmol).
C ALGBRC(86) is q_L_RM2_Gi in component environment (fmol).
C ALGBRC(90) is q_Gibetagamma in component environment (fmol).
C ALGBRC(94) is q_Gia_GDP in component environment (fmol).
C ALGBRC(98) is q_RM2_tag in component environment (fmol).
C ALGBRC(100) is q_L_RM2_tag in component environment (fmol).
C ALGBRC(102) is q_RM2_GRKArr in component environment (fmol).
C ALGBRC(104) is q_L_RM2_GRKArr in component environment (fmol).
C STATES(1) is q_ATP in component cAMP (fmol).
C STATES(2) is q_cAMP in component cAMP (fmol).
C STATES(3) is q_AC in component cAMP (fmol).
C STATES(4) is q_AC_ATP in component cAMP (fmol).
C STATES(5) is q_Gsa_GTP_AC in component cAMP (fmol).
C STATES(6) is q_Gsa_GTP_AC_ATP in component cAMP (fmol).
C STATES(7) is q_FSK_AC in component cAMP (fmol).
C STATES(8) is q_FSK_AC_ATP in component cAMP (fmol).
C STATES(9) is q_PDE in component cAMP (fmol).
C STATES(10) is q_PDE_cAMP in component cAMP (fmol).
C STATES(11) is q_five_AMP in component cAMP (fmol).
C STATES(12) is q_IBMX in component cAMP (fmol).
C STATES(13) is q_PDEinh in component cAMP (fmol).
C STATES(14) is q_Gsa_GTP in component cAMP (fmol).
C STATES(15) is q_FSK in component cAMP (fmol).
C STATES(16) is q_Gia_GTP in component cAMP (fmol).
C STATES(17) is q_ACinh in component cAMP (fmol).
C STATES(18) is q_PPi in component cAMP (fmol).
C STATES(19) is q_RB1_inactive in component GPCR_B1AR_reduced (fmol).
C STATES(20) is q_L_RB1_inactive in component GPCR_B1AR_reduced (fmol).
C STATES(21) is q_LB1 in component GPCR_B1AR_reduced (fmol).
C STATES(22) is q_RB1 in component GPCR_B1AR_reduced (fmol).
C STATES(23) is q_Gs in component GPCR_B1AR_reduced (fmol).
C STATES(24) is q_RB1_Gs in component GPCR_B1AR_reduced (fmol).
C STATES(25) is q_L_RB1 in component GPCR_B1AR_reduced (fmol).
C STATES(26) is q_L_RB1_Gs in component GPCR_B1AR_reduced (fmol).
C STATES(27) is q_Gsa_GTP in component GPCR_B1AR_reduced (fmol).
C STATES(28) is q_Gsbetagamma in component GPCR_B1AR_reduced (fmol).
C STATES(29) is q_Gsa_GDP in component GPCR_B1AR_reduced (fmol).
C STATES(30) is q_GTP in component GPCR_B1AR_reduced (fmol).
C STATES(31) is q_GDP in component GPCR_B1AR_reduced (fmol).
C STATES(32) is q_Pi in component GPCR_B1AR_reduced (fmol).
C STATES(33) is q_RB1_tag in component GPCR_B1AR_reduced (fmol).
C STATES(34) is q_L_RB1_tag in component GPCR_B1AR_reduced (fmol).
C STATES(35) is q_RB1_GRKArr in component GPCR_B1AR_reduced (fmol).
C STATES(36) is q_L_RB1_GRKArr in component GPCR_B1AR_reduced (fmol).
C STATES(37) is q_GRKArr in component GPCR_B1AR_reduced (fmol).
C STATES(38) is q_RM2_inactive in component GPCR_M2_reduced (fmol).
C STATES(39) is q_L_RM2_inactive in component GPCR_M2_reduced (fmol).
C STATES(40) is q_LM2 in component GPCR_M2_reduced (fmol).
C STATES(41) is q_RM2 in component GPCR_M2_reduced (fmol).
C STATES(42) is q_Gi in component GPCR_M2_reduced (fmol).
C STATES(43) is q_RM2_Gi in component GPCR_M2_reduced (fmol).
C STATES(44) is q_L_RM2 in component GPCR_M2_reduced (fmol).
C STATES(45) is q_L_RM2_Gi in component GPCR_M2_reduced (fmol).
C STATES(46) is q_Gia_GTP in component GPCR_M2_reduced (fmol).
C STATES(47) is q_Gibetagamma in component GPCR_M2_reduced (fmol).
C STATES(48) is q_Gia_GDP in component GPCR_M2_reduced (fmol).
C STATES(49) is q_GTP in component GPCR_M2_reduced (fmol).
C STATES(50) is q_GDP in component GPCR_M2_reduced (fmol).
C STATES(51) is q_Pi in component GPCR_M2_reduced (fmol).
C STATES(52) is q_RM2_tag in component GPCR_M2_reduced (fmol).
C STATES(53) is q_L_RM2_tag in component GPCR_M2_reduced (fmol).
C STATES(54) is q_RM2_GRKArr in component GPCR_M2_reduced (fmol).
C STATES(55) is q_L_RM2_GRKArr in component GPCR_M2_reduced (fmol).
C STATES(56) is q_GRKArr in component GPCR_M2_reduced (fmol).
C CONSTS(145) is R in component constants (J_per_K_per_mol).
C CONSTS(146) is T in component constants (kelvin).
C CONSTS(147) is F in component constants (C_per_mol).
C ALGBRC(70) is v1a in component cAMP (fmol_per_sec).
C ALGBRC(73) is v1b in component cAMP (fmol_per_sec).
C ALGBRC(76) is v2a in component cAMP (fmol_per_sec).
C ALGBRC(79) is v2b in component cAMP (fmol_per_sec).
C ALGBRC(82) is v3a in component cAMP (fmol_per_sec).
C ALGBRC(84) is v3b in component cAMP (fmol_per_sec).
C ALGBRC(87) is v4a in component cAMP (fmol_per_sec).
C ALGBRC(91) is v4b in component cAMP (fmol_per_sec).
C ALGBRC(95) is v5 in component cAMP (fmol_per_sec).
C ALGBRC(88) is v6 in component cAMP (fmol_per_sec).
C ALGBRC(92) is v7 in component cAMP (fmol_per_sec).
C ALGBRC(96) is vGiAC in component cAMP (fmol_per_sec).
C ALGBRC(29) is mu_ATP in component cAMP (J_per_mol).
C ALGBRC(33) is mu_AC in component cAMP (J_per_mol).
C ALGBRC(31) is mu_cAMP in component cAMP (J_per_mol).
C ALGBRC(35) is mu_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(58) is mu_FSK in component cAMP (J_per_mol).
C ALGBRC(41) is mu_FSK_AC in component cAMP (J_per_mol).
C ALGBRC(43) is mu_FSK_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(56) is mu_Gsa_GTP in component cAMP (J_per_mol).
C ALGBRC(37) is mu_Gsa_GTP_AC in component cAMP (J_per_mol).
C ALGBRC(39) is mu_Gsa_GTP_AC_ATP in component cAMP (J_per_mol).
C ALGBRC(46) is mu_PDE in component cAMP (J_per_mol).
C ALGBRC(54) is mu_PDEinh in component cAMP (J_per_mol).
C ALGBRC(48) is mu_PDE_cAMP in component cAMP (J_per_mol).
C ALGBRC(52) is mu_IBMX in component cAMP (J_per_mol).
C ALGBRC(50) is mu_five_AMP in component cAMP (J_per_mol).
C ALGBRC(62) is mu_Gia_GTP in component cAMP (J_per_mol).
C ALGBRC(64) is mu_ACinh in component cAMP (J_per_mol).
C ALGBRC(67) is mu_PPi in component cAMP (J_per_mol).
C CONSTS(148) is vol in component cAMP (pL).
C ALGBRC(12) is ATP_T in component cAMP (fmol).
C ALGBRC(13) is AC_T in component cAMP (fmol).
C ALGBRC(22) is Gs_T in component cAMP (fmol).
C ALGBRC(16) is cAMP_T in component cAMP (fmol).
C ALGBRC(65) is mu_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(68) is mu_L_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(71) is mu_LB1 in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(74) is mu_RB1 in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(77) is mu_Gs in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(80) is mu_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(83) is mu_L_RB1 in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(85) is mu_L_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(89) is mu_Gsa_GTP in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(93) is mu_Gsbetagamma in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(97) is mu_Gsa_GDP in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(99) is mu_GTP in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(101) is mu_GDP in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(103) is mu_Pi in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(105) is mu_RB1_tag in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(107) is mu_L_RB1_tag in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(109) is mu_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(111) is mu_L_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(113) is mu_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
C ALGBRC(115) is v_Rswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(117) is v_LRswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(119) is v_C_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(122) is v_R_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(120) is v_L_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(125) is v_Act1_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(127) is v_Act2_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(129) is v_Hyd_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(132) is v_Reassoc_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(130) is v_InternR_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(133) is v_InternLR_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(123) is v_Rsynthesis in component GPCR_B1AR_reduced (fmol_per_sec).
C ALGBRC(106) is mu_RM2_inactive in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(108) is mu_L_RM2_inactive in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(110) is mu_LM2 in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(112) is mu_RM2 in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(114) is mu_Gi in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(116) is mu_RM2_Gi in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(118) is mu_L_RM2 in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(121) is mu_L_RM2_Gi in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(124) is mu_Gia_GTP in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(126) is mu_Gibetagamma in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(128) is mu_Gia_GDP in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(131) is mu_GTP in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(134) is mu_GDP in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(135) is mu_Pi in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(136) is mu_RM2_tag in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(137) is mu_L_RM2_tag in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(138) is mu_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(139) is mu_L_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(140) is mu_GRKArr in component GPCR_M2_reduced (J_per_mol).
C ALGBRC(141) is v_Rswitch_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(142) is v_LRswitch_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(143) is v_C_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(144) is v_R_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(145) is v_L_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(146) is v_Act1_Gi in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(147) is v_Act2_Gi in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(149) is v_Hyd_Gi in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(151) is v_Reassoc_Gi in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(150) is v_InternR_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(152) is v_InternLR_M2 in component GPCR_M2_reduced (fmol_per_sec).
C ALGBRC(148) is v_Rsynthesis in component GPCR_M2_reduced (fmol_per_sec).
C RATES(1) is d/dt q_ATP in component cAMP (fmol).
C RATES(3) is d/dt q_AC in component cAMP (fmol).
C RATES(4) is d/dt q_AC_ATP in component cAMP (fmol).
C RATES(2) is d/dt q_cAMP in component cAMP (fmol).
C RATES(15) is d/dt q_FSK in component cAMP (fmol).
C RATES(7) is d/dt q_FSK_AC in component cAMP (fmol).
C RATES(8) is d/dt q_FSK_AC_ATP in component cAMP (fmol).
C RATES(14) is d/dt q_Gsa_GTP in component cAMP (fmol).
C RATES(5) is d/dt q_Gsa_GTP_AC in component cAMP (fmol).
C RATES(6) is d/dt q_Gsa_GTP_AC_ATP in component cAMP (fmol).
C RATES(10) is d/dt q_PDE_cAMP in component cAMP (fmol).
C RATES(9) is d/dt q_PDE in component cAMP (fmol).
C RATES(12) is d/dt q_IBMX in component cAMP (fmol).
C RATES(13) is d/dt q_PDEinh in component cAMP (fmol).
C RATES(11) is d/dt q_five_AMP in component cAMP (fmol).
C RATES(16) is d/dt q_Gia_GTP in component cAMP (fmol).
C RATES(17) is d/dt q_ACinh in component cAMP (fmol).
C RATES(18) is d/dt q_PPi in component cAMP (fmol).
C RATES(19) is d/dt q_RB1_inactive in component GPCR_B1AR_reduced (fmol).
C RATES(20) is d/dt q_L_RB1_inactive in component GPCR_B1AR_reduced (fmol).
C RATES(21) is d/dt q_LB1 in component GPCR_B1AR_reduced (fmol).
C RATES(22) is d/dt q_RB1 in component GPCR_B1AR_reduced (fmol).
C RATES(23) is d/dt q_Gs in component GPCR_B1AR_reduced (fmol).
C RATES(24) is d/dt q_RB1_Gs in component GPCR_B1AR_reduced (fmol).
C RATES(25) is d/dt q_L_RB1 in component GPCR_B1AR_reduced (fmol).
C RATES(26) is d/dt q_L_RB1_Gs in component GPCR_B1AR_reduced (fmol).
C RATES(27) is d/dt q_Gsa_GTP in component GPCR_B1AR_reduced (fmol).
C RATES(28) is d/dt q_Gsbetagamma in component GPCR_B1AR_reduced (fmol).
C RATES(29) is d/dt q_Gsa_GDP in component GPCR_B1AR_reduced (fmol).
C RATES(30) is d/dt q_GTP in component GPCR_B1AR_reduced (fmol).
C RATES(31) is d/dt q_GDP in component GPCR_B1AR_reduced (fmol).
C RATES(32) is d/dt q_Pi in component GPCR_B1AR_reduced (fmol).
C RATES(33) is d/dt q_RB1_tag in component GPCR_B1AR_reduced (fmol).
C RATES(34) is d/dt q_L_RB1_tag in component GPCR_B1AR_reduced (fmol).
C RATES(35) is d/dt q_RB1_GRKArr in component GPCR_B1AR_reduced (fmol).
C RATES(36) is d/dt q_L_RB1_GRKArr in component GPCR_B1AR_reduced (fmol).
C RATES(37) is d/dt q_GRKArr in component GPCR_B1AR_reduced (fmol).
C RATES(38) is d/dt q_RM2_inactive in component GPCR_M2_reduced (fmol).
C RATES(39) is d/dt q_L_RM2_inactive in component GPCR_M2_reduced (fmol).
C RATES(40) is d/dt q_LM2 in component GPCR_M2_reduced (fmol).
C RATES(41) is d/dt q_RM2 in component GPCR_M2_reduced (fmol).
C RATES(42) is d/dt q_Gi in component GPCR_M2_reduced (fmol).
C RATES(43) is d/dt q_RM2_Gi in component GPCR_M2_reduced (fmol).
C RATES(44) is d/dt q_L_RM2 in component GPCR_M2_reduced (fmol).
C RATES(45) is d/dt q_L_RM2_Gi in component GPCR_M2_reduced (fmol).
C RATES(46) is d/dt q_Gia_GTP in component GPCR_M2_reduced (fmol).
C RATES(47) is d/dt q_Gibetagamma in component GPCR_M2_reduced (fmol).
C RATES(48) is d/dt q_Gia_GDP in component GPCR_M2_reduced (fmol).
C RATES(49) is d/dt q_GTP in component GPCR_M2_reduced (fmol).
C RATES(50) is d/dt q_GDP in component GPCR_M2_reduced (fmol).
C RATES(51) is d/dt q_Pi in component GPCR_M2_reduced (fmol).
C RATES(52) is d/dt q_RM2_tag in component GPCR_M2_reduced (fmol).
C RATES(53) is d/dt q_L_RM2_tag in component GPCR_M2_reduced (fmol).
C RATES(54) is d/dt q_RM2_GRKArr in component GPCR_M2_reduced (fmol).
C RATES(55) is d/dt q_L_RM2_GRKArr in component GPCR_M2_reduced (fmol).
C RATES(56) is d/dt q_GRKArr in component GPCR_M2_reduced (fmol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 7.05076e+06
      CONSTS(2) = 0.00147213
      CONSTS(3) = 186465
      CONSTS(4) = 0.0503158
      CONSTS(5) = 2.50455e+08
      CONSTS(6) = 2.91227e-17
      CONSTS(7) = 30352.7
      CONSTS(8) = 0.116741
      CONSTS(9) = 833.759
      CONSTS(10) = 4236.98
      CONSTS(11) = 569102
      CONSTS(12) = 1052.47
      CONSTS(13) = 93.1496
      CONSTS(14) = 0.00202063
      CONSTS(15) = 1459.71
      CONSTS(16) = 31664.5
      CONSTS(17) = 0.00057588
      CONSTS(18) = 0.00326844
      CONSTS(19) = 37.7371
      CONSTS(20) = 0.0787034
      CONSTS(21) = 0.00156706
      CONSTS(22) = 0.00130092
      CONSTS(23) = 15.0203
      CONSTS(24) = 271.765
      CONSTS(25) = 0.0101296
      CONSTS(26) = 7036.75
      CONSTS(27) = 262283
      CONSTS(28) = 0.111426
      CONSTS(29) = 0.00270806
      CONSTS(30) = 1.32138
      CONSTS(31) = 0.48875
      CONSTS(32) = 0.00258927
      CONSTS(33) = 0.000652343
      CONSTS(34) = 0.318306
      CONSTS(35) = 6.60544e-05
      CONSTS(36) = 0.0166528
      CONSTS(37) = 1.68744
      CONSTS(38) = 3.94935
      CONSTS(39) = 6.86096
      CONSTS(40) = 4.91084
      CONSTS(41) = 0.0510801
      CONSTS(42) = 0.0998184
      CONSTS(43) = 1.67186
      CONSTS(44) = 1.24505
      CONSTS(45) = 0.0166528
      CONSTS(46) = 0.020208
      CONSTS(47) = 34.8659
      CONSTS(48) = 0.295487
      CONSTS(49) = 1.99992e-05
      CONSTS(50) = 0.0475823
      CONSTS(51) = 27.6205
      CONSTS(52) = 0.000115307
      CONSTS(53) = 3.12076
      CONSTS(54) = 143865
      CONSTS(55) = 4702.09
      CONSTS(56) = 0.0312076
      CONSTS(57) = 0.185505
      CONSTS(58) = 6.57187
      CONSTS(59) = 0.00143865
      CONSTS(60) = 0.000569194
      CONSTS(61) = 9.19376
      CONSTS(62) = 70.9724
      CONSTS(63) = 0.629469
      CONSTS(64) = 5.92518e-07
      CONSTS(65) = 1.51286e-06
      CONSTS(66) = 1.35735
      CONSTS(67) = 0.000117561
      CONSTS(68) = 0.00223456
      CONSTS(69) = 1.93537e-07
      CONSTS(70) = 4785.64
      CONSTS(71) = 1.06966
      CONSTS(72) = 28697.8
      CONSTS(73) = 70.9007
      CONSTS(74) = 0.0106966
      CONSTS(75) = 0.11227
      CONSTS(76) = 1.23934
      CONSTS(77) = 0.000286978
      CONSTS(78) = 0.00253994
      CONSTS(79) = 34.5537
      CONSTS(80) = 11.4287
      CONSTS(81) = 2.70687
      CONSTS(82) = 0.00554752
      CONSTS(83) = 0.00445621
      CONSTS(84) = 9.13266e-06
      CONSTS(85) = 34.4
      CONSTS(86) = 500
      CONSTS(87) = 3.5e-4
      CONSTS(88) = 0.25e-4
      CONSTS(89) = 1.8e-4
      CONSTS(90) = 1e1
      CONSTS(91) = 1e-5
      CONSTS(92) = 7e-4
      CONSTS(93) = 0.25e-4
      CONSTS(94) = 1.8e-4
      CONSTS(95) = 1e1
      CONSTS(96) = 1e-5
      CONSTS(97) = 190
      CONSTS(98) = 1.889E-03
      CONSTS(99) = 1e-18
      CONSTS(100) = 1e-18
      CONSTS(101) = 3.800E-05
      CONSTS(102) = 1e-18
      CONSTS(103) = 1e-18
      CONSTS(104) = 1e-18
      CONSTS(105) = 1e-18
      CONSTS(106) = 1e-18
      CONSTS(107) = 1.482E-03
      CONSTS(108) = 1e-18
      CONSTS(109) = 1e-18
      CONSTS(110) = 3.80E-02
      CONSTS(111) = 1e-18
      CONSTS(112) = 1e-18
      CONSTS(113) = 1e-18
      CONSTS(114) = 1e-18
      CONSTS(115) = 0.0004579000
      CONSTS(116) = 1e-18
      CONSTS(117) = 1e-18
      CONSTS(118) = 0.1455400000
      CONSTS(119) = 1e-18
      CONSTS(120) = 1e-18
      CONSTS(121) = 1e-18
      CONSTS(122) = 1e-18
      CONSTS(123) = 1e-18
      CONSTS(124) = 2.2
      CONSTS(125) = 1.1
      CONSTS(126) = 570
      CONSTS(127) = 1e-18
      CONSTS(128) = 1e-18
      CONSTS(129) = 1e-18
      CONSTS(130) = 1e-18
      CONSTS(131) = 1e-8
      CONSTS(132) = 0.00072
      CONSTS(133) = 1e-18
      CONSTS(134) = 1e-18
      CONSTS(135) = 0.00836
      CONSTS(136) = 1e-18
      CONSTS(137) = 1e-18
      CONSTS(138) = 1e-18
      CONSTS(139) = 1e-18
      CONSTS(140) = 1e-18
      CONSTS(141) = 1e-18
      CONSTS(142) = 1e-18
      CONSTS(143) = 1e-18
      CONSTS(144) = 1e-18
      STATES(1) = 1e-16
      STATES(2) = 1e-16
      STATES(3) = 1e-16
      STATES(4) = 1e-16
      STATES(5) = 1e-16
      STATES(6) = 1e-16
      STATES(7) = 1e-16
      STATES(8) = 1e-16
      STATES(9) = 1e-16
      STATES(10) = 1e-16
      STATES(11) = 1e-16
      STATES(12) = 1e-16
      STATES(13) = 1e-16
      STATES(14) = 1e-16
      STATES(15) = 1e-16
      STATES(16) = 1e-16
      STATES(17) = 1e-16
      STATES(18) = 1e-16
      STATES(19) = 1e-16
      STATES(20) = 1e-16
      STATES(21) = 1e-16
      STATES(22) = 1e-16
      STATES(23) = 1e-16
      STATES(24) = 1e-16
      STATES(25) = 1e-16
      STATES(26) = 1e-16
      STATES(27) = 1e-16
      STATES(28) = 1e-16
      STATES(29) = 1e-16
      STATES(30) = 1e-16
      STATES(31) = 1e-16
      STATES(32) = 1e-16
      STATES(33) = 1e-16
      STATES(34) = 1e-16
      STATES(35) = 1e-16
      STATES(36) = 1e-16
      STATES(37) = 1e-16
      STATES(38) = 1e-16
      STATES(39) = 1e-16
      STATES(40) = 1e-16
      STATES(41) = 1e-16
      STATES(42) = 1e-16
      STATES(43) = 1e-16
      STATES(44) = 1e-16
      STATES(45) = 1e-16
      STATES(46) = 1e-16
      STATES(47) = 1e-16
      STATES(48) = 1e-16
      STATES(49) = 1e-16
      STATES(50) = 1e-16
      STATES(51) = 1e-16
      STATES(52) = 1e-16
      STATES(53) = 1e-16
      STATES(54) = 1e-16
      STATES(55) = 1e-16
      STATES(56) = 1e-16
      CONSTS(145) = 8.31
      CONSTS(146) = 310
      CONSTS(147) = 96485
      CONSTS(148) = 38.0
      CONSTS(149) = CONSTS(90)/CONSTS(89)
      CONSTS(150) = CONSTS(95)/CONSTS(94)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = CONSTS(97)+STATES(1)
      ALGBRC(29) =  CONSTS(145)*CONSTS(146)*log( CONSTS(35)*ALGBRC(3))
      ALGBRC(5) = CONSTS(98)+STATES(3)
      ALGBRC(33) =  CONSTS(145)*CONSTS(146)*log( CONSTS(37)*ALGBRC(5))
      ALGBRC(6) = CONSTS(100)+STATES(4)
      ALGBRC(35) =  CONSTS(145)*CONSTS(146)*log( CONSTS(38)*ALGBRC(6))
      ALGBRC(70) =  CONSTS(1)*(EXP((ALGBRC(33)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(35)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(4) = CONSTS(99)+STATES(2)
      ALGBRC(31) =  CONSTS(145)*CONSTS(146)*log( CONSTS(36)*ALGBRC(4))
      ALGBRC(26) = CONSTS(114)+STATES(18)
      ALGBRC(67) =  CONSTS(145)*CONSTS(146)*log( CONSTS(52)*ALGBRC(26))
      ALGBRC(73) =  CONSTS(2)*(EXP(ALGBRC(35)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(33)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      RATES(4) = ALGBRC(70) - ALGBRC(73)
      ALGBRC(7) = CONSTS(105)+STATES(5)
      ALGBRC(37) =  CONSTS(145)*CONSTS(146)*log( CONSTS(39)*ALGBRC(7))
      ALGBRC(8) = CONSTS(106)+STATES(6)
      ALGBRC(39) =  CONSTS(145)*CONSTS(146)*log( CONSTS(40)*ALGBRC(8))
      ALGBRC(76) =  CONSTS(3)*(EXP((ALGBRC(37)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(39)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(79) =  CONSTS(4)*(EXP(ALGBRC(39)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(37)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      RATES(6) = ALGBRC(76) - ALGBRC(79)
      ALGBRC(9) = CONSTS(102)+STATES(7)
      ALGBRC(41) =  CONSTS(145)*CONSTS(146)*log( CONSTS(41)*ALGBRC(9))
      ALGBRC(10) = CONSTS(103)+STATES(8)
      ALGBRC(43) =  CONSTS(145)*CONSTS(146)*log( CONSTS(42)*ALGBRC(10))
      ALGBRC(82) =  CONSTS(5)*(EXP((ALGBRC(41)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(43)/( CONSTS(145)*CONSTS(146))))
      RATES(1) = (- ALGBRC(70) - ALGBRC(82)) - ALGBRC(76)
      ALGBRC(84) =  CONSTS(6)*(EXP(ALGBRC(43)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(41)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      RATES(8) = ALGBRC(82) - ALGBRC(84)
      ALGBRC(11) = CONSTS(107)+STATES(9)
      ALGBRC(46) =  CONSTS(145)*CONSTS(146)*log( CONSTS(43)*ALGBRC(11))
      ALGBRC(14) = CONSTS(109)+STATES(10)
      ALGBRC(48) =  CONSTS(145)*CONSTS(146)*log( CONSTS(44)*ALGBRC(14))
      ALGBRC(87) =  CONSTS(7)*(EXP((ALGBRC(46)+ALGBRC(31))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(48)/( CONSTS(145)*CONSTS(146))))
      RATES(2) = (ALGBRC(73)+ALGBRC(84)+ALGBRC(79)) - ALGBRC(87)
      ALGBRC(20) = CONSTS(104)+STATES(14)+STATES(27)
      ALGBRC(56) =  CONSTS(145)*CONSTS(146)*log( CONSTS(48)*ALGBRC(20))
      ALGBRC(88) =  CONSTS(10)*(EXP((ALGBRC(33)+ALGBRC(56))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(37)/( CONSTS(145)*CONSTS(146))))
      RATES(14) = - ALGBRC(88)
      RATES(5) = (ALGBRC(88) - ALGBRC(76))+ALGBRC(79)
      ALGBRC(21) = CONSTS(101)+STATES(15)
      ALGBRC(58) =  CONSTS(145)*CONSTS(146)*log( CONSTS(49)*ALGBRC(21))
      ALGBRC(92) =  CONSTS(11)*(EXP((ALGBRC(58)+ALGBRC(33))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(41)/( CONSTS(145)*CONSTS(146))))
      RATES(15) = - ALGBRC(92)
      RATES(7) = (ALGBRC(92)+ALGBRC(84)) - ALGBRC(82)
      ALGBRC(15) = CONSTS(111)+STATES(11)
      ALGBRC(50) =  CONSTS(145)*CONSTS(146)*log( CONSTS(45)*ALGBRC(15))
      ALGBRC(91) =  CONSTS(8)*(EXP(ALGBRC(48)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(46)+ALGBRC(50))/( CONSTS(145)*CONSTS(146))))
      RATES(10) = ALGBRC(87) - ALGBRC(91)
      RATES(11) = ALGBRC(91)
      ALGBRC(23) = CONSTS(112)+STATES(16)+STATES(46)
      ALGBRC(62) =  CONSTS(145)*CONSTS(146)*log( CONSTS(50)*ALGBRC(23))
      ALGBRC(24) = CONSTS(113)+STATES(17)
      ALGBRC(64) =  CONSTS(145)*CONSTS(146)*log( CONSTS(51)*ALGBRC(24))
      ALGBRC(96) =  CONSTS(12)*(EXP((ALGBRC(33)+ALGBRC(62))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(64)/( CONSTS(145)*CONSTS(146))))
      RATES(3) = (((ALGBRC(73) - ALGBRC(70)) - ALGBRC(88)) - ALGBRC(92)) - ALGBRC(96)
      ALGBRC(19) = CONSTS(108)+STATES(13)
      ALGBRC(54) =  CONSTS(145)*CONSTS(146)*log( CONSTS(47)*ALGBRC(19))
      ALGBRC(18) = CONSTS(110)+STATES(12)
      ALGBRC(52) =  CONSTS(145)*CONSTS(146)*log( CONSTS(46)*ALGBRC(18))
      ALGBRC(95) =  CONSTS(9)*(EXP((ALGBRC(46)+ALGBRC(52))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(54)/( CONSTS(145)*CONSTS(146))))
      RATES(9) = (ALGBRC(91) - ALGBRC(87)) - ALGBRC(95)
      RATES(12) = - ALGBRC(95)
      RATES(13) = ALGBRC(95)
      RATES(16) = - ALGBRC(96)
      RATES(17) = ALGBRC(96)
      RATES(18) = ALGBRC(96)
      ALGBRC(27) = CONSTS(116)+STATES(20)
      ALGBRC(68) =  CONSTS(145)*CONSTS(146)*log( CONSTS(54)*ALGBRC(27))
      ALGBRC(36) = CONSTS(120)+STATES(25)
      ALGBRC(83) =  CONSTS(145)*CONSTS(146)*log( CONSTS(59)*ALGBRC(36))
      ALGBRC(117) =  CONSTS(14)*(EXP(ALGBRC(68)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(83)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(25) = CONSTS(115)+STATES(19)
      ALGBRC(65) =  CONSTS(145)*CONSTS(146)*log( CONSTS(53)*ALGBRC(25))
      ALGBRC(1) = TERNRY(VOI.LT.CONSTS(87).AND.VOI.GT.CONSTS(87) - CONSTS(89), CONSTS(91)+ CONSTS(149)*((VOI - CONSTS(87))+CONSTS(89)), TERNRY(VOI.GE.CONSTS(87).AND.VOI.LT.CONSTS(87)+CONSTS(88), CONSTS(90)+CONSTS(91), TERNRY(VOI.LT.CONSTS(87)+CONSTS(89)+CONSTS(88).AND.VOI.GE.CONSTS(87)+CONSTS(88), CONSTS(91)+ - CONSTS(149)*(((VOI - CONSTS(87)) - CONSTS(89)) - CONSTS(88)), CONSTS(91))
      ALGBRC(28) = ALGBRC(1)+STATES(21)
      ALGBRC(71) =  CONSTS(145)*CONSTS(146)*log( CONSTS(55)*ALGBRC(28))
      ALGBRC(120) =  CONSTS(17)*(EXP((ALGBRC(65)+ALGBRC(71))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(68)/( CONSTS(145)*CONSTS(146))))
      RATES(20) = - ALGBRC(117)+ALGBRC(120)
      RATES(21) = - ALGBRC(120)
      ALGBRC(30) = CONSTS(117)+STATES(22)
      ALGBRC(74) =  CONSTS(145)*CONSTS(146)*log( CONSTS(56)*ALGBRC(30))
      ALGBRC(115) =  CONSTS(13)*(EXP(ALGBRC(65)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(74)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(32) = CONSTS(118)+STATES(23)
      ALGBRC(77) =  CONSTS(145)*CONSTS(146)*log( CONSTS(57)*ALGBRC(32))
      ALGBRC(34) = CONSTS(119)+STATES(24)
      ALGBRC(80) =  CONSTS(145)*CONSTS(146)*log( CONSTS(58)*ALGBRC(34))
      ALGBRC(119) =  CONSTS(15)*(EXP((ALGBRC(74)+ALGBRC(77))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(80)/( CONSTS(145)*CONSTS(146))))
      RATES(22) = ALGBRC(115) - ALGBRC(119)
      ALGBRC(123) = ALGBRC(115)+ALGBRC(120)
      RATES(19) = (- ALGBRC(115) - ALGBRC(120))+ 1.00000*ALGBRC(123)
      ALGBRC(38) = CONSTS(121)+STATES(26)
      ALGBRC(85) =  CONSTS(145)*CONSTS(146)*log( CONSTS(60)*ALGBRC(38))
      ALGBRC(122) =  CONSTS(16)*(EXP((ALGBRC(83)+ALGBRC(77))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(85)/( CONSTS(145)*CONSTS(146))))
      RATES(25) = ALGBRC(117) - ALGBRC(122)
      ALGBRC(89) =  CONSTS(145)*CONSTS(146)*log( CONSTS(48)*ALGBRC(20))
      ALGBRC(40) = CONSTS(122)+STATES(28)
      ALGBRC(93) =  CONSTS(145)*CONSTS(146)*log( CONSTS(61)*ALGBRC(40))
      ALGBRC(45) = CONSTS(124)+STATES(30)+STATES(49)
      ALGBRC(99) =  CONSTS(145)*CONSTS(146)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(47) = CONSTS(125)+STATES(31)+STATES(50)
      ALGBRC(101) =  CONSTS(145)*CONSTS(146)*log( CONSTS(64)*ALGBRC(47))
      ALGBRC(51) = CONSTS(127)+STATES(33)
      ALGBRC(105) =  CONSTS(145)*CONSTS(146)*log( CONSTS(66)*ALGBRC(51))
      ALGBRC(125) =  CONSTS(18)*(EXP((ALGBRC(80)+ALGBRC(99))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(89)+ALGBRC(93)+ALGBRC(105)+ALGBRC(101))/( CONSTS(145)*CONSTS(146))))
      RATES(24) = ALGBRC(119) - ALGBRC(125)
      ALGBRC(53) = CONSTS(128)+STATES(34)
      ALGBRC(107) =  CONSTS(145)*CONSTS(146)*log( CONSTS(67)*ALGBRC(53))
      ALGBRC(127) =  CONSTS(19)*(EXP((ALGBRC(85)+ALGBRC(99))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(89)+ALGBRC(93)+ALGBRC(107)+ALGBRC(101))/( CONSTS(145)*CONSTS(146))))
      RATES(26) = ALGBRC(122) - ALGBRC(127)
      RATES(30) = - ALGBRC(125) - ALGBRC(127)
      RATES(31) = ALGBRC(125)+ALGBRC(127)
      ALGBRC(42) = CONSTS(123)+STATES(29)
      ALGBRC(97) =  CONSTS(145)*CONSTS(146)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(49) = CONSTS(126)+STATES(32)+STATES(51)
      ALGBRC(103) =  CONSTS(145)*CONSTS(146)*log( CONSTS(65)*ALGBRC(49))
      ALGBRC(129) =  CONSTS(20)*(EXP(ALGBRC(89)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(97)+ALGBRC(103))/( CONSTS(145)*CONSTS(146))))
      RATES(27) = (ALGBRC(125)+ALGBRC(127)) - ALGBRC(129)
      RATES(32) = ALGBRC(129)
      ALGBRC(55) = CONSTS(129)+STATES(35)
      ALGBRC(109) =  CONSTS(145)*CONSTS(146)*log( CONSTS(68)*ALGBRC(55))
      ALGBRC(61) = CONSTS(131)+STATES(37)+STATES(56)
      ALGBRC(113) =  CONSTS(145)*CONSTS(146)*log( CONSTS(70)*ALGBRC(61))
      ALGBRC(130) =  CONSTS(22)*(EXP((ALGBRC(105)+ALGBRC(113))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(109)/( CONSTS(145)*CONSTS(146))))
      RATES(33) = ALGBRC(125) - ALGBRC(130)
      RATES(35) = ALGBRC(130)
      ALGBRC(132) =  CONSTS(21)*(EXP((ALGBRC(97)+ALGBRC(93))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(77)/( CONSTS(145)*CONSTS(146))))
      RATES(23) = (- ALGBRC(119) - ALGBRC(122))+ALGBRC(132)
      RATES(28) = (ALGBRC(125)+ALGBRC(127)) - ALGBRC(132)
      RATES(29) = ALGBRC(129) - ALGBRC(132)
      ALGBRC(57) = CONSTS(130)+STATES(36)
      ALGBRC(111) =  CONSTS(145)*CONSTS(146)*log( CONSTS(69)*ALGBRC(57))
      ALGBRC(133) =  CONSTS(23)*(EXP((ALGBRC(107)+ALGBRC(113))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(111)/( CONSTS(145)*CONSTS(146))))
      RATES(34) = ALGBRC(127) - ALGBRC(133)
      RATES(36) = ALGBRC(133)
      RATES(37) = - ALGBRC(130) - ALGBRC(133)
      ALGBRC(66) = CONSTS(133)+STATES(39)
      ALGBRC(108) =  CONSTS(145)*CONSTS(146)*log( CONSTS(72)*ALGBRC(66))
      ALGBRC(81) = CONSTS(137)+STATES(44)
      ALGBRC(118) =  CONSTS(145)*CONSTS(146)*log( CONSTS(77)*ALGBRC(81))
      ALGBRC(142) =  CONSTS(25)*(EXP(ALGBRC(108)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(118)/( CONSTS(145)*CONSTS(146))))
      RATES(39) = - ALGBRC(142)
      ALGBRC(63) = CONSTS(132)+STATES(38)
      ALGBRC(106) =  CONSTS(145)*CONSTS(146)*log( CONSTS(71)*ALGBRC(63))
      ALGBRC(72) = CONSTS(134)+STATES(41)
      ALGBRC(112) =  CONSTS(145)*CONSTS(146)*log( CONSTS(74)*ALGBRC(72))
      ALGBRC(141) =  CONSTS(24)*(EXP(ALGBRC(106)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(112)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(92).AND.VOI.GT.CONSTS(92) - CONSTS(94), CONSTS(96)+ CONSTS(150)*((VOI - CONSTS(92))+CONSTS(94)), TERNRY(VOI.GE.CONSTS(92).AND.VOI.LT.CONSTS(92)+CONSTS(93), CONSTS(95)+CONSTS(96), TERNRY(VOI.LT.CONSTS(92)+CONSTS(94)+CONSTS(93).AND.VOI.GE.CONSTS(92)+CONSTS(93), CONSTS(96)+ - CONSTS(150)*(((VOI - CONSTS(92)) - CONSTS(94)) - CONSTS(93)), CONSTS(96))
      ALGBRC(69) = ALGBRC(2)+STATES(40)
      ALGBRC(110) =  CONSTS(145)*CONSTS(146)*log( CONSTS(73)*ALGBRC(69))
      ALGBRC(145) =  CONSTS(28)*(EXP((ALGBRC(106)+ALGBRC(110))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(108)/( CONSTS(145)*CONSTS(146))))
      RATES(38) = - ALGBRC(141) - ALGBRC(145)
      RATES(40) = - ALGBRC(145)
      ALGBRC(75) = CONSTS(135)+STATES(42)
      ALGBRC(114) =  CONSTS(145)*CONSTS(146)*log( CONSTS(75)*ALGBRC(75))
      ALGBRC(78) = CONSTS(136)+STATES(43)
      ALGBRC(116) =  CONSTS(145)*CONSTS(146)*log( CONSTS(76)*ALGBRC(78))
      ALGBRC(143) =  CONSTS(26)*(EXP((ALGBRC(112)+ALGBRC(114))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(116)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(124) =  CONSTS(145)*CONSTS(146)*log( CONSTS(50)*ALGBRC(23))
      ALGBRC(90) = CONSTS(139)+STATES(47)
      ALGBRC(126) =  CONSTS(145)*CONSTS(146)*log( CONSTS(79)*ALGBRC(90))
      ALGBRC(131) =  CONSTS(145)*CONSTS(146)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(134) =  CONSTS(145)*CONSTS(146)*log( CONSTS(64)*ALGBRC(47))
      ALGBRC(98) = CONSTS(141)+STATES(52)
      ALGBRC(136) =  CONSTS(145)*CONSTS(146)*log( CONSTS(81)*ALGBRC(98))
      ALGBRC(146) =  CONSTS(29)*(EXP((ALGBRC(116)+ALGBRC(131))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(124)+ALGBRC(126)+ALGBRC(136)+ALGBRC(134))/( CONSTS(145)*CONSTS(146))))
      RATES(43) = ALGBRC(143) - ALGBRC(146)
      ALGBRC(86) = CONSTS(138)+STATES(45)
      ALGBRC(121) =  CONSTS(145)*CONSTS(146)*log( CONSTS(78)*ALGBRC(86))
      ALGBRC(144) =  CONSTS(27)*(EXP((ALGBRC(118)+ALGBRC(114))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(121)/( CONSTS(145)*CONSTS(146))))
      RATES(44) = (ALGBRC(142) - ALGBRC(144))+ALGBRC(145)
      ALGBRC(148) = (ALGBRC(143))+ALGBRC(145)
      RATES(41) = (ALGBRC(141) - ALGBRC(143))+ 1.00000*ALGBRC(148)
      ALGBRC(100) = CONSTS(142)+STATES(53)
      ALGBRC(137) =  CONSTS(145)*CONSTS(146)*log( CONSTS(82)*ALGBRC(100))
      ALGBRC(147) =  CONSTS(30)*(EXP((ALGBRC(121)+ALGBRC(131))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(124)+ALGBRC(126)+ALGBRC(137)+ALGBRC(134))/( CONSTS(145)*CONSTS(146))))
      RATES(45) = ALGBRC(144) - ALGBRC(147)
      RATES(49) = - ALGBRC(146) - ALGBRC(147)
      RATES(50) = ALGBRC(146)+ALGBRC(147)
      ALGBRC(94) = CONSTS(140)+STATES(48)
      ALGBRC(128) =  CONSTS(145)*CONSTS(146)*log( CONSTS(80)*ALGBRC(94))
      ALGBRC(135) =  CONSTS(145)*CONSTS(146)*log( CONSTS(65)*ALGBRC(49))
      ALGBRC(149) =  CONSTS(31)*(EXP(ALGBRC(124)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(128)+ALGBRC(135))/( CONSTS(145)*CONSTS(146))))
      RATES(46) = (ALGBRC(146)+ALGBRC(147)) - ALGBRC(149)
      RATES(51) = ALGBRC(149)
      ALGBRC(102) = CONSTS(143)+STATES(54)
      ALGBRC(138) =  CONSTS(145)*CONSTS(146)*log( CONSTS(83)*ALGBRC(102))
      ALGBRC(140) =  CONSTS(145)*CONSTS(146)*log( CONSTS(70)*ALGBRC(61))
      ALGBRC(150) =  CONSTS(33)*(EXP((ALGBRC(136)+ALGBRC(140))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(138)/( CONSTS(145)*CONSTS(146))))
      RATES(52) = ALGBRC(146) - ALGBRC(150)
      RATES(54) = ALGBRC(150)
      ALGBRC(151) =  CONSTS(32)*(EXP((ALGBRC(128)+ALGBRC(126))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(114)/( CONSTS(145)*CONSTS(146))))
      RATES(42) = (- ALGBRC(143) - ALGBRC(144))+ALGBRC(151)
      RATES(47) = (ALGBRC(146)+ALGBRC(147)) - ALGBRC(151)
      RATES(48) = ALGBRC(149) - ALGBRC(151)
      ALGBRC(104) = CONSTS(144)+STATES(55)
      ALGBRC(139) =  CONSTS(145)*CONSTS(146)*log( CONSTS(84)*ALGBRC(104))
      ALGBRC(152) =  CONSTS(34)*(EXP((ALGBRC(137)+ALGBRC(140))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(139)/( CONSTS(145)*CONSTS(146))))
      RATES(53) = ALGBRC(147) - ALGBRC(152)
      RATES(55) = ALGBRC(152)
      RATES(56) = - ALGBRC(150) - ALGBRC(152)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = CONSTS(97)+STATES(1)
      ALGBRC(29) =  CONSTS(145)*CONSTS(146)*log( CONSTS(35)*ALGBRC(3))
      ALGBRC(5) = CONSTS(98)+STATES(3)
      ALGBRC(33) =  CONSTS(145)*CONSTS(146)*log( CONSTS(37)*ALGBRC(5))
      ALGBRC(6) = CONSTS(100)+STATES(4)
      ALGBRC(35) =  CONSTS(145)*CONSTS(146)*log( CONSTS(38)*ALGBRC(6))
      ALGBRC(70) =  CONSTS(1)*(EXP((ALGBRC(33)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(35)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(4) = CONSTS(99)+STATES(2)
      ALGBRC(31) =  CONSTS(145)*CONSTS(146)*log( CONSTS(36)*ALGBRC(4))
      ALGBRC(26) = CONSTS(114)+STATES(18)
      ALGBRC(67) =  CONSTS(145)*CONSTS(146)*log( CONSTS(52)*ALGBRC(26))
      ALGBRC(73) =  CONSTS(2)*(EXP(ALGBRC(35)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(33)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(7) = CONSTS(105)+STATES(5)
      ALGBRC(37) =  CONSTS(145)*CONSTS(146)*log( CONSTS(39)*ALGBRC(7))
      ALGBRC(8) = CONSTS(106)+STATES(6)
      ALGBRC(39) =  CONSTS(145)*CONSTS(146)*log( CONSTS(40)*ALGBRC(8))
      ALGBRC(76) =  CONSTS(3)*(EXP((ALGBRC(37)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(39)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(79) =  CONSTS(4)*(EXP(ALGBRC(39)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(37)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(9) = CONSTS(102)+STATES(7)
      ALGBRC(41) =  CONSTS(145)*CONSTS(146)*log( CONSTS(41)*ALGBRC(9))
      ALGBRC(10) = CONSTS(103)+STATES(8)
      ALGBRC(43) =  CONSTS(145)*CONSTS(146)*log( CONSTS(42)*ALGBRC(10))
      ALGBRC(82) =  CONSTS(5)*(EXP((ALGBRC(41)+ALGBRC(29))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(43)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(84) =  CONSTS(6)*(EXP(ALGBRC(43)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(41)+ALGBRC(31)+ALGBRC(67))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(11) = CONSTS(107)+STATES(9)
      ALGBRC(46) =  CONSTS(145)*CONSTS(146)*log( CONSTS(43)*ALGBRC(11))
      ALGBRC(14) = CONSTS(109)+STATES(10)
      ALGBRC(48) =  CONSTS(145)*CONSTS(146)*log( CONSTS(44)*ALGBRC(14))
      ALGBRC(87) =  CONSTS(7)*(EXP((ALGBRC(46)+ALGBRC(31))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(48)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(20) = CONSTS(104)+STATES(14)+STATES(27)
      ALGBRC(56) =  CONSTS(145)*CONSTS(146)*log( CONSTS(48)*ALGBRC(20))
      ALGBRC(88) =  CONSTS(10)*(EXP((ALGBRC(33)+ALGBRC(56))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(37)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(21) = CONSTS(101)+STATES(15)
      ALGBRC(58) =  CONSTS(145)*CONSTS(146)*log( CONSTS(49)*ALGBRC(21))
      ALGBRC(92) =  CONSTS(11)*(EXP((ALGBRC(58)+ALGBRC(33))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(41)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(15) = CONSTS(111)+STATES(11)
      ALGBRC(50) =  CONSTS(145)*CONSTS(146)*log( CONSTS(45)*ALGBRC(15))
      ALGBRC(91) =  CONSTS(8)*(EXP(ALGBRC(48)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(46)+ALGBRC(50))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(23) = CONSTS(112)+STATES(16)+STATES(46)
      ALGBRC(62) =  CONSTS(145)*CONSTS(146)*log( CONSTS(50)*ALGBRC(23))
      ALGBRC(24) = CONSTS(113)+STATES(17)
      ALGBRC(64) =  CONSTS(145)*CONSTS(146)*log( CONSTS(51)*ALGBRC(24))
      ALGBRC(96) =  CONSTS(12)*(EXP((ALGBRC(33)+ALGBRC(62))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(64)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(19) = CONSTS(108)+STATES(13)
      ALGBRC(54) =  CONSTS(145)*CONSTS(146)*log( CONSTS(47)*ALGBRC(19))
      ALGBRC(18) = CONSTS(110)+STATES(12)
      ALGBRC(52) =  CONSTS(145)*CONSTS(146)*log( CONSTS(46)*ALGBRC(18))
      ALGBRC(95) =  CONSTS(9)*(EXP((ALGBRC(46)+ALGBRC(52))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(54)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(27) = CONSTS(116)+STATES(20)
      ALGBRC(68) =  CONSTS(145)*CONSTS(146)*log( CONSTS(54)*ALGBRC(27))
      ALGBRC(36) = CONSTS(120)+STATES(25)
      ALGBRC(83) =  CONSTS(145)*CONSTS(146)*log( CONSTS(59)*ALGBRC(36))
      ALGBRC(117) =  CONSTS(14)*(EXP(ALGBRC(68)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(83)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(25) = CONSTS(115)+STATES(19)
      ALGBRC(65) =  CONSTS(145)*CONSTS(146)*log( CONSTS(53)*ALGBRC(25))
      ALGBRC(1) = TERNRY(VOI.LT.CONSTS(87).AND.VOI.GT.CONSTS(87) - CONSTS(89), CONSTS(91)+ CONSTS(149)*((VOI - CONSTS(87))+CONSTS(89)), TERNRY(VOI.GE.CONSTS(87).AND.VOI.LT.CONSTS(87)+CONSTS(88), CONSTS(90)+CONSTS(91), TERNRY(VOI.LT.CONSTS(87)+CONSTS(89)+CONSTS(88).AND.VOI.GE.CONSTS(87)+CONSTS(88), CONSTS(91)+ - CONSTS(149)*(((VOI - CONSTS(87)) - CONSTS(89)) - CONSTS(88)), CONSTS(91))
      ALGBRC(28) = ALGBRC(1)+STATES(21)
      ALGBRC(71) =  CONSTS(145)*CONSTS(146)*log( CONSTS(55)*ALGBRC(28))
      ALGBRC(120) =  CONSTS(17)*(EXP((ALGBRC(65)+ALGBRC(71))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(68)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(30) = CONSTS(117)+STATES(22)
      ALGBRC(74) =  CONSTS(145)*CONSTS(146)*log( CONSTS(56)*ALGBRC(30))
      ALGBRC(115) =  CONSTS(13)*(EXP(ALGBRC(65)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(74)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(32) = CONSTS(118)+STATES(23)
      ALGBRC(77) =  CONSTS(145)*CONSTS(146)*log( CONSTS(57)*ALGBRC(32))
      ALGBRC(34) = CONSTS(119)+STATES(24)
      ALGBRC(80) =  CONSTS(145)*CONSTS(146)*log( CONSTS(58)*ALGBRC(34))
      ALGBRC(119) =  CONSTS(15)*(EXP((ALGBRC(74)+ALGBRC(77))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(80)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(123) = ALGBRC(115)+ALGBRC(120)
      ALGBRC(38) = CONSTS(121)+STATES(26)
      ALGBRC(85) =  CONSTS(145)*CONSTS(146)*log( CONSTS(60)*ALGBRC(38))
      ALGBRC(122) =  CONSTS(16)*(EXP((ALGBRC(83)+ALGBRC(77))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(85)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(89) =  CONSTS(145)*CONSTS(146)*log( CONSTS(48)*ALGBRC(20))
      ALGBRC(40) = CONSTS(122)+STATES(28)
      ALGBRC(93) =  CONSTS(145)*CONSTS(146)*log( CONSTS(61)*ALGBRC(40))
      ALGBRC(45) = CONSTS(124)+STATES(30)+STATES(49)
      ALGBRC(99) =  CONSTS(145)*CONSTS(146)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(47) = CONSTS(125)+STATES(31)+STATES(50)
      ALGBRC(101) =  CONSTS(145)*CONSTS(146)*log( CONSTS(64)*ALGBRC(47))
      ALGBRC(51) = CONSTS(127)+STATES(33)
      ALGBRC(105) =  CONSTS(145)*CONSTS(146)*log( CONSTS(66)*ALGBRC(51))
      ALGBRC(125) =  CONSTS(18)*(EXP((ALGBRC(80)+ALGBRC(99))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(89)+ALGBRC(93)+ALGBRC(105)+ALGBRC(101))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(53) = CONSTS(128)+STATES(34)
      ALGBRC(107) =  CONSTS(145)*CONSTS(146)*log( CONSTS(67)*ALGBRC(53))
      ALGBRC(127) =  CONSTS(19)*(EXP((ALGBRC(85)+ALGBRC(99))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(89)+ALGBRC(93)+ALGBRC(107)+ALGBRC(101))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(42) = CONSTS(123)+STATES(29)
      ALGBRC(97) =  CONSTS(145)*CONSTS(146)*log( CONSTS(62)*ALGBRC(42))
      ALGBRC(49) = CONSTS(126)+STATES(32)+STATES(51)
      ALGBRC(103) =  CONSTS(145)*CONSTS(146)*log( CONSTS(65)*ALGBRC(49))
      ALGBRC(129) =  CONSTS(20)*(EXP(ALGBRC(89)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(97)+ALGBRC(103))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(55) = CONSTS(129)+STATES(35)
      ALGBRC(109) =  CONSTS(145)*CONSTS(146)*log( CONSTS(68)*ALGBRC(55))
      ALGBRC(61) = CONSTS(131)+STATES(37)+STATES(56)
      ALGBRC(113) =  CONSTS(145)*CONSTS(146)*log( CONSTS(70)*ALGBRC(61))
      ALGBRC(130) =  CONSTS(22)*(EXP((ALGBRC(105)+ALGBRC(113))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(109)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(132) =  CONSTS(21)*(EXP((ALGBRC(97)+ALGBRC(93))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(77)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(57) = CONSTS(130)+STATES(36)
      ALGBRC(111) =  CONSTS(145)*CONSTS(146)*log( CONSTS(69)*ALGBRC(57))
      ALGBRC(133) =  CONSTS(23)*(EXP((ALGBRC(107)+ALGBRC(113))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(111)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(66) = CONSTS(133)+STATES(39)
      ALGBRC(108) =  CONSTS(145)*CONSTS(146)*log( CONSTS(72)*ALGBRC(66))
      ALGBRC(81) = CONSTS(137)+STATES(44)
      ALGBRC(118) =  CONSTS(145)*CONSTS(146)*log( CONSTS(77)*ALGBRC(81))
      ALGBRC(142) =  CONSTS(25)*(EXP(ALGBRC(108)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(118)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(63) = CONSTS(132)+STATES(38)
      ALGBRC(106) =  CONSTS(145)*CONSTS(146)*log( CONSTS(71)*ALGBRC(63))
      ALGBRC(72) = CONSTS(134)+STATES(41)
      ALGBRC(112) =  CONSTS(145)*CONSTS(146)*log( CONSTS(74)*ALGBRC(72))
      ALGBRC(141) =  CONSTS(24)*(EXP(ALGBRC(106)/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(112)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(92).AND.VOI.GT.CONSTS(92) - CONSTS(94), CONSTS(96)+ CONSTS(150)*((VOI - CONSTS(92))+CONSTS(94)), TERNRY(VOI.GE.CONSTS(92).AND.VOI.LT.CONSTS(92)+CONSTS(93), CONSTS(95)+CONSTS(96), TERNRY(VOI.LT.CONSTS(92)+CONSTS(94)+CONSTS(93).AND.VOI.GE.CONSTS(92)+CONSTS(93), CONSTS(96)+ - CONSTS(150)*(((VOI - CONSTS(92)) - CONSTS(94)) - CONSTS(93)), CONSTS(96))
      ALGBRC(69) = ALGBRC(2)+STATES(40)
      ALGBRC(110) =  CONSTS(145)*CONSTS(146)*log( CONSTS(73)*ALGBRC(69))
      ALGBRC(145) =  CONSTS(28)*(EXP((ALGBRC(106)+ALGBRC(110))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(108)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(75) = CONSTS(135)+STATES(42)
      ALGBRC(114) =  CONSTS(145)*CONSTS(146)*log( CONSTS(75)*ALGBRC(75))
      ALGBRC(78) = CONSTS(136)+STATES(43)
      ALGBRC(116) =  CONSTS(145)*CONSTS(146)*log( CONSTS(76)*ALGBRC(78))
      ALGBRC(143) =  CONSTS(26)*(EXP((ALGBRC(112)+ALGBRC(114))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(116)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(124) =  CONSTS(145)*CONSTS(146)*log( CONSTS(50)*ALGBRC(23))
      ALGBRC(90) = CONSTS(139)+STATES(47)
      ALGBRC(126) =  CONSTS(145)*CONSTS(146)*log( CONSTS(79)*ALGBRC(90))
      ALGBRC(131) =  CONSTS(145)*CONSTS(146)*log( CONSTS(63)*ALGBRC(45))
      ALGBRC(134) =  CONSTS(145)*CONSTS(146)*log( CONSTS(64)*ALGBRC(47))
      ALGBRC(98) = CONSTS(141)+STATES(52)
      ALGBRC(136) =  CONSTS(145)*CONSTS(146)*log( CONSTS(81)*ALGBRC(98))
      ALGBRC(146) =  CONSTS(29)*(EXP((ALGBRC(116)+ALGBRC(131))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(124)+ALGBRC(126)+ALGBRC(136)+ALGBRC(134))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(86) = CONSTS(138)+STATES(45)
      ALGBRC(121) =  CONSTS(145)*CONSTS(146)*log( CONSTS(78)*ALGBRC(86))
      ALGBRC(144) =  CONSTS(27)*(EXP((ALGBRC(118)+ALGBRC(114))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(121)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(148) = (ALGBRC(143))+ALGBRC(145)
      ALGBRC(100) = CONSTS(142)+STATES(53)
      ALGBRC(137) =  CONSTS(145)*CONSTS(146)*log( CONSTS(82)*ALGBRC(100))
      ALGBRC(147) =  CONSTS(30)*(EXP((ALGBRC(121)+ALGBRC(131))/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(124)+ALGBRC(126)+ALGBRC(137)+ALGBRC(134))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(94) = CONSTS(140)+STATES(48)
      ALGBRC(128) =  CONSTS(145)*CONSTS(146)*log( CONSTS(80)*ALGBRC(94))
      ALGBRC(135) =  CONSTS(145)*CONSTS(146)*log( CONSTS(65)*ALGBRC(49))
      ALGBRC(149) =  CONSTS(31)*(EXP(ALGBRC(124)/( CONSTS(145)*CONSTS(146))) - EXP((ALGBRC(128)+ALGBRC(135))/( CONSTS(145)*CONSTS(146))))
      ALGBRC(102) = CONSTS(143)+STATES(54)
      ALGBRC(138) =  CONSTS(145)*CONSTS(146)*log( CONSTS(83)*ALGBRC(102))
      ALGBRC(140) =  CONSTS(145)*CONSTS(146)*log( CONSTS(70)*ALGBRC(61))
      ALGBRC(150) =  CONSTS(33)*(EXP((ALGBRC(136)+ALGBRC(140))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(138)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(151) =  CONSTS(32)*(EXP((ALGBRC(128)+ALGBRC(126))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(114)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(104) = CONSTS(144)+STATES(55)
      ALGBRC(139) =  CONSTS(145)*CONSTS(146)*log( CONSTS(84)*ALGBRC(104))
      ALGBRC(152) =  CONSTS(34)*(EXP((ALGBRC(137)+ALGBRC(140))/( CONSTS(145)*CONSTS(146))) - EXP(ALGBRC(139)/( CONSTS(145)*CONSTS(146))))
      ALGBRC(12) = ALGBRC(3)+ALGBRC(6)+ALGBRC(10)+ALGBRC(8)
      ALGBRC(13) = ALGBRC(5)+ALGBRC(6)+ALGBRC(9)+ALGBRC(10)+ALGBRC(7)+ALGBRC(8)
      ALGBRC(16) = ALGBRC(4)+ALGBRC(14)+STATES(11)
      ALGBRC(17) = ALGBRC(4)+ALGBRC(14)+ALGBRC(15)+ALGBRC(3)+ALGBRC(6)+ALGBRC(8)+ALGBRC(10)
      ALGBRC(22) = ALGBRC(20)+ALGBRC(7)+ALGBRC(8)
      ALGBRC(44) = ALGBRC(32)+ALGBRC(34)+ALGBRC(38)+ALGBRC(20)+ALGBRC(42)
      ALGBRC(59) = ALGBRC(27)+ALGBRC(28)+ALGBRC(38)+ALGBRC(36)+ALGBRC(53)+ALGBRC(57)
      ALGBRC(60) = ALGBRC(25)+ALGBRC(27)+ALGBRC(30)+ALGBRC(34)+ALGBRC(36)+ALGBRC(38)+ALGBRC(51)+ALGBRC(53)+ALGBRC(55)+ALGBRC(57)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END