<?xml version='1.0' encoding='UTF-8'?>
<model name="Calvo_Gallego_2022" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<!-----------------------------------------------------------------
Units
--------------------------------------------------------------- -->
<units name="day">
<unit multiplier="86400.0" units="second"/>
</units>
<units name="pM">
<unit prefix="pico" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="flux">
<unit units="pM"/>
<unit exponent="-1" units="day"/>
</units>
<units name="first_order_rate_constant">
<unit exponent="-1" units="day"/>
</units>
<units name="second_order_rate_constant">
<unit exponent="-1" units="pM"/>
<unit exponent="-1" units="day"/>
</units>
<units name="MPa">
<unit prefix="mega" units="pascal"/>
</units>
<units name="per_MPa">
<unit exponent="-1" prefix="mega" units="pascal"/>
</units>
<units name="g_per_cm3">
<unit units="gram"/>
<unit exponent="-3" prefix="centi" units="metre"/>
</units>
<!-----------------------------------------------------------------
Components
---------------------------------------------------------------
Environment-->
<component name="environment">
<variable name="time" public_interface="out" units="day"/>
</component>
<!-- OBp-->
<component name="OBp">
<variable initial_value="1.1631869976e-03" name="OBp" public_interface="out" units="pM"/>
<variable initial_value="1.1631869976e-03" name="OBp_t0" public_interface="out" units="pM"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OBp_in" units="flux"/>
<variable name="OBp_out" public_interface="out" units="flux"/>
<variable name="OBu_t0" public_interface="in" units="pM"/>
<variable name="DifferOBpNormal" units="first_order_rate_constant"/>
<variable name="D_OBp_t0" public_interface="in" units="first_order_rate_constant"/>
<variable initial_value="0.0001" name="TGFbNormal" units="pM"/>
<variable name="KD_TGF_beta_repress" public_interface="in" units="pM"/>
<variable name="ProlifOBpNormal" units="first_order_rate_constant"/>
<variable initial_value="0.5" name="frac_prolifOBp_vs_differOBu" units="dimensionless"/>
<variable initial_value="0.001" name="OBpNormal" units="pM"/>
<variable initial_value="0.005" name="OBp_sat" units="pM"/>
<variable name="DifferOBu" units="first_order_rate_constant"/>
<variable name="ProlifOBp" units="first_order_rate_constant"/>
<variable name="D_OBu_t0" public_interface="in" units="first_order_rate_constant"/>
<variable name="pd_OBp_t0" public_interface="in" units="flux"/>
<variable name="Pi_TGFbeta_OBu_act_t0" public_interface="in" units="dimensionless"/>
<variable name="Pi_TGFbeta_OBp_rep" public_interface="in" units="dimensionless"/>
<variable name="Pi_WNT" public_interface="in" units="dimensionless"/>
<variable name="Pi_WNT_0" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>DifferOBpNormal</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>D_OBp_t0</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>TGFbNormal</ci>
<ci>KD_TGF_beta_repress</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ProlifOBpNormal</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>frac_prolifOBp_vs_differOBu</ci>
<ci>DifferOBpNormal</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>OBpNormal</ci>
<ci>OBp_sat</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>DifferOBu</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>frac_prolifOBp_vs_differOBu</ci>
</apply>
<ci>D_OBu_t0</ci>
<ci>Pi_TGFbeta_OBu_act_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>ProlifOBp</ci>
<apply>
<times/>
<ci>ProlifOBpNormal</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>OBp</ci>
<ci>OBp_sat</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>OBp_in</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>DifferOBu</ci>
<ci>OBu_t0</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>ProlifOBp</ci>
<ci>OBp</ci>
<ci>Pi_WNT</ci>
</apply>
<ci>Pi_WNT_0</ci>
</apply>
<ci>pd_OBp_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OBp_out</ci>
<apply>
<times/>
<ci>D_OBp_t0</ci>
<ci>Pi_TGFbeta_OBp_rep</ci>
<ci>OBp</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>OBp</ci>
</apply>
<apply>
<minus/>
<ci>OBp_in</ci>
<ci>OBp_out</ci>
</apply>
</apply>
</math>
</component>
<!-- OBa-->
<component name="OBa">
<variable initial_value="9.1880026722e-04" name="OBa" public_interface="out" units="pM"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OBa_in" units="flux"/>
<variable name="OBa_out" units="flux"/>
<variable name="OBp_out" public_interface="in" units="flux"/>
<variable name="pd_OBa_t0" public_interface="in" units="flux"/>
<variable name="D_OBa_t0" public_interface="in" units="first_order_rate_constant"/>
<variable name="A_OBa_t0" public_interface="in" units="first_order_rate_constant"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>OBa_in</ci>
<apply>
<plus/>
<ci>OBp_out</ci>
<ci>pd_OBa_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OBa_out</ci>
<apply>
<times/>
<apply>
<plus/>
<ci>D_OBa_t0</ci>
<ci>A_OBa_t0</ci>
</apply>
<ci>OBa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>OBa</ci>
</apply>
<apply>
<minus/>
<ci>OBa_in</ci>
<ci>OBa_out</ci>
</apply>
</apply>
</math>
</component>
<!-- OCp-->
<component name="OCp">
<variable initial_value="1.3556879825e-03" name="OCp" public_interface="out" units="pM"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OCp_in" units="flux"/>
<variable name="OCp_out" public_interface="out" units="flux"/>
<variable name="OCu_t0" public_interface="in" units="pM"/>
<variable name="Pi_RANKL_act_OCp" public_interface="in" units="dimensionless"/>
<variable name="Pi_RANKL_act_OCu" public_interface="in" units="dimensionless"/>
<variable name="D_OCu_t0" public_interface="in" units="first_order_rate_constant"/>
<variable name="pd_OCp_t0" public_interface="in" units="flux"/>
<variable name="D_OCp_t0" public_interface="in" units="first_order_rate_constant"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>OCp_in</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>D_OCu_t0</ci>
<ci>OCu_t0</ci>
<ci>Pi_RANKL_act_OCu</ci>
</apply>
<ci>pd_OCp_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OCp_out</ci>
<apply>
<times/>
<ci>D_OCp_t0</ci>
<ci>Pi_RANKL_act_OCp</ci>
<ci>OCp</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>OCp</ci>
</apply>
<apply>
<minus/>
<ci>OCp_in</ci>
<ci>OCp_out</ci>
</apply>
</apply>
</math>
</component>
<!-- OCa-->
<component name="OCa">
<variable initial_value="1.8376005344e-05" name="OCa" public_interface="out" units="pM"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OCa_in" units="flux"/>
<variable name="OCa_out" units="flux"/>
<variable name="OCp_out" public_interface="in" units="flux"/>
<variable name="pd_OCa_t0" public_interface="in" units="flux"/>
<variable name="A_OCa_t0" public_interface="in" units="first_order_rate_constant"/>
<variable name="Pi_TGFbeta_OCa_act" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>OCa_in</ci>
<apply>
<plus/>
<ci>OCp_out</ci>
<ci>pd_OCa_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OCa_out</ci>
<apply>
<times/>
<ci>A_OCa_t0</ci>
<ci>Pi_TGFbeta_OCa_act</ci>
<ci>OCa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>OCa</ci>
</apply>
<apply>
<minus/>
<ci>OCa_in</ci>
<ci>OCa_out</ci>
</apply>
</apply>
</math>
</component>
<!-- Bone marix fraction in [%]-->
<component name="fbm">
<variable initial_value="15" name="fbm" public_interface="out" units="dimensionless"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OBa" public_interface="in" units="pM"/>
<variable name="OCa" public_interface="in" units="pM"/>
<variable name="K_form" public_interface="in" units="second_order_rate_constant"/>
<variable name="K_res" public_interface="in" units="second_order_rate_constant"/>
<variable name="dfbmdt" public_interface="out" units="first_order_rate_constant"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>dfbmdt</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>K_form</ci>
<ci>OBa</ci>
</apply>
<apply>
<times/>
<ci>K_res</ci>
<ci>OCa</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>fbm</ci>
</apply>
<ci>dfbmdt</ci>
</apply>
</math>
</component>
<!-- Osteocytes-->
<component name="OCY">
<variable initial_value="1.2734398187e-02" name="OCY" public_interface="out" units="pM"/>
<variable name="OCY_act" public_interface="out" units="pM"/>
<variable name="OCY_act_0" public_interface="out" units="pM"/>
<variable name="fbm" public_interface="in" units="dimensionless"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="fact_0" public_interface="in" units="pM"/>
<variable name="dfbmdt" public_interface="in" units="first_order_rate_constant"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>OCY_act</ci>
<apply>
<times/>
<ci>fact_0</ci>
<ci>fbm</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OCY_act_0</ci>
<apply>
<times/>
<ci>fact_0</ci>
<cn cellml:units="dimensionless">20</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>OCY</ci>
</apply>
<apply>
<times/>
<ci>fact_0</ci>
<ci>dfbmdt</ci>
</apply>
</apply>
</math>
</component>
<!-- Mineralization-->
<component name="vm">
<variable initial_value="3.5893417161e-01" name="vm" public_interface="out" units="dimensionless"/>
<variable name="time" public_interface="in" units="day"/>
<variable name="OBa" public_interface="in" units="pM"/>
<variable name="fbm" public_interface="in" units="dimensionless"/>
<variable name="XKAPPA" public_interface="in" units="first_order_rate_constant"/>
<variable name="vmmax" public_interface="in" units="dimensionless"/>
<variable name="K_form" public_interface="in" units="second_order_rate_constant"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>vm</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>XKAPPA</ci>
<apply>
<minus/>
<ci>vmmax</ci>
<ci>vm</ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>vm</ci>
<ci>K_form</ci>
<ci>OBa</ci>
</apply>
<ci>fbm</ci>
</apply>
</apply>
</apply>
</math>
</component>
<!-- PTH-->
<component name="PTH">
<variable name="PTH_tot" units="pM"/>
<variable name="Pi_PTH_act" units="dimensionless"/>
<variable name="Pi_PTH_rep" public_interface="out" units="dimensionless"/>
<variable name="Beta_PTH" public_interface="in" units="flux"/>
<variable name="P_PTH_d" public_interface="in" units="flux"/>
<variable name="Deg_PTH" public_interface="in" units="first_order_rate_constant"/>
<variable name="KD_PTH_act" public_interface="in" units="pM"/>
<variable name="KD_PTH_rep" public_interface="in" units="pM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>PTH_tot</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>Beta_PTH</ci>
<ci>P_PTH_d</ci>
</apply>
<ci>Deg_PTH</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_PTH_act</ci>
<apply>
<divide/>
<ci>PTH_tot</ci>
<apply>
<plus/>
<ci>PTH_tot</ci>
<ci>KD_PTH_act</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_PTH_rep</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>PTH_tot</ci>
<ci>KD_PTH_rep</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- RANK-RANKL-OPG-->
<component name="RANK_RANKL_OPG">
<variable name="RANK" units="pM"/>
<variable name="RANKL" units="pM"/>
<variable name="OPG" units="pM"/>
<variable name="RANKL_tot" units="pM"/>
<variable name="RANKL_max" units="pM"/>
<variable name="P_RANKL1" units="flux"/>
<variable name="P_RANKL2" units="flux"/>
<variable name="P_RANKL" units="flux"/>
<variable name="P_RANKL_d" public_interface="in" units="flux"/>
<variable name="Pi_RANKL_act_OCp" public_interface="out" units="dimensionless"/>
<variable name="Pi_RANKL_act_OCu" public_interface="out" units="dimensionless"/>
<variable name="OBp" public_interface="in" units="pM"/>
<variable name="OBp_t0" public_interface="in" units="pM"/>
<variable name="OBa" public_interface="in" units="pM"/>
<variable name="OCp" public_interface="in" units="pM"/>
<variable name="OCY_act" public_interface="in" units="pM"/>
<variable name="Pi_PTH_rep" public_interface="in" units="dimensionless"/>
<variable name="Pi_NO_PTH_act_rep" public_interface="in" units="dimensionless"/>
<variable name="N_RANK_OCp" public_interface="in" units="dimensionless"/>
<variable name="N_RANKL_OBp_max" public_interface="in" units="dimensionless"/>
<variable name="K_RANK_RANKL" public_interface="in" units="pM"/>
<variable name="Beta_OPG" public_interface="in" units="first_order_rate_constant"/>
<variable name="OPG_max" public_interface="in" units="pM"/>
<variable name="Deg_OPG" public_interface="in" units="first_order_rate_constant"/>
<variable name="Deg_OPG_RANKL" public_interface="in" units="first_order_rate_constant"/>
<variable name="K_OPG_RANKL" public_interface="in" units="pM"/>
<variable name="Beta_RANKL_OCY" public_interface="in" units="first_order_rate_constant"/>
<variable name="Beta_RANKL_OBp" public_interface="in" units="first_order_rate_constant"/>
<variable name="Deg_RANKL" public_interface="in" units="first_order_rate_constant"/>
<variable name="Deg_RANK_RANKL" public_interface="in" units="first_order_rate_constant"/>
<variable name="KD_RANKL_act_OCp" public_interface="in" units="pM"/>
<variable name="KD_RANKL_act_OCu" public_interface="in" units="pM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>RANK</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>N_RANK_OCp</ci>
<ci>OCp</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>K_RANK_RANKL</ci>
</apply>
<ci>RANKL</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>OPG</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Beta_OPG</ci>
<ci>Pi_PTH_rep</ci>
<ci>OBa</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>Beta_OPG</ci>
<ci>Pi_PTH_rep</ci>
<ci>OBa</ci>
</apply>
<ci>OPG_max</ci>
</apply>
<ci>Deg_OPG</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>Deg_OPG_RANKL</ci>
<ci>K_OPG_RANKL</ci>
</apply>
<ci>RANKL</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RANKL_tot</ci>
<apply>
<times/>
<ci>RANKL</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>RANK</ci>
<ci>K_RANK_RANKL</ci>
</apply>
<apply>
<divide/>
<ci>OPG</ci>
<ci>K_OPG_RANKL</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RANKL_max</ci>
<apply>
<times/>
<ci>N_RANKL_OBp_max</ci>
<ci>OBp_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_RANKL1</ci>
<apply>
<times/>
<ci>Beta_RANKL_OCY</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>RANKL_tot</ci>
<ci>RANKL_max</ci>
</apply>
</apply>
<ci>OCY_act</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_RANKL2</ci>
<apply>
<times/>
<ci>Beta_RANKL_OBp</ci>
<ci>Pi_NO_PTH_act_rep</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>RANKL_tot</ci>
<ci>RANKL_max</ci>
</apply>
</apply>
<ci>OBp</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_RANKL</ci>
<apply>
<plus/>
<ci>P_RANKL1</ci>
<ci>P_RANKL2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>RANKL</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>P_RANKL</ci>
<ci>P_RANKL_d</ci>
</apply>
<apply>
<plus/>
<ci>Deg_RANKL</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>Deg_RANK_RANKL</ci>
<ci>K_RANK_RANKL</ci>
</apply>
<ci>RANK</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>Deg_OPG_RANKL</ci>
<ci>K_OPG_RANKL</ci>
</apply>
<ci>OPG</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_RANKL_act_OCp</ci>
<apply>
<divide/>
<ci>RANKL</ci>
<apply>
<plus/>
<ci>KD_RANKL_act_OCp</ci>
<ci>RANKL</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_RANKL_act_OCu</ci>
<apply>
<divide/>
<ci>RANKL</ci>
<apply>
<plus/>
<ci>KD_RANKL_act_OCu</ci>
<ci>RANKL</ci>
</apply>
</apply>
</apply>
</math>
</component>
<!-- TGF-beta-->
<component name="TGF_beta">
<variable name="TGF_beta" units="pM"/>
<variable name="TGF_beta_t0" units="pM"/>
<variable name="OCa" public_interface="in" units="pM"/>
<variable initial_value="0.0001" name="OCa_t0" units="pM"/>
<variable name="Alpha" public_interface="in" units="dimensionless"/>
<variable name="Pi_TGFbeta_OBu_act_t0" public_interface="out" units="dimensionless"/>
<variable name="KD_TGF_beta_activate" public_interface="in" units="pM"/>
<variable name="KD_TGF_beta_repress" public_interface="in" units="pM"/>
<variable name="Pi_TGFbeta_OBu_act" units="dimensionless"/>
<variable name="Pi_TGFbeta_OCa_act" public_interface="out" units="dimensionless"/>
<variable name="Pi_TGFbeta_OBp_rep" public_interface="out" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>TGF_beta_t0</ci>
<apply>
<times/>
<ci>Alpha</ci>
<ci>OCa_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_TGFbeta_OBu_act_t0</ci>
<apply>
<divide/>
<ci>TGF_beta_t0</ci>
<apply>
<plus/>
<ci>KD_TGF_beta_activate</ci>
<ci>TGF_beta_t0</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>TGF_beta</ci>
<apply>
<times/>
<ci>Alpha</ci>
<ci>OCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_TGFbeta_OBu_act</ci>
<apply>
<divide/>
<ci>TGF_beta</ci>
<apply>
<plus/>
<ci>KD_TGF_beta_activate</ci>
<ci>TGF_beta</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_TGFbeta_OCa_act</ci>
<apply>
<divide/>
<ci>TGF_beta</ci>
<apply>
<plus/>
<ci>KD_TGF_beta_activate</ci>
<ci>TGF_beta</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_TGFbeta_OBp_rep</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>TGF_beta</ci>
<ci>KD_TGF_beta_repress</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- Scl-->
<component name="Scl">
<variable name="Scl" public_interface="out" units="pM"/>
<variable name="Scl_0" public_interface="out" units="pM"/>
<variable name="P_Scl_b" units="flux"/>
<variable name="P_Scl_b_0" units="flux"/>
<variable name="A" units="second_order_rate_constant"/>
<variable name="A_0" units="second_order_rate_constant"/>
<variable name="B" units="first_order_rate_constant"/>
<variable name="B_0" units="first_order_rate_constant"/>
<variable name="C" units="flux"/>
<variable name="C_0" units="flux"/>
<variable name="OCYprod" units="flux"/>
<variable name="OCYprod_0" units="flux"/>
<variable name="OBp" public_interface="in" units="pM"/>
<variable name="OBp_t0" public_interface="in" units="pM"/>
<variable name="Wnt" public_interface="in" units="pM"/>
<variable name="Wnt_0" public_interface="in" units="pM"/>
<variable name="OCY_act" public_interface="in" units="pM"/>
<variable name="OCY_act_0" public_interface="in" units="pM"/>
<variable name="Pi_eps_rep" public_interface="in" units="dimensionless"/>
<variable name="Pi_eps_rep_stst" public_interface="in" units="dimensionless"/>
<variable name="Beta_Scl" units="first_order_rate_constant"/>
<variable name="Beta_Scl_0" public_interface="in" units="first_order_rate_constant"/>
<variable name="KD_SclLRP5" public_interface="in" units="pM"/>
<variable name="Deg_Scl" public_interface="in" units="first_order_rate_constant"/>
<variable name="Scl_max" public_interface="in" units="pM"/>
<variable name="KD_WntLRP5" public_interface="in" units="pM"/>
<variable name="Deg_SclLRP5" public_interface="in" units="first_order_rate_constant"/>
<variable name="LRP5perCell" public_interface="in" units="dimensionless"/>
<variable name="LRP5tot" units="pM"/>
<variable name="LRP5tot_0" units="pM"/>
<variable name="P_Scl_d" public_interface="in" units="flux"/>
<!-- Scl-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Beta_Scl</ci>
<ci>Beta_Scl_0</ci>
</apply>
<apply>
<eq/>
<ci>OCYprod</ci>
<apply>
<times/>
<ci>Beta_Scl</ci>
<ci>OCY_act</ci>
<ci>Pi_eps_rep</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_Scl_b</ci>
<apply>
<times/>
<ci>OCYprod</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Scl</ci>
<ci>Scl_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>LRP5tot</ci>
<apply>
<times/>
<ci>LRP5perCell</ci>
<ci>OBp</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>A</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>KD_SclLRP5</ci>
</apply>
<apply>
<plus/>
<ci>Deg_Scl</ci>
<apply>
<divide/>
<ci>OCYprod</ci>
<ci>Scl_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>B</ci>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>A</ci>
<ci>KD_SclLRP5</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt</ci>
<ci>KD_WntLRP5</ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>Deg_SclLRP5</ci>
<ci>LRP5tot</ci>
</apply>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<ci>P_Scl_d</ci>
<ci>OCYprod</ci>
</apply>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>C</ci>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<ci>P_Scl_d</ci>
<ci>OCYprod</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt</ci>
<ci>KD_WntLRP5</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Scl</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>B</ci>
</apply>
<apply>
<root/>
<apply>
<minus/>
<apply>
<power/>
<ci>B</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>A</ci>
<ci>C</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>A</ci>
</apply>
</apply>
</apply>
<!-- Scl_0-->
<apply>
<eq/>
<ci>OCYprod_0</ci>
<apply>
<times/>
<ci>Beta_Scl_0</ci>
<ci>OCY_act_0</ci>
<ci>Pi_eps_rep_stst</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_Scl_b_0</ci>
<apply>
<times/>
<ci>OCYprod_0</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Scl_0</ci>
<ci>Scl_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>LRP5tot_0</ci>
<apply>
<times/>
<ci>LRP5perCell</ci>
<ci>OBp_t0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>A_0</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>KD_SclLRP5</ci>
</apply>
<apply>
<plus/>
<ci>Deg_Scl</ci>
<apply>
<divide/>
<ci>OCYprod_0</ci>
<ci>Scl_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>B_0</ci>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>A_0</ci>
<ci>KD_SclLRP5</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt_0</ci>
<ci>KD_WntLRP5</ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>Deg_SclLRP5</ci>
<ci>LRP5tot_0</ci>
</apply>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<ci>P_Scl_d</ci>
<ci>OCYprod_0</ci>
</apply>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>C_0</ci>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<ci>P_Scl_d</ci>
<ci>OCYprod_0</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt_0</ci>
<ci>KD_WntLRP5</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Scl_0</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>B_0</ci>
</apply>
<apply>
<root/>
<apply>
<minus/>
<apply>
<power/>
<ci>B_0</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>A_0</ci>
<ci>C_0</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>A_0</ci>
</apply>
</apply>
</apply>
</math>
</component>
<!-- Wnt-->
<component name="Wnt">
<variable name="Wnt" public_interface="out" units="pM"/>
<variable name="Wnt_0" public_interface="in" units="pM"/>
<variable name="Scl" public_interface="in" units="pM"/>
<variable name="Scl_0" public_interface="in" units="pM"/>
<variable name="Pi_WNT" public_interface="out" units="dimensionless"/>
<variable name="Pi_WNT_0" public_interface="out" units="dimensionless"/>
<variable name="KD_WntLRP5" public_interface="in" units="pM"/>
<variable name="KD_SclLRP5" public_interface="in" units="pM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Wnt</ci>
<ci>Wnt_0</ci>
</apply>
<apply>
<eq/>
<ci>Pi_WNT</ci>
<apply>
<divide/>
<ci>Wnt</ci>
<apply>
<times/>
<ci>KD_WntLRP5</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt</ci>
<ci>KD_WntLRP5</ci>
</apply>
<apply>
<divide/>
<ci>Scl</ci>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_WNT_0</ci>
<apply>
<divide/>
<ci>Wnt_0</ci>
<apply>
<times/>
<ci>KD_WntLRP5</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Wnt_0</ci>
<ci>KD_WntLRP5</ci>
</apply>
<apply>
<divide/>
<ci>Scl_0</ci>
<ci>KD_SclLRP5</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- NO and PTH-->
<component name="NO_PTH">
<variable name="Pi_NO_PTH_act_rep" public_interface="out" units="dimensionless"/>
<variable name="PTH_tot" units="pM"/>
<variable name="Beta_PTH" public_interface="in" units="flux"/>
<variable name="P_PTH_d" public_interface="in" units="flux"/>
<variable name="Deg_PTH" public_interface="in" units="first_order_rate_constant"/>
<variable name="NO_tot" units="pM"/>
<variable name="NO_eq" units="pM"/>
<variable name="Beta_NO" units="first_order_rate_constant"/>
<variable name="Beta_NO_0" public_interface="in" units="first_order_rate_constant"/>
<variable name="P_NO_d" public_interface="in" units="flux"/>
<variable name="Deg_NO" public_interface="in" units="first_order_rate_constant"/>
<variable name="Pi_eps_act" public_interface="in" units="dimensionless"/>
<variable name="Pi_eps_act_stst" public_interface="in" units="dimensionless"/>
<variable name="OCY_act" public_interface="in" units="pM"/>
<variable name="OCY_act_0" public_interface="in" units="pM"/>
<variable name="NO_max" public_interface="in" units="pM"/>
<variable name="KD_rep_NO" units="pM"/>
<variable initial_value="2" name="aa" units="dimensionless"/>
<variable name="Pi_PTH_act" units="dimensionless"/>
<variable name="KD_PTH_act" public_interface="in" units="pM"/>
<variable name="Pi_NO_rep" units="dimensionless"/>
<variable name="lambda_s" public_interface="in" units="dimensionless"/>
<variable name="lambda_c" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>PTH_tot</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>Beta_PTH</ci>
<ci>P_PTH_d</ci>
</apply>
<ci>Deg_PTH</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Beta_NO</ci>
<ci>Beta_NO_0</ci>
</apply>
<apply>
<eq/>
<ci>NO_tot</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<ci>Beta_NO</ci>
<ci>Pi_eps_act</ci>
<ci>OCY_act</ci>
</apply>
<ci>P_NO_d</ci>
</apply>
<apply>
<plus/>
<ci>Deg_NO</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Beta_NO</ci>
<ci>Pi_eps_act</ci>
<ci>OCY_act</ci>
</apply>
<ci>NO_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>NO_eq</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Beta_NO_0</ci>
<ci>Pi_eps_act_stst</ci>
<ci>OCY_act_0</ci>
</apply>
<apply>
<plus/>
<ci>Deg_NO</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Beta_NO_0</ci>
<ci>Pi_eps_act_stst</ci>
<ci>OCY_act_0</ci>
</apply>
<ci>NO_max</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>KD_rep_NO</ci>
<apply>
<divide/>
<ci>NO_eq</ci>
<ci>aa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_PTH_act</ci>
<apply>
<divide/>
<ci>PTH_tot</ci>
<apply>
<plus/>
<ci>PTH_tot</ci>
<ci>KD_PTH_act</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_NO_rep</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>NO_tot</ci>
<ci>KD_rep_NO</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_NO_PTH_act_rep</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>lambda_s</ci>
<apply>
<plus/>
<ci>Pi_PTH_act</ci>
<ci>Pi_NO_rep</ci>
</apply>
</apply>
<apply>
<times/>
<ci>lambda_c</ci>
<ci>Pi_PTH_act</ci>
<ci>Pi_NO_rep</ci>
</apply>
</apply>
</apply>
</math>
</component>
<!-- Strain energy density-->
<component name="SED">
<variable name="SED" units="MPa"/>
<variable name="SED_bm" units="MPa"/>
<variable name="Pi_eps_act" public_interface="out" units="dimensionless"/>
<variable name="Pi_eps_rep" public_interface="out" units="dimensionless"/>
<variable name="Pi_eps_act_stst" public_interface="out" units="dimensionless"/>
<variable name="Pi_eps_rep_stst" public_interface="out" units="dimensionless"/>
<variable name="sig_macro_t0" public_interface="in" units="MPa"/>
<variable name="de_sig_macro" public_interface="in" units="MPa"/>
<variable name="sig_macro" units="MPa"/>
<variable name="sig_1" units="MPa"/>
<variable name="sig_2" units="MPa"/>
<variable name="sig_3" units="MPa"/>
<variable name="sig_4" units="MPa"/>
<variable name="sig_5" units="MPa"/>
<variable name="sig_6" units="MPa"/>
<variable name="eps_1" units="dimensionless"/>
<variable name="eps_2" units="dimensionless"/>
<variable name="eps_3" units="dimensionless"/>
<variable name="eps_4" units="dimensionless"/>
<variable name="eps_5" units="dimensionless"/>
<variable name="eps_6" units="dimensionless"/>
<variable name="dens_tis" units="g_per_cm3"/>
<variable name="vm" public_interface="in" units="dimensionless"/>
<variable name="vo" public_interface="in" units="dimensionless"/>
<variable name="rho_m" public_interface="in" units="g_per_cm3"/>
<variable name="rho_o" public_interface="in" units="g_per_cm3"/>
<variable name="fbm" public_interface="in" units="dimensionless"/>
<variable name="fvas" units="dimensionless"/>
<variable name="dens" units="g_per_cm3"/>
<variable name="dens_1" units="dimensionless"/>
<variable name="E_mod" units="MPa"/>
<variable name="nu" public_interface="in" units="dimensionless"/>
<variable name="S_11" units="per_MPa"/>
<variable name="S_12" units="per_MPa"/>
<variable name="S_44" units="per_MPa"/>
<variable name="omega" public_interface="in" units="dimensionless"/>
<variable name="valeq" units="dimensionless"/>
<variable name="taueq" public_interface="in" units="MPa"/>
<variable name="alphaAct" public_interface="in" units="dimensionless"/>
<variable name="alphaRep" public_interface="in" units="dimensionless"/>
<variable name="gammaAct" public_interface="in" units="dimensionless"/>
<variable name="gammaRep" public_interface="in" units="dimensionless"/>
<variable name="rhoAct" public_interface="in" units="dimensionless"/>
<variable name="rhoRep" public_interface="in" units="dimensionless"/>
<variable name="deltaAct" units="MPa"/>
<variable name="deltaRep" units="MPa"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>sig_1</ci>
<cn cellml:units="MPa">0.0</cn>
</apply>
<apply>
<eq/>
<ci>sig_2</ci>
<cn cellml:units="MPa">0.0</cn>
</apply>
<!--sig_3 = 0.0{MPa};-->
<apply>
<eq/>
<ci>sig_4</ci>
<cn cellml:units="MPa">0.0</cn>
</apply>
<apply>
<eq/>
<ci>sig_5</ci>
<cn cellml:units="MPa">0.0</cn>
</apply>
<apply>
<eq/>
<ci>sig_6</ci>
<cn cellml:units="MPa">0.0</cn>
</apply>
<apply>
<eq/>
<ci>sig_macro</ci>
<apply>
<plus/>
<ci>sig_macro_t0</ci>
<ci>de_sig_macro</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sig_3</ci>
<ci>sig_macro</ci>
</apply>
<apply>
<eq/>
<ci>dens_tis</ci>
<apply>
<plus/>
<cn cellml:units="g_per_cm3">1</cn>
<apply>
<times/>
<apply>
<minus/>
<ci>rho_m</ci>
<cn cellml:units="g_per_cm3">1</cn>
</apply>
<ci>vm</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>rho_o</ci>
<cn cellml:units="g_per_cm3">1</cn>
</apply>
<ci>vo</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>fvas</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>fbm</ci>
<cn cellml:units="dimensionless">100</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>dens</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>dens_tis</ci>
<ci>fbm</ci>
</apply>
<cn cellml:units="dimensionless">100</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>dens_1</ci>
<apply>
<divide/>
<ci>dens</ci>
<cn cellml:units="g_per_cm3">1</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>E_mod</ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="MPa">2014</cn>
<apply>
<power/>
<ci>dens_1</ci>
<cn cellml:units="dimensionless">2.5</cn>
</apply>
</apply>
<apply>
<lt/>
<ci>dens_1</ci>
<cn cellml:units="dimensionless">1.2</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<cn cellml:units="MPa">1763</cn>
<apply>
<power/>
<ci>dens_1</ci>
<cn cellml:units="dimensionless">3.2</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<!-- Compute the components of complaince tensor
! Bone is assumed as an isotropic material, so zero elements are ignored here,
and only the independent components are calculated-->
<apply>
<eq/>
<ci>S_11</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>E_mod</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>S_12</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>nu</ci>
</apply>
<ci>E_mod</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>S_44</ci>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="dimensionless">2</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>nu</ci>
</apply>
</apply>
<ci>E_mod</ci>
</apply>
</apply>
<!-- Calculate macroscopic strain using Hooke's law-->
<apply>
<eq/>
<ci>eps_1</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>S_11</ci>
<ci>sig_1</ci>
</apply>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_2</ci>
</apply>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>eps_2</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_1</ci>
</apply>
<apply>
<times/>
<ci>S_11</ci>
<ci>sig_2</ci>
</apply>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>eps_3</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_1</ci>
</apply>
<apply>
<times/>
<ci>S_12</ci>
<ci>sig_2</ci>
</apply>
<apply>
<times/>
<ci>S_11</ci>
<ci>sig_3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>eps_4</ci>
<apply>
<times/>
<ci>S_44</ci>
<ci>sig_4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>eps_5</ci>
<apply>
<times/>
<ci>S_44</ci>
<ci>sig_5</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>eps_6</ci>
<apply>
<times/>
<ci>S_44</ci>
<ci>sig_6</ci>
</apply>
</apply>
<!-- Calculate strain energy density at the continuum level-->
<apply>
<eq/>
<ci>SED</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<plus/>
<apply>
<times/>
<ci>sig_1</ci>
<ci>eps_1</ci>
</apply>
<apply>
<times/>
<ci>sig_2</ci>
<ci>eps_2</ci>
</apply>
<apply>
<times/>
<ci>sig_3</ci>
<ci>eps_3</ci>
</apply>
<apply>
<times/>
<ci>sig_4</ci>
<ci>eps_4</ci>
</apply>
<apply>
<times/>
<ci>sig_5</ci>
<ci>eps_5</ci>
</apply>
<apply>
<times/>
<ci>sig_6</ci>
<ci>eps_6</ci>
</apply>
</apply>
</apply>
</apply>
<!-- Calculate strain energy density in bone matrix, at the tissue level-->
<apply>
<eq/>
<ci>SED_bm</ci>
<apply>
<divide/>
<ci>SED</ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fvas</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<!-- Mechancial feedback activator function that accounts for the production of NO by osteocytes-->
<apply>
<eq/>
<ci>valeq</ci>
<ci>omega</ci>
</apply>
<apply>
<eq/>
<ci>deltaAct</ci>
<apply>
<times/>
<ci>taueq</ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>alphaAct</ci>
<ci>valeq</ci>
</apply>
<apply>
<minus/>
<ci>valeq</ci>
<ci>rhoAct</ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>gammaAct</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_eps_act</ci>
<apply>
<plus/>
<ci>rhoAct</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>alphaAct</ci>
<ci>rhoAct</ci>
</apply>
<apply>
<power/>
<ci>SED_bm</ci>
<ci>gammaAct</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>deltaAct</ci>
<ci>gammaAct</ci>
</apply>
<apply>
<power/>
<ci>SED_bm</ci>
<ci>gammaAct</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_eps_act_stst</ci>
<apply>
<plus/>
<ci>rhoAct</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>alphaAct</ci>
<ci>rhoAct</ci>
</apply>
<apply>
<power/>
<ci>taueq</ci>
<ci>gammaAct</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>deltaAct</ci>
<ci>gammaAct</ci>
</apply>
<apply>
<power/>
<ci>taueq</ci>
<ci>gammaAct</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>deltaRep</ci>
<apply>
<times/>
<ci>taueq</ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<minus/>
<ci>valeq</ci>
<ci>rhoRep</ci>
</apply>
<apply>
<minus/>
<ci>alphaRep</ci>
<ci>valeq</ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>gammaRep</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_eps_rep</ci>
<apply>
<minus/>
<ci>alphaRep</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>alphaRep</ci>
<ci>rhoRep</ci>
</apply>
<apply>
<power/>
<ci>SED_bm</ci>
<ci>gammaRep</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>deltaRep</ci>
<ci>gammaRep</ci>
</apply>
<apply>
<power/>
<ci>SED_bm</ci>
<ci>gammaRep</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pi_eps_rep_stst</ci>
<apply>
<minus/>
<ci>alphaRep</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>alphaRep</ci>
<ci>rhoRep</ci>
</apply>
<apply>
<power/>
<ci>taueq</ci>
<ci>gammaRep</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>deltaRep</ci>
<ci>gammaRep</ci>
</apply>
<apply>
<power/>
<ci>taueq</ci>
<ci>gammaRep</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-----------------------------------------------------------------
Model Parameters
--------------------------------------------------------------- -->
<component name="model_parameters">
<!-- Cell-->
<variable initial_value="0.01" name="OBu_t0" public_interface="out" units="pM"/>
<variable initial_value="0.01" name="OCu_t0" public_interface="out" units="pM"/>
<!-- OBp-->
<variable initial_value="0.185" name="D_OBp_t0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0.166" name="D_OBu_t0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0" name="pd_OBp_t0" public_interface="out" units="flux"/>
<!--OBa-->
<variable initial_value="0" name="pd_OBa_t0" public_interface="out" units="flux"/>
<variable initial_value="0.0212" name="D_OBa_t0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0.1908" name="A_OBa_t0" public_interface="out" units="first_order_rate_constant"/>
<!-- OCp-->
<variable initial_value="0.0219" name="D_OCu_t0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0" name="pd_OCp_t0" public_interface="out" units="flux"/>
<variable initial_value="0.01958" name="D_OCp_t0" public_interface="out" units="first_order_rate_constant"/>
<!-- OCa-->
<variable initial_value="0" name="pd_OCa_t0" public_interface="out" units="flux"/>
<variable initial_value="10" name="A_OCa_t0" public_interface="out" units="first_order_rate_constant"/>
<!-- Bone matrix fraction-->
<variable initial_value="50" name="K_form" public_interface="out" units="second_order_rate_constant"/>
<variable initial_value="2500" name="K_res" public_interface="out" units="second_order_rate_constant"/>
<!-- Osteocytes-->
<variable initial_value="0.000414138" name="fact_0" public_interface="out" units="pM"/>
<!-- Mineralization-->
<variable initial_value="0.007" name="XKAPPA" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0.516" name="vmmax" public_interface="out" units="dimensionless"/>
<!-- PTH-->
<variable initial_value="250" name="Beta_PTH" public_interface="out" units="flux"/>
<variable initial_value="0" name="P_PTH_d" public_interface="out" units="flux"/>
<variable initial_value="86" name="Deg_PTH" public_interface="out" units="first_order_rate_constant"/>
<!--var KD_PTH_act: pM {init: 0.899382, pub: out};-->
<variable initial_value="0.65" name="KD_PTH_act" public_interface="out" units="pM"/>
<variable initial_value="0.222581427709954" name="KD_PTH_rep" public_interface="out" units="pM"/>
<!-- NO_PTH -->
<variable initial_value="3440" name="Beta_NO_0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0" name="P_NO_d" public_interface="out" units="flux"/>
<variable initial_value="0.0021" name="Deg_NO" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="200000000" name="NO_max" public_interface="out" units="pM"/>
<variable initial_value="0.450454" name="lambda_s" public_interface="out" units="dimensionless"/>
<variable initial_value="0.900909" name="lambda_c" public_interface="out" units="dimensionless"/>
<!-- RANK-RANKL-OPG-->
<variable initial_value="2703476" name="N_RANKL_OBp_max" public_interface="out" units="dimensionless"/>
<variable initial_value="4160" name="N_RANK_OCp" public_interface="out" units="dimensionless"/>
<variable initial_value="10" name="K_RANK_RANKL" public_interface="out" units="pM"/>
<variable initial_value="162490033.783568" name="Beta_OPG" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="131.428571428571" name="OPG_max" public_interface="out" units="pM"/>
<variable initial_value="532608.695652174" name="Deg_OPG" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="10.132471014805" name="Deg_OPG_RANKL" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0.0151142857142857" name="K_OPG_RANKL" public_interface="out" units="pM"/>
<variable initial_value="5660" name="Beta_RANKL_OCY" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="23600" name="Beta_RANKL_OBp" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="10.132471014805" name="Deg_RANKL" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="10.132471014805" name="Deg_RANK_RANKL" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="3.34" name="KD_RANKL_act_OCp" public_interface="out" units="pM"/>
<variable initial_value="16.7" name="KD_RANKL_act_OCu" public_interface="out" units="pM"/>
<variable initial_value="0" name="P_RANKL_d" public_interface="out" units="flux"/>
<!-- TGF-beta-->
<variable initial_value="1" name="Alpha" public_interface="out" units="dimensionless"/>
<variable initial_value="0.000563278809675429" name="KD_TGF_beta_activate" public_interface="out" units="pM"/>
<variable initial_value="0.000175426051821094" name="KD_TGF_beta_repress" public_interface="out" units="pM"/>
<!-- Scl-->
<variable initial_value="24000" name="Beta_Scl_0" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="10" name="KD_SclLRP5" public_interface="out" units="pM"/>
<variable initial_value="1" name="Deg_Scl" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="70" name="Scl_max" public_interface="out" units="pM"/>
<variable initial_value="1000" name="KD_WntLRP5" public_interface="out" units="pM"/>
<variable initial_value="50" name="Deg_SclLRP5" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="5" name="LRP5perCell" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="P_Scl_d" public_interface="out" units="flux"/>
<!-- Wnt-->
<variable initial_value="170" name="Wnt_0" public_interface="out" units="pM"/>
<!-- Strain energy density-->
<variable initial_value="-0.1457" name="sig_macro_t0" public_interface="out" units="MPa"/>
<variable initial_value="0" name="de_sig_macro" public_interface="out" units="MPa"/>
<variable initial_value="0.428571428571428" name="vo" public_interface="out" units="dimensionless"/>
<variable initial_value="3.2" name="rho_m" public_interface="out" units="g_per_cm3"/>
<variable initial_value="1.41" name="rho_o" public_interface="out" units="g_per_cm3"/>
<variable initial_value="0.3" name="nu" public_interface="out" units="dimensionless"/>
<variable initial_value="0.95" name="omega" public_interface="out" units="dimensionless"/>
<variable initial_value="0.006652" name="taueq" public_interface="out" units="MPa"/>
<variable initial_value="1" name="alphaAct" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alphaRep" public_interface="out" units="dimensionless"/>
<variable initial_value="7" name="gammaAct" public_interface="out" units="dimensionless"/>
<variable initial_value="8.01559" name="gammaRep" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="rhoAct" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="rhoRep" public_interface="out" units="dimensionless"/>
</component>
<!-----------------------------------------------------------------
Mappings between components for variables that are shared
--------------------------------------------------------------- -->
<connection>
<map_components component_1="OBp" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="model_parameters"/>
<map_variables variable_1="OBu_t0" variable_2="OBu_t0"/>
<map_variables variable_1="D_OBp_t0" variable_2="D_OBp_t0"/>
<map_variables variable_1="D_OBu_t0" variable_2="D_OBu_t0"/>
<map_variables variable_1="pd_OBp_t0" variable_2="pd_OBp_t0"/>
<map_variables variable_1="KD_TGF_beta_repress" variable_2="KD_TGF_beta_repress"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="OBa"/>
<map_variables variable_1="OBp_out" variable_2="OBp_out"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="RANK_RANKL_OPG"/>
<map_variables variable_1="OBp" variable_2="OBp"/>
<map_variables variable_1="OBp_t0" variable_2="OBp_t0"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="TGF_beta"/>
<map_variables variable_1="Pi_TGFbeta_OBu_act_t0" variable_2="Pi_TGFbeta_OBu_act_t0"/>
<map_variables variable_1="Pi_TGFbeta_OBp_rep" variable_2="Pi_TGFbeta_OBp_rep"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="Wnt"/>
<map_variables variable_1="Pi_WNT" variable_2="Pi_WNT"/>
<map_variables variable_1="Pi_WNT_0" variable_2="Pi_WNT_0"/>
</connection>
<connection>
<map_components component_1="OBp" component_2="Scl"/>
<map_variables variable_1="OBp" variable_2="OBp"/>
<map_variables variable_1="OBp_t0" variable_2="OBp_t0"/>
</connection>
<connection>
<map_components component_1="OBa" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="OBa" component_2="model_parameters"/>
<map_variables variable_1="pd_OBa_t0" variable_2="pd_OBa_t0"/>
<map_variables variable_1="D_OBa_t0" variable_2="D_OBa_t0"/>
<map_variables variable_1="A_OBa_t0" variable_2="A_OBa_t0"/>
</connection>
<connection>
<map_components component_1="OBa" component_2="fbm"/>
<map_variables variable_1="OBa" variable_2="OBa"/>
</connection>
<connection>
<map_components component_1="OBa" component_2="vm"/>
<map_variables variable_1="OBa" variable_2="OBa"/>
</connection>
<connection>
<map_components component_1="OBa" component_2="RANK_RANKL_OPG"/>
<map_variables variable_1="OBa" variable_2="OBa"/>
</connection>
<connection>
<map_components component_1="OCp" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="OCp" component_2="model_parameters"/>
<map_variables variable_1="OCu_t0" variable_2="OCu_t0"/>
<map_variables variable_1="D_OCu_t0" variable_2="D_OCu_t0"/>
<map_variables variable_1="pd_OCp_t0" variable_2="pd_OCp_t0"/>
<map_variables variable_1="D_OCp_t0" variable_2="D_OCp_t0"/>
</connection>
<connection>
<map_components component_1="OCp" component_2="OCa"/>
<map_variables variable_1="OCp_out" variable_2="OCp_out"/>
</connection>
<connection>
<map_components component_1="OCp" component_2="RANK_RANKL_OPG"/>
<map_variables variable_1="OCp" variable_2="OCp"/>
<map_variables variable_1="Pi_RANKL_act_OCp" variable_2="Pi_RANKL_act_OCp"/>
<map_variables variable_1="Pi_RANKL_act_OCu" variable_2="Pi_RANKL_act_OCu"/>
</connection>
<connection>
<map_components component_1="OCa" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="OCa" component_2="model_parameters"/>
<map_variables variable_1="pd_OCa_t0" variable_2="pd_OCa_t0"/>
<map_variables variable_1="A_OCa_t0" variable_2="A_OCa_t0"/>
</connection>
<connection>
<map_components component_1="OCa" component_2="fbm"/>
<map_variables variable_1="OCa" variable_2="OCa"/>
</connection>
<connection>
<map_components component_1="OCa" component_2="TGF_beta"/>
<map_variables variable_1="OCa" variable_2="OCa"/>
<map_variables variable_1="Pi_TGFbeta_OCa_act" variable_2="Pi_TGFbeta_OCa_act"/>
</connection>
<connection>
<map_components component_1="fbm" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="fbm" component_2="model_parameters"/>
<map_variables variable_1="K_form" variable_2="K_form"/>
<map_variables variable_1="K_res" variable_2="K_res"/>
</connection>
<connection>
<map_components component_1="fbm" component_2="OCY"/>
<map_variables variable_1="fbm" variable_2="fbm"/>
<map_variables variable_1="dfbmdt" variable_2="dfbmdt"/>
</connection>
<connection>
<map_components component_1="fbm" component_2="vm"/>
<map_variables variable_1="fbm" variable_2="fbm"/>
</connection>
<connection>
<map_components component_1="fbm" component_2="SED"/>
<map_variables variable_1="fbm" variable_2="fbm"/>
</connection>
<connection>
<map_components component_1="OCY" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="OCY" component_2="model_parameters"/>
<map_variables variable_1="fact_0" variable_2="fact_0"/>
</connection>
<connection>
<map_components component_1="OCY" component_2="NO_PTH"/>
<map_variables variable_1="OCY_act" variable_2="OCY_act"/>
<map_variables variable_1="OCY_act_0" variable_2="OCY_act_0"/>
</connection>
<connection>
<map_components component_1="OCY" component_2="RANK_RANKL_OPG"/>
<map_variables variable_1="OCY_act" variable_2="OCY_act"/>
</connection>
<connection>
<map_components component_1="OCY" component_2="Scl"/>
<map_variables variable_1="OCY_act" variable_2="OCY_act"/>
<map_variables variable_1="OCY_act_0" variable_2="OCY_act_0"/>
</connection>
<connection>
<map_components component_1="vm" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="vm" component_2="model_parameters"/>
<map_variables variable_1="XKAPPA" variable_2="XKAPPA"/>
<map_variables variable_1="vmmax" variable_2="vmmax"/>
<map_variables variable_1="K_form" variable_2="K_form"/>
</connection>
<connection>
<map_components component_1="vm" component_2="SED"/>
<map_variables variable_1="vm" variable_2="vm"/>
</connection>
<connection>
<map_components component_1="PTH" component_2="model_parameters"/>
<map_variables variable_1="Beta_PTH" variable_2="Beta_PTH"/>
<map_variables variable_1="P_PTH_d" variable_2="P_PTH_d"/>
<map_variables variable_1="Deg_PTH" variable_2="Deg_PTH"/>
<map_variables variable_1="KD_PTH_act" variable_2="KD_PTH_act"/>
<map_variables variable_1="KD_PTH_rep" variable_2="KD_PTH_rep"/>
</connection>
<connection>
<map_components component_1="PTH" component_2="RANK_RANKL_OPG"/>
<map_variables variable_1="Pi_PTH_rep" variable_2="Pi_PTH_rep"/>
</connection>
<connection>
<map_components component_1="NO_PTH" component_2="SED"/>
<map_variables variable_1="Pi_eps_act" variable_2="Pi_eps_act"/>
<map_variables variable_1="Pi_eps_act_stst" variable_2="Pi_eps_act_stst"/>
</connection>
<connection>
<map_components component_1="RANK_RANKL_OPG" component_2="model_parameters"/>
<map_variables variable_1="N_RANKL_OBp_max" variable_2="N_RANKL_OBp_max"/>
<map_variables variable_1="N_RANK_OCp" variable_2="N_RANK_OCp"/>
<map_variables variable_1="K_RANK_RANKL" variable_2="K_RANK_RANKL"/>
<map_variables variable_1="Beta_OPG" variable_2="Beta_OPG"/>
<map_variables variable_1="OPG_max" variable_2="OPG_max"/>
<map_variables variable_1="Deg_OPG" variable_2="Deg_OPG"/>
<map_variables variable_1="Deg_OPG_RANKL" variable_2="Deg_OPG_RANKL"/>
<map_variables variable_1="K_OPG_RANKL" variable_2="K_OPG_RANKL"/>
<map_variables variable_1="Beta_RANKL_OCY" variable_2="Beta_RANKL_OCY"/>
<map_variables variable_1="Beta_RANKL_OBp" variable_2="Beta_RANKL_OBp"/>
<map_variables variable_1="Deg_RANKL" variable_2="Deg_RANKL"/>
<map_variables variable_1="Deg_RANK_RANKL" variable_2="Deg_RANK_RANKL"/>
<map_variables variable_1="KD_RANKL_act_OCp" variable_2="KD_RANKL_act_OCp"/>
<map_variables variable_1="KD_RANKL_act_OCu" variable_2="KD_RANKL_act_OCu"/>
<map_variables variable_1="P_RANKL_d" variable_2="P_RANKL_d"/>
</connection>
<connection>
<map_components component_1="RANK_RANKL_OPG" component_2="NO_PTH"/>
<map_variables variable_1="Pi_NO_PTH_act_rep" variable_2="Pi_NO_PTH_act_rep"/>
</connection>
<connection>
<map_components component_1="TGF_beta" component_2="model_parameters"/>
<map_variables variable_1="Alpha" variable_2="Alpha"/>
<map_variables variable_1="KD_TGF_beta_activate" variable_2="KD_TGF_beta_activate"/>
<map_variables variable_1="KD_TGF_beta_repress" variable_2="KD_TGF_beta_repress"/>
</connection>
<connection>
<map_components component_1="Scl" component_2="model_parameters"/>
<map_variables variable_1="Beta_Scl_0" variable_2="Beta_Scl_0"/>
<map_variables variable_1="KD_SclLRP5" variable_2="KD_SclLRP5"/>
<map_variables variable_1="Deg_Scl" variable_2="Deg_Scl"/>
<map_variables variable_1="Scl_max" variable_2="Scl_max"/>
<map_variables variable_1="KD_WntLRP5" variable_2="KD_WntLRP5"/>
<map_variables variable_1="Deg_SclLRP5" variable_2="Deg_SclLRP5"/>
<map_variables variable_1="LRP5perCell" variable_2="LRP5perCell"/>
<map_variables variable_1="P_Scl_d" variable_2="P_Scl_d"/>
<map_variables variable_1="Wnt_0" variable_2="Wnt_0"/>
</connection>
<connection>
<map_components component_1="Scl" component_2="Wnt"/>
<map_variables variable_1="Wnt" variable_2="Wnt"/>
<map_variables variable_1="Scl" variable_2="Scl"/>
<map_variables variable_1="Scl_0" variable_2="Scl_0"/>
</connection>
<connection>
<map_components component_1="Scl" component_2="SED"/>
<map_variables variable_1="Pi_eps_rep" variable_2="Pi_eps_rep"/>
<map_variables variable_1="Pi_eps_rep_stst" variable_2="Pi_eps_rep_stst"/>
</connection>
<connection>
<map_components component_1="Wnt" component_2="model_parameters"/>
<map_variables variable_1="Wnt_0" variable_2="Wnt_0"/>
<map_variables variable_1="KD_WntLRP5" variable_2="KD_WntLRP5"/>
<map_variables variable_1="KD_SclLRP5" variable_2="KD_SclLRP5"/>
</connection>
<connection>
<map_components component_1="NO_PTH" component_2="model_parameters"/>
<map_variables variable_1="Beta_PTH" variable_2="Beta_PTH"/>
<map_variables variable_1="Beta_NO_0" variable_2="Beta_NO_0"/>
<map_variables variable_1="P_PTH_d" variable_2="P_PTH_d"/>
<map_variables variable_1="Deg_PTH" variable_2="Deg_PTH"/>
<map_variables variable_1="P_NO_d" variable_2="P_NO_d"/>
<map_variables variable_1="Deg_NO" variable_2="Deg_NO"/>
<map_variables variable_1="NO_max" variable_2="NO_max"/>
<map_variables variable_1="KD_PTH_act" variable_2="KD_PTH_act"/>
<map_variables variable_1="lambda_s" variable_2="lambda_s"/>
<map_variables variable_1="lambda_c" variable_2="lambda_c"/>
</connection>
<connection>
<map_components component_1="SED" component_2="model_parameters"/>
<map_variables variable_1="sig_macro_t0" variable_2="sig_macro_t0"/>
<map_variables variable_1="de_sig_macro" variable_2="de_sig_macro"/>
<map_variables variable_1="vo" variable_2="vo"/>
<map_variables variable_1="rho_m" variable_2="rho_m"/>
<map_variables variable_1="rho_o" variable_2="rho_o"/>
<map_variables variable_1="nu" variable_2="nu"/>
<map_variables variable_1="omega" variable_2="omega"/>
<map_variables variable_1="taueq" variable_2="taueq"/>
<map_variables variable_1="alphaAct" variable_2="alphaAct"/>
<map_variables variable_1="alphaRep" variable_2="alphaRep"/>
<map_variables variable_1="gammaAct" variable_2="gammaAct"/>
<map_variables variable_1="gammaRep" variable_2="gammaRep"/>
<map_variables variable_1="rhoAct" variable_2="rhoAct"/>
<map_variables variable_1="rhoRep" variable_2="rhoRep"/>
</connection>
</model>
