We have been using a new MLD definition (delta_rho calculated from 10m depth, not the first model layer). This requires setting `HREF_FOR_MLD = 10.0`