Skip to content
Snippets Groups Projects
Commit af873adb authored by RIETTE Sébastien's avatar RIETTE Sébastien
Browse files

Sébastien Riette 14/02/2022 bit-reproduction PRFR option 1

Modification to reproduce old results with the PRFR option
parent 538ca2b7
No related branches found
No related tags found
No related merge requests found
...@@ -309,7 +309,11 @@ ELSE ...@@ -309,7 +309,11 @@ ELSE
CALL PRINT_MSG( NVERB_FATAL, 'GEN', 'ICE4_COMPUTE_PDF', 'wrong HSUBG_AUCV_RI case' ) CALL PRINT_MSG( NVERB_FATAL, 'GEN', 'ICE4_COMPUTE_PDF', 'wrong HSUBG_AUCV_RI case' )
ENDIF ENDIF
! !
#ifdef REPRO48
PRF=PHLC_HCF
#else
PRF=MAX(PHLC_HCF,PHLI_HCF) PRF=MAX(PHLC_HCF,PHLI_HCF)
#endif
! !
IF (LHOOK) CALL DR_HOOK('ICE4_COMPUTE_PDF', 1, ZHOOK_HANDLE) IF (LHOOK) CALL DR_HOOK('ICE4_COMPUTE_PDF', 1, ZHOOK_HANDLE)
END SUBROUTINE ICE4_COMPUTE_PDF END SUBROUTINE ICE4_COMPUTE_PDF
......
...@@ -328,8 +328,10 @@ IF (LLRFR) THEN ...@@ -328,8 +328,10 @@ IF (LLRFR) THEN
DO JL=1,KSIZE DO JL=1,KSIZE
PRAINFR(K1(JL), K2(JL), K3(JL)) = ZRAINFR(JL) PRAINFR(K1(JL), K2(JL), K3(JL)) = ZRAINFR(JL)
ZRRT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRR) ZRRT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRR)
#ifndef REPRO48
ZRST3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRS) ZRST3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRS)
ZRGT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRG) ZRGT3D (K1(JL), K2(JL), K3(JL)) = ZVART(JL,IRG)
#endif
END DO END DO
IF (KRR==7) THEN IF (KRR==7) THEN
DO JL=1,KSIZE DO JL=1,KSIZE
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment