Location: cellLib @ 7a0891202d18 / Protocols / X_clamp_protocol.cellml

Author:
WeiweiAi <wai484@aucklanduni.ac.nz>
Date:
2022-08-17 11:26:32+12:00
Desc:
Add readMatrices and writeBGparas; modify whichComp to be string array; update readme
Permanent Source URI:
https://models.cellml.org/workspace/6bc/rawfile/7a0891202d18fe6d01956b6768e3f10c3dc802c6/Protocols/X_clamp_protocol.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="X_clamp_protocol" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <!-- clamp the concentration for the activation behaviour-->
    <import xlink:href="../Components/units.cellml">
        <units name="mM" units_ref="mM"/>
    </import>
    <component name="X_clamp_protocol">
        <variable name="time" public_interface="in" units="second"/>
        <variable cmeta:id="steady_state_time" name="t_ss" public_interface="in" units="second"/>
        <variable cmeta:id="holding_concentration" name="X_Holding" public_interface="in" units="mM"/>
        <variable cmeta:id="activation_duration" name="t_act" public_interface="in" units="second"/>
        <variable cmeta:id="test_concentration" name="X_Test" public_interface="in" units="mM"/>
        <variable name="X" public_interface="out" units="mM"/>
        <!--get the model to steady state at the holding concentration.-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply id="concentration_clamping">
                <eq/>
                <ci>X</ci>
                <piecewise>
                    <piece>
                        <ci>X_Holding</ci>
                        <apply>
                            <lt/>
                            <ci>time</ci>
                            <ci>t_ss</ci>
                        </apply>
                    </piece>
                    <piece>
                        <ci>X_Test</ci>
                        <apply>
                            <and/>
                            <apply>
                                <geq/>
                                <ci>time</ci>
                                <ci>t_ss</ci>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <apply>
                                    <plus/>
                                    <ci>t_act</ci>
                                    <ci>t_ss</ci>
                                </apply>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <ci>X_Holding</ci>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
    </component>
</model>