The extended model code suite computes steady states of the spatially extended drainage model described in the accompanying paper,  
An analysis of instabilities and limit cycles in glacier-dammed reservoirs. The MATLAB code can be run from stability_grid_v2, which creates a grid of values of lake storage capacity V_p and water input q_in to the lake (defined as inputs q_vals and V_vals to stability_grid_v2). The remaining parameters of the model are contained in a structure 'parameters' supplied as an input to stability_grid_v2. For each combination of values V_p and q_in, stability_grid_v2 first computes a steady state (which is independent of V_p) using line_steady_v2, and then computes the spectrum of the linearized problem using network_stability_v5.
A parameter structure of the format used by stability_grid_v2 can be generated using parameter_generator_2018, which is called by default by stability_grid_v2 if the parameter structure is not supplied at input. parameter_generator_2018 produces a dimensionless parameter structure analogous to that output by parallel_parameter_generator, and the scales required to convert output back to dimensional variables are included in the output argument "scales" produced by parallel_parameter_generator.

The code for solving the time-dependent version of the spatially extended model can be found in the online supplementary material to Rada and Schoof (Channelized, distributed, and disconnected: subglacial drainage under a valley glacier in the Yukon, 2018, The Cryopshere)

(c) 2020 Christian Schoof

Distributed under the BSD-3 license, see https://opensource.org/licenses/BSD-3-Clause
