Investigate thermo conservation errors

[anton-seaice]

In the current esm1.6 preindustrial configuration, we get these conservation errors:

 --------------------------------
   CICE model diagnostic output  
 --------------------------------
  
 Thermo energy conservation error
 istep1, my_task, i, j:    13894481          10           9         294
 Flux error (W/m^2) =  1.111111111111111E-003
 Energy error (J) =   4.00000000000000     
 Initial energy =  -177658543.217718     
 Final energy   = -2.846865814515282E+016
 efinal - einit  = -2.846865796749428E+016
 fsurfn,flatn,fswint,fhocn, fsnow*Lfresh:
  -7907960546527.84       -1.75764449185750       0.000000000000000E+000
  2.526374171591467E-016  0.102906905954250     
 Input energy = -2.846865796749428E+016
 fbot(i,j),fcondbot(ij):
  0.000000000000000E+000  -7907960546527.84     
 fcondtop_solve(i,j), fcondtopn_extra(i,j):
  0.000000000000000E+000  0.000000000000000E+000
 enum(ij):  0.000000000000000E+000
 Intermediate energy =  -177658543.217718     
 efinal - einter = -2.846865796749428E+016
 einter - einit  =  0.000000000000000E+000
 Conduction Error =  0.000000000000000E+000
 Melt/Growth Error =   4.00000000000000     
 Advection Error =  0.000000000000000E+000
 dt*(fsurfn, flatn, fswint, fhocn, fsnow*Lfresh, fadvocn):
 -2.846865796750024E+016  -6327.52017068699       0.000000000000000E+000
  9.094947017729282E-013   370.464861435300       0.000000000000000E+000
 istep1, my_task, iblk =    13894481          10           1
 category n =            1
 Global block:          11
 Global i and j:         308         293
 Lat, Lon:   80.7662741933158        61.1715717767987     
 aice:  0.999987721910044     
 n:            1 aicen:   1.085986601001615E-011

This are allowed to continue because the limits are increased:

cice5/source/ice_therm_vertical.F90

Lines 2560 to 2565 in 5828e0c

	if (ferr > 10.0 * ferrmax) then
	l_stop = .true.
	write(nu_diag,*) 'BBB: TOO BAD --- CICE is to stop!'
	else
	write(nu_diag,*) 'BBB: WARNING -- too big flux error --'
	endif