Advanced 2D-Subduction cookbook with reactive fluid transport

[danieldouglas92]

Implement a cookbook which demonstrates reactive fluid transport in a 2D subduction model. The cookbook will have multiple different setups:

• A cookbook which keeps the temperature, and all compositions (except for the porosity) as static fields. The bound fluid will not be allowed to decrease, acting as a constant fluid source within a stationary slab that is initialized with the Geodynamic WorldBuilder. The fluid will advect according to the Darcy advection scheme. This will be a step above in complexity from the cookbook Simplified Subduction Model with Parameterized Solid-Fluid Reactions because there is no velocities being prescribed within the slab or mantle wedge.

• Repeat the above with viscous weakening from the water fugacity term implemented here to investigate the impact of this coupling

• Repeat the above two steps with the fluid being advected according to the fully coupled McKenzie equations.

• A dynamic 2D subduction model that begins with a small subducted segment and runs for ~10-20 Myr. The fluid would be advected with the Darcy velocity.

• Repeat the above step with viscous weakening from the water fugacity

• Repeat the above two steps with the fluid advected according to the fully coupled McKenzie equations.

[jdannberg]

Just some initial thoughts:
You made the checklist with all the cases you want to add; would it make sense to split this into two cookbooks? One being instantaneous and one being time-dependent? Just to make the files and the pull requests smaller.

I know this is WIP, but the images look really cool already!

[naliboff]

@danieldouglas92 -I second @jdannberg that although this is a WIP the initial results are really cool and this will make a great cookbook. Nice work! I've left some initial comments and suggestions to consider as you work on expanding out the cookbook.

[naliboff]

@danieldouglas92 - A few more minor comments from my end, but I think is otherwise very close. Please let me know when it is ready for a final review!

[danieldouglas92]

The text is ready for another look, I just need to update the figures (same figure ideas, just new model output)!

I was having issues with convergence with the fully coupled models, so my thoughts are to try to merge what is currently here, which is just the results and discussion of the darcy advection case, and then add the output and discussion for the mckenzie advection case in a follow up PR. Does that sound reasonable @jdannberg @naliboff?

[naliboff]

@danieldouglas92 - I just read through the cookbook again and it is in great shape.

A few questions and suggestions for how to proceed:

1. Would you like to merge this cookbook in pieces, or wait until all the components outlined in the first comment are ready?

2. If you prefer to merge it in pieces (that would be my preference), maybe make an associated issue that can keep track of individual components and associated PRs?

3. Again, if you prefer to merge in pieces, then this PR would just need an change log, ideally a test case, and perhaps squashing some of the commits. Other than that, this is in great shape and nice work!

[danieldouglas92]

@naliboff I think merging the different extensions of the cookbooks as follow up PR's makes the most sense. I can definitely make a tracking issue for that. I've squashed all the commits and added a change log, I'll add a test case shortly.

[danieldouglas92]

No rush, but this is ready for another look! @naliboff @jdannberg