From 31d4d0c52eedf53e24b65c69aa2fa27b59da6b05 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Riette?= <sebastien.riette@meteo.fr>
Date: Tue, 4 Oct 2022 14:33:04 +0200
Subject: [PATCH] S. Riette 4 oct 2022: bugfix for OCEAN key in
 mode_tke_eps_sources

---
 src/common/turb/mode_tke_eps_sources.F90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/common/turb/mode_tke_eps_sources.F90 b/src/common/turb/mode_tke_eps_sources.F90
index 2b9faf564..1e1e4bc40 100644
--- a/src/common/turb/mode_tke_eps_sources.F90
+++ b/src/common/turb/mode_tke_eps_sources.F90
@@ -266,7 +266,7 @@ END IF
 !*       2.2  Explicit TKE sources except horizontal turbulent transport 
 !
 ! extrapolate the dynamic production with a 1/Z law from its value at the 
-IF (LOCEAN) THEN
+IF (OOCEAN) THEN
   ! W(IKE) value stored in PDP(IKE) to the mass localization of tke(IKE)
   PDP(:,:,IKE) = PDP(:,:,IKE) * (1. + PDZZ(:,:,IKE)/PDZZ(:,:,IKE+1))
 ELSE
@@ -286,7 +286,7 @@ ZSOURCE(:,:,:) = ( PRTKES(:,:,:) +  PRTKEMS(:,:,:) ) / PRHODJ(:,:,:) &
 !
 ! To add here in ZSOURCE surface flux of TKE 
 !(assumed to be 0 for ATM, 
-IF (LOCEAN) THEN
+IF (OOCEAN) THEN
   !for ocean:wave breaking  simple/very rough param wE = 100 Ustar**3 where ustar is the Tau_atmi/rhocea  
   ZSOURCE (:,:,IKE)=ZSOURCE(:,:,IKE)-1.E2*((PSFUM(:,:)**2 + PSFVM(:,:)**2)**1.5) /PDZZ(:,:,IKE)
 END IF
-- 
GitLab