Olympic Mountains Glacier Flowline Model

This MATLAB program is designed to simulate the flow dynamics of glaciers within the Elwha and Quinault basins, situated in the Olympic mountains. Using parameters that affect glacial flow, such as the density of ice, the gravitational pull, and the coefficients for deformation and sliding of the glacier, this model provides an insight into the process of glacial evolution.

The model initializes by importing the desired glacier profile from two available .mat files, \"Elwha_Basin_Input.mat\" or \"Quinault_Basin_Input.mat,\" which provide comprehensive data about the unique conditions of the Elwha and Quinault basins, respectively. The geometric attributes of the glacier, including its bed elevation and width, are then meticulously defined and interpolated onto a calculated grid for accurate modelling.

Climate parameters integral to the functioning of glaciers are also incorporated into this model, encompassing annual precipitation rates, sea-level temperature, and the melt factor. These values, along with other parameters, feed into the critical computation of the mass balance of the glacier, which balances the accumulation against the melting of the glacier ice.

At the heart of the simulation is a time-stepping loop that updates and records the glacier's height, thickness, and the rate of change of thickness over the defined period. This iterative approach provides dynamic results that demonstrate the evolution of the glacier over time, providing outputs of both glacier height and thickness at each time step.

The generated outputs from this model provide an invaluable resource for understanding glacial behavior and the effects of varying climatic conditions on glaciers. This can be of great use in broader studies of climate change and its impact on glacial ecosystems, providing a data-driven foundation for further research and policy making.