Location: Guyton, Autonomics, 2008 @ 50ce90269e36 / autonomics.cellml

Author:
Tommy Yu <tommy.yu@auckland.ac.nz>
Date:
2010-08-25 18:18:49+12:00
Desc:
xml:base fix
Permanent Source URI:
http://models.cellml.org/workspace/guyton_autonomics_2008/rawfile/50ce90269e36dbc4ef48734f78b672330f5d19aa/autonomics.cellml

<?xml version="1.0" encoding="utf-8"?>
<model 
	name="autonomics_model"
	cmeta:id="autonomics_model"
	
	xmlns="http://www.cellml.org/cellml/1.1#"
	xmlns:cellml="http://www.cellml.org/cellml/1.1#"
	xmlns:cmeta="http://www.cellml.org/metadata/1.0#"
	xmlns:xlink="http://www.w3.org/1999/xlink">

	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#autonomics_model">
			<rdf:value>
			Autonomic Control Of The Circulation
			Autonomic control of the circulation primarily operates through the sympathetic
			system, though to a slight extent through parasympathetic signals to the heart.  
			These have been lumped together, and there are basically three separate feedback 
			mechanisms in this computational block.  These are:  (1)  feedback from the 
			baroreceptor control system; (2) feedback from the peripheral chemoreceptors in 
			the carotid and aortic bodies,; and (3) feedback control of the circulatory system 
			caused by central nervous system ischemia, that is, ischemia of the vasomotor center 
			of the brainstem.  Several other inputs that affect the autonomic nervous system are 
			also included.  These are: activation of the autonomic nervous system during exercise; 
			baroreceptor feedback effects from pulmonary artery pressure (PPA), left atrial 
			pressure (PLA), and an effect of low blood PO2 (PO2ART).

			Note:  Not shown in the diagram is also a variable (STA) that is normally zero.  
			When it is set to any level above zero, the value of the general autonomic multiplier (AU) 
			becomes fixed to the value of STA. 
			WHERE DO I PUT THIS NOTE?????
			</rdf:value>
		</rdf:Description>
	</rdf:RDF>

<!-- =============================================	Import the required units from the Units file	==================================================== -->
	<import xlink:href="units.cellml">
		<units name="minute" units_ref="minute"/>
		<units name="per_minute" units_ref="per_minute"/>
		<units name="beats_per_minute" units_ref="beats_per_minute"/>
		<units name="beats_per_minute_per_mmHg" units_ref="beats_per_minute_per_mmHg"/>
		<units name="minute_per_L" units_ref="minute_per_L"/>
		<units name="mmHg" units_ref="mmHg"/>
		<units name="per_mmHg" units_ref="per_mmHg"/>
		<units name="mmHg_per_mL" units_ref="mmHg_per_mL"/>
		<units name="mmHg_L" units_ref="mmHg_L"/>
		<units name="per_mmHg2" units_ref="per_mmHg2"/>
		<units name="mmHg3" units_ref="mmHg3"/>
		<units name="monovalent_mEq" units_ref="monovalent_mEq"/>
		<units name="monovalent_mEq_per_minute" units_ref="monovalent_mEq_per_minute"/>
		<units name="monovalent_mEq_per_litre" units_ref="monovalent_mEq_per_litre"/>
		<units name="mOsm" units_ref="mOsm"/>
		<units name="mOsm_per_litre" units_ref="mOsm_per_litre"/>
		<units name="mOsm_per_minute" units_ref="mOsm_per_minute"/>
		<units name="monovalent_mEq_per_litre_per_minute" units_ref="monovalent_mEq_per_litre_per_minute"/>
		<units name="litre2_per_monovalent_mEq_per_minute" units_ref="litre2_per_monovalent_mEq_per_minute"/>
		<units name="L_per_minute" units_ref="L_per_minute"/>
		<units name="per_mmHg_per_minute" units_ref="per_mmHg_per_minute"/>
		<units name="mL" units_ref="mL"/>
		<units name="gram_per_L" units_ref="gram_per_L"/>
		<units name="L_mmHg_per_gram" units_ref="L_mmHg_per_gram"/>
		<units name="L2_mmHg_per_gram2" units_ref="L2_mmHg_per_gram2"/>
		<units name="mmHg_minute_per_L" units_ref="mmHg_minute_per_L"/>
		<units name="mmHg_L_per_minute" units_ref="mmHg_L_per_minute"/>
		<units name="gram_per_minute" units_ref="gram_per_minute"/>
		<units name="mL_per_L" units_ref="mL_per_L"/>
		<units name="mL_per_L_per_mmHg" units_ref="mL_per_L_per_mmHg"/>
		<units name="mL_per_L_per_minute" units_ref="mL_per_L_per_minute"/>
		<units name="mL_per_minute_per_mmHg" units_ref="mL_per_minute_per_mmHg"/>
		<units name="L_mL_per_minute_per_mmHg" units_ref="L_mL_per_minute_per_mmHg"/>
		<units name="L_per_mL" units_ref="L_per_mL"/>
		<units name="L_per_mmHg" units_ref="L_per_mmHg"/>
		<units name="mL_per_minute" units_ref="mL_per_minute"/>
		<units name="L_per_minute_per_mmHg" units_ref="L_per_minute_per_mmHg"/>
		<units name="L_per_minute_per_mmHg2" units_ref="L_per_minute_per_mmHg2"/>
	</import>

<!-- =====================================	Import all Parameters and State Variables from the Parameter file	============================================== -->
	<import xlink:href="parameters.cellml">
		<component component_ref="parameter_values" name="parameter_values"/>
		<component component_ref="state_variables" name="state_variables"/>
	</import>


<!-- ========================================	AUTONOMICS TOP-LEVEL COMPONENT	============================================= -->
	<component 	name="autonomics" 
				cmeta:id="autonomics">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#autonomics">
				<rdf:value>
				Encapsulation grouping component containing all the components in the Autonomics Model. The inputs and 
				outputs of the Autonomics Model must be passed by this component.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>

		<variable name="time"	units="minute"	private_interface="out" public_interface="in"/>

