Thermodynamics

[skygering]

Equations between Subzero and Mukund's model are very different. Brandon said something about not really being able to call Subzero model's value a heatflux. Make sure that these refer to the same thing for Oceananigans/Subzero calculations for coupling.

[skygering]

Current thermodynamics are very simple and could be improved.

They assume both that the ocean is at a freezing temperature, and also that the atmosphere temperature is less than the ocean temperature as all grid points.

For now we will add in warnings if these conditions are not met, but in the future we should improve the thermodynamics to work in other situations.

[skygering]

See hflx_factor for both the ocean and individual floes to get a sense of where thermodynamics are right now. They are a direct copy from original Subzero code.