Location: NVC @ 4e32e1754f46 / coupling / Mnm1_p7.cellml

Author:
ssaf006 <ssaf006@aucklanduni.ac.nz>
Date:
2019-06-06 20:47:37+12:00
Desc:
fixed the nvc model
Permanent Source URI:
https://models.physiomeproject.org/workspace/4bd/rawfile/4e32e1754f465d7374b6a52d1371fba0a346c5c4/coupling/Mnm1_p7.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="Neurotransmitter" 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="Units.cellml">
        <units name="mM" units_ref="mM"/>
        <units name="per_s" units_ref="per_s"/>
        <units name="mM_per_s" units_ref="mM_per_s"/>
        <units name="per_mM" units_ref="per_mM"/>
        <units name="per_mM_s" units_ref="per_mM_s"/>
    </import>
    <import xlink:href="Neuron.cellml">
        <component component_ref="Neuron" name="Neuron"/>
    </import>
    <import xlink:href="Astrocyte.cellml">
        <component component_ref="Astrocyte" name="Astrocyte"/>
    </import>
    <import xlink:href="Blood.cellml">
        <component component_ref="Blood" name="Blood"/>
    </import>
    <component name="environment">
        <variable name="time" public_interface="out" units="second"/>
    </component>
    <component name="Mnm1">
        <variable name="t" private_interface="out" public_interface="in" units="second"/>
        <!-- State variables-->
        <variable name="q_stimulus" units="mM"/>
        <variable initial_value="0.0" name="q_inputDelay" private_interface="out" units="mM"/>
        <variable initial_value="0.0" name="q_Glu" private_interface="out" units="mM"/>
        <variable name="y" units="dimensionless"/>
        <variable name="q_oHb" private_interface="in" units="mM"/>
        <variable name="q_dHb" private_interface="in" units="mM"/>
        <variable name="q_O_2" private_interface="in" units="mM"/>
        <variable name="q_glucose" private_interface="in" units="mM"/>
        <!-- Parameters-->
        <variable name="k_flow" private_interface="in" units="per_s"/>
        <variable name="k_flow_basal" private_interface="in" units="per_s"/>
        <variable name="k_Stimulating" private_interface="in" units="per_s"/>
        <variable name="k_Inhibiting" private_interface="in" units="per_s"/>
        <variable initial_value="100.5114" name="k_y" units="dimensionless"/>
        <variable initial_value="22.2302" name="k_metabolic" units="per_s"/>
        <variable initial_value="107.8534" name="k_neurotrans" units="per_s"/>
        <variable initial_value="0.78" name="k_sink_Input" units="per_s"/>
        <variable initial_value="1.3363" name="k_sink_Glu" units="per_s"/>
        <!-- ODEs-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>q_stimulus</ci>
                <piecewise>
                    <piece>
                        <cn cellml:units="mM">1</cn>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>t</ci>
                                <cn cellml:units="second">0.5</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>t</ci>
                                <cn cellml:units="second">1.0</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <cn cellml:units="mM">0</cn>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_inputDelay</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <ci>k_metabolic</ci>
                        <ci>q_stimulus</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>k_sink_Input</ci>
                        <ci>q_inputDelay</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Glu</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <ci>k_neurotrans</ci>
                        <ci>q_stimulus</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>q_Glu</ci>
                        <ci>k_sink_Glu</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>y</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>k_y</ci>
                        <ci>q_oHb</ci>
                    </apply>
                    <ci>q_dHb</ci>
                </apply>
            </apply>
        </math>
    </component>
    <group>
        <relationship_ref relationship="encapsulation"/>
        <component_ref component="Mnm1">
            <component_ref component="Astrocyte"/>
            <component_ref component="Neuron"/>
            <component_ref component="Blood"/>
        </component_ref>
    </group>
    <connection>
        <map_components component_1="environment" component_2="Mnm1"/>
        <map_variables variable_1="time" variable_2="t"/>
    </connection>
    <connection>
        <map_components component_1="Mnm1" component_2="Blood"/>
        <map_variables variable_1="t" variable_2="t"/>
        <map_variables variable_1="q_inputDelay" variable_2="q_inputDelay"/>
        <map_variables variable_1="q_oHb" variable_2="q_oHb"/>
        <map_variables variable_1="q_dHb" variable_2="q_dHb"/>
        <map_variables variable_1="q_glucose" variable_2="q_glucose"/>
        <map_variables variable_1="q_O_2" variable_2="q_O_2"/>
        <map_variables variable_1="k_flow" variable_2="k_flow"/>
        <map_variables variable_1="k_flow_basal" variable_2="k_flow_basal"/>
        <map_variables variable_1="k_Stimulating" variable_2="k_Stimulating"/>
        <map_variables variable_1="k_Inhibiting" variable_2="k_Inhibiting"/>
    </connection>
    <connection>
        <map_components component_1="Mnm1" component_2="Astrocyte"/>
        <map_variables variable_1="t" variable_2="t"/>
        <map_variables variable_1="q_Glu" variable_2="q_Glu"/>
    </connection>
    <connection>
        <map_components component_1="Mnm1" component_2="Neuron"/>
        <map_variables variable_1="t" variable_2="t"/>
        <map_variables variable_1="q_Glu" variable_2="q_Glu"/>
    </connection>
    <connection>
        <map_components component_1="Astrocyte" component_2="Neuron"/>
        <map_variables variable_1="q_NO" variable_2="q_NO"/>
    </connection>
    <connection>
        <map_components component_1="Astrocyte" component_2="Blood"/>
        <map_variables variable_1="q_HETE" variable_2="q_HETE"/>
        <map_variables variable_1="q_EET" variable_2="q_EET"/>
        <map_variables variable_1="q_PG" variable_2="q_PG"/>
    </connection>
    <connection>
        <map_components component_1="Neuron" component_2="Blood"/>
        <map_variables variable_1="q_cGMP" variable_2="q_cGMP"/>
    </connection>
</model>