Setting the glacier rheology and flow properties

I used the NonNewtonian built in rheology in GALE. Right now, I have the temperature dependence turned off and everything set for temperate ice. Eventually, using the temperature dependence to set the flow parameter (A) could be easily implemented, I think.

<struct name="nonNewtonian">
<param name="Type">NonNewtonian</param> <!-- εdot = Aτ^n exp(−T0/T) -->
<param name="StrainRateInvariantField">StrainRateInvariantField</param>
<param name="TemperatureField">TemperatureField</param>
<param name="n">3</param>
<param name="T_0">0.0</param> <!-- having this be 0 should decouple the temperature -->
<param name="A">6.8e-24</param> <!-- Temperate ice deformation constant; see Physics of Ice -->
<param name="refStrainRate">0.01</param>
</struct>

<struct name="viscous">
<param name="Type">RheologyMaterial</param>
<param name="Shape">iceShape</param> <!-- shape defined for glacier -->
<param name="density">917</param> <!-- ice density -->

<param name="alpha">1.0</param> <!-- leaving this as default since temp is decoupled -->
<param name="diffusivity">1.0</param> <!-- leaving this as default since temp is decoupled -->
<list name="heatingElements"> <!-- leaving this as default since temp is decoupled -->
<struct>
<param name="Q">1.0</param>
<param name="lambda">1.0</param>
</struct>
</list>

<list name="Rheology">
<param>nonNewtonian</param>
<param>storeViscosity</param>
<param>storeStress</param>
</list>
</struct>