Location: Model of Human Jejunal Smooth Muscle Cell Electrophysiology @ b98a10f69ee7 / Experiments / Periodic_IStim_protocol.cellml

Author:
WeiweiAi <wai484@aucklanduni.ac.nz>
Date:
2021-07-06 11:22:35+12:00
Desc:
Change the figures size to fit the Physiome paper; Adjust the solver setting to remove the notch in Fig6. Minor editting in README.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/692/rawfile/b98a10f69ee76a70243efdb648c598db5a866a2c/Experiments/Periodic_IStim_protocol.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="Periodic_IStim_protocol" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <import xlink:href="../Components/units.cellml">
        <units name="ms" units_ref="ms"/>
        <units name="pA" units_ref="pA"/>
        <units name="mV" units_ref="mV"/>
        <units name="nS" units_ref="nS"/>
    </import>
    <component name="Periodic_IStim_protocol">
        <variable name="I_stim" public_interface="out" units="pA"/>
        <variable name="V_ICC" public_interface="out" units="mV"/>
        <variable name="time" public_interface="in" units="ms"/>
        <variable name="V" public_interface="in" units="mV"/>
        <variable name="t_period" public_interface="in" units="ms"/>
        <variable name="Gcouple" public_interface="in" units="nS"/>
        <variable name="V_ICCrest" public_interface="in" units="mV"/>
        <variable name="V_ICCamp" public_interface="in" units="mV"/>
        <variable name="t_ICCpeak" public_interface="in" units="ms"/>
        <variable name="t_slope" public_interface="in" units="ms"/>
        <variable name="f_1" public_interface="in" units="ms"/>
        <variable name="f_2" public_interface="in" units="ms"/>
        <variable name="t" units="ms"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_stim</ci>
                <apply>
                    <times/>
                    <ci>Gcouple</ci>
                    <apply>
                        <minus/>
                        <ci>V</ci>
                        <ci>V_ICC</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>t</ci>
                <apply>
                    <rem/>
                    <ci>time</ci>
                    <ci>t_period</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_ICC</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <plus/>
                            <ci>V_ICCrest</ci>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>V_ICCamp</ci>
                                    <ci>t</ci>
                                </apply>
                                <ci>f_2</ci>
                            </apply>
                        </apply>
                        <apply>
                            <lt/>
                            <ci>t</ci>
                            <ci>t_ICCpeak</ci>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <plus/>
                            <ci>V_ICCrest</ci>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>V_ICCamp</ci>
                                    <apply>
                                        <plus/>
                                        <cn cellml:units="dimensionless">1</cn>
                                        <apply>
                                            <exp/>
                                            <apply>
                                                <divide/>
                                                <apply>
                                                    <minus/>
                                                    <ci>f_1</ci>
                                                </apply>
                                                <apply>
                                                    <times/>
                                                    <cn cellml:units="dimensionless">2</cn>
                                                    <ci>t_slope</ci>
                                                </apply>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="dimensionless">1</cn>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <minus/>
                                                <apply>
                                                    <minus/>
                                                    <ci>t</ci>
                                                    <ci>f_2</ci>
                                                </apply>
                                                <apply>
                                                    <times/>
                                                    <cn cellml:units="dimensionless">0.5</cn>
                                                    <ci>f_1</ci>
                                                </apply>
                                            </apply>
                                            <ci>t_slope</ci>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
    </component>
</model>