Location: Liu, Hu, Peng, Liu, 1999 @ 9fdc78c888be / liu_hu_peng_liu_1999.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2010-07-01 09:31:19+12:00
Desc:
Added xul and session file
Permanent Source URI:
https://staging.physiomeproject.org/workspace/liu_hu_peng_liu_1999/rawfile/9fdc78c888be578144819f1db224ac36544369a6/liu_hu_peng_liu_1999.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 30/03/2010 at 12:44:59 at p.m. using:

COR (0.9.31.1371)
Copyright 2002-2010 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.cellml.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" cmeta:id="liu_model_1999" name="liu_model_1999">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Mathematical modeling and qualitative analysis of insulin therapies</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
				This CellML model runs in PCenv, COR and OpenCell to recreate the published results. The initial conditions were not explicitly defined in the paper, but were extrapolated from Table 1 and Figure 2. The model was then run until it reached a steady state, and new intial values were taken from the steady state. 
				</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT:  We propose a dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis. This model takes into account both the binding of hormone with proteins in the plasma and tissues and mutual interactions of the hormones in the system. The deductions from 
this model are in good agreement with experiment results.
</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="liu_1999.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the hormone flow and interactions between the different components of the hypothalamic-pituitary-adrenal (HPA) axis which are described by the mathematical model.  CRH represents corticotropin releasing hormone, ACTH represents corticotropin, and CBG is corticosteroid binding globulin.  Red arrows represent hormone disposal, while black arrows represent flow between the different glands and plasma. + represents a stimulatory effect while - represents an inhibitory effect.</caption>
				</informalfigure>
				
				<para>
The original paper reference is cited below:
</para>
				<para>
A dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis, Yi-Wei Liu, Zhi-Hong Hu, Jian-Hua Peng, and Bing-Zheng Liu, 1999, Mathematical and Computer Modelling, 29 (Issue 4), 103-110.  (No Pubmed ID)  
</para>
			</sect1>
		</article>
	</documentation>
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
				<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Family>Lloyd</vCard:Family>
					<vCard:Given>Catherine</vCard:Given>
					<vCard:Other>May</vCard:Other>
				</vCard:N>
				<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
				<dcterms:W3CDTF>2007-07-13</dcterms:W3CDTF>
			</dcterms:created>
			<dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">
        The University of Auckland, Auckland Bioengineering Institute
      </dc:publisher>
		</rdf:Description>
		<rdf:Description rdf:about="#liu_model_1999">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
        Liu et al's 1999 mathematical model for the pulsatile secretion of the hypothalamic-pituitary-adrenal axis.
      </dc:title>
			<cmeta:comment rdf:parseType="Resource">
				<rdf:value>
          This is a CellML description of Liu et al's 1999 mathematical model for the pulsatile secretion of the hypothalamic-pituitary-adrenal axis.
        </rdf:value>
				<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
					<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
				</dc:creator>
			</cmeta:comment>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
				<dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
					<bqs:subject_type>keyword</bqs:subject_type>
					<rdf:value>
						<rdf:Bag>
							<rdf:li>cortisol</rdf:li>
							<rdf:li>endocrine</rdf:li>
							<rdf:li>adrenal cortex</rdf:li>
							<rdf:li>HPA axis</rdf:li>
							<rdf:li>hypothalamus</rdf:li>
							<rdf:li>hpa axis</rdf:li>
							<rdf:li>anterior pituitary</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
				<bqs:JournalArticle rdf:parseType="Resource">
					<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
						<rdf:Seq>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Liu</vCard:Family>
										<vCard:Given>Yi-Wei</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Hu</vCard:Family>
										<vCard:Given>Zhi-Hong</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Peng</vCard:Family>
										<vCard:Given>Jian-Hua</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Liu</vCard:Family>
										<vCard:Given>Bing-Zheng</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
            A dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis
          </dc:title>
					<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
						<dcterms:W3CDTF>1999-02</dcterms:W3CDTF>
					</dcterms:issued>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
              Mathematical and Computer Modelling
            </dc:title>
					</bqs:Journal>
					<bqs:volume>29</bqs:volume>
					<bqs:issue>4</bqs:issue>
					<bqs:first_page>103</bqs:first_page>
					<bqs:last_page>110</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">CRH</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">corticotropin-releasing hormone</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">ACTH</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">corticotropin</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x3">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x3</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">free cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x4">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x4</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">CBG-bound cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x5">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x5</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">albumin-bound cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	
	<units name="minute">
      <unit multiplier="60" units="second"/>
   </units>
   <units name="microg_l">
      <unit prefix="micro" units="gram"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="per_microg_l">
      <unit exponent="-1" units="microg_l"/>
   </units>
   <units name="per_microg_l2">
      <unit exponent="-2" units="microg_l"/>
   </units>
   <units name="flux">
      <unit units="microg_l"/>
      <unit exponent="-1" units="minute"/>
   </units>
   <units name="first_order_rate_constant">
      <unit exponent="-1" units="minute"/>
   </units>
   <units name="second_order_rate_constant">
      <unit exponent="-1" units="microg_l"/>
      <unit exponent="-1" units="minute"/>
   </units>
   
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="minute"/>
   </component>
   
   <component name="x1">
      <variable cmeta:id="x1_x1" initial_value="0.01067" name="x1" public_interface="out" units="microg_l"/>
      <variable name="lambda_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x3" public_interface="in" units="microg_l"/>
      <variable name="a1" public_interface="in" units="flux"/>
      <variable name="a2" public_interface="in" units="flux"/>
      <variable name="a3" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a4" public_interface="in" units="second_order_rate_constant"/>
      <variable name="a5" public_interface="in" units="per_microg_l"/>
      <variable name="a6" public_interface="in" units="per_microg_l2"/>
      <variable name="a7" public_interface="in" units="per_microg_l"/>
      <variable name="a8" public_interface="in" units="per_microg_l2"/>
      <variable name="time" public_interface="in" units="minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>a1</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>a2</ci>
                        <apply>
                           <times/>
                           <ci>a3</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a4</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>a5</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a6</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a7</ci>
                           <ci>x3</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a8</ci>
                           <apply>
                              <power/>
                              <ci>x3</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_1</ci>
                  <ci>x1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="x2">
      <variable cmeta:id="x2_x2" initial_value="0.04665" name="x2" public_interface="out" units="microg_l"/>
      <variable name="lambda_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x1" public_interface="in" units="microg_l"/>
      <variable name="x3" public_interface="in" units="microg_l"/>
      <variable name="a9" public_interface="in" units="flux"/>
      <variable name="a10" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a11" public_interface="in" units="second_order_rate_constant"/>
      <variable name="a12" public_interface="in" units="per_microg_l"/>
      <variable name="a13" public_interface="in" units="per_microg_l2"/>
      <variable name="a14" public_interface="in" units="per_microg_l"/>
      <variable name="a15" public_interface="in" units="per_microg_l2"/>
      <variable name="time" public_interface="in" units="minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <plus/>
                     <ci>a9</ci>
                     <apply>
                        <times/>
                        <ci>a10</ci>
                        <ci>x1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a11</ci>
                        <apply>
                           <power/>
                           <ci>x1</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <ci>a12</ci>
                        <ci>x1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a13</ci>
                        <apply>
                           <power/>
                           <ci>x1</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a14</ci>
                        <ci>x3</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a15</ci>
                        <apply>
                           <power/>
                           <ci>x3</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_2</ci>
                  <ci>x2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="x3">
      <variable cmeta:id="x3_x3" initial_value="6.51" name="x3" public_interface="out" units="microg_l"/>
      <variable name="lambda_3_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x1" public_interface="in" units="microg_l"/>
      <variable name="x2" public_interface="in" units="microg_l"/>
      <variable name="x4" public_interface="in" units="microg_l"/>
      <variable name="x5" public_interface="in" units="microg_l"/>
      <variable name="a16" public_interface="in" units="flux"/>
      <variable name="a17" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a18" public_interface="in" units="second_order_rate_constant"/>
      <variable name="a19" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a20" public_interface="in" units="second_order_rate_constant"/>
      <variable name="a21" public_interface="in" units="per_microg_l"/>
      <variable name="a22" public_interface="in" units="per_microg_l2"/>
      <variable name="a23" public_interface="in" units="per_microg_l"/>
      <variable name="a24" public_interface="in" units="per_microg_l2"/>
      <variable name="a25" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a26" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>a16</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>a17</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a18</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a19</ci>
                           <ci>x2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a20</ci>
                           <apply>
                              <power/>
                              <ci>x2</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>a21</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a22</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a23</ci>
                           <ci>x2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a24</ci>
                           <apply>
                              <power/>
                              <ci>x2</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a25</ci>
                     <ci>x4</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a26</ci>
                     <ci>x5</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_3_</ci>
                  <ci>x3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="x4">
      <variable cmeta:id="x4_x4" initial_value="60.61" name="x4" public_interface="out" units="microg_l"/>
      <variable name="lambda_4_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x3" public_interface="in" units="microg_l"/>
      <variable name="a27" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x4</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>a27</ci>
                  <ci>x3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_4_</ci>
                  <ci>x4</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="x5">
      <variable cmeta:id="x5_x5" initial_value="12.61" name="x5" public_interface="out" units="microg_l"/>
      <variable name="lambda_5_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x3" public_interface="in" units="microg_l"/>
      <variable name="a28" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x5</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>a28</ci>
                  <ci>x3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_5_</ci>
                  <ci>x5</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="model_parameters">
      <variable initial_value="0.000017" name="a1" public_interface="out" units="flux"/>
      <variable initial_value="0.0023" name="a2" public_interface="out" units="flux"/>
      <variable initial_value="0.6" name="a3" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="45" name="a4" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="36" name="a5" public_interface="out" units="per_microg_l"/>
      <variable initial_value="216" name="a6" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="0.28" name="a7" public_interface="out" units="per_microg_l"/>
      <variable initial_value="0.36" name="a8" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="0.0003" name="a9" public_interface="out" units="flux"/>
      <variable initial_value="0.18" name="a10" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="150" name="a11" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="18" name="a12" public_interface="out" units="per_microg_l"/>
      <variable initial_value="460" name="a13" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="0.46" name="a14" public_interface="out" units="per_microg_l"/>
      <variable initial_value="0.1" name="a15" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="0.04" name="a16" public_interface="out" units="flux"/>
      <variable initial_value="150" name="a17" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="3800" name="a18" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="57" name="a19" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="2600" name="a20" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="200" name="a21" public_interface="out" units="per_microg_l"/>
      <variable initial_value="9400" name="a22" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="10" name="a23" public_interface="out" units="per_microg_l"/>
      <variable initial_value="320" name="a24" public_interface="out" units="per_microg_l2"/>
      <variable initial_value="0.04" name="a25" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.00097" name="a26" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.57" name="a27" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.0017" name="a28" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.059" name="lambda_1" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.028" name="lambda_2" public_interface="out" units="first_order_rate_constant"/>
      <variable name="lambda_3_" public_interface="out" units="first_order_rate_constant"/>
      <variable name="lambda_4_" public_interface="out" units="first_order_rate_constant"/>
      <variable name="lambda_5_" public_interface="out" units="first_order_rate_constant"/>
	  
	  <variable initial_value="0.0986" name="lambda_3" units="first_order_rate_constant"/>
      <variable initial_value="0.024" name="lambda_4" units="first_order_rate_constant"/>
      <variable initial_value="3e-5" name="lambda_5" units="first_order_rate_constant"/>
   
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_3_</ci>
            <apply>
               <plus/>
               <ci>lambda_3</ci>
               <ci>a27</ci>
               <ci>a28</ci>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_4_</ci>
            <apply>
               <plus/>
               <ci>lambda_4</ci>
               <ci>a25</ci>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_5_</ci>
            <apply>
               <plus/>
               <ci>lambda_5</ci>
               <ci>a26</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   
   <connection>
      <map_components component_1="x1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x3" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x4" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x5" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x1" component_2="x3"/>
      <map_variables variable_1="x1" variable_2="x1"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection>
      <map_components component_1="x1" component_2="x2"/>
      <map_variables variable_1="x1" variable_2="x1"/>
   </connection>
   <connection>
      <map_components component_1="x1" component_2="model_parameters"/>
      <map_variables variable_1="a1" variable_2="a1"/>
      <map_variables variable_1="a2" variable_2="a2"/>
      <map_variables variable_1="a3" variable_2="a3"/>
      <map_variables variable_1="a4" variable_2="a4"/>
      <map_variables variable_1="a5" variable_2="a5"/>
      <map_variables variable_1="a6" variable_2="a6"/>
      <map_variables variable_1="a7" variable_2="a7"/>
      <map_variables variable_1="a8" variable_2="a8"/>
      <map_variables variable_1="lambda_1" variable_2="lambda_1"/>
   </connection>
   <connection>
      <map_components component_1="x2" component_2="x3"/>
      <map_variables variable_1="x2" variable_2="x2"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection>
      <map_components component_1="x2" component_2="model_parameters"/>
      <map_variables variable_1="a9" variable_2="a9"/>
      <map_variables variable_1="a10" variable_2="a10"/>
      <map_variables variable_1="a11" variable_2="a11"/>
      <map_variables variable_1="a12" variable_2="a12"/>
      <map_variables variable_1="a13" variable_2="a13"/>
      <map_variables variable_1="a14" variable_2="a14"/>
      <map_variables variable_1="a15" variable_2="a15"/>
      <map_variables variable_1="lambda_2" variable_2="lambda_2"/>
   </connection>
   <connection>
      <map_components component_1="x4" component_2="x3"/>
      <map_variables variable_1="x4" variable_2="x4"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection>
      <map_components component_1="x5" component_2="x3"/>
      <map_variables variable_1="x5" variable_2="x5"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection>
      <map_components component_1="x3" component_2="model_parameters"/>
      <map_variables variable_1="a16" variable_2="a16"/>
      <map_variables variable_1="a17" variable_2="a17"/>
      <map_variables variable_1="a18" variable_2="a18"/>
      <map_variables variable_1="a19" variable_2="a19"/>
      <map_variables variable_1="a20" variable_2="a20"/>
      <map_variables variable_1="a21" variable_2="a21"/>
      <map_variables variable_1="a22" variable_2="a22"/>
      <map_variables variable_1="a23" variable_2="a23"/>
      <map_variables variable_1="a24" variable_2="a24"/>
      <map_variables variable_1="a25" variable_2="a25"/>
      <map_variables variable_1="a26" variable_2="a26"/>
      <map_variables variable_1="lambda_3_" variable_2="lambda_3_"/>
   </connection>
   <connection>
      <map_components component_1="x4" component_2="model_parameters"/>
      <map_variables variable_1="a27" variable_2="a27"/>
      <map_variables variable_1="lambda_4_" variable_2="lambda_4_"/>
   </connection>
   <connection>
      <map_components component_1="x5" component_2="model_parameters"/>
      <map_variables variable_1="a28" variable_2="a28"/>
      <map_variables variable_1="lambda_5_" variable_2="lambda_5_"/>
   </connection>
</model>