melt rate is assumed to depend quadratically on the thermal forcing</a>:</p><pre><code class="language-julia hljs">melt_rate = QuadraticMeltRate();</code></pre><p>Finally, build the ice sheet model and pass the melt rate model</p><pre><code class="language-julia hljs">model = Model(grid = grid, bed_elevation = bed_elevation, melt_rate = melt_rate)</code></pre><h2 id="Analytic-Melt-Rate-Parametrizations"><a class="docs-heading-anchor" href="#Analytic-Melt-Rate-Parametrizations">Analytic Melt Rate Parametriza