<!-- Inputs from components in other models -->
		<variable name="PA"		units="mmHg" private_interface="out" public_interface="in"/>
		<variable name="PO2ART"	units="mmHg" private_interface="out" public_interface="in"/>
		<variable name="PLA"	units="mmHg" private_interface="out" public_interface="in"/>
		<variable name="PRA"	units="mmHg" private_interface="out" public_interface="in"/>
		<variable name="PPA"	units="mmHg" private_interface="out" public_interface="in"/>

<!-- Outputs to components in other models -->
		<variable name="PA1" 	units="mmHg" private_interface="in" public_interface="out"/>
		<variable name="VVR" 	units="litre" private_interface="in" public_interface="out"/>
		<variable name="AUH" 	units="dimensionless" private_interface="in" public_interface="out"/>
		<variable name="AUR" 	units="dimensionless" private_interface="in" public_interface="out"/>
		<variable name="AUM" 	units="dimensionless" private_interface="in" public_interface="out"/>
		<variable name="AOM" 	units="dimensionless" private_interface="in" public_interface="out"/>
		<variable name="AVE" 	units="dimensionless" private_interface="in" public_interface="out"/>
		<variable name="AU" 	units="dimensionless" private_interface="in" public_interface="out"/>
	</component>

<!-- OUTPUT CONNECTIONS -->
	<connection>
		<map_components component_1="autonomics" component_2="pressure_driving_autonomic_receptors"/>
		<map_variables variable_1="PA1" variable_2="PA1"/>
		<map_variables variable_1="PA" variable_2="PA"/>
	</connection>
	<connection>
		<map_components component_1="autonomics" component_2="autonomic_drive_on_target_organs_and_tissues"/>
		<map_variables variable_1="VVR" variable_2="VVR"/>
		<map_variables variable_1="AUH" variable_2="AUH"/>
		<map_variables variable_1="AUR" variable_2="AUR"/>
		<map_variables variable_1="AUM" variable_2="AUM"/>
		<map_variables variable_1="AOM" variable_2="AOM"/>
		<map_variables variable_1="AVE" variable_2="AVE"/>
	</connection>
	<connection>
		<map_components component_1="autonomics" component_2="actual_autonomic_stimulation"/>
		<map_variables variable_1="AU" variable_2="AU"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>


<!-- ========================================	PRESSURE DRIVING AUTONOMIC RECEPTORS, PA1		============================================= -->
	<component 	name="pressure_driving_autonomic_receptors" 
				cmeta:id="pressure_driving_autonomic_receptors">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#pressure_driving_autonomic_receptors">
				<rdf:value>
				AU1:  
				Calculation of the effective systemic arterial pressure (PA1) by subtracting 
				any pressure drop (EXE) between the output point of the heart where the pressure 
				is equal to PA and the remainder of the aortic tree where the pressure will be 
				equal to PA1.  This block allows simulation of aortic constriction at its root.

				AU2:  
				A block that will allow one to set the effective systemic arterial pressure (PA1) 
				to any constant value desired by increasing the value CRRFLX to any value above zero.  
				As long as it remains at zero, there is no effect.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU1_and_AU2">
				<rdf:value>
				AU1:  
				Calculation of the effective systemic arterial pressure (PA1) by subtracting 
				any pressure drop (EXE) between the output point of the heart where the pressure 
				is equal to PA and the remainder of the aortic tree where the pressure will be 
				equal to PA1.  This block allows simulation of aortic constriction at its root.

				AU2:  
				A block that will allow one to set the effective systemic arterial pressure (PA1) 
				to any constant value desired by increasing the value CRRFLX to any value above zero.  
				As long as it remains at zero, there is no effect.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PA"	units="mmHg" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable 	cmeta:id="pressure_driving_autonomic_receptors_PA1"
					name="PA1" 	units="mmHg" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="CRRFLX"	units="mmHg" private_interface="none" public_interface="in"/>
		<variable name="EXE"	units="mmHg" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU1_and_AU2">
				<eq/>
				<ci>PA1</ci>
				<piecewise>
					<piece>
						<ci>CRRFLX</ci>
						<apply>
							<gt/>
							<ci>CRRFLX</ci>
							<cn cellml:units="mmHg">0.0000001</cn>
						</apply>
					</piece>
					<otherwise>
						<apply>
							<minus/>
							<ci>PA</ci>
							<ci>EXE</ci>
						</apply>
					</otherwise>
				</piecewise>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
<!-- see Autonomics connections -->

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="pressure_driving_autonomic_receptors" component_2="parameter_values"/>
		<map_variables variable_1="CRRFLX" variable_2="CRRFLX"/>
		<map_variables variable_1="EXE" variable_2="EXE"/>
	</connection>

<!-- ========================================	CHEMORECEPTORS - EFFECT	OF PA	============================================= -->
	<component 	name="chemoreceptors_effect_of_PA" 
				cmeta:id="chemoreceptors_effect_of_PA">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#chemoreceptors_effect_of_PA">
				<rdf:value>
				AU4:  
				Calculation of the nervous output from the chemoreceptors (AUC) at the 
				different systemic arterial pressure levels (PA1).  
				
				AU19:  
				Sensitivity control for increasing or decreasing the degree of response of 
				AUC to chemoreceptor nervous output.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU4_and_AU19">
				<rdf:value>
				AU4:  
				Calculation of the nervous output from the chemoreceptors (AUC) at the 
				different systemic arterial pressure levels (PA1).  
				
				AU19:  
				Sensitivity control for increasing or decreasing the degree of response of 
				AUC to chemoreceptor nervous output.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PA1"	units="mmHg" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AUC" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="AUC1"	units="dimensionless" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU4_and_AU19">
				<eq/>
				<ci>AUC</ci>
				<piecewise>
					<piece>
						<apply>
							<times/>
							<cn cellml:units="per_mmHg">0.005</cn>
							<apply>
								<minus/>
								<cn cellml:units="mmHg">80</cn>
								<ci>PA1</ci>
							</apply>
							<ci>AUC1</ci>
						</apply>
						<apply>
							<and/>
							<apply>
								<lt/>
								<ci>PA1</ci>
								<cn cellml:units="mmHg">80</cn>
							</apply>
							<apply>
								<geq/>
								<ci>PA1</ci>
								<cn cellml:units="mmHg">40</cn>
							</apply>
						</apply>
					</piece>
					<piece>
						<apply>
							<times/>
							<cn cellml:units="dimensionless">0.2</cn>
							<ci>AUC1</ci>
						</apply>
						<apply>
							<lt/>
							<ci>PA1</ci>
							<cn cellml:units="mmHg">40</cn>
						</apply>
					</piece>
					<otherwise>
						<cn cellml:units="dimensionless">0</cn>
					</otherwise>
				</piecewise>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="chemoreceptors_effect_of_PA" component_2="pressure_driving_autonomic_receptors"/>
		<map_variables variable_1="PA1" variable_2="PA1"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="chemoreceptors_effect_of_PA" component_2="parameter_values"/>
		<map_variables variable_1="AUC1" variable_2="AUC1"/>
	</connection>

