fix the .prm file

Author: danieldouglas92

cookbooks/2D_subduction_with_two_phase_flow/fixed_slab.prm

@@ -3,7 +3,8 @@
 # layer. The layers dehydrate as they are subducted, with the solid-fluid
 # reactions governed by an approximation published by Tian et al., 2018.
 
-set World builder file                         = fixed_slab.wb
+set Additional shared libraries                = /Users/danieldouglas/software/aspect/cookbooks/2D_subduction_with_two_phase_flow/plugin/lib2D_subduction_with_two_phase_flow.so
+set World builder file                         = /Users/danieldouglas/software/aspect/cookbooks/2D_subduction_with_two_phase_flow/fixed_slab.wb
 set Adiabatic surface temperature              = 1573
 set Nonlinear solver scheme                    = iterated Advection and Stokes
 set Nonlinear solver failure strategy          = continue with next timestep
@@ -58,7 +59,7 @@ subsection Temperature field
 end
 
 subsection Initial composition model
-  set List of model names = world builder
+  set List of model names = world builder disequilibrium
   subsection World builder
     set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment
   end
@@ -117,7 +118,7 @@ subsection Material model
     set Reference fluid density                          = 1000
     set Shear to bulk viscosity ratio                    = 0.1
     set Reference fluid viscosity                        = 1
-    set Reference permeability                           = 1e-7
+    set Reference permeability                           = 5e-8
     set Exponential fluid weakening factor               = 30
     set Fluid compressibility                            = 0
     set Fluid reaction time scale for operator splitting = 1e4
@@ -157,8 +158,6 @@ subsection Material model
                                                             MORB:        1.577e-21, \
                                                             sediment:    1.577e-21
 
-    set Query bound fluid for water fugacity              = true
-
     set Water fugacity exponents for diffusion creep      = 0.7
 
     set Stress exponents for diffusion creep              = 1.0
@@ -232,10 +231,10 @@ end
 subsection Mesh refinement
   set Coarsening fraction                      = 0.3
   set Refinement fraction                      = 0.5
-  set Initial adaptive refinement              = 0
-  set Initial global refinement                = 0
+  set Initial adaptive refinement              = 3
+  set Initial global refinement                = 5
   set Strategy                                 = temperature, viscosity, composition threshold, minimum refinement function
-  set Time steps between mesh refinement       = 5
+  set Time steps between mesh refinement       = 0
   set Refinement criteria scaling factors      = 2, 1.5, 1.5, 1
   set Refinement criteria merge operation      = max
   set Run postprocessors on initial refinement = false
@@ -261,11 +260,11 @@ subsection Boundary velocity model
 end
 
 subsection Melt settings
-  set Include melt transport                  = false
+  set Include melt transport = false
 end
 
 subsection Postprocess
-  set List of postprocessors = visualization, composition statistics, velocity statistics
+  set List of postprocessors          = visualization, composition statistics, velocity statistics
   subsection Visualization
     set List of output variables      = material properties, strain rate, named additional outputs, darcy velocity
     set Output format                 = vtu

cookbooks/2D_subduction_with_two_phase_flow/fixed_slab.wb

