public subroutine do_stim_winton(ice_cover, dz, dt, Tw, S, Ta, precip, Qsw, Qfluxes)
This SUBROUTINE updates the sea ice prognostic variables. The updated
variables are upper ice layer temperature (T1), lower ice layer temperature
(T2), snow thickness (hs), and ice thickness (hi).
The ice model performs this in two steps. First the temperatures are updated
and secondly the changes in ice and snow thickness are calculated.
Any surplus energy that is not used for melting is returned in tmelt and
bmelt.
Evaporation and bottom ablation formation are not included in
this version of the model. Furthermore we do not keep an explicit water
and salt budget for the sea ice and how that affects the water and salt
budgets in the ocean.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.