Philippe 23/06/2022: gamma: add support for MNH_BITREP_OMP

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(:) ) )