diff --git a/src/MNH/ice4_sedimentation_split.f90 b/src/MNH/ice4_sedimentation_split.f90
index ebd23fb6923ef2b023b899d7014d0252693ad348..1b19e954faba944f2152dfb1762f1e0cee8e023d 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
 !