- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2018-08-02 05:37:03+12:00
- Desc:
- fixing diffusion
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/436/rawfile/89cd53399e44e02d8db73780f55d1b1c8b85fe85/BK.cellml
<?xml version='1.0'?>
<model name="BK" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<component name="BK">
<variable name="q_K_c" public_interface="in" units="M"/>
<variable name="q_K_a" public_interface="in" units="M"/>
<variable name="q_Ca_c" public_interface="in" units="M"/>
<variable initial_value="0.0003" name="q_K_b_BK" public_interface="out" units="M"/>
<variable name="z_K" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_mol_K"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="V_a" public_interface="in" units="volt"/>
<variable name="V_A_BK" units="volt"/>
<variable name="J_BK" public_interface="out" units="umol_per_s"/>
<variable name="I_BK" public_interface="out" units="uA"/>
<variable initial_value="30" name="G_BK" units="siemens"/>
<variable name="P_BK" units="dimensionless"/>
<variable initial_value="1.7" name="eta" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<variable initial_value="0.1" name="alpha_BK" units="per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_BK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>q_K_b_BK</ci>
<ci>q_Ca_c</ci>
</apply>
<ci>eta</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_BK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>q_K_a</ci>
<ci>q_K_c</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_BK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
<ci>G_BK</ci>
<ci>alpha_BK</ci>
<ci>P_BK</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>V_a</ci>
<ci>V_A_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_BK</ci>
<apply>
<divide/>
<ci>I_BK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
</model>