From a43447a3c1f6357246e8e9148dcd8348e6373285 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Thu, 23 Jun 2022 12:30:27 +0000
Subject: [PATCH] Philippe 23/06/2022: gamma: add support for MNH_BITREP_OMP

---
 src/MNH/gamma.f90 | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/src/MNH/gamma.f90 b/src/MNH/gamma.f90
index 763042054..626aab326 100644
--- a/src/MNH/gamma.f90
+++ b/src/MNH/gamma.f90
@@ -72,7 +72,7 @@ END MODULE MODI_GAMMA
 !*       0. DECLARATIONS
 !           ------------
 !
-#ifdef MNH_BITREP
+#if defined(MNH_BITREP) || defined(MNH_BITREP_OMP)
 USE MODI_BITREP
 #endif
 !
@@ -119,7 +119,7 @@ ELSE
 END IF
 ZY = ZX
 ZTMP =  ZX + 5.5
-#ifndef MNH_BITREP
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
 ZTMP = (ZX + 0.5) * ALOG(ZTMP) - ZTMP
 #else
 ZTMP = (ZX + 0.5) * BR_LOG(ZTMP) - ZTMP
@@ -132,13 +132,13 @@ DO JJ = 1, 6
 END DO
 !
 IF (PX .LT. 0.) THEN
-#ifndef MNH_BITREP
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
   PGAMMA = ZPI / SIN(ZPI*PX) / EXP(ZTMP + ALOG(ZSTP*ZSER/ZX))
 #else
   PGAMMA = ZPI / SIN(ZPI*PX) / BR_EXP(ZTMP + BR_LOG(ZSTP*ZSER/ZX))
 #endif
 ELSE
-#ifndef MNH_BITREP
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
   PGAMMA = EXP(ZTMP + ALOG(ZSTP*ZSER/ZX))
 #else
   PGAMMA = BR_EXP(ZTMP + BR_LOG(ZSTP*ZSER/ZX))
@@ -197,7 +197,7 @@ END FUNCTION GAMMA_X0D
 !*       0. DECLARATIONS
 !           ------------
 !
-#ifdef MNH_BITREP
+#if defined(MNH_BITREP) || defined(MNH_BITREP_OMP)
 USE MODI_BITREP
 #endif
 !
@@ -235,7 +235,7 @@ WHERE ( PX(:)<0.0 )
 END WHERE
 ZY(:) = ZX(:)
 ZTMP(:) =  ZX(:) + 5.5
-#ifndef MNH_BITREP
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
 ZTMP(:) = (ZX(:) + 0.5)*ALOG(ZTMP(:)) - ZTMP(:)
 #else
 ZTMP(:) = (ZX(:) + 0.5)*BR_LOG(ZTMP(:)) - ZTMP(:)
@@ -247,7 +247,7 @@ DO JJ = 1 , 6
   ZSER(:) = ZSER(:) + ZCOEF(JJ)/ZY(:)
 END DO
 !
-#ifndef MNH_BITREP
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
 PGAMMA(:) = EXP( ZTMP(:) + ALOG( ZSTP*ZSER(:)/ZX(:) ) )
 #else
 PGAMMA(:) = BR_EXP( ZTMP(:) + BR_LOG( ZSTP*ZSER(:)/ZX(:) ) )
-- 
GitLab