From af873adbc15b60bd0920af8b12b1a0af11fba93a Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Riette?= <sebastien.riette@meteo.fr>
Date: Mon, 14 Feb 2022 09:49:01 +0100
Subject: [PATCH] =?UTF-8?q?S=C3=A9bastien=20Riette=2014/02/2022=20bit-repr?=
=?UTF-8?q?oduction=20PRFR=20option=201?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
Modification to reproduce old results with the PRFR option
---
src/common/micro/mode_ice4_compute_pdf.F90 | 4 ++++
src/common/micro/mode_ice4_tendencies.F90 | 2 ++
2 files changed, 6 insertions(+)
diff --git a/src/common/micro/mode_ice4_compute_pdf.F90 b/src/common/micro/mode_ice4_compute_pdf.F90
index 6fb091d6c..942c65c25 100644
--- a/src/common/micro/mode_ice4_compute_pdf.F90
+++ b/src/common/micro/mode_ice4_compute_pdf.F90
@@ -309,7 +309,11 @@ ELSE
CALL PRINT_MSG( NVERB_FATAL, 'GEN', 'ICE4_COMPUTE_PDF', 'wrong HSUBG_AUCV_RI case' )
ENDIF
!
+#ifdef REPRO48
+PRF=PHLC_HCF
+#else
PRF=MAX(PHLC_HCF,PHLI_HCF)
+#endif
!
IF (LHOOK) CALL DR_HOOK('ICE4_COMPUTE_PDF', 1, ZHOOK_HANDLE)
END SUBROUTINE ICE4_COMPUTE_PDF
diff --git a/src/common/micro/mode_ice4_tendencies.F90 b/src/common/micro/mode_ice4_tendencies.F90
index 766b7e4b0..f9af1458d 100644
--- a/src/common/micro/mode_ice4_tendencies.F90
+++ b/src/common/micro/mode_ice4_tendencies.F90
@@ -328,8 +328,10 @@ IF (LLRFR) THEN
DO JL=1,KSIZE
PRAINFR(K1(JL), K2(JL), K3(JL)) = ZRAINFR(JL)
ZRRT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRR)
+#ifndef REPRO48
ZRST3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRS)
ZRGT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRG)
+#endif
END DO
IF (KRR==7) THEN
DO JL=1,KSIZE
--
GitLab