compute heat flux at mixed layer depth for CICE standalone

[dpath2o]

https://github.com/dpath2o/AFIM/blob/87aaed86398d90489607aaadea6bb4ef1fb0739b/src/python/afim.py#L267

Comments from Siobhan: The stand alone CICE model is set up with a fixed depth mixed layer, ther nCAR folk ran it with their q-flux mixed layer version which did have variable mixed layers which came form their CESM runs, when we ran qflux run in CSIRO MK2 (and MK1 but we never called it back then) we had fixed depths.

In q-flux runs the mixed layers are constraining both the advective transport in the 1 d set up by also changing the vertical depths they NCAR approach was an attempt to allow for entrainment/detrainment as well. But they need all the components of a coupled model to set those data up and closing the heat budget is hard we tried to replicate unsuccessfully with ACCESS 1-0 data back in 2012/13.

The reason we talked about the Alek Petty-CPOM paper was that was a better approach for a mixed layer model under CICE but for Dan’s fast ice study a fixed mixed layer will be sufficient as the focus is on fast ice.

Comments from Dan: There is not a ‘default’ way to have ‘dynamic’ MLD in CICE standalone – i.e. an MLD that is fed to CICE per grid cell per time step. I do see how this can be done in Fortran code, but for Antarctic FI it is likely not as important as the atmosphere.

[dpath2o]

this issue is not going to be pursued within the context of my PhD

[dpath2o]

this issue is not going to be pursued within the context of my PhD