Rendering of the source text

<?xml version='1.0' encoding='UTF-8'?>
<model name="Periodic_stimulation" 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="mM" units_ref="mM"/>
        <units name="nS" units_ref="nS"/>
    </import>
    <import xlink:href="../Components/Membrane_potential.cellml">
        <component component_ref="Membrane_potential" name="Membrane_potential"/>
    </import>
    <import xlink:href="Periodic_IStim_protocol.cellml">
        <component component_ref="Periodic_IStim_protocol" name="Periodic_IStim_protocol"/>
    </import>
    <import xlink:href="../Components/time.cellml">
        <component component_ref="time" name="time"/>
    </import>
    <import xlink:href="../Components/default_values.cellml">
        <component component_ref="model_parameters" name="model_parameters"/>
        <component component_ref="initial_conditions" name="initial_conditions"/>
    </import>
    <!-- Connect the time source variable to the model and stimulus protocol -->
    <connection>
        <map_components component_1="time" component_2="Periodic_IStim_protocol"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="time" component_2="Membrane_potential"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <!-- Connect the parameters to the protocol definition -->
    <component name="Istim_parameters">
        <variable initial_value="10000" name="t_period" public_interface="out" units="ms"/>
        <variable initial_value="-57" name="V_ICCrest" public_interface="out" units="mV"/>
        <!--var V_ICCrest: mV {init: -58, pub: out};-->
        <variable initial_value="33.5" name="V_ICCamp" public_interface="out" units="mV"/>
        <!--var V_ICCamp: mV {init: 32, pub: out};-->
        <variable initial_value="300" name="t_ICCpeak" public_interface="out" units="ms"/>
        <variable initial_value="600" name="t_slope" public_interface="out" units="ms"/>
        <variable initial_value="12000" name="f_1" public_interface="out" units="ms"/>
        <variable initial_value="300" name="f_2" public_interface="out" units="ms"/>
        <variable initial_value="2.6" name="Gcouple" public_interface="out" units="nS"/>
    </component>
    <connection>
        <map_components component_1="Istim_parameters" component_2="Periodic_IStim_protocol"/>
        <map_variables variable_1="t_period" variable_2="t_period"/>
        <map_variables variable_1="V_ICCrest" variable_2="V_ICCrest"/>
        <map_variables variable_1="V_ICCamp" variable_2="V_ICCamp"/>
        <map_variables variable_1="t_ICCpeak" variable_2="t_ICCpeak"/>
        <map_variables variable_1="t_slope" variable_2="t_slope"/>
        <map_variables variable_1="f_1" variable_2="f_1"/>
        <map_variables variable_1="f_2" variable_2="f_2"/>
        <map_variables variable_1="Gcouple" variable_2="Gcouple"/>
    </connection>
    <!-- and connect the actual stimulus current to the model -->
    <connection>
        <map_components component_1="Periodic_IStim_protocol" component_2="Membrane_potential"/>
        <map_variables variable_1="I_stim" variable_2="I_stim"/>
        <map_variables variable_1="V" variable_2="V"/>
    </connection>
    <!-- SED-ML (and CSim simulation description above) is only able to address variables at the top level of the model hierarchy, so we need to define the variables we want SED-ML to be able to address here. -->
    <component name="outputs">
        <variable name="time" public_interface="in" units="ms"/>
        <variable name="I_stim" public_interface="in" units="pA"/>
        <variable name="V_ICC" public_interface="in" units="mV"/>
        <variable name="V" public_interface="in" units="mV"/>
        <variable name="I_ion" public_interface="in" units="pA"/>
        <variable name="I_CaL" public_interface="in" units="pA"/>
        <variable name="I_CaT" public_interface="in" units="pA"/>
        <variable name="I_Kv" public_interface="in" units="pA"/>
        <variable name="I_BK" public_interface="in" units="pA"/>
        <variable name="I_Na" public_interface="in" units="pA"/>
        <variable name="I_NCX" public_interface="in" units="pA"/>
        <variable name="I_NaK" public_interface="in" units="pA"/>
        <variable name="I_NsNa" public_interface="in" units="pA"/>
        <variable name="I_NsK" public_interface="in" units="pA"/>
        <variable name="Cai" public_interface="in" units="mM"/>
        <variable name="Nai" public_interface="in" units="mM"/>
        <variable name="Ki" public_interface="in" units="mM"/>
        <variable name="P_CaL" public_interface="in" units="dimensionless"/>
        <variable name="P_CaT" public_interface="in" units="dimensionless"/>
        <variable name="P_BK" public_interface="in" units="dimensionless"/>
        <variable name="P_Kv" public_interface="in" units="dimensionless"/>
        <variable name="P_Na" public_interface="in" units="dimensionless"/>
        <variable name="norm_CaL" public_interface="in" units="dimensionless"/>
        <variable name="norm_BK" public_interface="in" units="dimensionless"/>
        <variable name="norm_Na" public_interface="in" units="dimensionless"/>
        <variable name="d_CaT" public_interface="in" units="dimensionless"/>
        <variable name="f_CaT" public_interface="in" units="dimensionless"/>
        <variable name="x_Kv" public_interface="in" units="dimensionless"/>
        <variable name="y_Kv" public_interface="in" units="dimensionless"/>
    </component>
    <!-- and connect the output variable definitions to their actual source in the model -->
    <connection>
        <map_components component_1="Membrane_potential" component_2="outputs"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="I_ion" variable_2="I_ion"/>
        <map_variables variable_1="I_CaL" variable_2="I_CaL"/>
        <map_variables variable_1="I_CaT" variable_2="I_CaT"/>
        <map_variables variable_1="I_Kv" variable_2="I_Kv"/>
        <map_variables variable_1="I_BK" variable_2="I_BK"/>
        <map_variables variable_1="I_Na" variable_2="I_Na"/>
        <map_variables variable_1="I_NCX" variable_2="I_NCX"/>
        <map_variables variable_1="I_NaK" variable_2="I_NaK"/>
        <map_variables variable_1="I_NsNa" variable_2="I_NsNa"/>
        <map_variables variable_1="I_NsK" variable_2="I_NsK"/>
        <map_variables variable_1="Cai" variable_2="Cai"/>
        <map_variables variable_1="Nai" variable_2="Nai"/>
        <map_variables variable_1="Ki" variable_2="Ki"/>
        <map_variables variable_1="P_CaL" variable_2="P_CaL"/>
        <map_variables variable_1="P_CaT" variable_2="P_CaT"/>
        <map_variables variable_1="P_BK" variable_2="P_BK"/>
        <map_variables variable_1="P_Kv" variable_2="P_Kv"/>
        <map_variables variable_1="P_Na" variable_2="P_Na"/>
        <map_variables variable_1="norm_CaL" variable_2="norm_CaL"/>
        <map_variables variable_1="norm_BK" variable_2="norm_BK"/>
        <map_variables variable_1="norm_Na" variable_2="norm_Na"/>
        <map_variables variable_1="d_CaT" variable_2="d_CaT"/>
        <map_variables variable_1="f_CaT" variable_2="f_CaT"/>
        <map_variables variable_1="x_Kv" variable_2="x_Kv"/>
        <map_variables variable_1="y_Kv" variable_2="y_Kv"/>
    </connection>
    <connection>
        <map_components component_1="outputs" component_2="time"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <!-- and connect the time output to the source time variable -->
    <connection>
        <map_components component_1="Periodic_IStim_protocol" component_2="outputs"/>
        <map_variables variable_1="I_stim" variable_2="I_stim"/>
        <map_variables variable_1="V_ICC" variable_2="V_ICC"/>
    </connection>
    <!-- We change these parameters to perform sensitivity analysis-->
    <component name="g_parameters">
        <variable initial_value="1.44" name="g_CaL" public_interface="out" units="nS"/>
        <variable initial_value="1.0217" name="g_Kv" public_interface="out" units="nS"/>
        <variable initial_value="80" name="g_BK" public_interface="out" units="nS"/>
        <variable initial_value="25.1" name="g_Na" public_interface="out" units="nS"/>
    </component>
    <connection>
        <map_components component_1="g_parameters" component_2="Membrane_potential"/>
        <map_variables variable_1="g_CaL" variable_2="g_CaL"/>
        <map_variables variable_1="g_Kv" variable_2="g_Kv"/>
        <map_variables variable_1="g_BK" variable_2="g_BK"/>
        <map_variables variable_1="g_Na" variable_2="g_Na"/>
    </connection>
    <!-- We connect the default parameters through to the model-->
    <connection>
        <map_components component_1="model_parameters" component_2="Membrane_potential"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="Cm" variable_2="Cm"/>
        <map_variables variable_1="V_cell" variable_2="V_cell"/>
        <map_variables variable_1="Cao" variable_2="Cao"/>
        <map_variables variable_1="Ko" variable_2="Ko"/>
        <map_variables variable_1="Nao" variable_2="Nao"/>
        <map_variables variable_1="Q10Ca" variable_2="Q10Ca"/>
        <map_variables variable_1="Q10K" variable_2="Q10K"/>
        <map_variables variable_1="Q10Na" variable_2="Q10Na"/>
        <map_variables variable_1="T0CaL" variable_2="T0CaL"/>
        <map_variables variable_1="T0CaT" variable_2="T0CaT"/>
        <map_variables variable_1="T0K" variable_2="T0K"/>
        <map_variables variable_1="T0Na" variable_2="T0Na"/>
        <map_variables variable_1="n_CRT" variable_2="n_CRT"/>
        <map_variables variable_1="n_CaM" variable_2="n_CaM"/>
        <map_variables variable_1="CRT_total" variable_2="CRT_total"/>
        <map_variables variable_1="CaM_total" variable_2="CaM_total"/>
        <map_variables variable_1="K_D_CRT" variable_2="K_D_CRT"/>
        <map_variables variable_1="K_D_CaM" variable_2="K_D_CaM"/>
        <map_variables variable_1="g_CaT" variable_2="g_CaT"/>
        <map_variables variable_1="g_NsNa" variable_2="g_NsNa"/>
        <map_variables variable_1="g_NsK" variable_2="g_NsK"/>
        <map_variables variable_1="tau_dCaT" variable_2="tau_dCaT"/>
        <map_variables variable_1="tau_xKv" variable_2="tau_xKv"/>
        <map_variables variable_1="tau_yKv" variable_2="tau_yKv"/>
        <map_variables variable_1="k_on" variable_2="k_on"/>
        <map_variables variable_1="k_off_C" variable_2="k_off_C"/>
        <map_variables variable_1="k_off_O" variable_2="k_off_O"/>
        <map_variables variable_1="P_NCX" variable_2="P_NCX"/>
        <map_variables variable_1="K_mCa" variable_2="K_mCa"/>
        <map_variables variable_1="K_mNai" variable_2="K_mNai"/>
        <map_variables variable_1="k_sat" variable_2="k_sat"/>
        <map_variables variable_1="gamma" variable_2="gamma"/>
        <map_variables variable_1="P_NaK" variable_2="P_NaK"/>
        <map_variables variable_1="K_mK" variable_2="K_mK"/>
        <map_variables variable_1="K_mNa" variable_2="K_mNa"/>
    </connection>
    <!-- We connect the initial conditions through to the model-->
    <connection>
        <map_components component_1="initial_conditions" component_2="Membrane_potential"/>
        <map_variables variable_1="V_init" variable_2="V_init"/>
        <map_variables variable_1="Cai_total_init" variable_2="Cai_total_init"/>
        <map_variables variable_1="Cai_init" variable_2="Cai_init"/>
        <map_variables variable_1="Ki_init" variable_2="Ki_init"/>
        <map_variables variable_1="Nai_init" variable_2="Nai_init"/>
        <map_variables variable_1="O_CaL_init" variable_2="O_CaL_init"/>
        <map_variables variable_1="C0_init" variable_2="C0_init"/>
        <map_variables variable_1="C1_init" variable_2="C1_init"/>
        <map_variables variable_1="C2_init" variable_2="C2_init"/>
        <map_variables variable_1="C3_init" variable_2="C3_init"/>
        <map_variables variable_1="C0Ca_init" variable_2="C0Ca_init"/>
        <map_variables variable_1="C1Ca_init" variable_2="C1Ca_init"/>
        <map_variables variable_1="C2Ca_init" variable_2="C2Ca_init"/>
        <map_variables variable_1="C3Ca_init" variable_2="C3Ca_init"/>
        <map_variables variable_1="ICa_init" variable_2="ICa_init"/>
        <map_variables variable_1="IVs_init" variable_2="IVs_init"/>
        <map_variables variable_1="IVf_init" variable_2="IVf_init"/>
        <map_variables variable_1="IVfCa_init" variable_2="IVfCa_init"/>
        <map_variables variable_1="IVsCa_init" variable_2="IVsCa_init"/>
        <map_variables variable_1="d_CaT_init" variable_2="d_CaT_init"/>
        <map_variables variable_1="f_CaT_init" variable_2="f_CaT_init"/>
        <map_variables variable_1="x_Kv_init" variable_2="x_Kv_init"/>
        <map_variables variable_1="y_Kv_init" variable_2="y_Kv_init"/>
        <map_variables variable_1="BK_C0_init" variable_2="BK_C0_init"/>
        <map_variables variable_1="BK_C1_init" variable_2="BK_C1_init"/>
        <map_variables variable_1="BK_C2_init" variable_2="BK_C2_init"/>
        <map_variables variable_1="BK_C3_init" variable_2="BK_C3_init"/>
        <map_variables variable_1="BK_C4_init" variable_2="BK_C4_init"/>
        <map_variables variable_1="O0_init" variable_2="O0_init"/>
        <map_variables variable_1="O1_init" variable_2="O1_init"/>
        <map_variables variable_1="O2_init" variable_2="O2_init"/>
        <map_variables variable_1="O3_init" variable_2="O3_init"/>
        <map_variables variable_1="O4_init" variable_2="O4_init"/>
        <map_variables variable_1="Na_C1_init" variable_2="Na_C1_init"/>
        <map_variables variable_1="Na_C2_init" variable_2="Na_C2_init"/>
        <map_variables variable_1="Na_C3_init" variable_2="Na_C3_init"/>
        <map_variables variable_1="I1_init" variable_2="I1_init"/>
        <map_variables variable_1="I2_init" variable_2="I2_init"/>
        <map_variables variable_1="O_Na_init" variable_2="O_Na_init"/>
    </connection>
</model>