update figures

Author: danieldouglas92

cookbooks/2D_subduction_with_two_phase_flow/doc/2D_subduction_two_phase_flow.md

@@ -38,14 +38,20 @@ where $i$ denote the deformation mechanism (diffusion or dislocation creep), $n_
 
 In these models, the water density (1000 kg/m$^3$) is significantly lower than the solid density (3300 kg/m$^3$), so the fluid velocity is dominantly vertical due to the high buoyancy force experienced by the fluid. However, corner flow in the mantle wedge does impose a trench-ward horizontal component to the fluid velocity. As the fluid ascends through the hot peridotite mantle wedge, the PT conditions do not allow for the free fluid to be reabsorbed into the solid phase. However, when the free fluid starts ascending through the cooler overriding plate, the temperature is sufficiently low that hydration of the overriding plate begins to occur. This leads to a local reduction of the viscosity. Within the model where the fluid is advected with the Darcy velocity, two distinct bands of free water can be seen seeping out of the subducting plate. The first band produces a larger flux of water, is centered around a distance of ~175 km landward of the trench, and is sourced from the subducting peridotite layer. The second band is lower in flux magnitude, is centered around a distance of ~400 km landward from the trench, and is sourced from the gabbro and MORB layers.
 
-```{figure-md} fig:final-water
-<img src="final_water.png" />
+```{figure-md} fig:fluid-pathways
+<img src="fluid_figures.png" />
 
  The viscosity distribution at the end of the model run, with the free fluid overlain on top. White contour shows where the bound water makes up 0.1 wt% of the solid.
 ```
 
-```{figure-md} fig:initial-viscosity
-<img src="slab_viscosity.png" />
+```{figure-md} fig:viscosity-comparison
+<img src="viscosity_profile.png" />
+
+ The viscosity of the model within the immediate vicinity of the subduction zone. The vectors show the direction of and are scaled by the velocity of the solid phase.
+```
+
+```{figure-md} fig:velocity-temperature
+<img src="solid_velocity.png" />
 
  The viscosity of the model within the immediate vicinity of the subduction zone. The vectors show the direction of and are scaled by the velocity of the solid phase.
 ```

cookbooks/2D_subduction_with_two_phase_flow/doc/bound_water.png

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