@@ -35,10 +35,10 @@
       "max depth": 300e3,
       "composition models": 
       [
-        {"model": "tian water content mod", "compositions": [1], "min depth":0.0, "max depth": 10e3, "lithology": "sediment", "initial water content": 2, "cutoff pressure": 1},
-        {"model": "tian water content mod", "compositions": [1], "min depth":10e3, "max depth": 20e3, "lithology": "MORB", "initial water content": 1, "cutoff pressure": 16},
-        {"model": "tian water content mod", "compositions": [1], "min depth":20e3, "max depth": 30e3, "lithology": "gabbro", "initial water content": 0.5, "cutoff pressure": 26},
-        {"model": "tian water content mod", "compositions": [1], "min depth":30e3, "max depth": 120e3, "lithology": "peridotite", "initial water content": 1, "cutoff pressure": 10},
+        {"model": "tian water content", "compositions": [1], "min depth":0.0, "max depth": 10e3, "lithology": "sediment", "initial water content": 2, "cutoff pressure": 1},
+        {"model": "tian water content", "compositions": [1], "min depth":10e3, "max depth": 20e3, "lithology": "MORB", "initial water content": 1, "cutoff pressure": 16},
+        {"model": "tian water content", "compositions": [1], "min depth":20e3, "max depth": 30e3, "lithology": "gabbro", "initial water content": 0.5, "cutoff pressure": 26},
+        {"model": "tian water content", "compositions": [1], "min depth":30e3, "max depth": 120e3, "lithology": "peridotite", "initial water content": 1, "cutoff pressure": 10},
         {"model": "uniform", "compositions": [5], "min depth":0.0, "max depth": 10e3, "operation":"add"},
         {"model": "uniform", "compositions": [4], "min depth":10e3, "max depth": 20e3, "operation":"add"},
         {"model": "uniform", "compositions": [3], "min depth":20e3, "max depth": 30e3, "operation":"add"},
@@ -58,10 +58,10 @@
                         {"length":300e3,"thickness":[300e3],"top truncation":[-50e3],"angle":[45,45]}],
 
 	    "composition models":[
-        {"model":"tian water content mod", "compositions":[1], "min distance slab top":0, "max distance slab top":10e3, "lithology":"sediment", "initial water content": 2, "cutoff pressure": 1},
-        {"model":"tian water content mod", "compositions":[1], "min distance slab top":10e3, "max distance slab top":20e3, "lithology":"MORB", "initial water content": 1, "cutoff pressure": 16},
-        {"model":"tian water content mod", "compositions":[1], "min distance slab top":20e3, "max distance slab top":30e3, "lithology":"gabbro", "initial water content": 0.5, "cutoff pressure": 26},
-        {"model":"tian water content mod", "compositions":[1], "min distance slab top":30e3, "max distance slab top":120e3, "lithology":"peridotite", "initial water content": 1, "cutoff pressure": 10},
+        {"model":"tian water content", "compositions":[1], "min distance slab top":0, "max distance slab top":10e3, "lithology":"sediment", "initial water content": 2, "cutoff pressure": 1},
+        {"model":"tian water content", "compositions":[1], "min distance slab top":10e3, "max distance slab top":20e3, "lithology":"MORB", "initial water content": 1, "cutoff pressure": 16},
+        {"model":"tian water content", "compositions":[1], "min distance slab top":20e3, "max distance slab top":30e3, "lithology":"gabbro", "initial water content": 0.5, "cutoff pressure": 26},
+        {"model":"tian water content", "compositions":[1], "min distance slab top":30e3, "max distance slab top":120e3, "lithology":"peridotite", "initial water content": 1, "cutoff pressure": 10},
         {"model":"uniform", "compositions":[5], "min distance slab top":0.0, "max distance slab top":10e3, "operation":"add"},
         {"model":"uniform", "compositions":[4], "min distance slab top":10e3, "max distance slab top":20e3, "operation":"add"},
         {"model":"uniform", "compositions":[3], "min distance slab top":20e3, "max distance slab top":30e3, "operation":"add"},

cookbooks/2D_subduction_with_two_phase_flow/plugin/bound_fluid_source.cc

@@ -44,7 +44,7 @@ namespace aspect
              MaterialModel::MaterialModelOutputs<dim> &out) const
     {
       ReactiveFluidTransport<dim>::evaluate(in, out);
-      const unsigned int bound_fluid_idx = this->introspection().compositional_index_for_name("bound_fluid_idx");
+      const unsigned int bound_fluid_idx = this->introspection().compositional_index_for_name("bound_fluid");
       for (unsigned int q=0; q < in.n_evaluation_points(); ++q)
         if (out.reaction_terms[q][bound_fluid_idx] <= 0.0)
           out.reaction_terms[q][bound_fluid_idx] = 0.0;

cookbooks/2D_subduction_with_two_phase_flow/plugin/world_builder_disequilibrium_fluid.cc

@@ -128,7 +128,7 @@ namespace aspect
                                                          n_comp);
           if (n_comp == bound_fluid_idx &&
               this->get_parameters().compositional_field_methods[bound_fluid_idx] == Parameters<dim>::AdvectionFieldMethod::static_field)
-            comp_value *= 1.1;
+            comp_value *= 1.1/100;
           return comp_value;
         }