se it in the same role: to the set the mean melt rate on the shelf.</p><p>Finally, a binary file melt rate, in which the melt rate is read in from a binary file. First, we&#39;ll create such a file, which will set the melt rate to be uniform on the shelf.</p><pre><code class="language-julia hljs">isfloat = (h .&lt; -918.0/1028.0 .* bed.(grid.xxh, grid.yyh)) #indices of floating elements
m = zeros(nx,ny);
m[isfloat] .= 10.0; #set everywhere floating to 10m/a
folder = joinpath(@__DIR__, &quot;melt_rate_parame