obal initial_conditions = InitialConditions(initial_thickness = h);</code></pre><p>(We make this a global variable because we&#39;ll use it in a loop later.)</p><h2 id="Melt-Rate-Models"><a class="docs-heading-anchor" href="#Melt-Rate-Models">Melt Rate Models</a><a id="Melt-Rate-Models-1"></a><a class="docs-heading-anchor-permalink" href="#Melt-Rate-Models" title="Permalink"></a></h2><p>We&#39;ll loop over each of the melt rate models mentioned above. For each, we&#39;ll produce a map of the melt rate in th