do not correct zqsn in condensation step when ktherm=2 #3
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR fixed a bug of correcting snow enthalpy
zqsn(1)
prematurely whenktherm=2
.In 1440x1080 configuration, snow layer sometimes could get above freezing when
mushy layer
scheme is used. There is no check for such condition in the thermo solver. Instead,thickness_change
subroutine handles such case in an energy and mass conserving mannerIcepack/columnphysics/icepack_therm_vertical.F90
Lines 1314 to 1336 in 008f5f6
In GEOS coupling, when condensation happens over snow, there is a line to correct round-off error of updated top snow layer enthalpy
zqsn(1)
Icepack/columnphysics/icepack_therm_vertical.F90
Lines 1347 to 1348 in abbbd62
This correction is supposed to be for round-off, not for the above freezing case as seen in
mushy layer
scheme. Note: forBL99
, snow layer temperature never gets above freezing so the correction is valid here. Ifmushy layer
scheme applies this correction, energy is not conserved and the model will crash in conservation check. So anif
condition is added here to only apply correction forBL99
scheme. Inmushy layer
, the correction is done as above ⬆️.Zero-diff for AMIP and coupled runs using CICE4 sea ice model