diff --git a/src/MNH/p_abs.f90 b/src/MNH/p_abs.f90 index 077bdf2f435e100f14ad772be535b2d6efbcec17..16f888bc871cd2519f079ed9883dd8e8aa7fb27f 100644 --- a/src/MNH/p_abs.f90 +++ b/src/MNH/p_abs.f90 @@ -261,7 +261,11 @@ IF ( CEQNSYS=='DUR' .OR. CEQNSYS=='MAE' ) THEN DO JJ = IJB,IJE DO JI = IIB,IIE ZMASSGUESS_2D(JI,JJ) = ZMASSGUESS_2D(JI,JJ) + & +#ifndef MNH_BITREP (PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK))**ZCVD_O_RD & +#else + BR_POW((PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK)),ZCVD_O_RD) & +#endif * PRHODJ(JI,JJ,JK) / PRHODREF(JI,JJ,JK) & / PTHETAV(JI,JJ,JK) ZMASS_O_PI_2D(JI,JJ) = ZMASS_O_PI_2D(JI,JJ) + & @@ -291,7 +295,11 @@ IF ( CEQNSYS=='DUR' .OR. CEQNSYS=='MAE' ) THEN DO JJ = IJB,IJE DO JI = IIB,IIE ZMASSGUESS_2D(JI,JJ) = ZMASSGUESS_2D(JI,JJ) + & +#ifndef MNH_BITREP (PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK))**ZCVD_O_RD & +#else + BR_POW((PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK)),ZCVD_O_RD) & +#endif * ZWORK(JI,JJ,JK) / PTHETAV(JI,JJ,JK) END DO END DO @@ -301,7 +309,11 @@ IF ( CEQNSYS=='DUR' .OR. CEQNSYS=='MAE' ) THEN DO JJ = IJB,IJE DO JI = IIB,IIE ZMASSGUESS_2D(JI,JJ) = ZMASSGUESS_2D(JI,JJ) + & +#ifndef MNH_BITREP (PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK))**ZCVD_O_RD & +#else + BR_POW((PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK)),ZCVD_O_RD) & +#endif * PRHODJ(JI,JJ,JK) / PRHODREF(JI,JJ,JK) / PTHETAV(JI,JJ,JK) END DO END DO @@ -331,7 +343,11 @@ ELSEIF( CEQNSYS == 'LHE' ) THEN IF (LBOUSS) THEN ZRHOREF(:,:,:) = PRHODREF(:,:,:) ELSE +#ifndef MNH_BITREP ZRHOREF(:,:,:) = PEXNREF(:,:,:) ** ZCVD_O_RD & +#else + ZRHOREF(:,:,:) = BR_POW( PEXNREF(:,:,:), ZCVD_O_RD )& +#endif * XP00 / ( XRD * PTHVREF(:,:,:) ) ENDIF ! diff --git a/src/ZSOLVER/p_abs.f90 b/src/ZSOLVER/p_abs.f90 index e04c5db9b8d3c6d43c48086983af10f9e615b030..3ea1219e40ef239168cdb6e2a901fc637f64e7e9 100644 --- a/src/ZSOLVER/p_abs.f90 +++ b/src/ZSOLVER/p_abs.f90 @@ -310,7 +310,11 @@ IF ( CEQNSYS=='DUR' .OR. CEQNSYS=='MAE' ) THEN DO JJ = IJB,IJE DO JI = IIB,IIE ZMASSGUESS_2D(JI,JJ) = ZMASSGUESS_2D(JI,JJ) + & +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) (PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK))**ZCVD_O_RD & +#else + BR_POW((PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK)),ZCVD_O_RD) & +#endif * PRHODJ(JI,JJ,JK) / PRHODREF(JI,JJ,JK) & / PTHETAV(JI,JJ,JK) ZMASS_O_PI_2D(JI,JJ) = ZMASS_O_PI_2D(JI,JJ) + & @@ -361,7 +365,11 @@ IF ( CEQNSYS=='DUR' .OR. CEQNSYS=='MAE' ) THEN DO JJ = IJB,IJE DO JI = IIB,IIE ZMASSGUESS_2D(JI,JJ) = ZMASSGUESS_2D(JI,JJ) + & +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) (PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK))**ZCVD_O_RD & +#else + BR_POW((PEXNREF(JI,JJ,JK)+PPHIT(JI,JJ,JK)),ZCVD_O_RD) & +#endif * PRHODJ(JI,JJ,JK) / PRHODREF(JI,JJ,JK) / PTHETAV(JI,JJ,JK) END DO END DO @@ -393,7 +401,11 @@ ELSEIF( CEQNSYS == 'LHE' ) THEN IF (LBOUSS) THEN ZRHOREF(:,:,:) = PRHODREF(:,:,:) ELSE +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) ZRHOREF(:,:,:) = PEXNREF(:,:,:) ** ZCVD_O_RD & +#else + ZRHOREF(:,:,:) = BR_POW( PEXNREF(:,:,:), ZCVD_O_RD )& +#endif * XP00 / ( XRD * PTHVREF(:,:,:) ) ENDIF !