From 36536e7e15918aec25360edd286e34c19aa93ae4 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 17 Dec 2021 13:25:31 +0100
Subject: [PATCH] Philippe 17/12/2021: added forced STOP in TURB_FLUC (ZUSTAR
 used but not set)

---
 src/MNH/subl_blowsnow.f90 | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/MNH/subl_blowsnow.f90 b/src/MNH/subl_blowsnow.f90
index ad0cf278f..96c6b9be4 100644
--- a/src/MNH/subl_blowsnow.f90
+++ b/src/MNH/subl_blowsnow.f90
@@ -696,6 +696,7 @@ ZFCRI = MAX(ZFCRI1,ZFCRI2)
 !*      2    Calculate variances of the horizontal and vertical velocity components
 !
 ZS0   = ZFCRI*PZZ/PVMOD
+STOP 'Bug in TURB_FLUC: ZUSTAR used but not set'
 ZSIGU = 4.77 *ZUSTAR**2/ (1+33*ZS0)**0.66666
 ZSIGV = 2.76 *ZUSTAR**2/ (1+9.5*ZS0)**0.66666
 ZSIGW = 1.31 *ZUSTAR**2/ (1+3.12*ZS0)**0.66666
-- 
GitLab