Temperature-dependent flow-rate factor.

pf4d

Hi again,

I would like to apply a rheology law like that of equation (5) of Larour et al., (2012). (ideally including water content dependence).

It is confusing in this case why we need to set md.materials.rheology_B to something. I assume that If I set md.materials.rheology_law = "Arrhenius", that your thermo-echanical coupler will update the flow law based to temperature as suggested by

https://issm.jpl.nasa.gov/doxygen/chtml/Arrhenius_8cpp.html#a2a3966946116e1a72f79159bfd7a1f01

Can anyone explain how to properly set up an energy-dependent flow-rate factor?

mathieumorlighem

Hi Evan,

first, we gave up on Doxygen so it is not up to date (sorry).
If you run a transient with md.transient.isthermal = 1, then you can use a rheology law to update B. You can choose between:
- None
- BuddJacka
- Cuffey
- CuffeyTemperate
- Paterson
- Arrhenius
- LliboutryDuval

They are all implemented in shared/Elements/BuddJacka.cpp if you want to take a look. But that's only for the update. To initialize B you can use the same functions (lower case though) that are in src/m/materials.
Best,
Mathieu