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

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

Suppress error message in ice4_tendencies for PRFR option as
it is redundant with the message in rain_ice
parent bea2de71
No related branches found
No related tags found
No related merge requests found
......@@ -304,15 +304,6 @@ CALL ICE4_COMPUTE_PDF(KSIZE, HSUBG_AUCV_RC, HSUBG_AUCV_RI, HSUBG_PR_PDF,&
PHLI_HCF, PHLI_LCF, PHLI_HRI, PHLI_LRI, ZRAINFR)
LLRFR=HSUBG_RC_RR_ACCR=='PRFR' .OR. HSUBG_RR_EVAP=='PRFR'
IF (LLRFR) THEN
CALL PRINT_MSG(NVERB_FATAL, 'GEN', 'MODE_ICE4_TENDENCIES', 'LLRFR case broken by optimisation, see comments in mode_ice4_tendencies to knwon why (and how to reapir)....')
!Microphyscs was optimized by introducing chunks of KPROMA size
!Thus, in ice4_tendencies, the 1D array represent only a fraction of the points where microphisical species are present
!We cannot rebuild the entire 3D arrays here, so we cannot call ice4_rainfr_vert here
!A solution would be to suppress optimisation in this case by setting KPROMA=KSIZE in rain_ice
!Another solution would be to compute column by column?
!Another one would be to cut tendencies in 3 parts: before rainfr_vert, rainfr_vert, after rainfr_vert
!Diagnostic of precipitation fraction
PRAINFR(:,:,:) = 0.
ZRRT3D (:,:,:) = 0.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment