tice_high in init_vertical_profile for ktherm = 2

[zhaobin74]

Hi All,

This is more of a question. In init_vertical_profile, logical tice_high is set to true if a given ice layer k has zTin above Tmlt.

Icepack/columnphysics/icepack_therm_vertical.F90

Lines 903 to 904 in 86cae16

	if (zTin(k) > Tmax) then
	tice_high = .true.

If mushy layer thermodynamic scheme is used, zTinis corrected and a warning is issued.

Icepack/columnphysics/icepack_therm_vertical.F90

Lines 933 to 941 in 86cae16

	if (ktherm == 2) then
	zqin(k) = enthalpy_of_melting(zSin(k)) - c1
	zTin(k) = icepack_mushy_temperature_mush(zqin(k),zSin(k))
	write(warnstr,*) subname, 'Corrected quantities'
	call icepack_warnings_add(warnstr)
	write(warnstr,*) subname, 'zqin=',zqin(k)
	call icepack_warnings_add(warnstr)
	write(warnstr,*) subname, 'zTin=',zTin(k)
	call icepack_warnings_add(warnstr)

The current logic will also correct all layers below k, regardless of whether those layers have higher zTin than Tmlt, since tice_high stays True within the k-loop. My question is: is this correction of lower layers intended to be that way?

Thank you.

[apcraig]

You could also try asking this question here, https://bb.cgd.ucar.edu/cesm/forums/cice-consortium.146/

[dabail10]

Thanks for posting on the forum. I do think this might be a bug. Unless the intention was to reset the whole column regardless? I will investigate.

[eclare108213]

@akturner, in the current mushy code, when a layer temperature is too high, it is corrected but tice_high is not reset, so all temperatures below that layer are also recalculated. Is this intended, or should we be resetting tice_high=.false. for lower layers? Thanks.