From e4bd9cbc337960ff7c56e0302092cba17e64c51a Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Riette?= <sebastien.riette@meteo.fr>
Date: Wed, 4 Jan 2023 10:20:12 +0100
Subject: [PATCH] S. Riette 4 Jan 2022: bf for belenos bit-repro
---
src/arome/ext/aroini_micro.F90 | 3 ++-
src/common/micro/mode_ice4_fast_rs.F90 | 4 ++++
2 files changed, 6 insertions(+), 1 deletion(-)
diff --git a/src/arome/ext/aroini_micro.F90 b/src/arome/ext/aroini_micro.F90
index 347033c9c..77cc42030 100644
--- a/src/arome/ext/aroini_micro.F90
+++ b/src/arome/ext/aroini_micro.F90
@@ -146,7 +146,8 @@ XTSTEP_TS=PTSTEP_TS ! Approximative time step for time-splitting (0 for no time-
!
! 2. Set implicit default values for MODD_RAIN_ICE_DESCR
! et MODD_RAIN_ICE_PARAM
-
+LSNOW_T=.FALSE.
+LRED=CMICRO=='ICE3' .OR. CMICRO=='ICE4'
CALL INI_RAIN_ICE (KULOUT, PTSTEP, 20.,KSPLITR,CMICRO)
CALL INI_TIWMX
diff --git a/src/common/micro/mode_ice4_fast_rs.F90 b/src/common/micro/mode_ice4_fast_rs.F90
index d1a1afb4b..adc9f4715 100644
--- a/src/common/micro/mode_ice4_fast_rs.F90
+++ b/src/common/micro/mode_ice4_fast_rs.F90
@@ -158,7 +158,11 @@ IF(.NOT. LDSOFT) THEN
! 5.1.1 select the PLBDAS
!
DO JJ = 1, IGRIM
+#if defined(REPRO48) || defined(REPRO55)
+ ZVEC1(JJ) = PLBDAS(I1(JJ))
+#else
ZVEC1(JJ) = (PLBDAS(I1(JJ))**ICED%XALPHAS + ICED%XFVELOS**ICED%XALPHAS)**(1./ICED%XALPHAS)
+#endif
END DO
!
! 5.1.2 find the next lower indice for the PLBDAS in the geometrical
--
GitLab