From 38e71abb2aff76ce7b04bf9a3549b5cd7c1388bc Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Thu, 27 Jun 2024 12:06:39 +0200
Subject: [PATCH] Philippe 27/06/2024: PHYEX: mode_prandtl: ocean: comment out
 some unused computations

---
 src/PHYEX/turb/mode_prandtl.f90 | 48 ++++++++++++++++-----------------
 1 file changed, 24 insertions(+), 24 deletions(-)

diff --git a/src/PHYEX/turb/mode_prandtl.f90 b/src/PHYEX/turb/mode_prandtl.f90
index ec204cab3..cc1736374 100644
--- a/src/PHYEX/turb/mode_prandtl.f90
+++ b/src/PHYEX/turb/mode_prandtl.f90
@@ -530,18 +530,18 @@ DO JSV=1,KSV
   ELSE  IF (O2D) THEN ! 3D case in a 2D model
 !
     IF (OOCEAN) THEN
-      !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      ZWORK1(:,:) = (CST%XG *CST%XALPHAOC * PLM(:,:) * PLEPS(:,:) &
-                                       / PTKEM(:,:))**2
-      !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      CALL MZM_PHY(D,ZWORK1,ZWORK2)
-      IF (KRR /= 0) THEN
-        !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-        ZW1(:,:) = ZWORK2(:,:) * PETHETA(:,:)
-        !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      ELSE
-        ZW1 = ZWORK2
-      END IF
+      !!$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !ZWORK1(:,:) = (CST%XG *CST%XALPHAOC * PLM(:,:) * PLEPS(:,:) &
+      !                                 / PTKEM(:,:))**2
+      !!$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !CALL MZM_PHY(D,ZWORK1,ZWORK2)
+      !IF (KRR /= 0) THEN
+      !  !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !  ZW1(:,:) = ZWORK2(:,:) * PETHETA(:,:)
+      !  !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !ELSE
+      !  ZW1 = ZWORK2
+      !END IF
     ELSE
       !Compute only once and reuse in next JSV iterations
       IF ( JSV == 1 ) THEN
@@ -587,18 +587,18 @@ DO JSV=1,KSV
   ELSE ! 3D case in a 3D model
 !
     IF (OOCEAN) THEN
-      !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      ZWORK1(:,:) = (CST%XG *CST%XALPHAOC * PLM(:,:) * PLEPS(:,:) &
-                                       / PTKEM(:,:))**2
-      !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      CALL MZM_PHY(D,ZWORK1,ZWORK2)
-      IF (KRR /= 0) THEN
-        !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-        ZW1(:,:) = ZWORK2(:,:) * PETHETA(:,:)
-        !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
-      ELSE
-        ZW1 = ZWORK2
-      END IF
+      !!$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !ZWORK1(:,:) = (CST%XG *CST%XALPHAOC * PLM(:,:) * PLEPS(:,:) &
+      !                                 / PTKEM(:,:))**2
+      !!$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !CALL MZM_PHY(D,ZWORK1,ZWORK2)
+      !IF (KRR /= 0) THEN
+      !  !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !  ZW1(:,:) = ZWORK2(:,:) * PETHETA(:,:)
+      !  !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+      !ELSE
+      !  ZW1 = ZWORK2
+      !END IF
     ELSE
       !Compute only once and reuse in next JSV iterations
       IF ( JSV == 1 ) THEN
-- 
GitLab