Hi all,

I have been working on plume-lithosphere interaction in 3d models. We are using the Visco plastic material model and we consider both diffusion and dislocation creep in our models. However, the output/postprocessor can only plot the total viscosity in Paraview. We want to plot a diff/disl creep figure in Paraview shown below in the upper panel.

https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2019JB018560

I have compared the Postprocess set up in this paper where the figure comes from with mine:

**Setup in the 2019 JGR paper:**

```
subsection Postprocess
set List of postprocessors = visualization, velocity statistics, basic statistics, depth average
.
.
.
subsection Visualization
set List of output variables = viscosity, density, strain rate, gamma # default
```

**Below is my setup:**

```
subsection Postprocess
set List of postprocessors = visualization, temperature statistics, depth average, topography
.
.
.
subsection Visualization
set List of output variables = material properties, vertical heat flux, strain rate, stress, stress second invariant
subsection Material properties
set List of material properties = density, viscosity, thermal conductivity, thermal expansivity, specific heat, thermal diffusivity
end
```

The difference is that the 2019 JGR paper output the strain rate. However, even if I output the strain rate, I am not sure how to plot the dislocation creep viscosity in Paraview. Do I plot the dislocation creep first with the following equation and then plot the diff/disl using (η-η_disl)/η_disl?

η_disl=1/2*A_disl^(-1/n)*(ε ̇_disl)^((1-n)/n)*exp((E_disl*+pV_disl*)/(nRT))

I would greatly appreciate it if someone could help me figure out a way to plot the diff/disl in Paraview. Thanks!

Best,

Ziqi