<!-- ========================================	CHEMORECEPTORS - EFFECT	OF ARTERIAL PO2	============================================= -->
	<component 	name="chemoreceptors_effect_of_art_PO2" 
				cmeta:id="chemoreceptors_effect_of_art_PO2">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#chemoreceptors_effect_of_art_PO2">
				<rdf:value>
				AU20:  
				Calculation of a nervous factor (AUC2) for effecting autonomic control of 
				the circulation depending on the peripheral level of oxygen in the blood (PO2ART).  
				O2CHMO is a sensitivity controller.

				AU21:  
				Addition of the chemoreflex output (AUC) caused by activation of the chemoreceptors 
				by low arterial pressure plus the chemoreceptor output (AUC2) caused by reduced 
				arterial oxygen saturation (PO2ART).  The output of Block 21 is equal to AUC3.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU20">
				<rdf:value>
				AU20:  
				Calculation of a nervous factor (AUC2) for effecting autonomic control of 
				the circulation depending on the peripheral level of oxygen in the blood (PO2ART).  
				O2CHMO is a sensitivity controller.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU21">
				<rdf:value>
				AU21:  
				Addition of the chemoreflex output (AUC) caused by activation of the chemoreceptors 
				by low arterial pressure plus the chemoreceptor output (AUC2) caused by reduced 
				arterial oxygen saturation (PO2ART).  The output of Block 21 is equal to AUC3.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PO2ART"	units="mmHg" private_interface="none" public_interface="in"/>
		<variable name="AUC"	units="dimensionless" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AUC3" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="O2CHMO"	units="per_mmHg" private_interface="none" public_interface="in"/>

<!-- Internal variables -->
		<variable name="AUC2"	units="dimensionless" private_interface="none" public_interface="none"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU20">
				<eq/>
				<ci>AUC2</ci>
				<piecewise>
					<piece>
						<apply>
							<times/>
							<ci>O2CHMO</ci>
							<apply>
								<minus/>
								<cn cellml:units="mmHg">80</cn>
								<ci>PO2ART</ci>
							</apply>
						</apply>
						<apply>
							<and/>
							<apply>
								<lt/>
								<ci>PO2ART</ci>
								<cn cellml:units="mmHg">80</cn>
							</apply>
							<apply>
								<geq/>
								<ci>PO2ART</ci>
								<cn cellml:units="mmHg">40</cn>
							</apply>
						</apply>
					</piece>
					<piece>
						<apply>
							<times/>
							<ci>O2CHMO</ci>
							<cn cellml:units="mmHg">40</cn>
						</apply>
						<apply>
							<lt/>
							<ci>PO2ART</ci>
							<cn cellml:units="mmHg">40</cn>
						</apply>
					</piece>
					<otherwise>
						<cn cellml:units="dimensionless">0</cn>
					</otherwise>
				</piecewise>
			</apply>
			<apply id="AU21">
				<eq/>
				<ci>AUC3</ci>
				<apply>
					<plus/>
					<ci>AUC</ci>
					<ci>AUC2</ci>
				</apply>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="chemoreceptors_effect_of_art_PO2" component_2="autonomics"/>
		<map_variables variable_1="PO2ART" variable_2="PO2ART"/>
	</connection>
	<connection>
		<map_components component_1="chemoreceptors_effect_of_art_PO2" component_2="chemoreceptors_effect_of_PA"/>
		<map_variables variable_1="AUC" variable_2="AUC"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="chemoreceptors_effect_of_art_PO2" component_2="parameter_values"/>
		<map_variables variable_1="O2CHMO" variable_2="O2CHMO"/>
	</connection>

