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 325 entries in the algebraic variable array.
C There are a total of 43 entries in each of the rate and state variable arrays.
C There are a total of 210 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is MYOGRS in component temp_myogrs_and_pamk (dimensionless).
C CONSTS(2) is PAMK in component temp_myogrs_and_pamk (dimensionless).
C CONSTS(3) is PAMKRN in component temp_myogrs_and_pamk (dimensionless).
C ALGBRC(7) is ANM in component general_angiotensin_multiplier (dimensionless).
C ALGBRC(85) is CKE in component extracellular_K_concentration (monovalent_mEq_per_litre).
C ALGBRC(4) is AMK in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
C ALGBRC(6) is AMNA in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
C ALGBRC(8) is ANMAL in component angiotensin_control_of_aldosterone_secretion (dimensionless).
C CONSTS(4) is ANMALD in component aldosterone_parameter_values (dimensionless).
C ALGBRC(87) is OSMAL in component osmotic_control_of_aldosterone_secretion (dimensionless).
C ALGBRC(99) is AMR1 in component aldosterone_secretion (dimensionless).
C CONSTS(5) is AMKMUL in component aldosterone_parameter_values (dimensionless).
C CONSTS(6) is ALDINF in component aldosterone_parameter_values (dimensionless).
C CONSTS(7) is ALDKNS in component aldosterone_parameter_values (dimensionless).
C ALGBRC(91) is AMRBSC in component aldosterone_secretion (dimensionless).
C ALGBRC(94) is AMRT in component aldosterone_secretion (dimensionless).
C ALGBRC(97) is AMR in component aldosterone_secretion (dimensionless).
C STATES(1) is AMC in component aldosterone_concentration (dimensionless).
C CONSTS(8) is AMT in component aldosterone_parameter_values (minute).
C ALGBRC(2) is AM in component general_aldosterone_multiplier (dimensionless).
C CONSTS(9) is AM1UL in component aldosterone_parameter_values (dimensionless).
C CONSTS(10) is AM1LL in component aldosterone_parameter_values (dimensionless).
C CONSTS(11) is AMCSNS in component aldosterone_parameter_values (dimensionless).
C CONSTS(12) is ALDMM in component aldosterone_parameter_values (dimensionless).
C ALGBRC(1) is AM1 in component general_aldosterone_multiplier (dimensionless).
C CONSTS(13) is AMKM in component aldosterone_parameter_values (dimensionless).
C ALGBRC(3) is AMKT in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
C CONSTS(14) is AMNAM in component aldosterone_parameter_values (dimensionless).
C CONSTS(15) is AMNAUL in component aldosterone_parameter_values (dimensionless).
C CONSTS(16) is AMNALL in component aldosterone_parameter_values (dimensionless).
C ALGBRC(5) is AMNAT in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
C ALGBRC(233) is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(11) is ANUVN in component angiotensin_effect_on_venous_constriction (dimensionless).
C ALGBRC(10) is ANU in component angiotensin_effect_on_circulation (dimensionless).
C ALGBRC(269) is ANGSCR in component instantaneous_angiotensin_formation (dimensionless).
C ALGBRC(263) is MDFLW3 in component instantaneous_angiotensin_formation (L_per_minute).
C STATES(2) is ANX1 in component time_delayed_angiotensin_formation (dimensionless).
C CONSTS(17) is ANXM in component angiotensin_parameter_values (dimensionless).
C CONSTS(18) is ANV in component angiotensin_parameter_values (minute).
C ALGBRC(276) is ANX in component time_delayed_angiotensin_formation (dimensionless).
C ALGBRC(284) is ANPR in component total_angiotensin_formation (dimensionless).
C CONSTS(19) is REK in component angiotensin_parameter_values (dimensionless).
C ALGBRC(277) is ANPRT in component total_angiotensin_formation (dimensionless).
C ALGBRC(294) is ANPR1 in component artificial_angiotensin_formation (dimensionless).
C CONSTS(20) is ANGKNS in component angiotensin_parameter_values (dimensionless).
C CONSTS(21) is ANGINF in component angiotensin_parameter_values (dimensionless).
C STATES(3) is ANC in component angiotensin_concentration (dimensionless).
C CONSTS(22) is ANT in component angiotensin_parameter_values (minute).
C CONSTS(23) is ANMUL in component angiotensin_parameter_values (dimensionless).
C CONSTS(24) is ANMLL in component angiotensin_parameter_values (dimensionless).
C CONSTS(25) is ANCSNS in component angiotensin_parameter_values (dimensionless).
C CONSTS(26) is ANUM in component angiotensin_parameter_values (dimensionless).
C CONSTS(27) is ANULL in component angiotensin_parameter_values (dimensionless).
C ALGBRC(9) is ANU1 in component angiotensin_effect_on_circulation (dimensionless).
C CONSTS(28) is ANUVM in component angiotensin_parameter_values (dimensionless).
C CONSTS(29) is Z12 in component angiotensin_parameter_values (dimensionless).
C ALGBRC(83) is CNA in component extracellular_Na_concentration (monovalent_mEq_per_litre).
C ALGBRC(66) is PA1 in component pressure_driving_autonomic_receptors (mmHg).
C ALGBRC(15) is ADHMK in component ADH_effect_on_kidney (dimensionless).
C ALGBRC(13) is ADHMV in component ADH_effect_on_nonrenal_vascular_resistance (dimensionless).
C STATES(4) is ADHC in component ADH_in_blood (dimensionless).
C ALGBRC(88) is ADHNA in component osmotic_control_of_ADH_secretion (dimensionless).
C CONSTS(30) is CNR in component antidiuretic_hormone_parameter_values (monovalent_mEq_per_litre).
C ALGBRC(84) is ADHNA1 in component osmotic_control_of_ADH_secretion (dimensionless).
C ALGBRC(70) is ADHPR in component pressure_control_of_ADH_secretion (dimensionless).
C CONSTS(31) is ADHPUL in component antidiuretic_hormone_parameter_values (mmHg).
C CONSTS(32) is ADHPAM in component antidiuretic_hormone_parameter_values (per_mmHg2).
C ALGBRC(68) is ADHPA in component pressure_control_of_ADH_secretion (mmHg).
C ALGBRC(95) is ADH in component total_ADH_secretion (dimensionless).
C CONSTS(33) is ADHINF in component antidiuretic_hormone_parameter_values (dimensionless).
C ALGBRC(92) is ADH1 in component total_ADH_secretion (dimensionless).
C CONSTS(34) is ADHTC in component antidiuretic_hormone_parameter_values (minute).
C CONSTS(35) is ADHVUL in component antidiuretic_hormone_parameter_values (dimensionless).
C CONSTS(36) is ADHVLL in component antidiuretic_hormone_parameter_values (dimensionless).
C ALGBRC(12) is ADHMV1 in component ADH_effect_on_nonrenal_vascular_resistance (dimensionless).
C CONSTS(37) is ADHKLL in component antidiuretic_hormone_parameter_values (dimensionless).
C CONSTS(38) is ADHKUL in component antidiuretic_hormone_parameter_values (dimensionless).
C ALGBRC(14) is ADHMK1 in component ADH_effect_on_kidney (dimensionless).
C ALGBRC(48) is PLA in component left_atrial_pressure (mmHg).
C ALGBRC(35) is PRA in component right_atrial_pressure (mmHg).
C ALGBRC(17) is ANPX in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
C ALGBRC(51) is ANP in component total_ANP_secreted (dimensionless).
C ALGBRC(49) is ANPL in component total_ANP_secreted (dimensionless).
C ALGBRC(36) is ANPR2 in component total_ANP_secreted (dimensionless).
C ALGBRC(53) is ANP1 in component ANP_into_circulation (dimensionless).
C CONSTS(39) is ANPKNS in component atrial_natriuretic_peptide_parameter_values (dimensionless).
C CONSTS(40) is ANPINF in component atrial_natriuretic_peptide_parameter_values (dimensionless).
C STATES(5) is ANPC in component ANP_in_plasma (dimensionless).
C CONSTS(41) is ANPTC in component atrial_natriuretic_peptide_parameter_values (minute).
C CONSTS(42) is ANPXUL in component atrial_natriuretic_peptide_parameter_values (dimensionless).
C ALGBRC(16) is ANPX1 in component ANP_effect_on_renal_afferent_arteriolar_resistance (dimensionless).
C ALGBRC(65) is PA in component arterial_pressure_and_pressure_gradient (mmHg).
C ALGBRC(170) is PO2ART in component arterial_PO2 (mmHg).
C ALGBRC(41) is PPA in component pulmonary_vasculature_pressure (mmHg).
C ALGBRC(21) is VVR in component autonomic_drive_on_target_organs_and_tissues (litre).
C ALGBRC(23) is AUH in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(24) is AUR in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(26) is AUM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(25) is AOM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(27) is AVE in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(20) is AU in component actual_autonomic_stimulation (dimensionless).
C CONSTS(43) is CRRFLX in component autonomics_parameter_values (mmHg).
C CONSTS(44) is EXE in component autonomics_parameter_values (mmHg).
C ALGBRC(72) is AUC in component chemoreceptors_effect_of_PA (dimensionless).
C CONSTS(45) is AUC1 in component autonomics_parameter_values (dimensionless).
C ALGBRC(174) is AUC3 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
C CONSTS(46) is O2CHMO in component autonomics_parameter_values (per_mmHg).
C ALGBRC(172) is AUC2 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
C ALGBRC(18) is AU6C in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(47) is AUX in component autonomics_parameter_values (dimensionless).
C CONSTS(48) is AUK in component autonomics_parameter_values (per_minute).
C CONSTS(49) is BAROTC in component autonomics_parameter_values (minute).
C ALGBRC(67) is AUB in component arterial_baroreceptor_reflex (dimensionless).
C ALGBRC(69) is A1B in component arterial_baroreceptor_reflex (dimensionless).
C ALGBRC(71) is AU6A in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(50) is AU4 in component arterial_baroreceptor_reflex (dimensionless).
C STATES(6) is AU6 in component arterial_baroreceptor_reflex (dimensionless).
C ALGBRC(73) is AUN in component CNS_ischemic_reflex (dimensionless).
C CONSTS(51) is AUN1 in component autonomics_parameter_values (dimensionless).
C ALGBRC(50) is AULP in component autonomic_response_to_vasculature_pressure (dimensionless).
C CONSTS(52) is AULPM in component autonomics_parameter_values (dimensionless).
C CONSTS(207) is AUEX in component autonomic_response_to_exercise (dimensionless).
C CONSTS(53) is EXC in component autonomics_parameter_values (dimensionless).
C CONSTS(54) is EXCXP in component autonomics_parameter_values (dimensionless).
C ALGBRC(178) is AUTTL in component total_autonomic_stimulation (dimensionless).
C ALGBRC(176) is AUTTL1 in component total_autonomic_stimulation (dimensionless).
C CONSTS(55) is EXCML in component autonomics_parameter_values (dimensionless).
C CONSTS(56) is AUDMP in component autonomics_parameter_values (minute).
C CONSTS(57) is AUMAX in component autonomics_parameter_values (dimensionless).
C CONSTS(58) is AUMIN in component autonomics_parameter_values (dimensionless).
C CONSTS(59) is AUSLP in component autonomics_parameter_values (dimensionless).
C ALGBRC(181) is DAU in component actual_autonomic_stimulation (per_minute).
C STATES(7) is AU1 in component actual_autonomic_stimulation (dimensionless).
C ALGBRC(19) is AUT in component actual_autonomic_stimulation (dimensionless).
C CONSTS(60) is VV9 in component autonomics_parameter_values (litre).
C CONSTS(61) is AUL in component autonomics_parameter_values (litre).
C CONSTS(62) is AUV in component autonomics_parameter_values (dimensionless).
C CONSTS(63) is AUS in component autonomics_parameter_values (dimensionless).
C CONSTS(64) is O2A in component autonomics_parameter_values (dimensionless).
C CONSTS(65) is AUM1 in component autonomics_parameter_values (dimensionless).
C CONSTS(66) is AUM2 in component autonomics_parameter_values (dimensionless).
C CONSTS(67) is AUY in component autonomics_parameter_values (dimensionless).
C ALGBRC(22) is AUO in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C CONSTS(68) is MDMP in component autonomics_parameter_values (dimensionless).
C ALGBRC(82) is VEC in component extracellular_fluid_volume (litre).
C ALGBRC(153) is PPD in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
C ALGBRC(234) is RVS in component venous_resistance (mmHg_minute_per_L).
C ALGBRC(150) is DFP in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
C ALGBRC(124) is VPF in component pulmonary_interstitial_free_fluid_volume (litre).
C ALGBRC(235) is BFN in component blood_flow_through_NM_NR_tissues (L_per_minute).
C ALGBRC(210) is PVS in component venous_average_pressure (mmHg).
C ALGBRC(236) is PC in component capillary_pressure (mmHg).
C ALGBRC(28) is CPP in component plasma_protein_concentration (gram_per_L).
C ALGBRC(31) is PPC in component plasma_colloid_osmotic_pressure (mmHg).
C STATES(8) is VP in component plasma_volume (litre).
C ALGBRC(135) is PGH in component hydrostatic_pressure_of_tissue_gel (mmHg).
C ALGBRC(133) is PTC in component total_osmotic_pressure_of_tissue_gel (mmHg).
C ALGBRC(270) is VTCPL in component plasma_leakage (L_per_minute).
C ALGBRC(278) is VTC in component rate_of_fluid_out_of_capillaries (L_per_minute).
C CONSTS(69) is CFC in component capillary_dynamics_parameter_values (L_per_minute_per_mmHg).
C ALGBRC(140) is VTL in component lymph_flow (L_per_minute).
C CONSTS(70) is TRPL in component capillary_dynamics_parameter_values (L_per_minute).
C ALGBRC(285) is VPD in component plasma_volume (L_per_minute).
C STATES(9) is PRP in component total_plasma_protein (gram).
C ALGBRC(30) is DLP in component protein_destruction_and_formation (gram_per_minute).
C CONSTS(71) is CPR in component capillary_dynamics_parameter_values (gram_per_L).
C CONSTS(72) is LPPR in component capillary_dynamics_parameter_values (gram_per_minute).
C CONSTS(73) is LPDE in component capillary_dynamics_parameter_values (dimensionless).
C CONSTS(74) is LPK in component capillary_dynamics_parameter_values (L_per_minute).
C ALGBRC(29) is CPPD in component protein_destruction_and_formation (gram_per_L).
C CONSTS(75) is PCR in component capillary_dynamics_parameter_values (mmHg).
C CONSTS(76) is CPK in component capillary_dynamics_parameter_values (L_per_minute_per_mmHg).
C CONSTS(77) is PCE in component capillary_dynamics_parameter_values (dimensionless).
C ALGBRC(264) is PRCD in component plasma_leakage (mmHg).
C ALGBRC(128) is CPI in component interstitial_protein_concentration (gram_per_L).
C ALGBRC(279) is DPC in component protein_influx_into_interstitium (gram_per_minute).
C ALGBRC(141) is DPL in component lymph_protein_flow (gram_per_minute).
C ALGBRC(286) is DPP in component total_plasma_protein (gram_per_minute).
C ALGBRC(125) is VTS in component total_systemic_fluid_volume (litre).
C ALGBRC(126) is VTS1 in component interstitial_fluid_volume (litre).
C CONSTS(78) is TSSLML in component capillary_dynamics_parameter_values (dimensionless).
C CONSTS(79) is TSSLTC in component capillary_dynamics_parameter_values (per_minute).
C STATES(10) is VTS2 in component interstitial_fluid_volume (litre).
C STATES(11) is TSP in component total_interstitial_protein (gram).
C ALGBRC(287) is DPI in component total_interstitial_protein (gram_per_minute).
C ALGBRC(130) is PTCPR in component interstitial_colloid_osmotic_pressure (mmHg).
C ALGBRC(134) is PTT in component total_tissue_pressure (mmHg).
C ALGBRC(131) is CHY in component hydrostatic_pressure_of_tissue_gel (gram_per_L).
C CONSTS(80) is HYL in component capillary_dynamics_parameter_values (gram).
C CONSTS(81) is CMPTSS in component capillary_dynamics_parameter_values (dimensionless).
C CONSTS(82) is PGHF in component capillary_dynamics_parameter_values (L_mmHg_per_gram).
C ALGBRC(132) is POSHYL in component total_osmotic_pressure_of_tissue_gel (mmHg).
C CONSTS(83) is GCOPF in component capillary_dynamics_parameter_values (per_mmHg).
C CONSTS(84) is VTSF in component capillary_dynamics_parameter_values (litre).
C ALGBRC(136) is PIF in component interstial_free_fluid_pressure (mmHg).
C ALGBRC(138) is PTS in component interstitial_solid_tissue_pressure (mmHg).
C CONSTS(85) is PLDF in component capillary_dynamics_parameter_values (mmHg).
C ALGBRC(137) is PLD1 in component lymph_flow (mmHg).
C ALGBRC(139) is PLD in component lymph_flow (mmHg).
C ALGBRC(127) is VG in component interstitial_gel_volume (litre).
C ALGBRC(129) is VIF in component interstitial_free_fluid_volume (litre).
C ALGBRC(101) is AMM in component global_M_blood_flow_autoregulation_output (dimensionless).
C ALGBRC(113) is ARM in component global_NM_blood_flow_autoregulation_output (dimensionless).
C ALGBRC(199) is ATRRFB in component volume_effect_on_arteries (dimensionless).
C ALGBRC(206) is ATRVFB in component volume_effect_on_unstressed_venous_volume (litre).
C ALGBRC(89) is HMD in component heart_deterioration (dimensionless).
C STATES(12) is HPL in component left_ventricular_hypertrophy (dimensionless).
C STATES(13) is HPR in component right_ventricular_hypertrophy (dimensionless).
C ALGBRC(159) is OSA in component arterial_PO2 (dimensionless).
C ALGBRC(265) is RBF in component actual_renal_blood_flow (L_per_minute).
C ALGBRC(183) is VIM in component blood_viscosity (dimensionless).
C STATES(14) is VRC in component RBC_volume (litre).
C STATES(15) is VV6 in component long_term_stress_relaxation (litre).
C STATES(16) is VV7 in component short_term_stress_relaxation (litre).
C ALGBRC(216) is BFM in component blood_flow_through_M_tissues (L_per_minute).
C ALGBRC(280) is QAO in component systemic_blood_flow (L_per_minute).
C ALGBRC(169) is QRO in component right_ventricular_output (L_per_minute).
C ALGBRC(165) is QLO in component left_ventricular_output (L_per_minute).
C ALGBRC(45) is RPA in component pulmonary_arterial_resistance (mmHg_minute_per_L).
C ALGBRC(54) is RPV in component pulmonary_venous_resistance (mmHg_minute_per_L).
C ALGBRC(208) is VVE in component venous_excess_volume (litre).
C STATES(17) is VVS1 in component venous_blood_volume (litre).
C STATES(18) is VAS1 in component arterial_blood_volume (litre).
C STATES(19) is VLA1 in component left_atrial_blood_volume (litre).
C STATES(20) is VPA1 in component pulmonary_vasculature_blood_volume (litre).
C STATES(21) is VRA1 in component right_atrial_blood_volume (litre).
C ALGBRC(32) is VBD in component total_blood_volume_change (litre).
C ALGBRC(213) is QVO in component rate_of_blood_flow_from_veins_to_right_atrium (L_per_minute).
C ALGBRC(33) is VRA in component right_atrial_blood_volume (litre).
C ALGBRC(215) is DRA in component right_atrial_blood_volume (L_per_minute).
C ALGBRC(34) is VRE in component right_atrial_pressure (litre).
C ALGBRC(37) is PRA1 in component autonomic_stimulation_effect_on_right_atrial_pressure (mmHg).
C CONSTS(86) is HTAUML in component circulatory_dynamics_parameter_values (dimensionless).
C ALGBRC(161) is RVM in component pressure_effect_on_right_ventricular_pumping (dimensionless).
C ALGBRC(160) is PP2 in component pressure_effect_on_right_ventricular_pumping (mmHg).
C ALGBRC(59) is QLN in component left_ventricular_output (L_per_minute).
C ALGBRC(166) is HPEF in component pumping_effectiveness_of_right_ventricle (L_per_minute).
C CONSTS(87) is QRF in component circulatory_dynamics_parameter_values (L_per_minute).
C CONSTS(88) is HSR in component circulatory_dynamics_parameter_values (dimensionless).
C ALGBRC(38) is QRN in component right_ventricular_output (dimensionless).
C ALGBRC(57) is QPO in component rate_of_blood_flow_from_pulmonary_veins_to_left_atrium (L_per_minute).
C ALGBRC(39) is VPA in component pulmonary_vasculature_blood_volume (litre).
C ALGBRC(171) is DPA in component pulmonary_vasculature_blood_volume (L_per_minute).
C ALGBRC(40) is VPE in component pulmonary_vasculature_pressure (litre).
C ALGBRC(42) is PP1T in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(43) is PP1 in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(44) is CPA in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(52) is PL1 in component pulmonary_venous_resistance (mmHg).
C ALGBRC(55) is RPT in component total_pulmonary_vascular_resistance (mmHg_minute_per_L).
C ALGBRC(56) is PGL in component pressure_gradient_through_the_lungs (mmHg).
C ALGBRC(46) is VLA in component left_atrial_blood_volume (litre).
C ALGBRC(167) is DLA in component left_atrial_blood_volume (L_per_minute).
C ALGBRC(47) is VLE in component left_atrial_pressure (litre).
C ALGBRC(58) is PLA1 in component autonomic_stimulation_effect_on_left_atrial_pressure (mmHg).
C ALGBRC(163) is LVM in component pumping_effectiveness_of_left_ventricle (dimensionless).
C ALGBRC(162) is PA2 in component pumping_effectiveness_of_left_ventricle (mmHg).
C ALGBRC(164) is QLOT in component left_ventricular_output (L_per_minute).
C CONSTS(89) is HSL in component circulatory_dynamics_parameter_values (dimensionless).
C ALGBRC(74) is QLO1 in component left_ventricular_output (L_per_minute).
C ALGBRC(60) is VVS in component venous_blood_volume (litre).
C ALGBRC(288) is DVS in component venous_blood_volume (L_per_minute).
C ALGBRC(61) is VVA in component angiotensin_induced_venous_constriction (litre).
C CONSTS(90) is ANY in component circulatory_dynamics_parameter_values (litre).
C ALGBRC(207) is VVE1 in component venous_excess_volume (litre).
C CONSTS(91) is CV in component circulatory_dynamics_parameter_values (L_per_mmHg).
C ALGBRC(209) is PVS1 in component venous_average_pressure (mmHg).
C ALGBRC(62) is PR1 in component venous_outflow_pressure_into_heart (mmHg).
C CONSTS(92) is PR1LL in component circulatory_dynamics_parameter_values (mmHg).
C ALGBRC(211) is RVG in component resistance_from_veins_to_right_atrium (mmHg_minute_per_L).
C ALGBRC(212) is PGV in component rate_of_blood_flow_from_veins_to_right_atrium (mmHg).
C CONSTS(93) is CN7 in component circulatory_dynamics_parameter_values (dimensionless).
C CONSTS(94) is CN2 in component circulatory_dynamics_parameter_values (per_mmHg).
C CONSTS(95) is RVSM in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
C ALGBRC(237) is CN3 in component venous_resistance (dimensionless).
C ALGBRC(238) is RV1 in component venous_resistance (mmHg_minute_per_L).
C ALGBRC(239) is NNRVR in component NM_NR_venous_resistance (mmHg_minute_per_L).
C ALGBRC(63) is VAS in component arterial_blood_volume (litre).
C ALGBRC(289) is DAS in component arterial_blood_volume (L_per_minute).
C ALGBRC(75) is PAG in component arterial_pressure_and_pressure_gradient (mmHg).
C ALGBRC(64) is VAE in component arterial_pressure_and_pressure_gradient (litre).
C ALGBRC(76) is PAM in component pressure_effect_on_arterial_distention (dimensionless).
C CONSTS(96) is PAEX in component circulatory_dynamics_parameter_values (dimensionless).
C ALGBRC(201) is R1 in component non_renal_systemic_arterial_resistance_multiplier (dimensionless).
C ALGBRC(203) is NNRAR in component NM_NR_arterial_resistance (mmHg_minute_per_L).
C CONSTS(97) is RAR in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
C CONSTS(98) is RMULT1 in component circulatory_dynamics_parameter_values (dimensionless).
C ALGBRC(214) is PGS in component pressure_gradient_from_arteries_to_veins (mmHg).
C ALGBRC(205) is RSM in component M_systemic_resistance (mmHg_minute_per_L).
C CONSTS(99) is RAM in component circulatory_dynamics_parameter_values (mmHg_minute_per_L).
C ALGBRC(240) is RSN in component total_NM_NR_systemic_resistance (mmHg_minute_per_L).
C ALGBRC(77) is FISFLO in component blood_flow_through_AV_fistulas (L_per_minute).
C CONSTS(100) is FIS in component circulatory_dynamics_parameter_values (L_per_minute_per_mmHg).
C ALGBRC(271) is SYSFLO in component systemic_blood_flow (L_per_minute).
C ALGBRC(290) is RTP in component total_peripheral_resistance (mmHg_minute_per_L).
C STATES(22) is TVD in component rate_of_fluid_intake (L_per_minute).
C ALGBRC(312) is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(188) is STH in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
C ALGBRC(316) is KOD in component actual_K_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(325) is VUD in component actual_urine_volume (L_per_minute).
C STATES(23) is VTW in component total_body_water (litre).
C CONSTS(101) is NID in component electrolytes_parameter_values (monovalent_mEq_per_minute).
C CONSTS(102) is TRPL in component electrolytes_parameter_values (L_per_minute).
C ALGBRC(314) is NED in component extracellular_Na_concentration (monovalent_mEq_per_minute).
C STATES(24) is NAE in component extracellular_Na_concentration (monovalent_mEq).
C ALGBRC(78) is AMK1 in component aldosterone_effect_on_cellular_K_distribution (dimensionless).
C CONSTS(103) is ALCLK in component electrolytes_parameter_values (dimensionless).
C ALGBRC(79) is KE in component extracellular_K_concentration (monovalent_mEq).
C STATES(25) is KTOT in component extracellular_K_concentration (monovalent_mEq).
C CONSTS(104) is KID in component electrolytes_parameter_values (monovalent_mEq_per_minute).
C ALGBRC(318) is KTOTD in component extracellular_K_concentration (monovalent_mEq_per_minute).
C STATES(26) is VIC in component intracellular_fluid_volume (litre).
C ALGBRC(81) is CKI in component intracellular_K_concentration (monovalent_mEq_per_litre).
C ALGBRC(80) is KI in component intracellular_K_concentration (monovalent_mEq).
C ALGBRC(90) is VID in component intracellular_fluid_volume (L_per_minute).
C CONSTS(105) is VIDML in component electrolytes_parameter_values (litre2_per_monovalent_mEq_per_minute).
C ALGBRC(86) is CCD in component intracellular_fluid_volume (monovalent_mEq_per_litre).
C ALGBRC(115) is POT in component pressure_of_O2_in_NM_tissue_cells (mmHg).
C CONSTS(106) is HSL in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
C CONSTS(107) is Z13 in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
C CONSTS(108) is HSR in component heart_hypertrophy_or_deterioration_parameter_values (dimensionless).
C CONSTS(109) is DHDTR in component heart_hypertrophy_or_deterioration_parameter_values (per_mmHg_per_minute).
C STATES(27) is HMD1 in component heart_deterioration (dimensionless).
C ALGBRC(117) is DHM in component heart_deterioration (per_minute).
C ALGBRC(93) is AUHR in component effect_of_autonomic_stimulation_on_HR (beats_per_minute).
C CONSTS(208) is PRHR in component effect_of_PRA_on_HR (beats_per_minute).
C CONSTS(110) is PR1LL in component HR_and_SV_parameter_values (mmHg).
C ALGBRC(96) is HDHR in component effect_of_heart_deterioration_on_HR (dimensionless).
C ALGBRC(98) is HR in component heart_rate (beats_per_minute).
C ALGBRC(168) is SVO in component stroke_volume_output (litre).
C ALGBRC(103) is PMO in component pressure_of_O2_in_M_tissue_cells (mmHg).
C ALGBRC(105) is PDO in component M_autoregulatory_driving_force (mmHg).
C ALGBRC(107) is POE in component M_ST_sensitivity_control (mmHg).
C CONSTS(111) is POM in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C ALGBRC(100) is AMM1 in component M_ST_time_delay_and_limit (dimensionless).
C CONSTS(112) is A4K in component M_autoregulatory_local_blood_flow_parameter_values (minute).
C CONSTS(113) is AMM4 in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C STATES(28) is AMM1T in component M_ST_time_delay_and_limit (dimensionless).
C ALGBRC(108) is POF in component M_LT_sensitivity_control (mmHg).
C CONSTS(114) is POM2 in component M_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C STATES(29) is AMM2 in component M_LT_time_delay (dimensionless).
C CONSTS(115) is A4K2 in component M_autoregulatory_local_blood_flow_parameter_values (minute).
C STATES(30) is OVA in component O2_volume_of_arterial_blood (mL_per_L).
C ALGBRC(158) is HM in component hematocrit_fraction (dimensionless).
C ALGBRC(230) is RMO in component delivery_of_O2_to_M_tissues (mL_per_minute).
C ALGBRC(217) is O2ARTM in component M_O2_blood_supply (mL_per_minute).
C ALGBRC(231) is PVO in component M_venous_O2_content (mmHg).
C CONSTS(116) is EXC in component M_O2_delivery_parameter_values (dimensionless).
C CONSTS(117) is EXCXP2 in component M_O2_delivery_parameter_values (dimensionless).
C ALGBRC(232) is OVS in component M_venous_O2_content (dimensionless).
C ALGBRC(106) is MMO in component metabolic_O2_consumption_by_M_tissue (mL_per_minute).
C ALGBRC(104) is P2O in component metabolic_O2_consumption_by_M_tissue (mmHg).
C CONSTS(118) is OMM in component M_O2_delivery_parameter_values (mL_per_minute).
C CONSTS(119) is PM5 in component M_O2_delivery_parameter_values (per_mmHg).
C ALGBRC(102) is QOM in component volume_of_O2_in_M_tissue (mL).
C ALGBRC(262) is DO2M in component volume_of_O2_in_M_tissue (mL_per_minute).
C STATES(31) is QOM1 in component volume_of_O2_in_M_tissue (mL).
C CONSTS(120) is PK2 in component M_O2_delivery_parameter_values (mmHg_per_mL).
C ALGBRC(118) is POD in component NM_autoregulatory_driving_force (mmHg).
C CONSTS(121) is POR in component NM_autoregulatory_local_blood_flow_parameter_values (mmHg).
C ALGBRC(120) is POB in component NM_ST_sensitivity_control (mmHg).
C CONSTS(122) is POK in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C ALGBRC(109) is AR1 in component NM_ST_time_delay_and_damping (dimensionless).
C CONSTS(123) is A1K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
C STATES(32) is AR1T in component NM_ST_time_delay_and_damping (dimensionless).
C ALGBRC(121) is POA in component NM_I_sensitivity_control (mmHg).
C CONSTS(124) is PON in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C ALGBRC(110) is AR2 in component NM_I_time_delay_and_limit (dimensionless).
C CONSTS(125) is A2K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
C STATES(33) is AR2T in component NM_I_time_delay_and_limit (dimensionless).
C ALGBRC(122) is POC in component NM_LT_sensitivity_control (mmHg).
C CONSTS(126) is POZ in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C ALGBRC(111) is AR3 in component NM_LT_time_delay_and_limit (dimensionless).
C CONSTS(127) is A3K in component NM_autoregulatory_local_blood_flow_parameter_values (minute).
C STATES(34) is AR3T in component NM_LT_time_delay_and_limit (dimensionless).
C ALGBRC(112) is ARM1 in component total_NM_autoregulation (dimensionless).
C CONSTS(128) is AUTOSN in component NM_autoregulatory_local_blood_flow_parameter_values (dimensionless).
C ALGBRC(272) is DOB in component delivery_of_O2_to_NM_tissues (mL_per_minute).
C ALGBRC(266) is O2ARTN in component NM_O2_blood_supply (mL_per_minute).
C ALGBRC(273) is POV in component NM_venous_O2_content (mmHg).
C ALGBRC(274) is OSV in component NM_venous_O2_content (dimensionless).
C ALGBRC(119) is MO2 in component O2_consumption_by_NM_tissue (mL_per_minute).
C CONSTS(129) is O2M in component NM_O2_delivery_parameter_values (mL_per_minute).
C ALGBRC(116) is P1O in component O2_consumption_by_NM_tissue (mmHg).
C ALGBRC(114) is QO2 in component volume_of_O2_in_NM_tissue (mL).
C ALGBRC(291) is DO2N in component volume_of_O2_in_NM_tissue (mL_per_minute).
C ALGBRC(281) is DO2N1 in component volume_of_O2_in_NM_tissue (mL_per_minute).
C STATES(35) is QO2T in component volume_of_O2_in_NM_tissue (mL).
C ALGBRC(123) is PCP in component pulmonary_capillary_pressure (mmHg).
C ALGBRC(145) is POS in component colloid_osmotic_pressure_of_pulmonary_interstitium (mmHg).
C ALGBRC(142) is PPI in component pulmonary_interstitial_fluid_pressure (mmHg).
C ALGBRC(146) is PFI in component fluid_filtration_into_pulmonary_interstitium (L_per_minute).
C CONSTS(130) is CPF in component pulmonary_fluid_dynamics_parameter_values (L_per_minute_per_mmHg).
C ALGBRC(148) is PLF in component lung_lymphatic_protein_flow (L_per_minute).
C ALGBRC(149) is DFZ in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
C STATES(36) is VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
C ALGBRC(151) is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
C ALGBRC(147) is PPN in component protein_leakage_into_pulmonary_interstitium (gram_per_minute).
C ALGBRC(144) is CPN in component concentration_of_protein_in_pulmonary_interstitium (gram_per_L).
C ALGBRC(152) is PPZ in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
C STATES(37) is PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
C ALGBRC(143) is PPR in component concentration_of_protein_in_pulmonary_interstitium (gram).
C ALGBRC(282) is O2UTIL in component total_O2_utilization (mL_per_minute).
C ALGBRC(155) is O2VAD2 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C ALGBRC(177) is O2VTS2 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C ALGBRC(292) is ALVENT in component alveolar_ventilation (L_per_minute).
C CONSTS(131) is VNTSTM in component pulmonary_O2_uptake_parameter_values (dimensionless).
C ALGBRC(295) is PO2ALV in component alveolar_PO2 (mmHg).
C CONSTS(132) is PO2AMB in component pulmonary_O2_uptake_parameter_values (mmHg).
C ALGBRC(297) is O2DFS in component respiratory_O2_diffusion_into_capillaries (mL_per_minute).
C CONSTS(133) is PL2 in component pulmonary_O2_uptake_parameter_values (L_mL_per_minute_per_mmHg).
C CONSTS(134) is VPTISS in component pulmonary_O2_uptake_parameter_values (litre).
C ALGBRC(154) is RSPDFC in component respiratory_O2_diffusion_into_capillaries (mL_per_minute_per_mmHg).
C ALGBRC(299) is DOVA in component O2_volume_of_arterial_blood (mL_per_L_per_minute).
C ALGBRC(175) is O2VTST in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C ALGBRC(173) is O2VTST1 in component acute_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C ALGBRC(180) is DO2VAD in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (per_minute).
C STATES(38) is O2VAD1 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C ALGBRC(157) is HM1 in component hematocrit_fraction (dimensionless).
C ALGBRC(156) is VB in component hematocrit_fraction (litre).
C ALGBRC(179) is VIE in component viscosity_due_to_RBCs (dimensionless).
C CONSTS(135) is HMK in component red_cells_and_viscosity_parameter_values (dimensionless).
C CONSTS(136) is HKM in component red_cells_and_viscosity_parameter_values (dimensionless).
C ALGBRC(182) is VIB in component blood_viscosity (dimensionless).
C ALGBRC(189) is HM7 in component oxygen_stimulation (mmHg).
C CONSTS(137) is PO2AMB in component red_cells_and_viscosity_parameter_values (mmHg).
C CONSTS(138) is HM6 in component red_cells_and_viscosity_parameter_values (mmHg).
C CONSTS(209) is PO2AM1 in component oxygen_stimulation (mmHg).
C ALGBRC(185) is HM3 in component oxygen_stimulation (mmHg).
C CONSTS(210) is HM4 in component oxygen_stimulation (mmHg).
C ALGBRC(187) is HM5 in component oxygen_stimulation (mmHg).
C ALGBRC(191) is RC1 in component RBC_production (L_per_minute).
C CONSTS(139) is HM8 in component red_cells_and_viscosity_parameter_values (L_per_minute_per_mmHg).
C CONSTS(140) is REK in component red_cells_and_viscosity_parameter_values (dimensionless).
C ALGBRC(194) is RC2 in component RBC_destruction (L_per_minute).
C CONSTS(141) is RKC in component red_cells_and_viscosity_parameter_values (per_minute).
C CONSTS(142) is TRRBC in component red_cells_and_viscosity_parameter_values (L_per_minute).
C ALGBRC(197) is RCD in component RBC_volume (L_per_minute).
C CONSTS(143) is SR in component stress_relaxation_parameter_values (dimensionless).
C CONSTS(144) is SRK in component stress_relaxation_parameter_values (minute).
C CONSTS(145) is SR2 in component stress_relaxation_parameter_values (dimensionless).
C CONSTS(146) is SRK2 in component stress_relaxation_parameter_values (minute).
C CONSTS(147) is ANMSLT in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
C CONSTS(148) is Z10 in component thirst_drinking_and_salt_appetite_parameter_values (mmHg).
C CONSTS(149) is Z11 in component thirst_drinking_and_salt_appetite_parameter_values (per_mmHg2).
C ALGBRC(184) is ANMSML in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
C ALGBRC(186) is STH1 in component effect_of_salt_appetite_stimulation_on_thirst (dimensionless).
C ALGBRC(192) is AHCM in component effect_of_antidiuretic_hormone_on_thirst (dimensionless).
C CONSTS(150) is AHTHM in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
C ALGBRC(195) is ANMTH in component effect_of_angiotensin_on_thirst (dimensionless).
C CONSTS(151) is ANMTM in component thirst_drinking_and_salt_appetite_parameter_values (dimensionless).
C CONSTS(152) is DR in component thirst_drinking_and_salt_appetite_parameter_values (L_per_minute).
C CONSTS(153) is TVDDL in component thirst_drinking_and_salt_appetite_parameter_values (minute).
C ALGBRC(200) is AHTH in component rate_of_fluid_intake (dimensionless).
C ALGBRC(198) is AHTH1 in component rate_of_fluid_intake (dimensionless).
C ALGBRC(204) is TVZ in component rate_of_fluid_intake (L_per_minute).
C ALGBRC(202) is TVZ1 in component rate_of_fluid_intake (L_per_minute).
C ALGBRC(193) is AHZ in component effect_of_pressure_on_volume_receptors (dimensionless).
C CONSTS(154) is AH10 in component volume_receptors_parameter_values (dimensionless).
C CONSTS(155) is AH9 in component volume_receptors_parameter_values (per_mmHg).
C ALGBRC(190) is AHZ1 in component effect_of_pressure_on_volume_receptors (dimensionless).
C STATES(39) is AHY in component time_dependent_volume_receptor_adaptation (dimensionless).
C CONSTS(156) is AH11 in component volume_receptors_parameter_values (minute).
C ALGBRC(196) is AH7 in component total_volume_nervous_feedback (dimensionless).
C CONSTS(157) is ATRFBM in component volume_receptors_parameter_values (dimensionless).
C CONSTS(158) is ATRVM in component volume_receptors_parameter_values (litre).
C ALGBRC(218) is PAR in component perfusion_pressure (mmHg).
C CONSTS(159) is GBL in component kidney_parameter_values (mmHg).
C CONSTS(160) is RAPRSP in component kidney_parameter_values (mmHg).
C CONSTS(161) is RFCDFT in component kidney_parameter_values (dimensionless).
C CONSTS(162) is RCDFPC in component kidney_parameter_values (dimensionless).
C CONSTS(163) is RCDFDP in component kidney_parameter_values (minute).
C STATES(40) is PAR1 in component perfusion_pressure (mmHg).
C ALGBRC(241) is RNAUG2 in component renal_autoregulatory_feedback_factor (dimensionless).
C CONSTS(164) is RNAUGN in component kidney_parameter_values (minute_per_L).
C CONSTS(165) is RNAULL in component kidney_parameter_values (dimensionless).
C CONSTS(166) is RNAUUL in component kidney_parameter_values (dimensionless).
C CONSTS(167) is RNAUAD in component kidney_parameter_values (per_minute).
C ALGBRC(242) is RNAUG1 in component renal_autoregulatory_feedback_factor (dimensionless).
C ALGBRC(243) is RNAUG1T in component renal_autoregulatory_feedback_factor (dimensionless).
C STATES(41) is RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C ALGBRC(220) is AUMK in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(168) is ARF in component kidney_parameter_values (dimensionless).
C ALGBRC(219) is AUMKT in component autonomic_effect_on_AAR (dimensionless).
C ALGBRC(222) is ANMAR in component angiotensin_effect_on_AAR (dimensionless).
C CONSTS(169) is ANMAM in component kidney_parameter_values (dimensionless).
C CONSTS(170) is ANMARL in component kidney_parameter_values (dimensionless).
C ALGBRC(221) is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
C ALGBRC(244) is AAR1 in component AAR_calculation (mmHg_minute_per_L).
C CONSTS(171) is AARK in component kidney_parameter_values (mmHg_minute_per_L).
C ALGBRC(245) is AAR in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(172) is ANPXAF in component kidney_parameter_values (mmHg_minute_per_L).
C CONSTS(173) is AARLL in component kidney_parameter_values (mmHg_minute_per_L).
C ALGBRC(246) is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C ALGBRC(223) is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
C CONSTS(174) is AUMK1 in component kidney_parameter_values (dimensionless).
C ALGBRC(224) is ANMER in component angiotensin_effect_on_EAR (dimensionless).
C CONSTS(175) is ANMEM in component kidney_parameter_values (dimensionless).
C ALGBRC(247) is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
C CONSTS(176) is EFAFR in component kidney_parameter_values (dimensionless).
C ALGBRC(248) is EAR in component EAR_calculation (mmHg_minute_per_L).
C CONSTS(177) is EARK in component kidney_parameter_values (mmHg_minute_per_L).
C CONSTS(178) is EARLL in component kidney_parameter_values (mmHg_minute_per_L).
C ALGBRC(249) is EAR1 in component EAR_calculation (mmHg_minute_per_L).
C ALGBRC(250) is RR in component total_renal_resistance (mmHg_minute_per_L).
C ALGBRC(251) is RFN in component normal_renal_blood_flow (L_per_minute).
C CONSTS(179) is REK in component kidney_parameter_values (dimensionless).
C ALGBRC(252) is GFN in component glomerular_filtration_rate (L_per_minute).
C ALGBRC(253) is GLPC in component glomerular_colloid_osmotic_pressure (mmHg).
C CONSTS(180) is GPPD in component kidney_parameter_values (dimensionless).
C CONSTS(181) is GLPCA in component kidney_parameter_values (mmHg).
C ALGBRC(254) is EFAFPR in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(255) is EFAFPR1 in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(256) is GLP in component glomerular_pressure (mmHg).
C ALGBRC(257) is APD in component glomerular_pressure (mmHg).
C ALGBRC(267) is GFR in component glomerular_filtration_rate (L_per_minute).
C CONSTS(182) is PXTP in component kidney_parameter_values (mmHg).
C CONSTS(183) is GFLC in component kidney_parameter_values (L_per_minute_per_mmHg).
C CONSTS(184) is GFNLL in component kidney_parameter_values (L_per_minute).
C ALGBRC(258) is PFL in component glomerular_filtration_rate (mmHg).
C ALGBRC(259) is GFN1 in component glomerular_filtration_rate (L_per_minute).
C CONSTS(185) is MDFL1 in component kidney_parameter_values (dimensionless).
C ALGBRC(260) is PTFL in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(261) is MDFLWT in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(275) is RTSPPC in component renal_tissue_osmotic_pressure (mmHg).
C CONSTS(186) is RTPPR in component kidney_parameter_values (dimensionless).
C CONSTS(187) is RTPPRS in component kidney_parameter_values (mmHg).
C ALGBRC(268) is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
C ALGBRC(319) is UROD in component actual_urea_excretion_rate (mOsm_per_minute).
C STATES(42) is PLUR in component glomerular_urea_concentration (mOsm).
C CONSTS(188) is URFORM in component kidney_parameter_values (mOsm_per_minute).
C ALGBRC(225) is PLURC in component plasma_urea_concentration (mOsm_per_litre).
C ALGBRC(283) is RCPRS in component peritubular_capillary_pressure (mmHg).
C CONSTS(189) is RFABX in component kidney_parameter_values (dimensionless).
C CONSTS(190) is RVRS in component kidney_parameter_values (mmHg_minute_per_L).
C ALGBRC(301) is RFABD in component peritubular_capillary_reabsorption_factor (dimensionless).
C CONSTS(191) is RTSPRS in component kidney_parameter_values (mmHg).
C CONSTS(192) is RABSC in component kidney_parameter_values (per_mmHg).
C CONSTS(193) is RFABDP in component kidney_parameter_values (dimensionless).
C CONSTS(194) is RFABDM in component kidney_parameter_values (dimensionless).
C ALGBRC(293) is RABSPR in component peritubular_capillary_reabsorption_factor (mmHg).
C ALGBRC(296) is RFAB1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(298) is RFAB in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(300) is RFABD1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(302) is DTNAI in component distal_tubular_Na_delivery (monovalent_mEq_per_minute).
C ALGBRC(304) is DTNARA in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(195) is DTNAR in component kidney_parameter_values (monovalent_mEq_per_minute).
C CONSTS(196) is DIURET in component kidney_parameter_values (dimensionless).
C CONSTS(197) is AHMNAR in component kidney_parameter_values (dimensionless).
C CONSTS(198) is DTNARL in component kidney_parameter_values (monovalent_mEq_per_minute).
C ALGBRC(303) is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(227) is DTNANG in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(199) is ANMNAM in component kidney_parameter_values (dimensionless).
C ALGBRC(226) is DTNANG1 in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(305) is DTKI in component distal_tubular_K_delivery (monovalent_mEq_per_minute).
C ALGBRC(306) is RFABK in component effect_of_physical_forces_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C CONSTS(200) is RFABKM in component kidney_parameter_values (monovalent_mEq_per_minute).
C ALGBRC(308) is MDFLK in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C CONSTS(201) is MDFLKM in component kidney_parameter_values (monovalent_mEq_per_litre).
C ALGBRC(307) is MDFLK1 in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C ALGBRC(315) is KODN in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(324) is VUDN in component normal_urine_volume (L_per_minute).
C STATES(43) is DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(309) is DTKSC in component K_secretion_from_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(202) is ANMKEM in component kidney_parameter_values (dimensionless).
C CONSTS(203) is ANMKEL in component kidney_parameter_values (dimensionless).
C CONSTS(204) is CKEEX in component kidney_parameter_values (dimensionless).
C ALGBRC(228) is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
C ALGBRC(229) is ANMKE in component K_secretion_from_distal_tubules (dimensionless).
C ALGBRC(311) is NODN in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(310) is NODN1 in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(313) is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(317) is DTURI in component normal_urea_excretion (mOsm_per_minute).
C ALGBRC(320) is OSMOPN1 in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(321) is OSMOPN in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(322) is OSMOP1T in component normal_urine_volume (mOsm_per_minute).
C ALGBRC(323) is OSMOP1 in component normal_urine_volume (mOsm_per_minute).
C CONSTS(205) is RNAGTC in component kidney_parameter_values (minute).
C CONSTS(206) is GFNDMP in component kidney_parameter_values (dimensionless).
C RATES(1) is d/dt AMC in component aldosterone_concentration (dimensionless).
C RATES(2) is d/dt ANX1 in component time_delayed_angiotensin_formation (dimensionless).
C RATES(3) is d/dt ANC in component angiotensin_concentration (dimensionless).
C RATES(4) is d/dt ADHC in component ADH_in_blood (dimensionless).
C RATES(5) is d/dt ANPC in component ANP_in_plasma (dimensionless).
C RATES(6) is d/dt AU6 in component arterial_baroreceptor_reflex (dimensionless).
C RATES(7) is d/dt AU1 in component actual_autonomic_stimulation (dimensionless).
C RATES(8) is d/dt VP in component plasma_volume (litre).
C RATES(9) is d/dt PRP in component total_plasma_protein (gram).
C RATES(10) is d/dt VTS2 in component interstitial_fluid_volume (litre).
C RATES(11) is d/dt TSP in component total_interstitial_protein (gram).
C RATES(21) is d/dt VRA1 in component right_atrial_blood_volume (litre).
C RATES(20) is d/dt VPA1 in component pulmonary_vasculature_blood_volume (litre).
C RATES(19) is d/dt VLA1 in component left_atrial_blood_volume (litre).
C RATES(17) is d/dt VVS1 in component venous_blood_volume (litre).
C RATES(18) is d/dt VAS1 in component arterial_blood_volume (litre).
C RATES(24) is d/dt NAE in component extracellular_Na_concentration (monovalent_mEq).
C RATES(25) is d/dt KTOT in component extracellular_K_concentration (monovalent_mEq).
C RATES(26) is d/dt VIC in component intracellular_fluid_volume (litre).
C RATES(23) is d/dt VTW in component total_body_water (litre).
C RATES(12) is d/dt HPL in component left_ventricular_hypertrophy (dimensionless).
C RATES(13) is d/dt HPR in component right_ventricular_hypertrophy (dimensionless).
C RATES(27) is d/dt HMD1 in component heart_deterioration (dimensionless).
C RATES(28) is d/dt AMM1T in component M_ST_time_delay_and_limit (dimensionless).
C RATES(29) is d/dt AMM2 in component M_LT_time_delay (dimensionless).
C RATES(31) is d/dt QOM1 in component volume_of_O2_in_M_tissue (mL).
C RATES(32) is d/dt AR1T in component NM_ST_time_delay_and_damping (dimensionless).
C RATES(33) is d/dt AR2T in component NM_I_time_delay_and_limit (dimensionless).
C RATES(34) is d/dt AR3T in component NM_LT_time_delay_and_limit (dimensionless).
C RATES(35) is d/dt QO2T in component volume_of_O2_in_NM_tissue (mL).
C RATES(36) is d/dt VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
C RATES(37) is d/dt PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
C RATES(30) is d/dt OVA in component O2_volume_of_arterial_blood (mL_per_L).
C RATES(38) is d/dt O2VAD1 in component progressive_chemoreceptor_adaptation_of_alveolar_ventilation (dimensionless).
C RATES(14) is d/dt VRC in component RBC_volume (litre).
C RATES(16) is d/dt VV7 in component short_term_stress_relaxation (litre).
C RATES(15) is d/dt VV6 in component long_term_stress_relaxation (litre).
C RATES(22) is d/dt TVD in component rate_of_fluid_intake (L_per_minute).
C RATES(39) is d/dt AHY in component time_dependent_volume_receptor_adaptation (dimensionless).
C RATES(40) is d/dt PAR1 in component perfusion_pressure (mmHg).
C RATES(41) is d/dt RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C RATES(42) is d/dt PLUR in component glomerular_urea_concentration (mOsm).
C RATES(43) is d/dt DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 1.0
CONSTS(2) = 1.0
CONSTS(3) = 1.0
CONSTS(4) = 2.5
CONSTS(5) = 12
CONSTS(6) = 0
CONSTS(7) = 0
STATES(1) = 1.0
CONSTS(8) = 60
CONSTS(9) = 5
CONSTS(10) = 0
CONSTS(11) = 0.65
CONSTS(12) = 2.5
CONSTS(13) = 0.5
CONSTS(14) = 0.8
CONSTS(15) = 15
CONSTS(16) = 0.04
STATES(2) = 0.0
CONSTS(17) = 0
CONSTS(18) = 5000
CONSTS(19) = 1
CONSTS(20) = 0
CONSTS(21) = 0
STATES(3) = 0.859476
CONSTS(22) = 12
CONSTS(23) = 1.8
CONSTS(24) = 0.7
CONSTS(25) = 0.4
CONSTS(26) = 6
CONSTS(27) = 0.8
CONSTS(28) = 0
CONSTS(29) = 5
STATES(4) = 1.0
CONSTS(30) = 139
CONSTS(31) = 85
CONSTS(32) = 0.3
CONSTS(33) = 0
CONSTS(34) = 15
CONSTS(35) = 2.5
CONSTS(36) = 0.93617
CONSTS(37) = 0.2
CONSTS(38) = 5
CONSTS(39) = 0
CONSTS(40) = 0
STATES(5) = 1.0
CONSTS(41) = 4
CONSTS(42) = 10
CONSTS(43) = 0
CONSTS(44) = 0
CONSTS(45) = 0.3
CONSTS(46) = 0.01
CONSTS(47) = 1
CONSTS(48) = 0.004
CONSTS(49) = 0.16
CONSTS(50) = -0.060024
STATES(6) = 1.00132
CONSTS(51) = 0.5
CONSTS(52) = 0
CONSTS(53) = 1
CONSTS(54) = 1.0
CONSTS(55) = 0.01
CONSTS(56) = 0.3
CONSTS(57) = 5.0
CONSTS(58) = 0.4
CONSTS(59) = 1.5
STATES(7) = 1.00007
CONSTS(60) = 2.51
CONSTS(61) = 1.5
CONSTS(62) = 0.55
CONSTS(63) = 1
CONSTS(64) = 0.1
CONSTS(65) = 3
CONSTS(66) = 1
CONSTS(67) = 0
CONSTS(68) = 0
STATES(8) = 3.00449
CONSTS(69) = 0.01167
CONSTS(70) = 0
STATES(9) = 216.243
CONSTS(71) = 40
CONSTS(72) = 0.03
CONSTS(73) = 8
CONSTS(74) = 2.728e-14
CONSTS(75) = 15
CONSTS(76) = 0.000253
CONSTS(77) = 1
CONSTS(78) = 0.15
CONSTS(79) = 0.005
STATES(10) = 0.0
STATES(11) = 279.945
CONSTS(80) = 60
CONSTS(81) = 2
CONSTS(82) = -2
CONSTS(83) = 0.8092
CONSTS(84) = 6
CONSTS(85) = 4.2
STATES(12) = 1.00163
STATES(13) = 1.00237
STATES(14) = 2.00439
STATES(15) = 0.0101913
STATES(16) = 0.00366525
STATES(17) = 3.28246
STATES(18) = 0.862514
STATES(19) = 0.379883
STATES(20) = 0.38131
STATES(21) = 0.100043
CONSTS(86) = 0.4
CONSTS(87) = 0.15
CONSTS(88) = 1
CONSTS(89) = 1
CONSTS(90) = -0.2
CONSTS(91) = 0.1
CONSTS(92) = 0
CONSTS(93) = 0.2
CONSTS(94) = 0.0212
CONSTS(95) = 1
CONSTS(96) = 2
CONSTS(97) = 30.52
CONSTS(98) = 1
CONSTS(99) = 96.3
CONSTS(100) = 0
STATES(22) = 0.000980838
STATES(23) = 39.8952
CONSTS(101) = 0.1
CONSTS(102) = 0
STATES(24) = 2109.91
CONSTS(103) = 0.3
STATES(25) = 3622.54
CONSTS(104) = 0.08
STATES(26) = 25.0404
CONSTS(105) = 0.01
CONSTS(106) = 1
CONSTS(107) = 0.625
CONSTS(108) = 1
CONSTS(109) = 0.05
STATES(27) = 1.0
CONSTS(110) = 0
CONSTS(111) = 0.04
CONSTS(112) = 0.1
CONSTS(113) = 0.005
STATES(28) = 1.00269
CONSTS(114) = 2
STATES(29) = 1.09071
CONSTS(115) = 40000
STATES(30) = 204.497
CONSTS(116) = 1
CONSTS(117) = 0.17
CONSTS(118) = 57.1
CONSTS(119) = 30
STATES(31) = 48.0839
CONSTS(120) = 0.79167
CONSTS(121) = 35
CONSTS(122) = 0.1
CONSTS(123) = 0.5
STATES(32) = 1.02127
CONSTS(124) = 0.1
CONSTS(125) = 60
STATES(33) = 1.01179
CONSTS(126) = 2
CONSTS(127) = 40000
STATES(34) = 1.1448
CONSTS(128) = 0.9
CONSTS(129) = 164
STATES(35) = 72.2362
CONSTS(130) = 0.0003
STATES(36) = 0.0123238
STATES(37) = 0.419998
CONSTS(131) = 1
CONSTS(132) = 150
CONSTS(133) = 1.8
CONSTS(134) = 0.0175
STATES(38) = 2.368e-07
CONSTS(135) = 90
CONSTS(136) = 0.53333
CONSTS(137) = 150
CONSTS(138) = 1850
CONSTS(139) = 4.714e-08
CONSTS(140) = 1
CONSTS(141) = 5.8e-06
CONSTS(142) = 0
CONSTS(143) = 1
CONSTS(144) = 5
CONSTS(145) = 1
CONSTS(146) = 10000
CONSTS(147) = 2
CONSTS(148) = 45
CONSTS(149) = 0.01
CONSTS(150) = 2
CONSTS(151) = 1.5
CONSTS(152) = 0
CONSTS(153) = 30
CONSTS(154) = 0.333
CONSTS(155) = 1
STATES(39) = 0.301963
CONSTS(156) = 1000
CONSTS(157) = 0
CONSTS(158) = 0
CONSTS(159) = 0
CONSTS(160) = 0
CONSTS(161) = 0
CONSTS(162) = 0
CONSTS(163) = 2000
STATES(40) = 103.525
CONSTS(164) = 0.6
CONSTS(165) = 0.3
CONSTS(166) = 10
CONSTS(167) = 0
STATES(41) = 0.0
CONSTS(168) = 0.5
CONSTS(169) = 1.4
CONSTS(170) = 0.86
CONSTS(171) = 1
CONSTS(172) = 1.5
CONSTS(173) = 4
CONSTS(174) = 0.3
CONSTS(175) = 1.6
CONSTS(176) = 0
CONSTS(177) = 1
CONSTS(178) = 24
CONSTS(179) = 1
CONSTS(180) = 1.0
CONSTS(181) = 1.0
CONSTS(182) = 8
CONSTS(183) = 0.0208333
CONSTS(184) = 0.001
CONSTS(185) = 10
CONSTS(186) = 0.9
CONSTS(187) = 15.2
STATES(42) = 159.549
CONSTS(188) = 0.24
CONSTS(189) = 0.8
CONSTS(190) = 19.167
CONSTS(191) = 6
CONSTS(192) = 0.5
CONSTS(193) = 1
CONSTS(194) = 0.3
CONSTS(195) = 0.675
CONSTS(196) = 1
CONSTS(197) = 0.3
CONSTS(198) = 1e-06
CONSTS(199) = 1
CONSTS(200) = 0.03
CONSTS(201) = 0.667
STATES(43) = 0.0367573
CONSTS(202) = 2
CONSTS(203) = 0.3
CONSTS(204) = 4
CONSTS(205) = 15
CONSTS(206) = 3
CONSTS(207) = CONSTS(53) ** CONSTS(54)
CONSTS(208) = CONSTS(110) ** 0.500000*5.00000
CONSTS(209) = TERNRY(CONSTS(137).GT.80.0000, 80.0000, CONSTS(137))
CONSTS(210) = CONSTS(137) - 40.0000
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(32) = ((((((STATES(8)+STATES(14)) - STATES(17)) - STATES(18)) - STATES(19)) - STATES(20)) - STATES(21))/2.00000
ALGBRC(46) = STATES(19)+ ALGBRC(32)*0.128000
ALGBRC(47) = ALGBRC(46) - 0.380000
ALGBRC(48) = ALGBRC(47)/0.0100000
ALGBRC(49) = TERNRY( (ALGBRC(48) - 1.00000)*1.00000.LT.0.00000, 0.00000, (ALGBRC(48) - 1.00000)*1.00000)
ALGBRC(33) = STATES(21)+ ALGBRC(32)*0.0574000
ALGBRC(34) = ALGBRC(33) - 0.100000
ALGBRC(35) = ALGBRC(34)/0.00500000
ALGBRC(36) = TERNRY( (ALGBRC(35)+1.00000)*2.00000.LT.0.00000, 0.00000, (ALGBRC(35)+1.00000)*2.00000)
ALGBRC(51) = (ALGBRC(49)+ALGBRC(36))/3.00000
ALGBRC(53) = TERNRY(CONSTS(39).GT.0.00000, CONSTS(39), ALGBRC(51)+CONSTS(40))
RATES(5) = (ALGBRC(53) - STATES(5))/CONSTS(41)
ALGBRC(63) = STATES(18)+ ALGBRC(32)*0.261000
ALGBRC(64) = ALGBRC(63) - 0.495000
ALGBRC(65) = ALGBRC(64)/0.00355000
RATES(40) = ((100.000+ (ALGBRC(65) - 100.000)*CONSTS(162)) - STATES(40))/CONSTS(163)
ALGBRC(66) = TERNRY(CONSTS(43).GT.1.00000e-07, CONSTS(43), ALGBRC(65) - CONSTS(44))
ALGBRC(67) = TERNRY(ALGBRC(66).LT.160.000.AND.ALGBRC(66).GE.80.0000, 0.0166670*(160.000 - ALGBRC(66)), TERNRY(ALGBRC(66).LT.80.0000, 1.33360, 0.00000)
ALGBRC(69) = (ALGBRC(67) - 1.00000)*CONSTS(47)+1.00000
ALGBRC(71) = ALGBRC(69) - CONSTS(50)
RATES(6) = (ALGBRC(71) - STATES(6))/CONSTS(49)
ALGBRC(82) = STATES(23) - STATES(26)
ALGBRC(83) = STATES(24)/ALGBRC(82)
ALGBRC(1) = CONSTS(9) - (CONSTS(9) - 1.00000)/( ((CONSTS(10) - 1.00000)/(CONSTS(10) - CONSTS(9)))*(STATES(1) - 1.00000)*CONSTS(11)+1.00000)
ALGBRC(2) = (ALGBRC(1) - 1.00000)*CONSTS(12)+1.00000
ALGBRC(3) = (ALGBRC(2) - 1.00000)*CONSTS(13)+1.00000
ALGBRC(4) = TERNRY(ALGBRC(3).LT.0.200000, 0.200000, ALGBRC(3))
ALGBRC(78) = (ALGBRC(4) - 1.00000)*CONSTS(103)+1.00000
ALGBRC(79) = (STATES(25) - 3000.00)/( ALGBRC(78)*9.33330)
ALGBRC(80) = STATES(25) - ALGBRC(79)
ALGBRC(81) = ALGBRC(80)/STATES(26)
ALGBRC(86) = ALGBRC(81) - ALGBRC(83)
ALGBRC(90) = ALGBRC(86)*CONSTS(105)
RATES(26) = ALGBRC(90)
ALGBRC(84) = (ALGBRC(83) - CONSTS(30))/(142.000 - CONSTS(30))
ALGBRC(88) = TERNRY(ALGBRC(84).LT.0.00000, 0.00000, ALGBRC(84))
ALGBRC(68) = TERNRY(ALGBRC(66).GT.CONSTS(31), CONSTS(31), ALGBRC(66))
ALGBRC(70) = CONSTS(31) - ALGBRC(68) ** 2.00000*CONSTS(32)
ALGBRC(92) = ALGBRC(88)+ALGBRC(70)+CONSTS(33)
ALGBRC(95) = TERNRY(ALGBRC(92).LT.0.00000, 0.00000, ALGBRC(92))
RATES(4) = (ALGBRC(95) - STATES(4))/CONSTS(34)
ALGBRC(7) = CONSTS(23) - (CONSTS(23) - 1.00000)/( ((CONSTS(24) - 1.00000)/(CONSTS(24) - CONSTS(23)))*(STATES(3) - 1.00000)*CONSTS(25)+1.00000)
ALGBRC(8) = (ALGBRC(7) - 1.00000)*CONSTS(4)+1.00000
ALGBRC(85) = ALGBRC(79)/ALGBRC(82)
ALGBRC(87) = (ALGBRC(85) - 3.30000)/1.00000
ALGBRC(91) = ALGBRC(8)*0.909000*ALGBRC(87)
ALGBRC(94) = (ALGBRC(91) - 1.00000)*CONSTS(5)+1.00000
ALGBRC(97) = TERNRY(ALGBRC(94).LT.0.00000, 0.00000, ALGBRC(94))
ALGBRC(99) = TERNRY(CONSTS(7).GT.0.00000, CONSTS(7), ALGBRC(97)+CONSTS(6))
RATES(1) = (ALGBRC(99) - STATES(1))/CONSTS(8)
ALGBRC(102) = TERNRY(STATES(31).LT.0.000100000, 0.000100000, STATES(31))
ALGBRC(103) = CONSTS(120)*ALGBRC(102)
ALGBRC(105) = ALGBRC(103) - 38.0000
ALGBRC(107) = ALGBRC(105)*CONSTS(111)+1.00000
RATES(28) = ( ALGBRC(107)*1.00000 - STATES(28))/CONSTS(112)
ALGBRC(108) = CONSTS(114)*ALGBRC(105)+1.00000
RATES(29) = ( ALGBRC(108)*1.00000 - STATES(29))/CONSTS(115)
ALGBRC(114) = TERNRY(STATES(35).LT.0.00000, 0.00000, STATES(35))
ALGBRC(115) = ALGBRC(114)*0.486110
ALGBRC(117) = (ALGBRC(115) - 10.0000)*CONSTS(109)
RATES(27) = ALGBRC(117)
ALGBRC(118) = ALGBRC(115) - CONSTS(121)
ALGBRC(120) = ALGBRC(118)*CONSTS(122)+1.00000
RATES(32) = ( ALGBRC(120)*1.00000 - STATES(32))/CONSTS(123)
ALGBRC(121) = CONSTS(124)*ALGBRC(118)+1.00000
RATES(33) = ( ALGBRC(121)*1.00000 - STATES(33))/CONSTS(125)
ALGBRC(122) = CONSTS(126)*ALGBRC(118)+1.00000
RATES(34) = ( ALGBRC(122)*1.00000 - STATES(34))/CONSTS(127)
ALGBRC(124) = TERNRY(STATES(36).LT.0.00100000, 0.00100000, STATES(36))
ALGBRC(125) = (ALGBRC(82) - STATES(8)) - ALGBRC(124)
RATES(10) = ( (ALGBRC(125) - 12.0000)*CONSTS(78) - STATES(10))*CONSTS(79)
ALGBRC(28) = STATES(9)/STATES(8)
ALGBRC(31) = 0.280000*ALGBRC(28)+ 0.00190000*ALGBRC(28) ** 2.00000
ALGBRC(39) = STATES(20)+ ALGBRC(32)*0.155000
ALGBRC(40) = ALGBRC(39) - 0.306250
ALGBRC(41) = ALGBRC(40)/0.00480000
ALGBRC(42) = 0.0260000*ALGBRC(41)
ALGBRC(43) = TERNRY(ALGBRC(42).LT.1.00000e-05, 1.00000e-05, ALGBRC(42))
ALGBRC(44) = ALGBRC(43) ** 0.500000
ALGBRC(45) = 1.00000/ALGBRC(44)
ALGBRC(52) = ALGBRC(48)+18.0000
ALGBRC(54) = 1.00000/( ALGBRC(52)*0.0357000)
ALGBRC(123) = ( (ALGBRC(41) - ALGBRC(48))*ALGBRC(54))/(ALGBRC(54)+ALGBRC(45))+ALGBRC(48)
ALGBRC(143) = TERNRY(STATES(37).LT.0.0250000, 0.0250000, STATES(37))
ALGBRC(144) = ALGBRC(143)/ALGBRC(124)
ALGBRC(145) = ALGBRC(144)*0.400000
ALGBRC(142) = 2.00000 - 0.150000/ALGBRC(124)
ALGBRC(146) = (((ALGBRC(123) - ALGBRC(142))+ALGBRC(145)) - ALGBRC(31))*CONSTS(130)
ALGBRC(148) = (ALGBRC(142)+11.0000)*0.000300000
ALGBRC(149) = ALGBRC(146) - ALGBRC(148)
ALGBRC(150) = ALGBRC(149)
RATES(36) = ALGBRC(150)
ALGBRC(151) = ALGBRC(148)*ALGBRC(144)
ALGBRC(147) = (ALGBRC(28) - ALGBRC(144))*0.000225000
ALGBRC(152) = ALGBRC(147) - ALGBRC(151)
ALGBRC(153) = ALGBRC(152)
RATES(37) = ALGBRC(153)
ALGBRC(19) = CONSTS(57) - (CONSTS(57) - 1.00000)/EXP( CONSTS(59)*(STATES(7) - 1.00000))
ALGBRC(20) = TERNRY(ALGBRC(19).LT.CONSTS(58), CONSTS(58), ALGBRC(19))
ALGBRC(22) = ALGBRC(20) - 1.00000
ALGBRC(23) = ALGBRC(22)*CONSTS(62)+1.00000
ALGBRC(89) = TERNRY(STATES(27).GT.1.00000, 1.00000, STATES(27))
ALGBRC(58) = (ALGBRC(48)+4.00000)*( CONSTS(86)*(ALGBRC(20) - 1.00000)+1.00000) - 4.00000
ALGBRC(59) = TERNRY(ALGBRC(58).LE.- 2.00000, 0.0100000, TERNRY(ALGBRC(58).GT.- 2.00000.AND.ALGBRC(58).LE.1.00000, 0.0100000+( (3.60000 - 0.0100000)*(ALGBRC(58) - - 2.00000))/(1.00000 - - 2.00000), TERNRY(ALGBRC(58).GT.1.00000.AND.ALGBRC(58).LE.5.00000, 3.60000+( (9.40000 - 3.60000)*(ALGBRC(58) - 1.00000))/(5.00000 - 1.00000), TERNRY(ALGBRC(58).GT.5.00000.AND.ALGBRC(58).LE.8.00000, 9.40000+( (11.6000 - 9.40000)*(ALGBRC(58) - 5.00000))/(8.00000 - 5.00000), TERNRY(ALGBRC(58).GT.8.00000.AND.ALGBRC(58).LE.12.0000, 11.6000+( (13.5000 - 11.6000)*(ALGBRC(58) - 8.00000))/(12.0000 - 8.00000), 13.5000)
ALGBRC(156) = STATES(8)+STATES(14)
ALGBRC(157) = STATES(14)/ALGBRC(156)
ALGBRC(158) = 100.000*ALGBRC(157)
ALGBRC(159) = (STATES(30)/ALGBRC(158))/5.25000
ALGBRC(162) = ALGBRC(65)/( ALGBRC(23)*ALGBRC(159))
ALGBRC(163) = TERNRY(ALGBRC(162).LE.0.00000, 1.04000, TERNRY(ALGBRC(162).GT.0.00000.AND.ALGBRC(162).LE.60.0000, 1.04000+( (1.02500 - 1.04000)*(ALGBRC(162) - 0.00000))/(60.0000 - 0.00000), TERNRY(ALGBRC(162).GT.60.0000.AND.ALGBRC(162).LE.125.000, 1.02500+( (0.970000 - 1.02500)*(ALGBRC(162) - 60.0000))/(125.000 - 60.0000), TERNRY(ALGBRC(162).GT.125.000.AND.ALGBRC(162).LE.160.000, 0.970000+( (0.880000 - 0.970000)*(ALGBRC(162) - 125.000))/(160.000 - 125.000), TERNRY(ALGBRC(162).GT.160.000.AND.ALGBRC(162).LE.200.000, 0.880000+( (0.590000 - 0.880000)*(ALGBRC(162) - 160.000))/(200.000 - 160.000), TERNRY(ALGBRC(162).GT.200.000.AND.ALGBRC(162).LE.240.000, 0.590000+( (0.00000 - 0.590000)*(ALGBRC(162) - 200.000))/(240.000 - 200.000), 0.00000)
ALGBRC(164) = ALGBRC(163)*ALGBRC(59)*ALGBRC(23)*CONSTS(89)*ALGBRC(89)*STATES(12)
ALGBRC(74) = (ALGBRC(48) - ALGBRC(65))/3.00000
ALGBRC(165) = TERNRY(ALGBRC(74).GT.0.00000, ALGBRC(164)+ALGBRC(74), ALGBRC(164))
ALGBRC(55) = ALGBRC(54)+ALGBRC(45)
ALGBRC(56) = ALGBRC(41) - ALGBRC(48)
ALGBRC(57) = ALGBRC(56)/ALGBRC(55)
ALGBRC(167) = ALGBRC(57) - ALGBRC(165)
RATES(19) = ALGBRC(167)
ALGBRC(160) = (ALGBRC(41)/ALGBRC(23))/ALGBRC(159)
ALGBRC(161) = TERNRY(ALGBRC(160).LE.0.00000, 1.06000, TERNRY(ALGBRC(160).GT.0.00000.AND.ALGBRC(160).LE.32.0000, 1.06000+( (0.970000 - 1.06000)*(ALGBRC(160) - 0.00000))/(32.0000 - 0.00000), TERNRY(ALGBRC(160).GT.32.0000.AND.ALGBRC(160).LE.38.4000, 0.970000+( (0.930000 - 0.970000)*(ALGBRC(160) - 32.0000))/(38.4000 - 32.0000), TERNRY(ALGBRC(160).GT.38.4000.AND.ALGBRC(160).LE.48.0000, 0.930000+( (0.800000 - 0.930000)*(ALGBRC(160) - 38.4000))/(48.0000 - 38.4000), TERNRY(ALGBRC(160).GT.48.0000.AND.ALGBRC(160).LE.60.8000, 0.800000+( (0.460000 - 0.800000)*(ALGBRC(160) - 48.0000))/(60.8000 - 48.0000), TERNRY(ALGBRC(160).GT.60.8000.AND.ALGBRC(160).LE.72.0000, 0.460000+( (0.00000 - 0.460000)*(ALGBRC(160) - 60.8000))/(72.0000 - 60.8000), 0.00000)
ALGBRC(166) = (1.00000 - CONSTS(87))*ALGBRC(23)*ALGBRC(161)*CONSTS(88)*ALGBRC(89)*STATES(13)+( CONSTS(87)*ALGBRC(165))/ALGBRC(59)
ALGBRC(37) = (ALGBRC(35)+8.00000)*( CONSTS(86)*(ALGBRC(20) - 1.00000)+1.00000) - 8.00000
ALGBRC(38) = TERNRY(ALGBRC(37).LE.- 8.00000, 0.00000, TERNRY(ALGBRC(37).GT.- 8.00000.AND.ALGBRC(37).LE.- 6.00000, 0.00000+( (0.750000 - 0.00000)*(ALGBRC(37) - - 8.00000))/(- 6.00000 - - 8.00000), TERNRY(ALGBRC(37).GT.- 6.00000.AND.ALGBRC(37).LE.- 2.00000, 0.750000+( (2.60000 - 0.750000)*(ALGBRC(37) - - 6.00000))/(- 2.00000 - - 6.00000), TERNRY(ALGBRC(37).GT.- 2.00000.AND.ALGBRC(37).LE.4.00000, 2.60000+( (9.80000 - 2.60000)*(ALGBRC(37) - - 2.00000))/(4.00000 - - 2.00000), TERNRY(ALGBRC(37).GT.4.00000.AND.ALGBRC(37).LE.12.0000, 9.80000+( (13.5000 - 9.80000)*(ALGBRC(37) - 4.00000))/(12.0000 - 4.00000), 13.5000)
ALGBRC(169) = ALGBRC(38)*ALGBRC(166)
ALGBRC(171) = ALGBRC(169) - ALGBRC(57)
RATES(20) = ALGBRC(171)
ALGBRC(170) = TERNRY(ALGBRC(159).GT.1.00000, 114.000+ (ALGBRC(159) - 1.00000)*6667.00, TERNRY(ALGBRC(159).GT.0.936000.AND.ALGBRC(159).LE.1.00000, 74.0000+ (ALGBRC(159) - 0.936000)*625.000, TERNRY(ALGBRC(159).GT.0.800000.AND.ALGBRC(159).LE.0.936000, 46.0000+ (ALGBRC(159) - 0.800000)*205.882, ALGBRC(159)*57.5000)
ALGBRC(173) = (ALGBRC(170) - 67.0000)/30.0000
ALGBRC(175) = TERNRY(ALGBRC(173).GT.1.00000, 1.00000, TERNRY(ALGBRC(173).LT.0.600000, 0.600000, ALGBRC(173))
ALGBRC(177) = 1.00000/ALGBRC(175)
ALGBRC(180) = ( (ALGBRC(177) - 1.00000)*3.00000 - STATES(38))*0.000500000
RATES(38) = ALGBRC(180)
ALGBRC(72) = TERNRY(ALGBRC(66).LT.80.0000.AND.ALGBRC(66).GE.40.0000, 0.00500000*(80.0000 - ALGBRC(66))*CONSTS(45), TERNRY(ALGBRC(66).LT.40.0000, 0.200000*CONSTS(45), 0.00000)
ALGBRC(172) = TERNRY(ALGBRC(170).LT.80.0000.AND.ALGBRC(170).GE.40.0000, CONSTS(46)*(80.0000 - ALGBRC(170)), TERNRY(ALGBRC(170).LT.40.0000, CONSTS(46)*40.0000, 0.00000)
ALGBRC(174) = ALGBRC(72)+ALGBRC(172)
ALGBRC(18) = STATES(6)
ALGBRC(73) = TERNRY(ALGBRC(66).LT.40.0000, 0.0400000*(40.0000 - ALGBRC(66))*CONSTS(51), 0.00000)
ALGBRC(50) = (15.0000/(ALGBRC(48)+ALGBRC(35)+ALGBRC(41)) - 1.00000)*CONSTS(52)+1.00000
ALGBRC(176) = ( CONSTS(207)*ALGBRC(50)*(ALGBRC(174)+ALGBRC(18)+ALGBRC(73)) - 1.00000)*CONSTS(55)+1.00000
ALGBRC(178) = TERNRY(ALGBRC(176).LT.0.00000, 0.00000, ALGBRC(176))
ALGBRC(181) = (ALGBRC(178) - STATES(7))/CONSTS(56)
RATES(7) = ALGBRC(181)
ALGBRC(190) = ABS(ALGBRC(35)) ** CONSTS(154)*CONSTS(155)
ALGBRC(193) = TERNRY(ALGBRC(35).LT.0.00000, - ALGBRC(190), ALGBRC(190))
RATES(39) = (ALGBRC(193) - STATES(39))/CONSTS(156)
ALGBRC(185) = (CONSTS(209) - 40.0000)*ALGBRC(158)
ALGBRC(187) = TERNRY(ALGBRC(185)+CONSTS(210).LT.0.00000, 0.00000, ALGBRC(185)+CONSTS(210))
ALGBRC(189) = CONSTS(138) - ALGBRC(187)
ALGBRC(191) = TERNRY( ALGBRC(189)*CONSTS(139)*CONSTS(140)+5.00000e-06.LT.0.00000, 0.00000, ALGBRC(189)*CONSTS(139)*CONSTS(140)+5.00000e-06)
ALGBRC(179) = ALGBRC(158)/( (CONSTS(135) - ALGBRC(158))*CONSTS(136))
ALGBRC(182) = ALGBRC(179)+1.50000
ALGBRC(183) = 0.333300*ALGBRC(182)
ALGBRC(194) = STATES(14)*CONSTS(141)*ALGBRC(183)
ALGBRC(197) = (ALGBRC(191) - ALGBRC(194))+CONSTS(142)
RATES(14) = ALGBRC(197)
ALGBRC(195) = (ALGBRC(7) - 1.00000)*CONSTS(151)*0.00100000
ALGBRC(184) = (ALGBRC(7) - 1.00000)*CONSTS(147)+1.00000
ALGBRC(186) = CONSTS(148) - ALGBRC(115) ** 2.00000*CONSTS(149)*ALGBRC(184)
ALGBRC(188) = TERNRY(ALGBRC(186).LT.0.800000, 0.800000, TERNRY(ALGBRC(186).GT.8.00000, 8.00000, ALGBRC(186))
ALGBRC(192) = (STATES(4) - 1.00000)*CONSTS(150)+1.00000
ALGBRC(198) = ALGBRC(192)*ALGBRC(188)*0.00100000
ALGBRC(200) = TERNRY(ALGBRC(198).LT.0.00000, 0.00000, ALGBRC(198))
ALGBRC(202) = (ALGBRC(195)+ALGBRC(200))*1.00000
ALGBRC(204) = TERNRY(ALGBRC(202).LT.0.00000, 0.00000, ALGBRC(202))
RATES(22) = ((ALGBRC(204)+CONSTS(152)) - STATES(22))/CONSTS(153)
ALGBRC(21) = (CONSTS(60) - ALGBRC(20)*CONSTS(61))+CONSTS(61)
ALGBRC(196) = ALGBRC(193) - STATES(39)
ALGBRC(206) = ALGBRC(196)*CONSTS(158)
ALGBRC(60) = STATES(17)+ ALGBRC(32)*0.398600
ALGBRC(9) = (ALGBRC(7) - 1.00000)*CONSTS(26)+1.00000
ALGBRC(10) = TERNRY(ALGBRC(9).LT.CONSTS(27), CONSTS(27), ALGBRC(9))
ALGBRC(61) = (ALGBRC(10) - 1.00000)*CONSTS(90)
ALGBRC(207) = ((((ALGBRC(60) - ALGBRC(21)) - ALGBRC(61)) - STATES(16)) - STATES(15)) - ALGBRC(206)
ALGBRC(208) = TERNRY(ALGBRC(207).LT.0.000100000, 0.000100000, ALGBRC(207))
RATES(16) = ( (ALGBRC(208) - 0.740000)*CONSTS(143) - STATES(16))/CONSTS(144)
RATES(15) = ( (ALGBRC(208) - 0.740000)*CONSTS(145) - STATES(15))/CONSTS(146)
ALGBRC(209) = 3.70000+(ALGBRC(208) - 0.740000)/CONSTS(91)
ALGBRC(210) = TERNRY(ALGBRC(209).LT.0.000100000, 0.000100000, ALGBRC(209))
ALGBRC(211) = 0.740000/ALGBRC(210)/( ALGBRC(183)*3.70000) ** 0.500000
ALGBRC(62) = TERNRY(ALGBRC(35).LT.CONSTS(92), CONSTS(92), ALGBRC(35))
ALGBRC(212) = ALGBRC(210) - ALGBRC(62)
ALGBRC(213) = ALGBRC(212)/ALGBRC(211)
ALGBRC(215) = ALGBRC(213) - ALGBRC(169)
RATES(21) = ALGBRC(215)
ALGBRC(214) = ALGBRC(65) - ALGBRC(210)
ALGBRC(100) = TERNRY(STATES(28).LT.CONSTS(113), CONSTS(113), STATES(28))
ALGBRC(101) = ALGBRC(100)*STATES(29)
ALGBRC(12) = CONSTS(35) - (CONSTS(35) - 1.00000)/( ((CONSTS(36) - 1.00000)/(CONSTS(36) - CONSTS(35)))*(STATES(4) - 1.00000)+1.00000)
ALGBRC(13) = TERNRY(ALGBRC(12).LT.CONSTS(36), CONSTS(36), ALGBRC(12))
ALGBRC(26) = ALGBRC(22)*CONSTS(65)+1.00000 ** CONSTS(66)
ALGBRC(199) = ALGBRC(196)*CONSTS(157)+1.00000
ALGBRC(76) = ALGBRC(65)/100.000 ** CONSTS(96)
ALGBRC(201) = (( ALGBRC(10)*ALGBRC(13)*ALGBRC(26)*ALGBRC(183)*CONSTS(2))/ALGBRC(76))/ALGBRC(199)
ALGBRC(205) = CONSTS(99)*ALGBRC(101)*ALGBRC(201)*CONSTS(1)*CONSTS(98)
ALGBRC(216) = ALGBRC(214)/ALGBRC(205)
ALGBRC(217) = STATES(30)*ALGBRC(216)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
ALGBRC(25) = ALGBRC(22)*CONSTS(64)+1.00000
ALGBRC(104) = TERNRY(ALGBRC(103).GT.38.0000, 38.0000, ALGBRC(103))
ALGBRC(106) = ALGBRC(25)*CONSTS(118)*CONSTS(116)*(1.00000 - 38.0001 - ALGBRC(104) ** 3.00000/54872.0)
ALGBRC(262) = ALGBRC(230) - ALGBRC(106)
RATES(31) = ALGBRC(262)
ALGBRC(11) = (ALGBRC(10) - 1.00000)*CONSTS(28)+1.00000
ALGBRC(16) = CONSTS(42) - CONSTS(42)/( 0.555556*(1.00000+STATES(5)))
ALGBRC(17) = TERNRY(ALGBRC(16).LT.- 1.00000, - 1.00000, ALGBRC(16))
ALGBRC(27) = ALGBRC(22)*CONSTS(67)+1.00000
ALGBRC(109) = TERNRY(STATES(32).LT.0.500000, 0.500000, STATES(32))
ALGBRC(110) = TERNRY(STATES(33).LT.0.500000, 0.500000, STATES(33))
ALGBRC(111) = TERNRY(STATES(34).LT.0.300000, 0.300000, STATES(34))
ALGBRC(112) = ALGBRC(109)*ALGBRC(110)*ALGBRC(111)
ALGBRC(113) = (ALGBRC(112) - 1.00000)*CONSTS(128)+1.00000
ALGBRC(203) = CONSTS(97)*ALGBRC(113)*ALGBRC(201)*CONSTS(1)*CONSTS(98)
ALGBRC(218) = TERNRY(CONSTS(160).GT.0.00000.AND.CONSTS(161).LE.0.00000, CONSTS(160), TERNRY(CONSTS(161).GT.0.00000, STATES(40), ALGBRC(65) - CONSTS(159))
ALGBRC(219) = (ALGBRC(26) - 1.00000)*CONSTS(168)+1.00000
ALGBRC(220) = TERNRY(ALGBRC(219).LT.0.800000, 0.800000, ALGBRC(219))
ALGBRC(221) = (ALGBRC(7) - 1.00000)*CONSTS(169)+1.00000
ALGBRC(222) = TERNRY(ALGBRC(221).LT.CONSTS(170), CONSTS(170), ALGBRC(221))
ALGBRC(223) = (ALGBRC(220) - 1.00000)*CONSTS(174)+1.00000
ALGBRC(224) = (ALGBRC(7) - 1.00000)*CONSTS(175)+1.00000
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
RATES(41) = (ALGBRC(241) - 1.00000)*CONSTS(167)
ALGBRC(263) = ALGBRC(233)
ALGBRC(269) = TERNRY(ALGBRC(263).GT.1.00000, 1.00000/(1.00000+ (ALGBRC(263) - 1.00000)*72.0000), 10.0000 - 9.00000/(1.00000+ (1.00000 - ALGBRC(263))*8.00000))
ALGBRC(276) = (ALGBRC(269) - 1.00000)*CONSTS(17)
RATES(2) = (ALGBRC(276) - STATES(2))/CONSTS(18)
ALGBRC(75) = ALGBRC(65) - ALGBRC(35)
ALGBRC(77) = ALGBRC(75)*CONSTS(100)
ALGBRC(265) = CONSTS(179)*ALGBRC(251)
ALGBRC(271) = ALGBRC(216)+ALGBRC(235)+ALGBRC(265)
ALGBRC(280) = ALGBRC(271)+ALGBRC(77)
RATES(12) = (( ALGBRC(65)*ALGBRC(280))/( 500.000*CONSTS(106)) ** CONSTS(107) - STATES(12))/57600.0
RATES(13) = (( ALGBRC(41)*ALGBRC(280))/( 75.0000*CONSTS(108)) ** CONSTS(107) - STATES(13))/57600.0
ALGBRC(126) = ALGBRC(125) - STATES(10)
ALGBRC(134) = (ALGBRC(126) - CONSTS(84))/CONSTS(84) ** 2.00000*1.00000
ALGBRC(131) = (CONSTS(80)/ALGBRC(125))/5.00000 ** CONSTS(81)
ALGBRC(135) = ALGBRC(131)*CONSTS(82)+ALGBRC(134)
ALGBRC(128) = STATES(11)/ALGBRC(125)
ALGBRC(130) = 0.280000*ALGBRC(128)+ 0.00190000*ALGBRC(128) ** 2.00000
ALGBRC(132) = ALGBRC(131)*2.00000
ALGBRC(133) = ALGBRC(132)*ALGBRC(130)*CONSTS(83)
ALGBRC(264) = TERNRY(ALGBRC(236) - CONSTS(75).LT.0.00000, 0.00000, ALGBRC(236) - CONSTS(75))
ALGBRC(270) = ALGBRC(264)*CONSTS(76) ** CONSTS(77)
ALGBRC(278) = (((ALGBRC(236) - ALGBRC(31)) - ALGBRC(135))+ALGBRC(133))*CONSTS(69)+ALGBRC(270)
ALGBRC(136) = ALGBRC(135) - ALGBRC(132)
ALGBRC(137) = (ALGBRC(136)+CONSTS(85)) - ALGBRC(134)
ALGBRC(139) = TERNRY(ALGBRC(137).GT.7.00000, 7.00000, ALGBRC(137))
ALGBRC(140) = TERNRY(ALGBRC(139).LT.0.00000, 0.00000, ALGBRC(139)*0.0200000)
ALGBRC(285) = ((ALGBRC(140) - ALGBRC(278)) - ALGBRC(150))+CONSTS(70)
RATES(8) = ALGBRC(285)
ALGBRC(29) = TERNRY(ALGBRC(28) - CONSTS(71).LT.0.00000, 0.00000, ALGBRC(28) - CONSTS(71))
ALGBRC(30) = CONSTS(72) - ALGBRC(29) ** CONSTS(73)*CONSTS(74)
ALGBRC(279) = ALGBRC(270)*ALGBRC(28)+ (ALGBRC(28) - ALGBRC(128))*0.00104000
ALGBRC(141) = ALGBRC(128)*ALGBRC(140)
ALGBRC(286) = (((ALGBRC(30)+ALGBRC(141)) - ALGBRC(279)) - ALGBRC(153))+ CONSTS(70)*72.0000
RATES(9) = ALGBRC(286)
ALGBRC(287) = ALGBRC(279) - ALGBRC(141)
RATES(11) = ALGBRC(287)
ALGBRC(288) = ALGBRC(280) - ALGBRC(213)
RATES(17) = ALGBRC(288)
ALGBRC(289) = ALGBRC(165) - ALGBRC(280)
RATES(18) = ALGBRC(289)
ALGBRC(266) = STATES(30)*ALGBRC(235)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
ALGBRC(116) = TERNRY(ALGBRC(115).GT.35.0000, 35.0000, ALGBRC(115))
ALGBRC(119) = ALGBRC(25)*CONSTS(129)*(1.00000 - 35.0001 - ALGBRC(116) ** 3.00000/42875.0)
ALGBRC(281) = ALGBRC(272) - ALGBRC(119)
ALGBRC(291) = TERNRY(ALGBRC(114).LT.6.00000.AND.ALGBRC(281).LT.0.00000, ALGBRC(281)*0.100000, ALGBRC(281))
RATES(35) = ALGBRC(291)
ALGBRC(277) = (ALGBRC(269)+STATES(2))*CONSTS(19)
ALGBRC(284) = TERNRY(ALGBRC(277).LT.1.00000e-05, 1.00000e-05, ALGBRC(277))
ALGBRC(294) = TERNRY(CONSTS(20).GT.0.00000, CONSTS(20), ALGBRC(284)+CONSTS(21))
RATES(3) = (ALGBRC(294) - STATES(3))/CONSTS(22)
ALGBRC(282) = ALGBRC(272)+ALGBRC(230)
ALGBRC(155) = STATES(38)+1.00000
ALGBRC(292) = ALGBRC(282)*CONSTS(131)*0.0266670*ALGBRC(177)*ALGBRC(155)
ALGBRC(295) = CONSTS(132) - (ALGBRC(282)/ALGBRC(292))/0.761000
ALGBRC(154) = CONSTS(133)/(CONSTS(134)+ALGBRC(124))
ALGBRC(297) = (ALGBRC(295) - ALGBRC(170))*ALGBRC(154)
ALGBRC(299) = (ALGBRC(297) - ALGBRC(282))/( ALGBRC(169)*1.00000)
RATES(30) = ALGBRC(299)
ALGBRC(302) = ALGBRC(233)*ALGBRC(83)*0.00616190
ALGBRC(5) = (ALGBRC(2) - 1.00000)*CONSTS(14)+1.00000
ALGBRC(6) = TERNRY(ALGBRC(5).LT.CONSTS(16), CONSTS(16), TERNRY(ALGBRC(5).GT.CONSTS(15), CONSTS(15), ALGBRC(5))
ALGBRC(14) = CONSTS(38) - (CONSTS(38) - 1.00000)/( ((CONSTS(37) - 1.00000)/(CONSTS(37) - CONSTS(38)))*(STATES(4) - 1.00000)+1.00000)
ALGBRC(15) = TERNRY(ALGBRC(14).LT.CONSTS(37), CONSTS(37), ALGBRC(14))
ALGBRC(268) = ALGBRC(253)*CONSTS(186) - CONSTS(187)
ALGBRC(275) = TERNRY(ALGBRC(268).LT.1.00000, 1.00000, ALGBRC(268))
ALGBRC(283) = ( (ALGBRC(251) - 1.20000)*CONSTS(189)+1.20000)*CONSTS(190)
ALGBRC(293) = ((ALGBRC(253)+CONSTS(191)) - ALGBRC(283)) - ALGBRC(275)
ALGBRC(296) = ALGBRC(293)*CONSTS(192)
ALGBRC(298) = ALGBRC(296)
ALGBRC(300) = (ALGBRC(298) - 1.00000)*CONSTS(194)+1.00000
ALGBRC(301) = TERNRY(ALGBRC(300).LT.0.000100000, 0.000100000, ALGBRC(300))
ALGBRC(303) = (( ALGBRC(6)*ALGBRC(301)*CONSTS(195))/CONSTS(196))*( (ALGBRC(15) - 1.00000)*CONSTS(197)+1.00000)
ALGBRC(304) = TERNRY(ALGBRC(303).LT.CONSTS(198), CONSTS(198), ALGBRC(303))
ALGBRC(226) = ( (ALGBRC(7) - 1.00000)*CONSTS(199)+1.00000)*0.100000
ALGBRC(227) = TERNRY(ALGBRC(226).LT.0.00000, 0.00000, ALGBRC(226))
ALGBRC(310) = (ALGBRC(302) - ALGBRC(304)) - ALGBRC(227)
ALGBRC(311) = TERNRY(ALGBRC(310).LT.1.00000e-08, 1.00000e-08, ALGBRC(310))
ALGBRC(312) = ALGBRC(311)*CONSTS(179)
ALGBRC(314) = ( CONSTS(101)*ALGBRC(188) - ALGBRC(312))+ CONSTS(102)*142.000
RATES(24) = ALGBRC(314)
ALGBRC(305) = ( ALGBRC(302)*ALGBRC(85))/ALGBRC(83)
ALGBRC(306) = (ALGBRC(301) - 1.00000)*CONSTS(200)
ALGBRC(307) = (ALGBRC(233) - 1.00000)*CONSTS(201)+1.00000
ALGBRC(308) = TERNRY(ALGBRC(307).LT.0.100000, 0.100000, ALGBRC(307))
ALGBRC(228) = (ALGBRC(7) - 1.00000)*CONSTS(202)+1.00000
ALGBRC(229) = TERNRY(ALGBRC(228).LT.CONSTS(203), CONSTS(203), ALGBRC(228))
ALGBRC(309) = ( ALGBRC(85)/4.40000 ** CONSTS(204)*ALGBRC(4)*0.0800000*ALGBRC(308))/ALGBRC(229)
ALGBRC(313) = ((ALGBRC(305)+ALGBRC(309)) - STATES(43)) - ALGBRC(306)
ALGBRC(315) = TERNRY(ALGBRC(313).LT.0.00000, 0.00000, ALGBRC(313))
ALGBRC(316) = ALGBRC(315)*CONSTS(179)
ALGBRC(318) = CONSTS(104) - ALGBRC(316)
RATES(25) = ALGBRC(318)
ALGBRC(225) = STATES(42)/STATES(23)
ALGBRC(317) = ALGBRC(252) ** 2.00000*ALGBRC(225)*3.84000
ALGBRC(319) = ALGBRC(317)*CONSTS(179)
RATES(42) = CONSTS(188) - ALGBRC(319)
ALGBRC(320) = ALGBRC(317)+ 2.00000*(ALGBRC(311)+ALGBRC(315))
ALGBRC(321) = TERNRY(ALGBRC(320).GT.0.600000, 0.600000, ALGBRC(320))
ALGBRC(322) = ALGBRC(320) - 0.600000
ALGBRC(323) = TERNRY(ALGBRC(322).LT.0.00000, 0.00000, ALGBRC(322))
ALGBRC(324) = ALGBRC(321)/( 600.000*ALGBRC(15))+ALGBRC(323)/360.000
RATES(43) = ( (ALGBRC(315)/ALGBRC(324))*0.000451800 - STATES(43))*1.00000
ALGBRC(325) = ALGBRC(324)*CONSTS(179)
RATES(23) = STATES(22) - ALGBRC(325)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(32) = ((((((STATES(8)+STATES(14)) - STATES(17)) - STATES(18)) - STATES(19)) - STATES(20)) - STATES(21))/2.00000
ALGBRC(46) = STATES(19)+ ALGBRC(32)*0.128000
ALGBRC(47) = ALGBRC(46) - 0.380000
ALGBRC(48) = ALGBRC(47)/0.0100000
ALGBRC(49) = TERNRY( (ALGBRC(48) - 1.00000)*1.00000.LT.0.00000, 0.00000, (ALGBRC(48) - 1.00000)*1.00000)
ALGBRC(33) = STATES(21)+ ALGBRC(32)*0.0574000
ALGBRC(34) = ALGBRC(33) - 0.100000
ALGBRC(35) = ALGBRC(34)/0.00500000
ALGBRC(36) = TERNRY( (ALGBRC(35)+1.00000)*2.00000.LT.0.00000, 0.00000, (ALGBRC(35)+1.00000)*2.00000)
ALGBRC(51) = (ALGBRC(49)+ALGBRC(36))/3.00000
ALGBRC(53) = TERNRY(CONSTS(39).GT.0.00000, CONSTS(39), ALGBRC(51)+CONSTS(40))
ALGBRC(63) = STATES(18)+ ALGBRC(32)*0.261000
ALGBRC(64) = ALGBRC(63) - 0.495000
ALGBRC(65) = ALGBRC(64)/0.00355000
ALGBRC(66) = TERNRY(CONSTS(43).GT.1.00000e-07, CONSTS(43), ALGBRC(65) - CONSTS(44))
ALGBRC(67) = TERNRY(ALGBRC(66).LT.160.000.AND.ALGBRC(66).GE.80.0000, 0.0166670*(160.000 - ALGBRC(66)), TERNRY(ALGBRC(66).LT.80.0000, 1.33360, 0.00000)
ALGBRC(69) = (ALGBRC(67) - 1.00000)*CONSTS(47)+1.00000
ALGBRC(71) = ALGBRC(69) - CONSTS(50)
ALGBRC(82) = STATES(23) - STATES(26)
ALGBRC(83) = STATES(24)/ALGBRC(82)
ALGBRC(1) = CONSTS(9) - (CONSTS(9) - 1.00000)/( ((CONSTS(10) - 1.00000)/(CONSTS(10) - CONSTS(9)))*(STATES(1) - 1.00000)*CONSTS(11)+1.00000)
ALGBRC(2) = (ALGBRC(1) - 1.00000)*CONSTS(12)+1.00000
ALGBRC(3) = (ALGBRC(2) - 1.00000)*CONSTS(13)+1.00000
ALGBRC(4) = TERNRY(ALGBRC(3).LT.0.200000, 0.200000, ALGBRC(3))
ALGBRC(78) = (ALGBRC(4) - 1.00000)*CONSTS(103)+1.00000
ALGBRC(79) = (STATES(25) - 3000.00)/( ALGBRC(78)*9.33330)
ALGBRC(80) = STATES(25) - ALGBRC(79)
ALGBRC(81) = ALGBRC(80)/STATES(26)
ALGBRC(86) = ALGBRC(81) - ALGBRC(83)
ALGBRC(90) = ALGBRC(86)*CONSTS(105)
ALGBRC(84) = (ALGBRC(83) - CONSTS(30))/(142.000 - CONSTS(30))
ALGBRC(88) = TERNRY(ALGBRC(84).LT.0.00000, 0.00000, ALGBRC(84))
ALGBRC(68) = TERNRY(ALGBRC(66).GT.CONSTS(31), CONSTS(31), ALGBRC(66))
ALGBRC(70) = CONSTS(31) - ALGBRC(68) ** 2.00000*CONSTS(32)
ALGBRC(92) = ALGBRC(88)+ALGBRC(70)+CONSTS(33)
ALGBRC(95) = TERNRY(ALGBRC(92).LT.0.00000, 0.00000, ALGBRC(92))
ALGBRC(7) = CONSTS(23) - (CONSTS(23) - 1.00000)/( ((CONSTS(24) - 1.00000)/(CONSTS(24) - CONSTS(23)))*(STATES(3) - 1.00000)*CONSTS(25)+1.00000)
ALGBRC(8) = (ALGBRC(7) - 1.00000)*CONSTS(4)+1.00000
ALGBRC(85) = ALGBRC(79)/ALGBRC(82)
ALGBRC(87) = (ALGBRC(85) - 3.30000)/1.00000
ALGBRC(91) = ALGBRC(8)*0.909000*ALGBRC(87)
ALGBRC(94) = (ALGBRC(91) - 1.00000)*CONSTS(5)+1.00000
ALGBRC(97) = TERNRY(ALGBRC(94).LT.0.00000, 0.00000, ALGBRC(94))
ALGBRC(99) = TERNRY(CONSTS(7).GT.0.00000, CONSTS(7), ALGBRC(97)+CONSTS(6))
ALGBRC(102) = TERNRY(STATES(31).LT.0.000100000, 0.000100000, STATES(31))
ALGBRC(103) = CONSTS(120)*ALGBRC(102)
ALGBRC(105) = ALGBRC(103) - 38.0000
ALGBRC(107) = ALGBRC(105)*CONSTS(111)+1.00000
ALGBRC(108) = CONSTS(114)*ALGBRC(105)+1.00000
ALGBRC(114) = TERNRY(STATES(35).LT.0.00000, 0.00000, STATES(35))
ALGBRC(115) = ALGBRC(114)*0.486110
ALGBRC(117) = (ALGBRC(115) - 10.0000)*CONSTS(109)
ALGBRC(118) = ALGBRC(115) - CONSTS(121)
ALGBRC(120) = ALGBRC(118)*CONSTS(122)+1.00000
ALGBRC(121) = CONSTS(124)*ALGBRC(118)+1.00000
ALGBRC(122) = CONSTS(126)*ALGBRC(118)+1.00000
ALGBRC(124) = TERNRY(STATES(36).LT.0.00100000, 0.00100000, STATES(36))
ALGBRC(125) = (ALGBRC(82) - STATES(8)) - ALGBRC(124)
ALGBRC(28) = STATES(9)/STATES(8)
ALGBRC(31) = 0.280000*ALGBRC(28)+ 0.00190000*ALGBRC(28) ** 2.00000
ALGBRC(39) = STATES(20)+ ALGBRC(32)*0.155000
ALGBRC(40) = ALGBRC(39) - 0.306250
ALGBRC(41) = ALGBRC(40)/0.00480000
ALGBRC(42) = 0.0260000*ALGBRC(41)
ALGBRC(43) = TERNRY(ALGBRC(42).LT.1.00000e-05, 1.00000e-05, ALGBRC(42))
ALGBRC(44) = ALGBRC(43) ** 0.500000
ALGBRC(45) = 1.00000/ALGBRC(44)
ALGBRC(52) = ALGBRC(48)+18.0000
ALGBRC(54) = 1.00000/( ALGBRC(52)*0.0357000)
ALGBRC(123) = ( (ALGBRC(41) - ALGBRC(48))*ALGBRC(54))/(ALGBRC(54)+ALGBRC(45))+ALGBRC(48)
ALGBRC(143) = TERNRY(STATES(37).LT.0.0250000, 0.0250000, STATES(37))
ALGBRC(144) = ALGBRC(143)/ALGBRC(124)
ALGBRC(145) = ALGBRC(144)*0.400000
ALGBRC(142) = 2.00000 - 0.150000/ALGBRC(124)
ALGBRC(146) = (((ALGBRC(123) - ALGBRC(142))+ALGBRC(145)) - ALGBRC(31))*CONSTS(130)
ALGBRC(148) = (ALGBRC(142)+11.0000)*0.000300000
ALGBRC(149) = ALGBRC(146) - ALGBRC(148)
ALGBRC(150) = ALGBRC(149)
ALGBRC(151) = ALGBRC(148)*ALGBRC(144)
ALGBRC(147) = (ALGBRC(28) - ALGBRC(144))*0.000225000
ALGBRC(152) = ALGBRC(147) - ALGBRC(151)
ALGBRC(153) = ALGBRC(152)
ALGBRC(19) = CONSTS(57) - (CONSTS(57) - 1.00000)/EXP( CONSTS(59)*(STATES(7) - 1.00000))
ALGBRC(20) = TERNRY(ALGBRC(19).LT.CONSTS(58), CONSTS(58), ALGBRC(19))
ALGBRC(22) = ALGBRC(20) - 1.00000
ALGBRC(23) = ALGBRC(22)*CONSTS(62)+1.00000
ALGBRC(89) = TERNRY(STATES(27).GT.1.00000, 1.00000, STATES(27))
ALGBRC(58) = (ALGBRC(48)+4.00000)*( CONSTS(86)*(ALGBRC(20) - 1.00000)+1.00000) - 4.00000
ALGBRC(59) = TERNRY(ALGBRC(58).LE.- 2.00000, 0.0100000, TERNRY(ALGBRC(58).GT.- 2.00000.AND.ALGBRC(58).LE.1.00000, 0.0100000+( (3.60000 - 0.0100000)*(ALGBRC(58) - - 2.00000))/(1.00000 - - 2.00000), TERNRY(ALGBRC(58).GT.1.00000.AND.ALGBRC(58).LE.5.00000, 3.60000+( (9.40000 - 3.60000)*(ALGBRC(58) - 1.00000))/(5.00000 - 1.00000), TERNRY(ALGBRC(58).GT.5.00000.AND.ALGBRC(58).LE.8.00000, 9.40000+( (11.6000 - 9.40000)*(ALGBRC(58) - 5.00000))/(8.00000 - 5.00000), TERNRY(ALGBRC(58).GT.8.00000.AND.ALGBRC(58).LE.12.0000, 11.6000+( (13.5000 - 11.6000)*(ALGBRC(58) - 8.00000))/(12.0000 - 8.00000), 13.5000)
ALGBRC(156) = STATES(8)+STATES(14)
ALGBRC(157) = STATES(14)/ALGBRC(156)
ALGBRC(158) = 100.000*ALGBRC(157)
ALGBRC(159) = (STATES(30)/ALGBRC(158))/5.25000
ALGBRC(162) = ALGBRC(65)/( ALGBRC(23)*ALGBRC(159))
ALGBRC(163) = TERNRY(ALGBRC(162).LE.0.00000, 1.04000, TERNRY(ALGBRC(162).GT.0.00000.AND.ALGBRC(162).LE.60.0000, 1.04000+( (1.02500 - 1.04000)*(ALGBRC(162) - 0.00000))/(60.0000 - 0.00000), TERNRY(ALGBRC(162).GT.60.0000.AND.ALGBRC(162).LE.125.000, 1.02500+( (0.970000 - 1.02500)*(ALGBRC(162) - 60.0000))/(125.000 - 60.0000), TERNRY(ALGBRC(162).GT.125.000.AND.ALGBRC(162).LE.160.000, 0.970000+( (0.880000 - 0.970000)*(ALGBRC(162) - 125.000))/(160.000 - 125.000), TERNRY(ALGBRC(162).GT.160.000.AND.ALGBRC(162).LE.200.000, 0.880000+( (0.590000 - 0.880000)*(ALGBRC(162) - 160.000))/(200.000 - 160.000), TERNRY(ALGBRC(162).GT.200.000.AND.ALGBRC(162).LE.240.000, 0.590000+( (0.00000 - 0.590000)*(ALGBRC(162) - 200.000))/(240.000 - 200.000), 0.00000)
ALGBRC(164) = ALGBRC(163)*ALGBRC(59)*ALGBRC(23)*CONSTS(89)*ALGBRC(89)*STATES(12)
ALGBRC(74) = (ALGBRC(48) - ALGBRC(65))/3.00000
ALGBRC(165) = TERNRY(ALGBRC(74).GT.0.00000, ALGBRC(164)+ALGBRC(74), ALGBRC(164))
ALGBRC(55) = ALGBRC(54)+ALGBRC(45)
ALGBRC(56) = ALGBRC(41) - ALGBRC(48)
ALGBRC(57) = ALGBRC(56)/ALGBRC(55)
ALGBRC(167) = ALGBRC(57) - ALGBRC(165)
ALGBRC(160) = (ALGBRC(41)/ALGBRC(23))/ALGBRC(159)
ALGBRC(161) = TERNRY(ALGBRC(160).LE.0.00000, 1.06000, TERNRY(ALGBRC(160).GT.0.00000.AND.ALGBRC(160).LE.32.0000, 1.06000+( (0.970000 - 1.06000)*(ALGBRC(160) - 0.00000))/(32.0000 - 0.00000), TERNRY(ALGBRC(160).GT.32.0000.AND.ALGBRC(160).LE.38.4000, 0.970000+( (0.930000 - 0.970000)*(ALGBRC(160) - 32.0000))/(38.4000 - 32.0000), TERNRY(ALGBRC(160).GT.38.4000.AND.ALGBRC(160).LE.48.0000, 0.930000+( (0.800000 - 0.930000)*(ALGBRC(160) - 38.4000))/(48.0000 - 38.4000), TERNRY(ALGBRC(160).GT.48.0000.AND.ALGBRC(160).LE.60.8000, 0.800000+( (0.460000 - 0.800000)*(ALGBRC(160) - 48.0000))/(60.8000 - 48.0000), TERNRY(ALGBRC(160).GT.60.8000.AND.ALGBRC(160).LE.72.0000, 0.460000+( (0.00000 - 0.460000)*(ALGBRC(160) - 60.8000))/(72.0000 - 60.8000), 0.00000)
ALGBRC(166) = (1.00000 - CONSTS(87))*ALGBRC(23)*ALGBRC(161)*CONSTS(88)*ALGBRC(89)*STATES(13)+( CONSTS(87)*ALGBRC(165))/ALGBRC(59)
ALGBRC(37) = (ALGBRC(35)+8.00000)*( CONSTS(86)*(ALGBRC(20) - 1.00000)+1.00000) - 8.00000
ALGBRC(38) = TERNRY(ALGBRC(37).LE.- 8.00000, 0.00000, TERNRY(ALGBRC(37).GT.- 8.00000.AND.ALGBRC(37).LE.- 6.00000, 0.00000+( (0.750000 - 0.00000)*(ALGBRC(37) - - 8.00000))/(- 6.00000 - - 8.00000), TERNRY(ALGBRC(37).GT.- 6.00000.AND.ALGBRC(37).LE.- 2.00000, 0.750000+( (2.60000 - 0.750000)*(ALGBRC(37) - - 6.00000))/(- 2.00000 - - 6.00000), TERNRY(ALGBRC(37).GT.- 2.00000.AND.ALGBRC(37).LE.4.00000, 2.60000+( (9.80000 - 2.60000)*(ALGBRC(37) - - 2.00000))/(4.00000 - - 2.00000), TERNRY(ALGBRC(37).GT.4.00000.AND.ALGBRC(37).LE.12.0000, 9.80000+( (13.5000 - 9.80000)*(ALGBRC(37) - 4.00000))/(12.0000 - 4.00000), 13.5000)
ALGBRC(169) = ALGBRC(38)*ALGBRC(166)
ALGBRC(171) = ALGBRC(169) - ALGBRC(57)
ALGBRC(170) = TERNRY(ALGBRC(159).GT.1.00000, 114.000+ (ALGBRC(159) - 1.00000)*6667.00, TERNRY(ALGBRC(159).GT.0.936000.AND.ALGBRC(159).LE.1.00000, 74.0000+ (ALGBRC(159) - 0.936000)*625.000, TERNRY(ALGBRC(159).GT.0.800000.AND.ALGBRC(159).LE.0.936000, 46.0000+ (ALGBRC(159) - 0.800000)*205.882, ALGBRC(159)*57.5000)
ALGBRC(173) = (ALGBRC(170) - 67.0000)/30.0000
ALGBRC(175) = TERNRY(ALGBRC(173).GT.1.00000, 1.00000, TERNRY(ALGBRC(173).LT.0.600000, 0.600000, ALGBRC(173))
ALGBRC(177) = 1.00000/ALGBRC(175)
ALGBRC(180) = ( (ALGBRC(177) - 1.00000)*3.00000 - STATES(38))*0.000500000
ALGBRC(72) = TERNRY(ALGBRC(66).LT.80.0000.AND.ALGBRC(66).GE.40.0000, 0.00500000*(80.0000 - ALGBRC(66))*CONSTS(45), TERNRY(ALGBRC(66).LT.40.0000, 0.200000*CONSTS(45), 0.00000)
ALGBRC(172) = TERNRY(ALGBRC(170).LT.80.0000.AND.ALGBRC(170).GE.40.0000, CONSTS(46)*(80.0000 - ALGBRC(170)), TERNRY(ALGBRC(170).LT.40.0000, CONSTS(46)*40.0000, 0.00000)
ALGBRC(174) = ALGBRC(72)+ALGBRC(172)
ALGBRC(18) = STATES(6)
ALGBRC(73) = TERNRY(ALGBRC(66).LT.40.0000, 0.0400000*(40.0000 - ALGBRC(66))*CONSTS(51), 0.00000)
ALGBRC(50) = (15.0000/(ALGBRC(48)+ALGBRC(35)+ALGBRC(41)) - 1.00000)*CONSTS(52)+1.00000
ALGBRC(176) = ( CONSTS(207)*ALGBRC(50)*(ALGBRC(174)+ALGBRC(18)+ALGBRC(73)) - 1.00000)*CONSTS(55)+1.00000
ALGBRC(178) = TERNRY(ALGBRC(176).LT.0.00000, 0.00000, ALGBRC(176))
ALGBRC(181) = (ALGBRC(178) - STATES(7))/CONSTS(56)
ALGBRC(190) = ABS(ALGBRC(35)) ** CONSTS(154)*CONSTS(155)
ALGBRC(193) = TERNRY(ALGBRC(35).LT.0.00000, - ALGBRC(190), ALGBRC(190))
ALGBRC(185) = (CONSTS(209) - 40.0000)*ALGBRC(158)
ALGBRC(187) = TERNRY(ALGBRC(185)+CONSTS(210).LT.0.00000, 0.00000, ALGBRC(185)+CONSTS(210))
ALGBRC(189) = CONSTS(138) - ALGBRC(187)
ALGBRC(191) = TERNRY( ALGBRC(189)*CONSTS(139)*CONSTS(140)+5.00000e-06.LT.0.00000, 0.00000, ALGBRC(189)*CONSTS(139)*CONSTS(140)+5.00000e-06)
ALGBRC(179) = ALGBRC(158)/( (CONSTS(135) - ALGBRC(158))*CONSTS(136))
ALGBRC(182) = ALGBRC(179)+1.50000
ALGBRC(183) = 0.333300*ALGBRC(182)
ALGBRC(194) = STATES(14)*CONSTS(141)*ALGBRC(183)
ALGBRC(197) = (ALGBRC(191) - ALGBRC(194))+CONSTS(142)
ALGBRC(195) = (ALGBRC(7) - 1.00000)*CONSTS(151)*0.00100000
ALGBRC(184) = (ALGBRC(7) - 1.00000)*CONSTS(147)+1.00000
ALGBRC(186) = CONSTS(148) - ALGBRC(115) ** 2.00000*CONSTS(149)*ALGBRC(184)
ALGBRC(188) = TERNRY(ALGBRC(186).LT.0.800000, 0.800000, TERNRY(ALGBRC(186).GT.8.00000, 8.00000, ALGBRC(186))
ALGBRC(192) = (STATES(4) - 1.00000)*CONSTS(150)+1.00000
ALGBRC(198) = ALGBRC(192)*ALGBRC(188)*0.00100000
ALGBRC(200) = TERNRY(ALGBRC(198).LT.0.00000, 0.00000, ALGBRC(198))
ALGBRC(202) = (ALGBRC(195)+ALGBRC(200))*1.00000
ALGBRC(204) = TERNRY(ALGBRC(202).LT.0.00000, 0.00000, ALGBRC(202))
ALGBRC(21) = (CONSTS(60) - ALGBRC(20)*CONSTS(61))+CONSTS(61)
ALGBRC(196) = ALGBRC(193) - STATES(39)
ALGBRC(206) = ALGBRC(196)*CONSTS(158)
ALGBRC(60) = STATES(17)+ ALGBRC(32)*0.398600
ALGBRC(9) = (ALGBRC(7) - 1.00000)*CONSTS(26)+1.00000
ALGBRC(10) = TERNRY(ALGBRC(9).LT.CONSTS(27), CONSTS(27), ALGBRC(9))
ALGBRC(61) = (ALGBRC(10) - 1.00000)*CONSTS(90)
ALGBRC(207) = ((((ALGBRC(60) - ALGBRC(21)) - ALGBRC(61)) - STATES(16)) - STATES(15)) - ALGBRC(206)
ALGBRC(208) = TERNRY(ALGBRC(207).LT.0.000100000, 0.000100000, ALGBRC(207))
ALGBRC(209) = 3.70000+(ALGBRC(208) - 0.740000)/CONSTS(91)
ALGBRC(210) = TERNRY(ALGBRC(209).LT.0.000100000, 0.000100000, ALGBRC(209))
ALGBRC(211) = 0.740000/ALGBRC(210)/( ALGBRC(183)*3.70000) ** 0.500000
ALGBRC(62) = TERNRY(ALGBRC(35).LT.CONSTS(92), CONSTS(92), ALGBRC(35))
ALGBRC(212) = ALGBRC(210) - ALGBRC(62)
ALGBRC(213) = ALGBRC(212)/ALGBRC(211)
ALGBRC(215) = ALGBRC(213) - ALGBRC(169)
ALGBRC(214) = ALGBRC(65) - ALGBRC(210)
ALGBRC(100) = TERNRY(STATES(28).LT.CONSTS(113), CONSTS(113), STATES(28))
ALGBRC(101) = ALGBRC(100)*STATES(29)
ALGBRC(12) = CONSTS(35) - (CONSTS(35) - 1.00000)/( ((CONSTS(36) - 1.00000)/(CONSTS(36) - CONSTS(35)))*(STATES(4) - 1.00000)+1.00000)
ALGBRC(13) = TERNRY(ALGBRC(12).LT.CONSTS(36), CONSTS(36), ALGBRC(12))
ALGBRC(26) = ALGBRC(22)*CONSTS(65)+1.00000 ** CONSTS(66)
ALGBRC(199) = ALGBRC(196)*CONSTS(157)+1.00000
ALGBRC(76) = ALGBRC(65)/100.000 ** CONSTS(96)
ALGBRC(201) = (( ALGBRC(10)*ALGBRC(13)*ALGBRC(26)*ALGBRC(183)*CONSTS(2))/ALGBRC(76))/ALGBRC(199)
ALGBRC(205) = CONSTS(99)*ALGBRC(101)*ALGBRC(201)*CONSTS(1)*CONSTS(98)
ALGBRC(216) = ALGBRC(214)/ALGBRC(205)
ALGBRC(217) = STATES(30)*ALGBRC(216)
ALGBRC(25) = ALGBRC(22)*CONSTS(64)+1.00000
ALGBRC(104) = TERNRY(ALGBRC(103).GT.38.0000, 38.0000, ALGBRC(103))
ALGBRC(106) = ALGBRC(25)*CONSTS(118)*CONSTS(116)*(1.00000 - 38.0001 - ALGBRC(104) ** 3.00000/54872.0)
ALGBRC(262) = ALGBRC(230) - ALGBRC(106)
ALGBRC(11) = (ALGBRC(10) - 1.00000)*CONSTS(28)+1.00000
ALGBRC(16) = CONSTS(42) - CONSTS(42)/( 0.555556*(1.00000+STATES(5)))
ALGBRC(17) = TERNRY(ALGBRC(16).LT.- 1.00000, - 1.00000, ALGBRC(16))
ALGBRC(27) = ALGBRC(22)*CONSTS(67)+1.00000
ALGBRC(109) = TERNRY(STATES(32).LT.0.500000, 0.500000, STATES(32))
ALGBRC(110) = TERNRY(STATES(33).LT.0.500000, 0.500000, STATES(33))
ALGBRC(111) = TERNRY(STATES(34).LT.0.300000, 0.300000, STATES(34))
ALGBRC(112) = ALGBRC(109)*ALGBRC(110)*ALGBRC(111)
ALGBRC(113) = (ALGBRC(112) - 1.00000)*CONSTS(128)+1.00000
ALGBRC(203) = CONSTS(97)*ALGBRC(113)*ALGBRC(201)*CONSTS(1)*CONSTS(98)
ALGBRC(218) = TERNRY(CONSTS(160).GT.0.00000.AND.CONSTS(161).LE.0.00000, CONSTS(160), TERNRY(CONSTS(161).GT.0.00000, STATES(40), ALGBRC(65) - CONSTS(159))
ALGBRC(219) = (ALGBRC(26) - 1.00000)*CONSTS(168)+1.00000
ALGBRC(220) = TERNRY(ALGBRC(219).LT.0.800000, 0.800000, ALGBRC(219))
ALGBRC(221) = (ALGBRC(7) - 1.00000)*CONSTS(169)+1.00000
ALGBRC(222) = TERNRY(ALGBRC(221).LT.CONSTS(170), CONSTS(170), ALGBRC(221))
ALGBRC(223) = (ALGBRC(220) - 1.00000)*CONSTS(174)+1.00000
ALGBRC(224) = (ALGBRC(7) - 1.00000)*CONSTS(175)+1.00000
ALGBRC(263) = ALGBRC(233)
ALGBRC(269) = TERNRY(ALGBRC(263).GT.1.00000, 1.00000/(1.00000+ (ALGBRC(263) - 1.00000)*72.0000), 10.0000 - 9.00000/(1.00000+ (1.00000 - ALGBRC(263))*8.00000))
ALGBRC(276) = (ALGBRC(269) - 1.00000)*CONSTS(17)
ALGBRC(75) = ALGBRC(65) - ALGBRC(35)
ALGBRC(77) = ALGBRC(75)*CONSTS(100)
ALGBRC(265) = CONSTS(179)*ALGBRC(251)
ALGBRC(271) = ALGBRC(216)+ALGBRC(235)+ALGBRC(265)
ALGBRC(280) = ALGBRC(271)+ALGBRC(77)
ALGBRC(126) = ALGBRC(125) - STATES(10)
ALGBRC(134) = (ALGBRC(126) - CONSTS(84))/CONSTS(84) ** 2.00000*1.00000
ALGBRC(131) = (CONSTS(80)/ALGBRC(125))/5.00000 ** CONSTS(81)
ALGBRC(135) = ALGBRC(131)*CONSTS(82)+ALGBRC(134)
ALGBRC(128) = STATES(11)/ALGBRC(125)
ALGBRC(130) = 0.280000*ALGBRC(128)+ 0.00190000*ALGBRC(128) ** 2.00000
ALGBRC(132) = ALGBRC(131)*2.00000
ALGBRC(133) = ALGBRC(132)*ALGBRC(130)*CONSTS(83)
ALGBRC(264) = TERNRY(ALGBRC(236) - CONSTS(75).LT.0.00000, 0.00000, ALGBRC(236) - CONSTS(75))
ALGBRC(270) = ALGBRC(264)*CONSTS(76) ** CONSTS(77)
ALGBRC(278) = (((ALGBRC(236) - ALGBRC(31)) - ALGBRC(135))+ALGBRC(133))*CONSTS(69)+ALGBRC(270)
ALGBRC(136) = ALGBRC(135) - ALGBRC(132)
ALGBRC(137) = (ALGBRC(136)+CONSTS(85)) - ALGBRC(134)
ALGBRC(139) = TERNRY(ALGBRC(137).GT.7.00000, 7.00000, ALGBRC(137))
ALGBRC(140) = TERNRY(ALGBRC(139).LT.0.00000, 0.00000, ALGBRC(139)*0.0200000)
ALGBRC(285) = ((ALGBRC(140) - ALGBRC(278)) - ALGBRC(150))+CONSTS(70)
ALGBRC(29) = TERNRY(ALGBRC(28) - CONSTS(71).LT.0.00000, 0.00000, ALGBRC(28) - CONSTS(71))
ALGBRC(30) = CONSTS(72) - ALGBRC(29) ** CONSTS(73)*CONSTS(74)
ALGBRC(279) = ALGBRC(270)*ALGBRC(28)+ (ALGBRC(28) - ALGBRC(128))*0.00104000
ALGBRC(141) = ALGBRC(128)*ALGBRC(140)
ALGBRC(286) = (((ALGBRC(30)+ALGBRC(141)) - ALGBRC(279)) - ALGBRC(153))+ CONSTS(70)*72.0000
ALGBRC(287) = ALGBRC(279) - ALGBRC(141)
ALGBRC(288) = ALGBRC(280) - ALGBRC(213)
ALGBRC(289) = ALGBRC(165) - ALGBRC(280)
ALGBRC(266) = STATES(30)*ALGBRC(235)
ALGBRC(116) = TERNRY(ALGBRC(115).GT.35.0000, 35.0000, ALGBRC(115))
ALGBRC(119) = ALGBRC(25)*CONSTS(129)*(1.00000 - 35.0001 - ALGBRC(116) ** 3.00000/42875.0)
ALGBRC(281) = ALGBRC(272) - ALGBRC(119)
ALGBRC(291) = TERNRY(ALGBRC(114).LT.6.00000.AND.ALGBRC(281).LT.0.00000, ALGBRC(281)*0.100000, ALGBRC(281))
ALGBRC(277) = (ALGBRC(269)+STATES(2))*CONSTS(19)
ALGBRC(284) = TERNRY(ALGBRC(277).LT.1.00000e-05, 1.00000e-05, ALGBRC(277))
ALGBRC(294) = TERNRY(CONSTS(20).GT.0.00000, CONSTS(20), ALGBRC(284)+CONSTS(21))
ALGBRC(282) = ALGBRC(272)+ALGBRC(230)
ALGBRC(155) = STATES(38)+1.00000
ALGBRC(292) = ALGBRC(282)*CONSTS(131)*0.0266670*ALGBRC(177)*ALGBRC(155)
ALGBRC(295) = CONSTS(132) - (ALGBRC(282)/ALGBRC(292))/0.761000
ALGBRC(154) = CONSTS(133)/(CONSTS(134)+ALGBRC(124))
ALGBRC(297) = (ALGBRC(295) - ALGBRC(170))*ALGBRC(154)
ALGBRC(299) = (ALGBRC(297) - ALGBRC(282))/( ALGBRC(169)*1.00000)
ALGBRC(302) = ALGBRC(233)*ALGBRC(83)*0.00616190
ALGBRC(5) = (ALGBRC(2) - 1.00000)*CONSTS(14)+1.00000
ALGBRC(6) = TERNRY(ALGBRC(5).LT.CONSTS(16), CONSTS(16), TERNRY(ALGBRC(5).GT.CONSTS(15), CONSTS(15), ALGBRC(5))
ALGBRC(14) = CONSTS(38) - (CONSTS(38) - 1.00000)/( ((CONSTS(37) - 1.00000)/(CONSTS(37) - CONSTS(38)))*(STATES(4) - 1.00000)+1.00000)
ALGBRC(15) = TERNRY(ALGBRC(14).LT.CONSTS(37), CONSTS(37), ALGBRC(14))
ALGBRC(268) = ALGBRC(253)*CONSTS(186) - CONSTS(187)
ALGBRC(275) = TERNRY(ALGBRC(268).LT.1.00000, 1.00000, ALGBRC(268))
ALGBRC(283) = ( (ALGBRC(251) - 1.20000)*CONSTS(189)+1.20000)*CONSTS(190)
ALGBRC(293) = ((ALGBRC(253)+CONSTS(191)) - ALGBRC(283)) - ALGBRC(275)
ALGBRC(296) = ALGBRC(293)*CONSTS(192)
ALGBRC(298) = ALGBRC(296)
ALGBRC(300) = (ALGBRC(298) - 1.00000)*CONSTS(194)+1.00000
ALGBRC(301) = TERNRY(ALGBRC(300).LT.0.000100000, 0.000100000, ALGBRC(300))
ALGBRC(303) = (( ALGBRC(6)*ALGBRC(301)*CONSTS(195))/CONSTS(196))*( (ALGBRC(15) - 1.00000)*CONSTS(197)+1.00000)
ALGBRC(304) = TERNRY(ALGBRC(303).LT.CONSTS(198), CONSTS(198), ALGBRC(303))
ALGBRC(226) = ( (ALGBRC(7) - 1.00000)*CONSTS(199)+1.00000)*0.100000
ALGBRC(227) = TERNRY(ALGBRC(226).LT.0.00000, 0.00000, ALGBRC(226))
ALGBRC(310) = (ALGBRC(302) - ALGBRC(304)) - ALGBRC(227)
ALGBRC(311) = TERNRY(ALGBRC(310).LT.1.00000e-08, 1.00000e-08, ALGBRC(310))
ALGBRC(312) = ALGBRC(311)*CONSTS(179)
ALGBRC(314) = ( CONSTS(101)*ALGBRC(188) - ALGBRC(312))+ CONSTS(102)*142.000
ALGBRC(305) = ( ALGBRC(302)*ALGBRC(85))/ALGBRC(83)
ALGBRC(306) = (ALGBRC(301) - 1.00000)*CONSTS(200)
ALGBRC(307) = (ALGBRC(233) - 1.00000)*CONSTS(201)+1.00000
ALGBRC(308) = TERNRY(ALGBRC(307).LT.0.100000, 0.100000, ALGBRC(307))
ALGBRC(228) = (ALGBRC(7) - 1.00000)*CONSTS(202)+1.00000
ALGBRC(229) = TERNRY(ALGBRC(228).LT.CONSTS(203), CONSTS(203), ALGBRC(228))
ALGBRC(309) = ( ALGBRC(85)/4.40000 ** CONSTS(204)*ALGBRC(4)*0.0800000*ALGBRC(308))/ALGBRC(229)
ALGBRC(313) = ((ALGBRC(305)+ALGBRC(309)) - STATES(43)) - ALGBRC(306)
ALGBRC(315) = TERNRY(ALGBRC(313).LT.0.00000, 0.00000, ALGBRC(313))
ALGBRC(316) = ALGBRC(315)*CONSTS(179)
ALGBRC(318) = CONSTS(104) - ALGBRC(316)
ALGBRC(225) = STATES(42)/STATES(23)
ALGBRC(317) = ALGBRC(252) ** 2.00000*ALGBRC(225)*3.84000
ALGBRC(319) = ALGBRC(317)*CONSTS(179)
ALGBRC(320) = ALGBRC(317)+ 2.00000*(ALGBRC(311)+ALGBRC(315))
ALGBRC(321) = TERNRY(ALGBRC(320).GT.0.600000, 0.600000, ALGBRC(320))
ALGBRC(322) = ALGBRC(320) - 0.600000
ALGBRC(323) = TERNRY(ALGBRC(322).LT.0.00000, 0.00000, ALGBRC(322))
ALGBRC(324) = ALGBRC(321)/( 600.000*ALGBRC(15))+ALGBRC(323)/360.000
ALGBRC(325) = ALGBRC(324)*CONSTS(179)
ALGBRC(24) = ALGBRC(22)*CONSTS(63)+1.00000
ALGBRC(93) = 72.0000*ALGBRC(24)
ALGBRC(96) = (ALGBRC(89) - 1.00000)*0.500000+1.00000
ALGBRC(98) = (ALGBRC(93)+CONSTS(208))*ALGBRC(96)
ALGBRC(127) = TERNRY(ALGBRC(125).LE.0.00000, 0.00000, TERNRY(ALGBRC(125).GT.0.00000.AND.ALGBRC(125).LE.12.0000, 0.00000+( (11.4000 - 0.00000)*(ALGBRC(125) - 0.00000))/(12.0000 - 0.00000), TERNRY(ALGBRC(125).GT.12.0000.AND.ALGBRC(125).LE.15.0000, 11.4000+( (14.0000 - 11.4000)*(ALGBRC(125) - 12.0000))/(15.0000 - 12.0000), TERNRY(ALGBRC(125).GT.15.0000.AND.ALGBRC(125).LE.18.0000, 14.0000+( (16.0000 - 14.0000)*(ALGBRC(125) - 15.0000))/(18.0000 - 15.0000), TERNRY(ALGBRC(125).GT.18.0000.AND.ALGBRC(125).LE.21.0000, 16.0000+( (17.3000 - 16.0000)*(ALGBRC(125) - 18.0000))/(21.0000 - 18.0000), TERNRY(ALGBRC(125).GT.21.0000.AND.ALGBRC(125).LE.24.0000, 17.3000+( (18.0000 - 17.3000)*(ALGBRC(125) - 21.0000))/(24.0000 - 21.0000), 18.0000)
ALGBRC(129) = ALGBRC(125) - ALGBRC(127)
ALGBRC(138) = ALGBRC(134) - ALGBRC(136)
ALGBRC(168) = ALGBRC(165)/ALGBRC(98)
ALGBRC(267) = ALGBRC(252)*CONSTS(179)
ALGBRC(290) = ALGBRC(75)/ALGBRC(280)
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0[1] = abs(ALGBRC(232) - (ALGBRC(217) - ALGBRC(230))/( ALGBRC(158)*5.25000*ALGBRC(216)))
minfunc_0[2] = abs(ALGBRC(231) - 57.1400*ALGBRC(232)*CONSTS(116) ** CONSTS(117))
minfunc_0[3] = abs(ALGBRC(230) - (ALGBRC(231) - ALGBRC(103))*CONSTS(119)*ALGBRC(216))
RETURN
END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_1[1] = abs(ALGBRC(236) - ( ALGBRC(234)*1.70000*ALGBRC(235)+ALGBRC(210)))
minfunc_1[2] = abs(ALGBRC(237) - ( (ALGBRC(236) - 17.0000)*CONSTS(93)+17.0000)*CONSTS(94))
minfunc_1[3] = abs(ALGBRC(238) - CONSTS(95)/ALGBRC(237))
minfunc_1[4] = abs(ALGBRC(234) - ALGBRC(27)*ALGBRC(238)*ALGBRC(183)*ALGBRC(11))
minfunc_1[5] = abs(ALGBRC(239) - ALGBRC(234)*1.79000)
minfunc_1[6] = abs(ALGBRC(240) - (ALGBRC(203)+ALGBRC(239)))
minfunc_1[7] = abs(ALGBRC(235) - ALGBRC(214)/ALGBRC(240))
minfunc_1[8] = abs(ALGBRC(243) - ( (ALGBRC(233) - 1.00000)*CONSTS(164)+1.00000))
minfunc_1[9] = abs(ALGBRC(242) - (TERNRY(ALGBRC(243).LT.CONSTS(165), CONSTS(165), TERNRY(ALGBRC(243).GT.CONSTS(166), CONSTS(166), ALGBRC(243))))
minfunc_1[10] = abs(ALGBRC(241) - (ALGBRC(242) - STATES(41)))
minfunc_1[11] = abs(ALGBRC(244) - CONSTS(171)*CONSTS(3)*ALGBRC(220)*ALGBRC(241)*ALGBRC(222)*40.0000*CONSTS(1))
minfunc_1[12] = abs(ALGBRC(246) - ((ALGBRC(244) - ALGBRC(17)*CONSTS(172))+CONSTS(172)))
minfunc_1[13] = abs(ALGBRC(245) - (TERNRY(ALGBRC(246).LT.CONSTS(173), CONSTS(173), ALGBRC(246))))
minfunc_1[14] = abs(ALGBRC(247) - ( (ALGBRC(241) - 1.00000)*CONSTS(176)+1.00000))
minfunc_1[15] = abs(ALGBRC(249) - 43.3330*CONSTS(177)*ALGBRC(224)*ALGBRC(247)*CONSTS(1)*ALGBRC(223))
minfunc_1[16] = abs(ALGBRC(248) - (TERNRY(ALGBRC(249).LT.CONSTS(178), CONSTS(178), ALGBRC(249))))
minfunc_1[17] = abs(ALGBRC(250) - (ALGBRC(245)+ALGBRC(248)))
minfunc_1[18] = abs(ALGBRC(251) - ALGBRC(218)/ALGBRC(250))
minfunc_1[19] = abs(ALGBRC(255) - ( ALGBRC(251)*(1.00000 - ALGBRC(157)))/( ALGBRC(251)*(1.00000 - ALGBRC(157)) - ALGBRC(252)))
minfunc_1[20] = abs(ALGBRC(254) - (TERNRY(ALGBRC(255).LT.1.00000, 1.00000, ALGBRC(255))))
minfunc_1[21] = abs(ALGBRC(253) - (TERNRY(CONSTS(181).GT.0.00000, ALGBRC(254) ** 1.35000*ALGBRC(31)*0.980000, ALGBRC(31)+4.00000)))
minfunc_1[22] = abs(ALGBRC(257) - ALGBRC(245)*ALGBRC(251))
minfunc_1[23] = abs(ALGBRC(256) - (ALGBRC(218) - ALGBRC(257)))
minfunc_1[24] = abs(ALGBRC(258) - ((ALGBRC(256) - ALGBRC(253)) - CONSTS(182)))
minfunc_1[25] = abs(ALGBRC(259) - ALGBRC(258)*CONSTS(183))
minfunc_1[26] = abs(ALGBRC(252) - (TERNRY(ALGBRC(259).LT.CONSTS(184), CONSTS(184), ALGBRC(259))))
minfunc_1[27] = abs(ALGBRC(260) - ALGBRC(252)*8.00000)
minfunc_1[28] = abs(ALGBRC(261) - ( (ALGBRC(260) - 1.00000)*CONSTS(185)+1.00000))
minfunc_1[29] = abs(ALGBRC(233) - (TERNRY(ALGBRC(261).LT.0.00000, 0.00000, ALGBRC(261))))
RETURN
END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_2[1] = abs(ALGBRC(274) - (ALGBRC(266) - ALGBRC(272))/( ALGBRC(158)*5.25000*ALGBRC(235)))
minfunc_2[2] = abs(ALGBRC(273) - ALGBRC(274)*57.1400)
minfunc_2[3] = abs(ALGBRC(272) - (ALGBRC(273) - ALGBRC(115))*12.8570*ALGBRC(235))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END