Currently we have a simple `ice_plot` function. Would like to enhance to change scaling, allow multiple plots, compare models etc