This is a link to the ipython notebook *.ipynb
Geometry and model parm summary
Set up a vertical channel and made the lower crust very weak (low visc) and hot (asth temp) so that it would be sure to detach.
Vertical & horizontal detachment
Vertical channel is set to vc=1.95 (so that is .05 wide or 132 km wide!)
No plastic rheology. Only viscosity variation with depth.
# Ideal continental crust: Jelly sandwich
AsthViscosity = 0.01
lowerMantleViscosity = 100.0
UCrustViscosity = 100.0
LCrustViscosity = 0.01
LithViscosity = 10.0
CUCViscosity = 100.0
CLCViscosity = 100.0
CLithViscosity = 100.0
Temp of lower crust set to 1. Top = 0 base =1.
I used a conductive lithosphere as a temperature initial condition defined by a variable heatflux between the craton (45mW/m2) and cordillera (60 mW/m2). The Cordillera Lithosphere is set to be 1300C at 120km, while the craton is set to be 1300C at 200km.
Density of the lith is greater than the asthenosphere. It wants to sink. Temperature dependent density is applied like in Currie et al.
- Implement Plastic rheology
- Implement strain weakening like shear stress example
- Figure out how to extract topography from models (see below)
- Investigate if need air layer in models
- Figure out min model parms to product delam and dripping like Adams stuff.
- Write up equations using mac (for the formulas)