<!-- ========================================	BARO REFLEX / ARTERIAL BARORECEPTORS	============================================= -->
<!-- ========================================	FIX!!!! - MISSING AU17 and AU14-16!!!	============================================= -->
	<component 	name="arterial_baroreceptor_reflex" 
				cmeta:id="arterial_baroreceptor_reflex">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#arterial_baroreceptor_reflex">
				<rdf:value>
				AU3:  
				Calculation of the nervous output from the baroreceptors (AUB) at the different 
				systemic arterial pressure levels (PA1).

				AU6, AU7, and AU8:  
				Sensitivity control of the effect of baroreceptor reflex output (AUB) on the 
				autonomic nervous system.  The output of this sensitivity controller is AU6A, 
				and the degree of sensitivity control is equal to AUX.

				AU9, AU10, and AU11:  
				Time delay in the buildup of sympathetic output (AU6) that results from changes 
				in baroreceptor reflex nervous output (AU6A).  The time constant of this delay 
				circuit is equal to BAROTC. 

				AU14, AU15, and AU16:  
				Calculation of adaptation of the baroreceptor feedback mechanism (AU4) over a 
				period of hours.  The time constant of this adaptation is equal to AUK.
				THIS IS COMMENTED OUT BECAUSE IT IS NOT IN THE MODSIM CODE

				AU17:  
				This block sets AU6 equal to AU6A irrespective of the time constant for buildup 
				of the nervous effect of the baroreceptor reflex when the iteration interval for 
				solution of the equations is greater than .16666.  This prevents some instability 
				when the equations are being calculated for long-term instead of short-term changes.

				AU18:  
				Damping of baroreceptor autonomic feedback output (AU6) when long-term solutions are 
				being calculated, to prevent oscillation in the circuit.  The output after the damping 
				is AU6C.  MDMP sets the degree of damping.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU3">
				<rdf:value>
				AU3:  
				Calculation of the nervous output from the baroreceptors (AUB) at the different 
				systemic arterial pressure levels (PA1).
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU6_AU7_and_part_AU8">
				<rdf:value>
				AU6, AU7, and AU8:  
				Sensitivity control of the effect of baroreceptor reflex output (AUB) on the 
				autonomic nervous system.  The output of this sensitivity controller is AU6A, 
				and the degree of sensitivity control is equal to AUX.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#rest_AU8">
				<rdf:value>
				AU6, AU7, and AU8:  
				Sensitivity control of the effect of baroreceptor reflex output (AUB) on the 
				autonomic nervous system.  The output of this sensitivity controller is AU6A, 
				and the degree of sensitivity control is equal to AUX.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU9_to_AU11">
				<rdf:value>
				AU9, AU10, and AU11:  
				Time delay in the buildup of sympathetic output (AU6) that results from changes 
				in baroreceptor reflex nervous output (AU6A).  The time constant of this delay 
				circuit is equal to BAROTC. 
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU14_to_AU16">
				<rdf:value>
				AU14, AU15, and AU16:  
				Calculation of adaptation of the baroreceptor feedback mechanism (AU4) over a 
				period of hours.  The time constant of this adaptation is equal to AUK.
				THIS IS COMMENTED OUT BECAUSE IT IS NOT IN THE MODSIM CODE
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU18">
				<rdf:value>
				AU18:  
				Damping of baroreceptor autonomic feedback output (AU6) when long-term solutions are 
				being calculated, to prevent oscillation in the circuit.  The output after the damping 
				is AU6C.  MDMP sets the degree of damping.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PA1"	units="mmHg" private_interface="none" public_interface="in"/>
		<variable name="AU6_init"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AU4_init"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="time"	units="minute" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AU6C" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="AUX"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUK"	units="per_minute" private_interface="none" public_interface="in"/>
		<variable name="BAROTC"	units="minute" private_interface="none" public_interface="in"/>

<!-- Internal variables -->
		<variable name="AUB"	units="dimensionless" private_interface="none" public_interface="none"/>
		<variable name="A1B"	units="dimensionless" private_interface="none" public_interface="none"/>
		<variable name="AU6A"	units="dimensionless" private_interface="none" public_interface="none"/>
<!--		<variable name="AU8"	units="per_minute" private_interface="none" public_interface="none"/>
		<variable name="AU2"	units="dimensionless" private_interface="none" public_interface="none"/> -->
		<variable name="AU4"	initial_value="AU4_init"	units="dimensionless" private_interface="none" public_interface="none"/>
		<variable name="AU6"	initial_value="AU6_init"	units="dimensionless" private_interface="none" public_interface="none"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU3">
				<eq/>
				<ci>AUB</ci>
				<piecewise>
					<piece>
						<apply>
							<times/>
							<cn cellml:units="per_mmHg">0.016667</cn>
							<apply>
								<minus/>
								<cn cellml:units="mmHg">160</cn>
								<ci>PA1</ci>
							</apply>
						</apply>
						<apply>
							<and/>
							<apply>
								<lt/>
								<ci>PA1</ci>
								<cn cellml:units="mmHg">160</cn>
							</apply>
							<apply>
								<geq/>
								<ci>PA1</ci>
								<cn cellml:units="mmHg">80</cn>
							</apply>
						</apply>
					</piece>
					<piece>
						<cn cellml:units="dimensionless">1.3336</cn>
						<apply>
							<lt/>
							<ci>PA1</ci>
							<cn cellml:units="mmHg">80</cn>
						</apply>
					</piece>
					<otherwise>
						<cn cellml:units="dimensionless">0</cn>
					</otherwise>
				</piecewise>
			</apply>
			<apply id="AU6_AU7_and_part_AU8">
				<eq/>
				<ci>A1B</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<apply>
							<minus/>
							<ci>AUB</ci>
							<cn cellml:units="dimensionless">1</cn>
						</apply>
						<ci>AUX</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="rest_AU8">
				<eq/>
				<ci>AU6A</ci>
				<apply>
					<minus/>
					<ci>A1B</ci>
					<ci>AU4</ci>
				</apply>
			</apply>
			<apply id="AU9_to_AU11">
				<eq/>
				<apply>
					<diff/>
					<bvar>
						<ci>time</ci>
					</bvar>
					<ci>AU6</ci>
				</apply>
				<apply>
					<divide/>
					<apply>
						<minus/>
						<ci>AU6A</ci>
						<ci>AU6</ci>
					</apply>
					<ci>BAROTC</ci>
				</apply>
			</apply>
<!--			<apply id="AU14">
                <eq/>
                <ci>AU2</ci>
                <apply>
                    <minus/>
                    <ci>AU6</ci>
                    <cn cellml:units="dimensionless">1</cn>
                </apply>
            </apply>
			<apply id="AU15">
                <eq/>
                <ci>AU8</ci>
                <apply>
                    <times/>
                    <ci>AUK</ci>
                    <ci>AU2</ci>
                </apply>
            </apply>
			<apply id="AU16">
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>AU4</ci>
                </apply>
                <ci>AU8</ci>
            </apply> -->
			<apply id="AU18">
				<eq/>
				<ci>AU6C</ci>
				<ci>AU6</ci>
			</apply>
		</math>
	</component>

