 <OpenGeoSysProjectDiff base_file="A2.prj">

   <add sel="/*/processes/process">

     <subtype>StressSaturation_StrainPressureTemperature</subtype>

   </add>

   <remove sel="/*/processes/process[1]/initial_stress" />

   <add sel="/*/processes/process[1]">

     <initial_stress type="effective">sigma0</initial_stress>

   </add>

   <remove sel="/*/processes/process/constitutive_relation[*]" />

   <add sel="/*/processes/process">

     <constitutive_relation id="0,1">

       <type>MFront</type>

       <behaviour>ThermoPoroElasticity</behaviour>

       <material_properties>

         <material_property name="YoungModulus" parameter="E" />

         <material_property name="PoissonRatio" parameter="nu" />

         <material_property name="ThermalExpansion" parameter="alpha" />

         <material_property name="BiotCoefficient" parameter="alpha_B" />

         <material_property name="BishopsExponent" parameter="m_chi" />

         <material_property name="ResidualLiquidSaturation" parameter="S_L_res" />

         <material_property name="ResidualGasSaturation" parameter="S_G_res" />

         <material_property name="BubblePressure" parameter="p_b" />

         <material_property name="VanGenuchtenExponent_m" parameter="m_S" />

       </material_properties>

     </constitutive_relation>

   </add>

   <add sel="/*/parameters">

     <parameter>

       <name>alpha</name>

       <type>Constant</type>

       <value>1e-5</value>

     </parameter>

     <parameter>

       <name>alpha_B</name>

       <type>Constant</type>

       <value>0.6</value>

     </parameter>

     <parameter>

       <name>m_chi</name>

       <type>Constant</type>

       <value>1</value>

     </parameter>

     <parameter>

       <name>S_L_res</name>

       <type>Constant</type>

       <value>0</value>

     </parameter>

     <parameter>

       <name>S_G_res</name>

       <type>Constant</type>

       <value>0</value>

     </parameter>

     <parameter>

       <name>p_b</name>

       <type>Constant</type>

       <value>42e6</value>

     </parameter>

     <parameter>

       <name>m_S</name>

       <type>Constant</type>

       <value>0.4</value>

     </parameter>

   </add>

   <remove sel="/*/time_loop/output/variables/variable[text()="sigma"]" />

   <add sel="/*/time_loop/output/variables">

     <variable>sigma_total</variable>

   </add>

   <remove sel="/*/processes/process/secondary_variables/secondary_variable[@name="sigma"]" />

   <add sel="/*/processes/process/secondary_variables">

     <secondary_variable name="sigma_total" />

   </add>

   <replace msel="/*/time_loop/output/prefix/text()">A2_effective_stess0_test</replace>

   <remove sel="/*/test_definition" />

 </OpenGeoSysProjectDiff>