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