From 32782e38ef9cf6d770672d1e316a85c992a76c97 Mon Sep 17 00:00:00 2001
From: Brian Groenke <brian.groenke@awi.de>
Date: Thu, 14 Nov 2024 19:13:15 +0100
Subject: [PATCH] Minor adjustment to advective heat flux

---
 src/Physics/Heat/heat_water.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/Physics/Heat/heat_water.jl b/src/Physics/Heat/heat_water.jl
index 5c513182..a187f4f9 100644
--- a/src/Physics/Heat/heat_water.jl
+++ b/src/Physics/Heat/heat_water.jl
@@ -10,7 +10,7 @@ WaterHeatBC(waterbc::WaterBC, heatbc::HeatBC) = Coupled(waterbc, heatbc)
 Computes the advective energy flux between grid cells with temperatures `T₁` and `T₂`
 given the heat capacity of water `cw` and latent heat of fusion `L`.
 """
-advectiveflux(jw, T₁, T₂, cw, L) = jw*(cw*T₁*(jw > zero(jw)) + cw*T₂*(jw < zero(jw)) + L)
+advectiveflux(jw, T₁, T₂, cw, L) = jw*(cw*(T₁-T₂) + L)
 
 """
     water_energy_advection!(::SubSurface, ::Coupled(WaterBalance, HeatBalance), state)
-- 
GitLab