diff --git a/src/MNH/fast_terms.f90 b/src/MNH/fast_terms.f90
index f618e6aeaa9380f771fa1a377b51a4ffdc5dd076..b52417e04a42ca5acf830747b3333db1dc572383 100644
--- a/src/MNH/fast_terms.f90
+++ b/src/MNH/fast_terms.f90
@@ -337,7 +337,11 @@ END WHERE
 !
 !*       2.2    estimate the Exner function at t+1
 !
-ZEXNS(:,:,:) = ( PPABST(:,:,:)  / XP00 ) ** (XRD/XCPD)  
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
+ZEXNS(:,:,:) = ( PPABST(:,:,:) / XP00 ) ** (XRD/XCPD)
+#else
+ZEXNS(:,:,:) = BR_POW( PPABST(:,:,:) / XP00, XRD / XCPD )
+#endif
 !$acc end kernels
 !
 !    beginning of the iterative loop
@@ -372,7 +376,11 @@ DO JITER =1,ITERMAX
   CALL GET_HALO_D(ZT)
 !$acc kernels
 #endif
+#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
   ZW1(:,:,:) = EXP( XALPW - XBETAW/ZT(:,:,:) - XGAMW*ALOG( ZT(:,:,:) ) )
+#else
+  ZW1(:,:,:) = BR_EXP( XALPW - XBETAW/ZT(:,:,:) - XGAMW * BR_LOG( ZT(:,:,:) ) )
+#endif
 !
 !*       2.7    compute the saturation mixing ratio at t+1
 !