From 90d3f90e88e7183265951454cbca2aeda593e5c6 Mon Sep 17 00:00:00 2001 From: Carolyn Begeman Date: Wed, 8 Jan 2025 10:26:11 -0600 Subject: [PATCH] Add temperature tendency due to mixing of frazil meltwater with seawater --- components/mpas-ocean/src/shared/mpas_ocn_frazil_forcing.F | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/components/mpas-ocean/src/shared/mpas_ocn_frazil_forcing.F b/components/mpas-ocean/src/shared/mpas_ocn_frazil_forcing.F index 5ca55901e3eb..5e363824c53d 100644 --- a/components/mpas-ocean/src/shared/mpas_ocn_frazil_forcing.F +++ b/components/mpas-ocean/src/shared/mpas_ocn_frazil_forcing.F @@ -599,6 +599,10 @@ subroutine ocn_frazil_forcing_build_arrays(domain, meshPool, forcingPool, stateP * config_frazil_ice_density ) / (config_specific_heat_sea_water * rho_sw) & / dt + ! ocean fluid temperature is modified due to mixing of existing water with meltwater at the freezing point + frazilTemperatureTendency(k,iCell) = frazilTemperatureTendency(k,iCell) + & + meltedFrazilIceThickness * oceanFreezingTemperature / dt + ! keep track of new frazil ice sumNewThicknessWeightedSaltContent = sumNewThicknessWeightedSaltContent - meltedThicknessWeightedSaltContent sumNewFrazilIceThickness = sumNewFrazilIceThickness - meltedFrazilIceThickness