Location: Sten2020-Neurovascular Coupling Model @ 6c4ba7884bad / sten2020_parameters_Fig3D.cellml

Author:
schdemp <sdem348@aucklanduni.ac.nz>
Date:
2022-06-16 18:11:58-04:00
Desc:
Change parameter file names and subsequent loading name strings
Permanent Source URI:
https://models.physiomeproject.org/workspace/8a2/rawfile/6c4ba7884bad10ef792565cfc3d30eb9fec6748c/sten2020_parameters_Fig3D.cellml

<?xml version='1.0'?>
<model name="sten2020_base_parameters" 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="sten2020_units.cellml">
        <units name="per_sec" units_ref="per_sec"/>
    </import>
    <component name="parameters">
        <variable name="tstart" public_interface="out" units="second"/>
        <variable name="tend" public_interface="out" units="second"/>
        <variable name="toggle" units="dimensionless"/>
        <variable name="k_u_NO" public_interface="out" units="per_sec"/>
        <variable name="k_u_NPY" public_interface="out" units="per_sec"/>
        <variable name="k_u_Pyr" public_interface="out" units="per_sec"/>
        <variable name="initCA_NO" public_interface="out" units="dimensionless"/>
        <variable name="initCA_NPY" public_interface="out" units="dimensionless"/>
        <variable name="initCA_Pyr" public_interface="out" units="dimensionless"/>
        <variable name="initNO" public_interface="out" units="dimensionless"/>
        <variable name="initNO_vsm" public_interface="out" units="dimensionless"/>
        <variable name="initNPY" public_interface="out" units="dimensionless"/>
        <variable name="initNPY_vsm" public_interface="out" units="dimensionless"/>
        <variable name="initAA" public_interface="out" units="dimensionless"/>
        <variable name="initPGE2" public_interface="out" units="dimensionless"/>
        <variable name="initPGE2_vsm" public_interface="out" units="dimensionless"/>
        <variable name="initS_o" public_interface="out" units="dimensionless"/>
        <variable initial_value="10" name="k_Ca" public_interface="out" units="per_sec"/>
        <variable name="k_pf1" public_interface="out" units="per_sec"/>
        <variable name="k_pf2" public_interface="out" units="per_sec"/>
        <variable name="k_in1" public_interface="out" units="per_sec"/>
        <variable name="k_in2" public_interface="out" units="per_sec"/>
        <variable name="k_inf1" public_interface="out" units="per_sec"/>
        <variable name="k_inf2" public_interface="out" units="per_sec"/>
        <variable name="sink_N_NO" public_interface="out" units="per_sec"/>
        <variable name="sink_N_NPY" public_interface="out" units="per_sec"/>
        <variable name="sink_N_Pyr" public_interface="out" units="per_sec"/>
        <variable name="sink_Ca_NO" public_interface="out" units="per_sec"/>
        <variable name="sink_Ca_NPY" public_interface="out" units="per_sec"/>
        <variable name="sink_Ca_Pyr" public_interface="out" units="per_sec"/>
        <variable name="k_PL" public_interface="out" units="per_sec"/>
        <variable name="k_COX" public_interface="out" units="per_sec"/>
        <variable name="k_PGE2" public_interface="out" units="per_sec"/>
        <variable name="sink_PGE2" public_interface="out" units="per_sec"/>
        <variable name="k_NOS" public_interface="out" units="per_sec"/>
        <variable name="k_NO" public_interface="out" units="per_sec"/>
        <variable name="sink_NO" public_interface="out" units="per_sec"/>
        <variable name="k_NPY" public_interface="out" units="per_sec"/>
        <variable name="V_max" public_interface="out" units="per_sec"/>
        <variable name="K_M" public_interface="out" units="dimensionless"/>
        <variable name="sink_NPY" public_interface="out" units="per_sec"/>
        <variable name="k_y1" public_interface="out" units="dimensionless"/>
        <variable name="k_y2" public_interface="out" units="dimensionless"/>
        <variable name="k_y3" public_interface="out" units="dimensionless"/>
        <!-- Variables that change between cases-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>tstart</ci>
                <cn cellml:units="second">0.55</cn>
            </apply>
            <apply>
                <eq/>
                <ci>tend</ci>
                <cn cellml:units="second">1</cn>
            </apply>
            <apply>
                <eq/>
                <ci>toggle</ci>
                <cn cellml:units="dimensionless">0</cn>
            </apply>
            <!--1 for awake, 0 for anesthesia-->
            <apply>
                <eq/>
                <ci>k_u_NO</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.0892</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_u_NPY</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.5481</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_u_Pyr</ci>
                <cn cellml:units="per_sec">0.0</cn>
            </apply>
            <!-- initial values-->
            <apply>
                <eq/>
                <ci>initCA_NO</ci>
                <cn cellml:units="dimensionless">7.12853031505862</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initCA_NPY</ci>
                <cn cellml:units="dimensionless">1.98655228530129</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initCA_Pyr</ci>
                <cn cellml:units="dimensionless">0.113005608843197</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initNO</ci>
                <cn cellml:units="dimensionless">2944.42163942016</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initNO_vsm</ci>
                <cn cellml:units="dimensionless">352.127546590394</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initNPY</ci>
                <cn cellml:units="dimensionless">0.00234534285336482</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initNPY_vsm</ci>
                <cn cellml:units="dimensionless">0.512153329031501</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initAA</ci>
                <cn cellml:units="dimensionless">6.22443590826843</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initPGE2</ci>
                <cn cellml:units="dimensionless">0.0405508534942328</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initPGE2_vsm</ci>
                <cn cellml:units="dimensionless">0.0840427006050188</cn>
            </apply>
            <apply>
                <eq/>
                <ci>initS_o</ci>
                <cn cellml:units="dimensionless">81.9425413806638</cn>
            </apply>
            <!--Awake or Anesthetic coefficients -->
            <apply>
                <eq/>
                <ci>k_pf1</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">0.4986</cn>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">0.4779</cn>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>k_pf2</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">1.3580</cn>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.1161</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>k_in1</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.8595</cn>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">0.5473</cn>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>k_in2</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">3.5</cn>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.2593</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>k_inf1</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.0240</cn>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">3.5000</cn>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>k_inf2</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">3.5</cn>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.2028</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>sink_N_NO</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.3303</cn>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.1079</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>sink_N_NPY</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.1175</cn>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.2502</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>sink_N_Pyr</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.1905</cn>
                            </apply>
                        </apply>
                        <apply>
                            <gt/>
                            <ci>toggle</ci>
                            <cn cellml:units="dimensionless">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="per_sec">1</cn>
                            <apply>
                                <power/>
                                <cn cellml:units="dimensionless">10</cn>
                                <cn cellml:units="dimensionless">0.7872</cn>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <!-- Coefficients that remain constant-->
            <apply>
                <eq/>
                <ci>sink_Ca_NO</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.1470</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sink_Ca_NPY</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.7019</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sink_Ca_Pyr</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">1.9469</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_PL</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.9056</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_COX</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.8354</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_PGE2</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">1.3507</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sink_PGE2</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">1.0342</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_NOS</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">2.1490</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_NO</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.4670</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sink_NO</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <cn cellml:units="dimensionless">0.4553</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_NPY</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.9595</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_max</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.6541</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K_M</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">4.4006</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sink_NPY</ci>
                <apply>
                    <times/>
                    <cn cellml:units="per_sec">1</cn>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.3708</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_y1</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">0.2763</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_y2</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <cn cellml:units="dimensionless">2.7055</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_y3</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <cn cellml:units="dimensionless">2.4582</cn>
                </apply>
            </apply>
        </math>
    </component>
</model>