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 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