From ea56b9327d1174cef68fa4ee3b206ca4891b678e Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Wed, 22 Jun 2022 12:54:52 +0000
Subject: [PATCH] Philippe 22/06/2022: use gamma intrinsics if CCE for
 performance problems

---
 src/MNH/ice4_sedimentation_split.f90 | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/MNH/ice4_sedimentation_split.f90 b/src/MNH/ice4_sedimentation_split.f90
index ebd23fb69..1b19e954f 100644
--- a/src/MNH/ice4_sedimentation_split.f90
+++ b/src/MNH/ice4_sedimentation_split.f90
@@ -87,7 +87,11 @@ USE MODE_MNH_ZWORK,      ONLY: MNH_MEM_GET, MNH_MEM_POSITION_PIN, MNH_MEM_RELEAS
 #endif
 USE MODE_MPPDB
 !
+#ifndef MNH_COMPILER_CCE
+!Do not use internal GAMMA function but the Fortran intrinsics one with Cray compiler
+!Warning: intrinsics gamma does not give same results between CPU and GPU with NVHPC (tested with 22.2 version)
 USE MODI_GAMMA
+#endif
 !
 IMPLICIT NONE
 !
-- 
GitLab