Location: Acid-Base Physiology @ bc71d1ee7761 / BG / StateDiagram.cellml
- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2017-01-23 18:49:42+13:00
- Desc:
- fix units
- Permanent Source URI:
- https://models.cellml.org/workspace/297/rawfile/bc71d1ee7761d7e7de319ebdaa317371ef939670/BG/StateDiagram.cellml
<?xml version='1.0'?>
<model name="StateDiagram" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<component name="state">
<variable name="t" public_interface="in" units="second"/>
<variable name="RT" public_interface="in" units="J_per_mol"/>
<variable name="z" public_interface="in" units="dimensionless"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="-0.06" name="u_E" units="J_per_C"/>
<variable name="e1" units="dimensionless"/>
<variable initial_value="1" name="rho_1" units="per_m"/>
<variable initial_value="1" name="rho_2" units="per_m"/>
<!-- Parameters-->
<variable name="K_q1" public_interface="in" units="m3_per_mol"/>
<variable name="K_q2" public_interface="in" units="m3_per_mol"/>
<variable name="K_q3" public_interface="in" units="m3_per_mol"/>
<variable name="K_q4" public_interface="in" units="m3_per_mol"/>
<variable name="K_q5" public_interface="in" units="m3_per_mol"/>
<variable name="K_q6" public_interface="in" units="m3_per_mol"/>
<variable name="K_q7" public_interface="in" units="m3_per_mol"/>
<variable name="K_q8" public_interface="in" units="m3_per_mol"/>
<variable name="K_q9" public_interface="in" units="m3_per_mol"/>
<variable name="K_q10" public_interface="in" units="m3_per_mol"/>
<variable name="K_Re1" public_interface="in" units="m_per_s"/>
<variable name="K_Re2" public_interface="in" units="mol_per_m3_s"/>
<variable name="K_Re3" public_interface="in" units="m_per_s"/>
<variable name="K_Re4" public_interface="in" units="m_per_s"/>
<variable name="K_Re5" public_interface="in" units="mol_per_m3_s"/>
<variable name="K_Re6" public_interface="in" units="m_per_s"/>
<variable name="K_Re7" public_interface="in" units="m_per_s"/>
<!-- State variables-->
<variable name="q1_init" public_interface="in" units="mol_per_m3"/>
<variable name="q2_init" public_interface="in" units="mol_per_m3"/>
<variable name="q3_init" public_interface="in" units="mol_per_m3"/>
<variable name="q4_init" public_interface="in" units="mol_per_m3"/>
<variable name="q5_init" public_interface="in" units="mol_per_m3"/>
<variable name="q6_init" public_interface="in" units="mol_per_m3"/>
<variable name="q7_init" public_interface="in" units="mol_per_m3"/>
<variable name="q8_init" public_interface="in" units="mol_per_m3"/>
<variable name="q9_init" public_interface="in" units="mol_per_m3"/>
<variable name="q10_init" public_interface="in" units="mol_per_m3"/>
<variable initial_value="0.0" name="q1_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q2_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q3_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q4_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q5_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q6_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q7_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q8_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q9_aux" units="mol_per_m3"/>
<variable initial_value="0.0" name="q10_aux" units="mol_per_m3"/>
<variable name="q1" public_interface="out" units="mol_per_m3"/>
<variable name="q2" public_interface="out" units="mol_per_m3"/>
<variable name="q3" public_interface="out" units="mol_per_m3"/>
<variable name="q4" public_interface="out" units="mol_per_m3"/>
<variable name="q5" public_interface="out" units="mol_per_m3"/>
<variable name="q6" public_interface="out" units="mol_per_m3"/>
<variable name="q7" public_interface="out" units="mol_per_m3"/>
<variable name="q8" public_interface="out" units="mol_per_m3"/>
<variable name="q9" public_interface="out" units="mol_per_m3"/>
<variable name="q10" public_interface="out" units="mol_per_m3"/>
<variable name="v_Re1" public_interface="out" units="mol_per_m2_s"/>
<variable name="v_Re2" public_interface="out" units="mol_per_m3_s"/>
<variable name="v_Re3" public_interface="out" units="mol_per_m2_s"/>
<variable name="v_Re4" public_interface="out" units="mol_per_m2_s"/>
<variable name="v_Re5" public_interface="out" units="mol_per_m3_s"/>
<variable name="v_Re6" public_interface="out" units="mol_per_m2_s"/>
<variable name="v_Re7" public_interface="out" units="mol_per_m2_s"/>
<!-- Conservation laws-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q1_aux</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>rho_2</ci>
</apply>
<ci>v_Re1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q2_aux</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>rho_1</ci>
<ci>v_Re1</ci>
</apply>
<ci>v_Re2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q3_aux</ci>
</apply>
<apply>
<minus/>
<ci>v_Re2</ci>
<apply>
<times/>
<ci>rho_1</ci>
<ci>v_Re3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q4_aux</ci>
</apply>
<apply>
<times/>
<ci>rho_2</ci>
<ci>v_Re3</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q5_aux</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>rho_2</ci>
</apply>
<ci>v_Re4</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q6_aux</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>rho_1</ci>
<ci>v_Re4</ci>
</apply>
<ci>v_Re5</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q7_aux</ci>
</apply>
<apply>
<minus/>
<ci>v_Re5</ci>
<apply>
<times/>
<ci>rho_1</ci>
<ci>v_Re6</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q8_aux</ci>
</apply>
<apply>
<times/>
<ci>rho_2</ci>
<ci>v_Re6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q9_aux</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_Re2</ci>
<ci>v_Re5</ci>
</apply>
<apply>
<times/>
<ci>rho_1</ci>
<ci>v_Re7</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q10_aux</ci>
</apply>
<apply>
<times/>
<ci>rho_2</ci>
<ci>v_Re7</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q1</ci>
<apply>
<plus/>
<ci>q1_aux</ci>
<ci>q1_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q2</ci>
<apply>
<plus/>
<ci>q2_aux</ci>
<ci>q2_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q3</ci>
<apply>
<plus/>
<ci>q3_aux</ci>
<ci>q3_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q4</ci>
<apply>
<plus/>
<ci>q4_aux</ci>
<ci>q4_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q5</ci>
<apply>
<plus/>
<ci>q5_aux</ci>
<ci>q5_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q6</ci>
<apply>
<plus/>
<ci>q6_aux</ci>
<ci>q6_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q7</ci>
<apply>
<plus/>
<ci>q7_aux</ci>
<ci>q7_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q8</ci>
<apply>
<plus/>
<ci>q8_aux</ci>
<ci>q8_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q9</ci>
<apply>
<plus/>
<ci>q9_aux</ci>
<ci>q9_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q10</ci>
<apply>
<plus/>
<ci>q10_aux</ci>
<ci>q10_init</ci>
</apply>
</apply>
<!-- Constitutive relations-->
<apply>
<eq/>
<ci>e1</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>u_E</ci>
<ci>RT</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re1</ci>
<apply>
<times/>
<ci>K_Re1</ci>
<apply>
<minus/>
<ci>q1</ci>
<ci>q2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re2</ci>
<apply>
<times/>
<ci>K_Re2</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>K_q3</ci>
<ci>q3</ci>
<ci>K_q9</ci>
<ci>q9</ci>
</apply>
<apply>
<times/>
<ci>K_q2</ci>
<ci>q2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re3</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_Re3</ci>
<ci>e1</ci>
<apply>
<minus/>
<ci>q3</ci>
<apply>
<times/>
<ci>q4</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re4</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_Re4</ci>
<ci>e1</ci>
<apply>
<minus/>
<ci>q5</ci>
<apply>
<times/>
<ci>q6</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re5</ci>
<apply>
<times/>
<ci>K_Re5</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>K_q7</ci>
<ci>q7</ci>
<ci>K_q9</ci>
<ci>q9</ci>
</apply>
<apply>
<times/>
<ci>K_q6</ci>
<ci>q6</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re6</ci>
<apply>
<times/>
<ci>K_Re6</ci>
<apply>
<minus/>
<ci>q7</ci>
<ci>q8</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Re7</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_Re7</ci>
<ci>e1</ci>
<apply>
<minus/>
<ci>q9</ci>
<apply>
<times/>
<ci>q10</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>e1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
</model>