  ylabel = &quot;y (km)&quot;,
                        colorbar_title = &quot;melt rate (m/yr)&quot;,
                        title = key,
                        framestyle = &quot;box&quot;)
    xlims!((420, 640))
    plot!(size = (500,300))

    # save the figure
    savefig
    display(plt) #uncomment to show in (e.g.) VSCode
end</code></pre><center><img src="https://raw.githubusercontent.com/RJArthern/WAVI.jl/build-docs/docs/src/assets/example-plots//melt_parametrizations//quadratic.png" alt="" title="