ve; in this example, we don&#39;t care about the velocity solve being completely accurate, so we&#39;ll do the minimum number of iterations in this process!)</p><pre><code class="language-julia hljs">for (key, melt) in melt_rates
    #instantiate the model with the appropriate melt rate
    local model = Model(grid = grid,
            bed_elevation = bed, 
            initial_conditions = initial_conditions,
            solver_params = SolverParams(maxiter_picard=1),
            melt_rate = melt);

    #bri