<!-- STATE VARIABLES INITIAL VALUE -->
	<connection>
		<map_components component_1="arterial_baroreceptor_reflex" component_2="state_variables"/>
		<map_variables variable_1="AU4_init" variable_2="AU4_init"/> 
		<map_variables variable_1="AU6_init" variable_2="AU6_init"/> 
	</connection>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="arterial_baroreceptor_reflex" component_2="pressure_driving_autonomic_receptors"/>
		<map_variables variable_1="PA1" variable_2="PA1"/>
	</connection>
	<connection>
		<map_components component_1="arterial_baroreceptor_reflex" component_2="autonomics"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="arterial_baroreceptor_reflex" component_2="parameter_values"/>
		<map_variables variable_1="AUX" variable_2="AUX"/>
		<map_variables variable_1="AUK" variable_2="AUK"/>
		<map_variables variable_1="BAROTC" variable_2="BAROTC"/>
	</connection>

<!-- ========================================	CNS ISCHEMIC REFLEX	============================================= -->
	<component 	name="CNS_ischemic_reflex" 
				cmeta:id="CNS_ischemic_reflex">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#CNS_ischemic_reflex">
				<rdf:value>
				AU5:  
				Calculation of the nervous output (AUN) caused by activation of the central 
				nervous system ischemic reflex, resulting from reduced systemic arterial 
				pressure (PA1).

				AU22:  
				Control of the sensitivity of the CNS ischemic reflex output by the sensitivity 
				controller AUN1.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU5_and_AU22">
				<rdf:value>
				AU5:  
				Calculation of the nervous output (AUN) caused by activation of the central 
				nervous system ischemic reflex, resulting from reduced systemic arterial 
				pressure (PA1).

				AU22:  
				Control of the sensitivity of the CNS ischemic reflex output by the sensitivity 
				controller AUN1.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PA1"	units="mmHg" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AUN" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="AUN1"	units="dimensionless" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU5_and_AU22">
				<eq/>
				<ci>AUN</ci>
				<piecewise>
					<piece>
						<apply>
							<times/>
							<cn cellml:units="per_mmHg">0.04</cn>
							<apply>
								<minus/>
								<cn cellml:units="mmHg">40</cn>
								<ci>PA1</ci>
							</apply>
							<ci>AUN1</ci>
						</apply>
						<apply>
							<lt/>
							<ci>PA1</ci>
							<cn cellml:units="mmHg">40</cn>
						</apply>
					</piece>
					<otherwise>
						<cn cellml:units="dimensionless">0</cn>
					</otherwise>
				</piecewise>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="CNS_ischemic_reflex" component_2="pressure_driving_autonomic_receptors"/>
		<map_variables variable_1="PA1" variable_2="PA1"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="CNS_ischemic_reflex" component_2="parameter_values"/>
		<map_variables variable_1="AUN1" variable_2="AUN1"/>
	</connection>

<!-- ==============================	STRETCH RECEPTORS - RESPONSE TO PRESSURE IN PULMONARY VASCULATURE	==================================== -->
	<component 	name="autonomic_response_to_vasculature_pressure" 
				cmeta:id="autonomic_response_to_vasculature_pressure">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#autonomic_response_to_vasculature_pressure">
				<rdf:value>
				AU24, AU25, AU26, AU27, and AU28:  
				Calculation of an additional factor (AULP) that affects the total autonomic 
				response, caused by stretch receptors in the pulmonary vasculature.  These 
				are in response to left atrial pressure (PLA), right atrial pressure (PRA), 
				and pulmonary arterial pressure (PPA).  The sensitivity controller for these 
				effects is AULPM in Block 27.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU24_to_AU28">
				<rdf:value>
				AU24, AU25, AU26, AU27, and AU28:  
				Calculation of an additional factor (AULP) that affects the total autonomic 
				response, caused by stretch receptors in the pulmonary vasculature.  These 
				are in response to left atrial pressure (PLA), right atrial pressure (PRA), 
				and pulmonary arterial pressure (PPA).  The sensitivity controller for these 
				effects is AULPM in Block 27.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="PLA"	units="mmHg" private_interface="none" public_interface="in"/>
		<variable name="PRA"	units="mmHg" private_interface="none" public_interface="in"/>
		<variable name="PPA"	units="mmHg" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AULP" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="AULPM"	units="dimensionless" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU24_to_AU28">
				<eq/>
				<ci>AULP</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<apply>
							<minus/>
							<apply>
								<divide/>
								<cn cellml:units="mmHg">15</cn>
								<apply>
									<plus/>
									<ci>PLA</ci>
									<ci>PRA</ci>
									<ci>PPA</ci>
								</apply>
							</apply>
							<cn cellml:units="dimensionless">1</cn>
						</apply>
						<ci>AULPM</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="autonomic_response_to_vasculature_pressure" component_2="autonomics"/>
		<map_variables variable_1="PLA" variable_2="PLA"/>
		<map_variables variable_1="PRA" variable_2="PRA"/>
		<map_variables variable_1="PPA" variable_2="PPA"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="autonomic_response_to_vasculature_pressure" component_2="parameter_values"/>
		<map_variables variable_1="AULPM" variable_2="AULPM"/>
	</connection>

<!-- ========================================	EFFECT OF EXERCISE ON AUTONOMIC OUTPUT	============================================= -->
	<component 	name="autonomic_response_to_exercise" 
				cmeta:id="autonomic_response_to_exercise">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#autonomic_response_to_exercise">
				<rdf:value>
				AU29, AU30, AU31, and AU32:  
				Effect of the exercise nervous signal (EXC) on autonomic output.  The 
				exponent (EXCXP) in Block 29 provides a curve fitting effect of exercise on the 
				autonomic output, and the factor (EXCML) is a sensitivity multiplier effect.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU29_to_AU32">
				<rdf:value>
				AU29, AU30, AU31, and AU32:  
				Effect of the exercise nervous signal (EXC) on autonomic output.  The 
				exponent (EXCXP) in Block 29 provides a curve fitting effect of exercise on the 
				autonomic output, and the factor (EXCML) is a sensitivity multiplier effect.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Outputs to other components -->
		<variable name="AUEX" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="EXC"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="EXCXP"	units="dimensionless" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU29_to_AU32">
				<eq/>
				<ci>AUEX</ci>
				<apply>
					<power/>
					<ci>EXC</ci>
					<ci>EXCXP</ci>
				</apply>
			</apply>
		</math>
	</component>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="autonomic_response_to_exercise" component_2="parameter_values"/>
		<map_variables variable_1="EXC" variable_2="EXC"/>
		<map_variables variable_1="EXCXP" variable_2="EXCXP"/>
	</connection>

