 Models&quot; section for information on all keyword parameters, here and for the below models.)</p><p>Next up: PICO. For this parametrization, we first have to make a mask defining where the ice front is, and then pass it as a keyword argument when we instantiate the <code>PICO</code> melt rate object.</p><pre><code class="language-julia hljs">ice_front_mask = zeros(nx,ny);
ice_front_mask[end,:] .= 1; #this tells julia that the ice front is at the downstream end of the domain
melt_PICO = PICO(ice_front_mas