Location: Wodarz, Sierro, Klenerman, 2007 @ 91cf00619cbc / wodarz_sierro_klenerman_2007_c.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2010-07-13 11:49:43+12:00
Desc:
Added xul and clickable sessions for all three models
Permanent Source URI:
https://models.cellml.org/workspace/wodarz_sierro_klenerman_2007/rawfile/91cf00619cbc2d3d0c654ba06e8c3226985e786e/wodarz_sierro_klenerman_2007_c.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 8/01/2010 at 1:28:35 at p.m. using:

COR (0.9.31.1333)
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="wodarz_2007c" name="wodarz_2007c">

      <documentation xmlns="http://cellml.org/tmp-documentation">
            <article>
                  <articleinfo>
                        <title>Dynamics of killer T cell inflation in viral infections</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 represents the third model in the original paper (this builds on the first 'core' model which describes the basic dynamics of the cytotoxic T lymphocyte (CTL) infection, and the second model which includes a description of the CTL response to the virus infection, by incorporating a natural killer (NK) cell response to the viral infection).  The model runs in OpenCell to recreate the results in the published paper.  The units have been checked and they are consistent.  The model also runs in COR but due to the long duration of the simulation is not ideal for use in COR.
          </para>
			</section>
			        <sect1 id="sec_structure">
                        <title>Model Structure</title>
                        <para>ABSTRACT: Upon acute viral infection, a typical cytotoxic T lymphocyte (CTL) response is characterized by a phase of expansion and contraction after which it settles at a relatively stable memory level. Recently, experimental data from mice infected with murine cytomegalovirus (MCMV) showed different and unusual dynamics. After acute infection had resolved, some antigen specific CTL started to expand over time despite the fact that no replicative virus was detectable. This phenomenon has been termed as "CTL memory inflation". In order to examine the dynamics of this system further, we developed a mathematical model analysing the impact of innate and adaptive immune responses. According to this model, a potentially important contributor to CTL inflation is competition between the specific CTL response and an innate natural killer (NK) cell response. Inflation occurs most readily if the NK cell response is more efficient than the CTL at reducing virus load during acute infection, but thereafter maintains a chronic virus load which is sufficient to induce CTL proliferation. The model further suggests that weaker NK cell mediated protection can correlate with more pronounced CTL inflation dynamics over time. We present experimental data from mice infected with MCMV which are consistent with the theoretical predictions. This model provides valuable information and may help to explain the inflation of CMV specific CD8+T cells seen in humans as they age.</para>
						
						<para>The original paper reference is cited below:</para>
                        <para>
                              Dynamics of killer T cell inflation in viral infections, Dominik Wodarz, Sophie Sierro and Paul Klenerman, 2007, <emphasis>Journal of the Royal Society, Interface</emphasis>, volume 4, issue 14, 533-543. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/17251133">PubMed ID: 17251133</ulink>
                        </para>
                        
                        <informalfigure float="0" id="fig_reaction_diagram">
                              <mediaobject>
                                    <imageobject>
                                          <objectinfo>
                                                <title>model diagram</title>
                                          </objectinfo>
                                          <imagedata fileref="wodarz_2007c.png"/>
                                    </imageobject>
                              </mediaobject>
                              <caption>Schematic diagram of the different components and variables included in the mathematical model of the dynamics of killer T cell inflation.  The core model comprises five components: susceptible host cells (x), free virus particles (v), early-infected cells (y0), late-infected cells (y1), and latently-infected cells (L).  When susceptible host cells become infected with virus, the product cells are either productively infected or are latently infected.  Productively infected cells can be divided into two subsets; those which express early-gene products and those which express late-gene products.  The latter produce new free virus particles and the replication cycle is completed.  Latently infected cells are initially silent but can become activated, giving rise to more productively infected cells.  Virus-specific CTL (za) divide on antigenic stimulation (ca) and kill virus-infected host cells through lysis (pa).  NK cells (zi) can also divide upon antigenic stimulation (ci) and are able to kill virus-infected host cells (pi).</caption>
                        </informalfigure>
                     
                  </sect1>
            </article>
      </documentation>
				
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#x">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">susceptible host cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#y0">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">y0</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">infected cells expressing early-viral gene products</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#y1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">y1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">infected cells expressing late-viral gene products</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#L">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">L</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">latently infected cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#v">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">v</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">free virus particles</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#za">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">za</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">CMV specific effector CTL</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#zi">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">zi</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">NK cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#R0">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R0</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">basic reproductive ratio of the virus</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Bag rdf:about="rdf:#36fa7a2a-6880-4817-a9a0-2a365badeba8">
			<rdf:li>immunology</rdf:li>
			<rdf:li>viral dynamics</rdf:li>
			<rdf:li>murine cytomegalovirus</rdf:li>
			<rdf:li>CD8+ T cell</rdf:li>
			<rdf:li>CTL</rdf:li>
			<rdf:li>NK cell</rdf:li>
		</rdf:Bag>
		<rdf:Seq rdf:about="rdf:#28c0b158-a4e5-4f9f-bbd1-06cb4105f6cc">
			<rdf:li rdf:resource="rdf:#27c5b7d5-640b-4564-af51-db589d2195df"/>
			<rdf:li rdf:resource="rdf:#c04ec10a-78e3-4995-9f3d-8fdcbeaa32ae"/>
			<rdf:li rdf:resource="rdf:#9b16f48b-738f-4231-bbb5-ce1be6974ee8"/>
		</rdf:Seq>
		<rdf:Description rdf:about="">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Dynamics of killer T cell inflation in viral infections (Core Model + CTL and NK Response)</dc:title>
			<dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
			<cmeta:comment rdf:resource="rdf:#83f872ff-dd58-4fe9-87f4-66c3a0a25a70"/>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#dd5dd081-33af-420b-95e7-4a3eeec7ed71"/>
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#b022109f-9bb8-41cf-9063-b04a50f600e3"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9b16f48b-738f-4231-bbb5-ce1be6974ee8">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#8a73e1ef-125e-408a-ba79-17d20857a07d"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#3c5e9454-560e-4a23-a5ab-454bf4edc15f">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Sophie</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Sierro</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#83f872ff-dd58-4fe9-87f4-66c3a0a25a70">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#6c2c9ce3-8a57-4211-a42d-a938aadb8980"/>
			<rdf:value/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9d647835-a41d-4246-b3b8-13035a7ac105">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lloyd</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#295c79ef-4e47-492a-8dcb-8a8b5b133549">
			<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
			<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="#wodarz_2007c">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
        Wodarz, Sierro and Klenerman's 2007 mathematical model of the dynamics of killer T cell inflation in viral infections.
      </dc:title>
			<cmeta:comment rdf:resource="rdf:#ab1a9998-754d-4b5a-b5f6-1e9f659d0f3d"/>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#a6c14a32-6c68-4a46-9ca7-2bc7477d7476"/>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#ca607a98-19b5-4bbf-b30d-6ee480467ed1"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8a73e1ef-125e-408a-ba79-17d20857a07d">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Paul</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Klenerman</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#79acfbcd-ee38-4e42-a432-7204fb44cf65">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of the Royal Society, Interface</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#6c2c9ce3-8a57-4211-a42d-a938aadb8980">
			<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#dd5dd081-33af-420b-95e7-4a3eeec7ed71">
			<dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-07-17T00:00:00+00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b022109f-9bb8-41cf-9063-b04a50f600e3">
			<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#e308f8cf-b70f-405b-8f71-2978007ba7a7"/>
			<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#295c79ef-4e47-492a-8dcb-8a8b5b133549"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#9d647835-a41d-4246-b3b8-13035a7ac105"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1b8d23ee-3fc4-4192-87ad-d1507eb09db7">
			<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#27c5b7d5-640b-4564-af51-db589d2195df">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#2bfe5905-6272-4166-abc7-ee002b0ea32a"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#660e895b-ff7c-4911-8eb9-e67561385883">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#28c0b158-a4e5-4f9f-bbd1-06cb4105f6cc"/>
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Dynamics of killer T cell inflation in viral infections</dc:title>
			<bqs:volume xmlns:bqs="http://www.cellml.org/bqs/1.0#">4</bqs:volume>
			<bqs:first_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">533</bqs:first_page>
			<bqs:Journal xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#79acfbcd-ee38-4e42-a432-7204fb44cf65"/>
			<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#e3324fc8-7011-4989-9e3f-9e8f4b33a44b"/>
			<bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">543</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a6c14a32-6c68-4a46-9ca7-2bc7477d7476">
			<bqs:Pubmed_id xmlns:bqs="http://www.cellml.org/bqs/1.0#">17251133</bqs:Pubmed_id>
			<bqs:JournalArticle xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#660e895b-ff7c-4911-8eb9-e67561385883"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2bfe5905-6272-4166-abc7-ee002b0ea32a">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Dominik</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Wodarz</vCard:Family>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a78e90ae-f8df-44ea-bad3-fe9ffc126049">
			<bqs:subject_type xmlns:bqs="http://www.cellml.org/bqs/1.0#">keyword</bqs:subject_type>
			<rdf:value rdf:resource="rdf:#36fa7a2a-6880-4817-a9a0-2a365badeba8"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#ca607a98-19b5-4bbf-b30d-6ee480467ed1">
			<dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#a78e90ae-f8df-44ea-bad3-fe9ffc126049"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e308f8cf-b70f-405b-8f71-2978007ba7a7">
			<vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The University of Auckland</vCard:Orgname>
			<vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The Bioengineering Institute</vCard:Orgunit>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#ab1a9998-754d-4b5a-b5f6-1e9f659d0f3d">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#1b8d23ee-3fc4-4192-87ad-d1507eb09db7"/>
			<rdf:value>This is a CellML description of Wodarz, Sierro and Klenerman's 2007 mathematical model of the dynamics of killer T cell inflation in viral infections.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#c04ec10a-78e3-4995-9f3d-8fdcbeaa32ae">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#3c5e9454-560e-4a23-a5ab-454bf4edc15f"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e3324fc8-7011-4989-9e3f-9e8f4b33a44b">
			<dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-22</dcterms:W3CDTF>
		</rdf:Description>
	</rdf:RDF>   
	
   <units name="hour">
      <unit multiplier="3600" units="second"/>
   </units>
   <units name="first_order_rate_constant">
      <unit exponent="-1" units="hour"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="hour"/>
   </component>
   <component name="x">
      <variable cmeta:id="x_x" initial_value="1" name="x" public_interface="out" units="dimensionless"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="gamma" public_interface="in" units="first_order_rate_constant"/>
      <variable name="lambda" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="v" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x</ci>
            </apply>
            <apply>
               <minus/>
               <ci>lambda</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>d</ci>
                     <ci>x</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta</ci>
                     <ci>x</ci>
                     <ci>v</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <ci>x</ci>
                     <ci>v</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="y0">
      <variable cmeta:id="y0_y0" initial_value="0" name="y0" public_interface="out" units="dimensionless"/>
      <variable name="beta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x" public_interface="in" units="dimensionless"/>
      <variable name="a0" public_interface="in" units="first_order_rate_constant"/>
      <variable name="eta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="phi" public_interface="in" units="first_order_rate_constant"/>
      <variable name="L" public_interface="in" units="dimensionless"/>
      <variable name="v" public_interface="in" units="dimensionless"/>
      <variable name="pa" public_interface="in" units="first_order_rate_constant"/>
      <variable name="za" public_interface="in" units="dimensionless"/>
      <variable name="p_i" public_interface="in" units="first_order_rate_constant"/>
      <variable name="zi" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y0</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta</ci>
                     <ci>x</ci>
                     <ci>v</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>a0</ci>
                        <ci>y0</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>eta</ci>
                        <ci>y0</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>pa</ci>
                        <ci>y0</ci>
                        <ci>za</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>p_i</ci>
                        <ci>y0</ci>
                        <ci>zi</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>phi</ci>
                  <ci>L</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="y1">
      <variable cmeta:id="y1_y1" initial_value="0" name="y1" public_interface="out" units="dimensionless"/>
      <variable name="a1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="eta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="y0" public_interface="in" units="dimensionless"/>
      <variable name="pa" public_interface="in" units="first_order_rate_constant"/>
      <variable name="za" public_interface="in" units="dimensionless"/>
      <variable name="p_i" public_interface="in" units="first_order_rate_constant"/>
      <variable name="zi" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>eta</ci>
                  <ci>y0</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>a1</ci>
                     <ci>y1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>pa</ci>
                     <ci>y1</ci>
                     <ci>za</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>p_i</ci>
                     <ci>y1</ci>
                     <ci>zi</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L">
      <variable cmeta:id="L_L" initial_value="0" name="L" public_interface="out" units="dimensionless"/>
      <variable name="gamma" public_interface="in" units="first_order_rate_constant"/>
      <variable name="x" public_interface="in" units="dimensionless"/>
      <variable name="v" public_interface="in" units="dimensionless"/>
      <variable name="phi" public_interface="in" units="first_order_rate_constant"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>L</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>gamma</ci>
                  <ci>x</ci>
                  <ci>v</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>phi</ci>
                     <ci>L</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>d</ci>
                     <ci>L</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="v">
      <variable cmeta:id="v_v" initial_value="1" name="v" public_interface="out" units="dimensionless"/>
      <variable name="k" public_interface="in" units="first_order_rate_constant"/>
      <variable name="y1" public_interface="in" units="dimensionless"/>
      <variable name="u" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>v</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k</ci>
                  <ci>y1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>u</ci>
                  <ci>v</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_8">
      <variable initial_value="0" name="m_8" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_7" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_8</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_7</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_8</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_7">
      <variable initial_value="0" name="m_7" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_6" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_7</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_6</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_7</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_6">
      <variable initial_value="0" name="m_6" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_5" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_6</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_5</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_6</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_5">
      <variable initial_value="0" name="m_5" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_4" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_5</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_4</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_5</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_4">
      <variable initial_value="0" name="m_4" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_3" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_4</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_4</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_3">
      <variable initial_value="0" name="m_3" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_2" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_2</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_2">
      <variable initial_value="0" name="m_2" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_1" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_1">
      <variable initial_value="0" name="m_1" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_0" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>r</ci>
                  <ci>m_0</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r</ci>
                  <ci>m_1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_0">
      <variable initial_value="0" name="m_0" public_interface="out" units="dimensionless"/>
      <variable name="r" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m_0</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>r</ci>
               </apply>
               <ci>m_0</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="za">
      <variable cmeta:id="za_za" initial_value="1" name="za" public_interface="out" units="dimensionless"/>
      <variable name="alpha" public_interface="in" units="first_order_rate_constant"/>
      <variable name="m_8" public_interface="in" units="dimensionless"/>
      <variable name="ca" public_interface="in" units="first_order_rate_constant"/>
      <variable name="y0" public_interface="in" units="dimensionless"/>
      <variable name="y1" public_interface="in" units="dimensionless"/>
      <variable name="ba" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>za</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>alpha</ci>
                     <ci>m_8</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>ca</ci>
                     <apply>
                        <plus/>
                        <ci>y0</ci>
                        <ci>y1</ci>
                     </apply>
                     <ci>za</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>ba</ci>
                  <ci>za</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="zi">
      <variable cmeta:id="zi_zi" initial_value="0.1" name="zi" public_interface="out" units="dimensionless"/>
      <variable name="xi" public_interface="in" units="first_order_rate_constant"/>
      <variable name="ci" public_interface="in" units="first_order_rate_constant"/>
      <variable name="y0" public_interface="in" units="dimensionless"/>
      <variable name="y1" public_interface="in" units="dimensionless"/>
      <variable name="bi" public_interface="in" units="first_order_rate_constant"/>
      <variable name="time" public_interface="in" units="hour"/>
      <variable name="log_zi" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>zi</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>xi</ci>
                  <apply>
                     <times/>
                     <ci>ci</ci>
                     <apply>
                        <plus/>
                        <ci>y0</ci>
                        <ci>y1</ci>
                     </apply>
                     <ci>zi</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>bi</ci>
                  <ci>zi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>log_zi</ci>
            <apply>
               <log/>
               <ci>zi</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="R0">
      <variable name="R0" units="dimensionless"/>
      <variable name="gamma" public_interface="in" units="first_order_rate_constant"/>
      <variable name="lambda" public_interface="in" units="first_order_rate_constant"/>
      <variable name="eta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="a0" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta" public_interface="in" units="first_order_rate_constant"/>
      <variable name="phi" public_interface="in" units="first_order_rate_constant"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>R0</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>lambda</ci>
                     <ci>eta</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>d</ci>
                     <ci>a1</ci>
                     <apply>
                        <plus/>
                        <ci>a0</ci>
                        <ci>eta</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>beta</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>gamma</ci>
                        <ci>phi</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>phi</ci>
                        <ci>d</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="model_parameters">
      <variable initial_value="0.1" name="a0" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.2" name="a1" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="k" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="u" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="10" name="lambda" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="d" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="beta" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.5" name="gamma" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.2" name="alpha" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="phi" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.01" name="eta" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.000001" name="pa" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="15.5" name="ca" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="r" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="ba" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="p_i" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="12" name="ci" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="bi" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.01" name="xi" public_interface="out" units="first_order_rate_constant"/>
   </component>
   <connection>
      <map_components component_1="x" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="y1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="v" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_8" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_7" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_6" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_5" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_4" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_3" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="m_0" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="za" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="zi" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="x" component_2="v"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="x" component_2="model_parameters"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="gamma" variable_2="gamma"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
      <map_variables variable_1="beta" variable_2="beta"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="v"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="x"/>
      <map_variables variable_1="x" variable_2="x"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="L"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="za"/>
      <map_variables variable_1="za" variable_2="za"/>
      <map_variables variable_1="y0" variable_2="y0"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="zi"/>
      <map_variables variable_1="zi" variable_2="zi"/>
      <map_variables variable_1="y0" variable_2="y0"/>
   </connection>
   <connection>
      <map_components component_1="y0" component_2="model_parameters"/>
      <map_variables variable_1="a0" variable_2="a0"/>
      <map_variables variable_1="pa" variable_2="pa"/>
      <map_variables variable_1="p_i" variable_2="p_i"/>
      <map_variables variable_1="beta" variable_2="beta"/>
      <map_variables variable_1="eta" variable_2="eta"/>
      <map_variables variable_1="phi" variable_2="phi"/>
   </connection>
   <connection>
      <map_components component_1="y1" component_2="y0"/>
      <map_variables variable_1="y0" variable_2="y0"/>
   </connection>
   <connection>
      <map_components component_1="y1" component_2="model_parameters"/>
      <map_variables variable_1="a1" variable_2="a1"/>
      <map_variables variable_1="pa" variable_2="pa"/>
      <map_variables variable_1="p_i" variable_2="p_i"/>
      <map_variables variable_1="eta" variable_2="eta"/>
   </connection>
   <connection>
      <map_components component_1="y1" component_2="za"/>
      <map_variables variable_1="za" variable_2="za"/>
      <map_variables variable_1="y1" variable_2="y1"/>
   </connection>
   <connection>
      <map_components component_1="y1" component_2="zi"/>
      <map_variables variable_1="zi" variable_2="zi"/>
      <map_variables variable_1="y1" variable_2="y1"/>
   </connection>
   <connection>
      <map_components component_1="L" component_2="v"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="L" component_2="x"/>
      <map_variables variable_1="x" variable_2="x"/>
   </connection>
   <connection>
      <map_components component_1="L" component_2="model_parameters"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="gamma" variable_2="gamma"/>
      <map_variables variable_1="phi" variable_2="phi"/>
   </connection>
   <connection>
      <map_components component_1="v" component_2="y1"/>
      <map_variables variable_1="y1" variable_2="y1"/>
   </connection>
   <connection>
      <map_components component_1="v" component_2="model_parameters"/>
      <map_variables variable_1="k" variable_2="k"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="m_8" component_2="m_7"/>
      <map_variables variable_1="m_7" variable_2="m_7"/>
   </connection>
   <connection>
      <map_components component_1="m_8" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_7" component_2="m_6"/>
      <map_variables variable_1="m_6" variable_2="m_6"/>
   </connection>
   <connection>
      <map_components component_1="m_7" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_6" component_2="m_5"/>
      <map_variables variable_1="m_5" variable_2="m_5"/>
   </connection>
   <connection>
      <map_components component_1="m_6" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_5" component_2="m_4"/>
      <map_variables variable_1="m_4" variable_2="m_4"/>
   </connection>
   <connection>
      <map_components component_1="m_5" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_4" component_2="m_3"/>
      <map_variables variable_1="m_3" variable_2="m_3"/>
   </connection>
   <connection>
      <map_components component_1="m_4" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_3" component_2="m_2"/>
      <map_variables variable_1="m_2" variable_2="m_2"/>
   </connection>
   <connection>
      <map_components component_1="m_3" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_2" component_2="m_1"/>
      <map_variables variable_1="m_1" variable_2="m_1"/>
   </connection>
   <connection>
      <map_components component_1="m_2" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_1" component_2="m_0"/>
      <map_variables variable_1="m_0" variable_2="m_0"/>
   </connection>
   <connection>
      <map_components component_1="m_1" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="m_0" component_2="model_parameters"/>
      <map_variables variable_1="r" variable_2="r"/>
   </connection>
   <connection>
      <map_components component_1="za" component_2="m_8"/>
      <map_variables variable_1="m_8" variable_2="m_8"/>
   </connection>
   <connection>
      <map_components component_1="za" component_2="model_parameters"/>
      <map_variables variable_1="alpha" variable_2="alpha"/>
      <map_variables variable_1="ca" variable_2="ca"/>
      <map_variables variable_1="ba" variable_2="ba"/>
   </connection>
   <connection>
      <map_components component_1="zi" component_2="model_parameters"/>
      <map_variables variable_1="xi" variable_2="xi"/>
      <map_variables variable_1="ci" variable_2="ci"/>
      <map_variables variable_1="bi" variable_2="bi"/>
   </connection>
   <connection>
      <map_components component_1="R0" component_2="model_parameters"/>
      <map_variables variable_1="a0" variable_2="a0"/>
      <map_variables variable_1="gamma" variable_2="gamma"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
      <map_variables variable_1="phi" variable_2="phi"/>
      <map_variables variable_1="a1" variable_2="a1"/>
      <map_variables variable_1="beta" variable_2="beta"/>
      <map_variables variable_1="eta" variable_2="eta"/>
      <map_variables variable_1="d" variable_2="d"/>
   </connection>
</model>