<!-- ========================================	TOTAL AUTONOMIC STIMULATION	============================================= -->
	<component 	name="total_autonomic_stimulation" 
				cmeta:id="total_autonomic_stimulation">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#total_autonomic_stimulation">
				<rdf:value>
				AU23:  
				Summation of the different nervous output controls of autonomic stimulation, AUC3 
				from the chemoreceptor component, AU6C from the arterial baroreceptor reflex component, 
				and the output of Block 22 from the CNS ischemic reflex component.

				AU33:  
				Summation of all of the different factors affecting autonomic stimulation 
				of the circulation, giving a total output of AUTTL. 

				AU34:  
				Limitation of the lower level of autonomic stimulation of the circulatory system 
				so that this cannot fall below the level of zero.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU23_and_AU33">
				<rdf:value>
				AU23:  
				Summation of the different nervous output controls of autonomic stimulation, AUC3 
				from the chemoreceptor component, AU6C from the arterial baroreceptor reflex component, 
				and the output of Block 22 from the CNS ischemic reflex component.

				AU33:  
				Summation of all of the different factors affecting autonomic stimulation 
				of the circulation, giving a total output of AUTTL. 
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU34">
				<rdf:value>
				AU34:  
				Limitation of the lower level of autonomic stimulation of the circulatory system 
				so that this cannot fall below the level of zero.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="AUC3"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AU6C"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUN"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AULP"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUEX"	units="dimensionless" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable name="AUTTL" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Internal variables -->
		<variable name="AUTTL1"	units="dimensionless" private_interface="none" public_interface="none"/>

<!-- Parameters from parameter_file -->
		<variable name="EXCML"	units="dimensionless" private_interface="none" public_interface="in"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU23_and_AU33">
				<eq/>
				<ci>AUTTL1</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<apply>
							<minus/>
							<apply>
								<times/>
								<ci>AUEX</ci>
								<ci>AULP</ci>
								<apply>
									<plus/>
									<ci>AUC3</ci>
									<ci>AU6C</ci>
									<ci>AUN</ci>
								</apply>
							</apply>
							<cn cellml:units="dimensionless">1</cn>
						</apply>
						<ci>EXCML</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="AU34">
				<eq/>
				<ci>AUTTL</ci>
				<piecewise>
					<piece>
						<cn cellml:units="dimensionless">0</cn>
						<apply>
							<lt/>
							<ci>AUTTL1</ci>
							<cn cellml:units="dimensionless">0</cn>
						</apply>
					</piece>
					<otherwise>
						<ci>AUTTL1</ci>
					</otherwise>
				</piecewise>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="chemoreceptors_effect_of_art_PO2"/>
		<map_variables variable_1="AUC3" variable_2="AUC3"/>
	</connection>
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="arterial_baroreceptor_reflex"/>
		<map_variables variable_1="AU6C" variable_2="AU6C"/>
	</connection>
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="CNS_ischemic_reflex"/>
		<map_variables variable_1="AUN" variable_2="AUN"/>
	</connection>
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="autonomic_response_to_vasculature_pressure"/>
		<map_variables variable_1="AULP" variable_2="AULP"/>
	</connection>
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="autonomic_response_to_exercise"/>
		<map_variables variable_1="AUEX" variable_2="AUEX"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="total_autonomic_stimulation" component_2="parameter_values"/>
		<map_variables variable_1="EXCML" variable_2="EXCML"/>
	</connection>

<!-- ========================================	ACTUAL AUTONOMIC STIMULATION		============================================= -->
	<component 	name="actual_autonomic_stimulation" 
				cmeta:id="actual_autonomic_stimulation">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#actual_autonomic_stimulation">
				<rdf:value>
				AU35, AU36, and AU37:  
				This is a time-delay circuit to delay the peripheral changes that occur in the 
				circulatory system for a fraction of a minute after changes in the nervous component 
				take place.  This results from the need to build up autonomic transmitter substance 
				and for the different organs to respond.  The output after this delay circuit is AU1.  
				The time constant of the delay is AUDMP.

				AU38:  
				This is a curve fitting step to fit the output strength of functional reaction to 
				sympathetic stimulation (AU) to the input level of nervous stimulation (AU1).  
				The maximum level of AU is set by the equation in this block to equal AUMAX.  
				AUSLPC determines the slope of the relationship.

				AU39:  
				This sets the minimum level of AU (the output functional reaction) equal to a 
				minimum value of AUMIN.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU35">
				<rdf:value>
				AU35, AU36, and AU37:  
				This is a time-delay circuit to delay the peripheral changes that occur in the 
				circulatory system for a fraction of a minute after changes in the nervous component 
				take place.  This results from the need to build up autonomic transmitter substance 
				and for the different organs to respond.  The output after this delay circuit is AU1.  
				The time constant of the delay is AUDMP.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU36_and_AU37">
				<rdf:value>
				AU35, AU36, and AU37:  
				This is a time-delay circuit to delay the peripheral changes that occur in the 
				circulatory system for a fraction of a minute after changes in the nervous component 
				take place.  This results from the need to build up autonomic transmitter substance 
				and for the different organs to respond.  The output after this delay circuit is AU1.  
				The time constant of the delay is AUDMP.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU38">
				<rdf:value>
				AU38:  
				This is a curve fitting step to fit the output strength of functional reaction to 
				sympathetic stimulation (AU) to the input level of nervous stimulation (AU1).  
				The maximum level of AU is set by the equation in this block to equal AUMAX.  
				AUSLPC determines the slope of the relationship.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU39">
				<rdf:value>
				AU39:  
				This sets the minimum level of AU (the output functional reaction) equal to a 
				minimum value of AUMIN.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="AUTTL"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AU1_init"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="time"	units="minute" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable 	cmeta:id="actual_autonomic_stimulation_AU"
					name="AU" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="AUDMP"	units="minute" private_interface="none" public_interface="in"/>
		<variable name="AUMAX"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUMIN"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUSLP"	units="dimensionless" private_interface="none" public_interface="in"/>

<!-- Internal variables -->
		<variable name="DAU"	units="per_minute" private_interface="none" public_interface="none"/>
		<variable name="AU1"	initial_value="AU1_init"	units="dimensionless" private_interface="none" public_interface="none"/>
		<variable name="AUT"	units="dimensionless" private_interface="none" public_interface="none"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU35">
				<eq/>
				<ci>DAU</ci>
				<apply>
					<divide/>
					<apply>
						<minus/>
						<ci>AUTTL</ci>
						<ci>AU1</ci>
					</apply>
					<ci>AUDMP</ci>
				</apply>
			</apply>
			<apply id="AU36_and_AU37">
				<eq/>
				<apply>
					<diff/>
					<bvar>
						<ci>time</ci>
					</bvar>
					<ci>AU1</ci>
				</apply>
				<ci>DAU</ci>
			</apply>
			<apply id="AU38">
				<eq/>
				<ci>AUT</ci>
				<apply>
					<minus/>
					<ci>AUMAX</ci>
					<apply>
						<divide/>
						<apply>
							<minus/>
							<ci>AUMAX</ci>
							<cn cellml:units="dimensionless">1</cn>
						</apply>
						<apply>
							<exp/>
							<apply>
								<times/>
								<ci>AUSLP</ci>
								<apply>
									<minus/>
									<ci>AU1</ci>
									<cn cellml:units="dimensionless">1</cn>
								</apply>
							</apply>
						</apply>
					</apply>
				</apply>
			</apply>
			<apply id="AU39">
				<eq/>
				<ci>AU</ci>
				<piecewise>
					<piece>
						<ci>AUMIN</ci>
						<apply>
							<lt/>
							<ci>AUT</ci>
							<ci>AUMIN</ci>
						</apply>
					</piece>
					<otherwise>
						<ci>AUT</ci>
					</otherwise>
				</piecewise>
			</apply>
		</math>
	</component>

<!-- STATE VARIABLES INITIAL VALUE -->
	<connection>
		<map_components component_1="actual_autonomic_stimulation" component_2="state_variables"/>
		<map_variables variable_1="AU1_init" variable_2="AU1_init"/>
	</connection>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="actual_autonomic_stimulation" component_2="total_autonomic_stimulation"/>
		<map_variables variable_1="AUTTL" variable_2="AUTTL"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="actual_autonomic_stimulation" component_2="parameter_values"/>
		<map_variables variable_1="AUDMP" variable_2="AUDMP"/>
		<map_variables variable_1="AUMAX" variable_2="AUMAX"/>
		<map_variables variable_1="AUMIN" variable_2="AUMIN"/>
		<map_variables variable_1="AUSLP" variable_2="AUSLP"/>
	</connection>

<!-- ==============================	AUTONOMIC DRIVE ON THE VARIOUS TARGET ORGANS AND TISSUES	================================= -->
<!-- ==============================	FIX!!! - AU47 AND AU48 NOT IN ORIGINAL CODE BUT ARE ON DIAGRAM!!!!	================================= -->
	<component 	name="autonomic_drive_on_target_organs_and_tissues" 
				cmeta:id="autonomic_drive_on_target_organs_and_tissues">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#autonomic_drive_on_target_organs_and_tissues">
				<rdf:value>
				AU40 and AU41:  
				Calculation of the effect on venous vascular resistance (VVR) of different 
				levels of autonomic functional reaction (AU).  The sensitivity control is AUL, 
				and VV9 determines the range.

				AU42:  
				A step to reduce the output effect of normal autonomic reaction (AU) equal to 
				zero (AUO) so that differences from control levels can be activated in 
				Blocks 43, 45, 47, 50, and 53.

				AU43 and AU44:  
				Sensitivity control for the autonomic effect on the heart (AUH).  The sensitivity 
				is controlled by AUV.

				AU45 and AU46:  
				Sensitivity control of the autonomic effect on heart rate (AUR).  The sensitivity 
				is controlled by AUS.

				AU47 and AU48:
				Calculation of the autonomic effect on muscle metabolism (AOM).  The sensitivity 
				control for this variable is O2A.

				AU50, AU51, and AU52:  
				Calculation of an autonomic multiplier effect that is used at multiple points in 
				the circulatory system (AUM).  The values AUM1 and AUM2 are curve fitting controls.

				AU53 and AU54:  
				Sensitivity control for adjusting the autonomic multiplier effect on the 
				venous tree (AVE).  The variable (AUY) controls the sensitivity.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU40_and_AU41">
				<rdf:value>
				AU40 and AU41:  
				Calculation of the effect on venous vascular resistance (VVR) of different 
				levels of autonomic functional reaction (AU).  The sensitivity control is AUL, 
				and VV9 determines the range.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU42">
				<rdf:value>
				AU42:  
				A step to reduce the output effect of normal autonomic reaction (AU) equal to 
				zero (AUO) so that differences from control levels can be activated in 
				Blocks 43, 45, 47, 50, and 53.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU43_and_AU44">
				<rdf:value>
				AU43 and AU44:  
				Sensitivity control for the autonomic effect on the heart (AUH).  The sensitivity 
				is controlled by AUV.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU45_and_AU46">
				<rdf:value>
				AU45 and AU46:  
				Sensitivity control of the autonomic effect on heart rate (AUR).  The sensitivity 
				is controlled by AUS.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU47_and_AU48">
				<rdf:value>
				AU47 and AU48:
				Calculation of the autonomic effect on muscle metabolism (AOM).  The sensitivity 
				control for this variable is O2A.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU50_to_AU52">
				<rdf:value>
				AU50, AU51, and AU52:  
				Calculation of an autonomic multiplier effect that is used at multiple points in 
				the circulatory system (AUM).  The values AUM1 and AUM2 are curve fitting controls.
				</rdf:value>
			</rdf:Description>
			<rdf:Description rdf:about="#AU53_and_AU54">
				<rdf:value>
				AU53 and AU54:  
				Sensitivity control for adjusting the autonomic multiplier effect on the 
				venous tree (AVE).  The variable (AUY) controls the sensitivity.
				</rdf:value>
			</rdf:Description>
		</rdf:RDF>
		
