From 583d42004bd3a03e74c7deee15f6c5cec37d7d7f Mon Sep 17 00:00:00 2001
From: Victor Stepanenko <stepanen@srcc.msu.ru>
Date: Fri, 26 Feb 2021 16:09:11 +0300
Subject: [PATCH] Outflow discharge ensured positive in advection calculations
---
source/model/trib.f90 | 1 +
1 file changed, 1 insertion(+)
diff --git a/source/model/trib.f90 b/source/model/trib.f90
index 0516339..0e76961 100644
--- a/source/model/trib.f90
+++ b/source/model/trib.f90
@@ -278,6 +278,7 @@ SUBROUTINE TRIBTEMP(time,dt,h1,dhwtrib,z_full,area_int,area_half,gsp,gas,wst,spi
wst%wArea(i-1) = wst%wArea(i) - gsp%ddz05(i-1)*h1*(inflow - outflow) !Solution of horizontally averaged continuity
!print*, i, wst%wArea(i-1), wst%wArea(i), gsp%ddz05(i),h1,inflow , outflow !Solution of horizontally averaged continuity
outflow = outflow/area_int(i)
+ outflow = max(outflow,0.) !Ensuring positivity of outflow discharge, used for advection only
!print*, 'outflow', outflow, U_tribin_(1,i),T_tribin_(1,i),width_tribin_(1,i), area_int(i), invdt
meanTin = 0.
--
GitLab