SI3 modular 1D-2D interface for thermodynamics

Development description

This is a large (and long) development that is required to 1) use SI3 thermodynamics in other models (e.g. nextsim) and 2) be compatible with GPU architecture

Code implementation

• Make icethd.F90 fully modular 1D-2D

  • extract DO loops from each of the subroutines
  • make a big DO loop inside icethd
  • extract the big DO loop from icethd and put it in icestp.F90

• Build an interface for passing the fields between thermo and the rest of the code

  • Make a list of input and output fields
  • [ ]
  • [ ]