<!-- Inputs from other components -->
		<variable name="AU"	units="dimensionless" private_interface="none" public_interface="in"/>

<!-- Outputs to other components -->
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_VVR"
					name="VVR" 	units="litre" private_interface="none" public_interface="out"/>
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_AUH"
					name="AUH" 	units="dimensionless" private_interface="none" public_interface="out"/>
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_AUR"
					name="AUR" 	units="dimensionless" private_interface="none" public_interface="out"/>
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_AOM"
					name="AOM" 	units="dimensionless" private_interface="none" public_interface="out"/>
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_AUM"
					name="AUM" 	units="dimensionless" private_interface="none" public_interface="out"/>
		<variable 	cmeta:id="autonomic_drive_on_target_organs_and_tissues_AVE"
					name="AVE" 	units="dimensionless" private_interface="none" public_interface="out"/>

<!-- Parameters from parameter_file -->
		<variable name="VV9"	units="litre" private_interface="none" public_interface="in"/>
		<variable name="AUL"	units="litre" private_interface="none" public_interface="in"/>
		<variable name="AUV"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUS"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="O2A"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUM1"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUM2"	units="dimensionless" private_interface="none" public_interface="in"/>
		<variable name="AUY"	units="dimensionless" private_interface="none" public_interface="in"/>

<!-- Internal variables -->
		<variable name="AUO"	units="dimensionless" private_interface="none" public_interface="none"/>

		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply id="AU40_and_AU41">
				<eq/>
				<ci>VVR</ci>
				<apply>
					<plus/>
					<apply>
						<minus/>
						<ci>VV9</ci>
						<apply>
							<times/>
							<ci>AU</ci>
							<ci>AUL</ci>
						</apply>
					</apply>
					<ci>AUL</ci>
				</apply>
			</apply>
			<apply id="AU42">
				<eq/>
				<ci>AUO</ci>
				<apply>
					<minus/>
					<ci>AU</ci>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="AU43_and_AU44">
				<eq/>
				<ci>AUH</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<ci>AUO</ci>
						<ci>AUV</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="AU45_and_AU46">
				<eq/>
				<ci>AUR</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<ci>AUO</ci>
						<ci>AUS</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="AU47_and_AU48">
				<eq/>
				<ci>AOM</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<ci>AUO</ci>
						<ci>O2A</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
			<apply id="AU50_to_AU52">
				<eq/>
				<ci>AUM</ci>
				<apply>
					<power/>
					<apply>
						<plus/>
						<apply>
							<times/>
							<ci>AUO</ci>
							<ci>AUM1</ci>
						</apply>
						<cn cellml:units="dimensionless">1</cn>
					</apply>
					<ci>AUM2</ci>
				</apply>
			</apply>
			<apply id="AU53_and_AU54">
				<eq/>
				<ci>AVE</ci>
				<apply>
					<plus/>
					<apply>
						<times/>
						<ci>AUO</ci>
						<ci>AUY</ci>
					</apply>
					<cn cellml:units="dimensionless">1</cn>
				</apply>
			</apply>
		</math>
	</component>

<!-- INPUT CONNECTIONS -->
	<connection>
		<map_components component_1="autonomic_drive_on_target_organs_and_tissues" component_2="actual_autonomic_stimulation"/>
		<map_variables variable_1="AU" variable_2="AU"/>
	</connection>

<!-- PARAMETER CONNECTIONS -->
	<connection>
		<map_components component_1="autonomic_drive_on_target_organs_and_tissues" component_2="parameter_values"/>
		<map_variables variable_1="VV9" variable_2="VV9"/>
		<map_variables variable_1="AUL" variable_2="AUL"/>
		<map_variables variable_1="AUV" variable_2="AUV"/>
		<map_variables variable_1="AUS" variable_2="AUS"/>
		<map_variables variable_1="O2A" variable_2="O2A"/>
		<map_variables variable_1="AUM1" variable_2="AUM1"/>
		<map_variables variable_1="AUM2" variable_2="AUM2"/>
		<map_variables variable_1="AUY" variable_2="AUY"/>
	</connection>


	
<!-- ============================================================		GROUPING		=============================================================== -->
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="autonomics">
			<component_ref component="pressure_driving_autonomic_receptors"/>
			<component_ref component="chemoreceptors_effect_of_PA"/>
			<component_ref component="chemoreceptors_effect_of_art_PO2"/>
			<component_ref component="arterial_baroreceptor_reflex"/>
			<component_ref component="CNS_ischemic_reflex"/>
			<component_ref component="autonomic_response_to_vasculature_pressure"/>
			<component_ref component="autonomic_response_to_exercise"/>
			<component_ref component="total_autonomic_stimulation"/>
			<component_ref component="actual_autonomic_stimulation"/>
			<component_ref component="autonomic_drive_on_target_organs_and_tissues"/>
      </component_ref>
   </group>
      
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="autonomics">
			<component_ref component="parameter_values"/>
			<component_ref component="state_variables"/>
			<component_ref component="pressure_driving_autonomic_receptors"/>
			<component_ref component="chemoreceptors_effect_of_PA"/>
			<component_ref component="chemoreceptors_effect_of_art_PO2"/>
			<component_ref component="arterial_baroreceptor_reflex"/>
			<component_ref component="CNS_ischemic_reflex"/>
			<component_ref component="autonomic_response_to_vasculature_pressure"/>
			<component_ref component="autonomic_response_to_exercise"/>
			<component_ref component="total_autonomic_stimulation"/>
			<component_ref component="actual_autonomic_stimulation"/>
			<component_ref component="autonomic_drive_on_target_organs_and_tissues"/>
      </component_ref>
   </group>
   
</model>