Skip to content
Snippets Groups Projects
Commit 4e8401e9 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 25/10/2021: bugfix RESTART reproductibility. In the case of no online...

Quentin 25/10/2021: bugfix RESTART reproductibility. In the case of no online chemistry (LORILAM=F) or LDUST/LSALT=T processes changing the aerosols (XAER)
at each time-step, XAER is computed once at the START segment in ini_ecmwf_radiations based on current pressure and not updated during MESONH except if LAERO_FT=T and CAER='TEGE'. At RESTART, the aerosols climatology is updated but the
pressure has changed. For reproductibility, one must use LAERO_FT=T and CAER='TEGE'. XAER_CLIM was not updated in the process done by CALL AEROZON in phys_param
parent 74ed2142
No related branches found
No related tags found
No related merge requests found
...@@ -680,6 +680,7 @@ IF (.NOT. OCLOUD_ONLY .AND. KTCOUNT /= 1) THEN ...@@ -680,6 +680,7 @@ IF (.NOT. OCLOUD_ONLY .AND. KTCOUNT /= 1) THEN
NDLON,NFLEV,CAER,NAER,NSTATM, & NDLON,NFLEV,CAER,NAER,NSTATM, &
XSINDEL,XCOSDEL,XTSIDER,XCORSOL, & XSINDEL,XCOSDEL,XTSIDER,XCORSOL, &
XSTATM,XOZON, XAER) XSTATM,XOZON, XAER)
XAER_CLIM = XAER
END IF END IF
END IF END IF
! !
......
...@@ -2149,6 +2149,13 @@ IF ( CRAD /= 'NONE' .AND. CPROGRAM=='MESONH' ) THEN ...@@ -2149,6 +2149,13 @@ IF ( CRAD /= 'NONE' .AND. CPROGRAM=='MESONH' ) THEN
IF(CCONF=='START') THEN IF(CCONF=='START') THEN
CGETRAD='INIT' CGETRAD='INIT'
END IF END IF
IF(CCONF=='RESTA' .AND. (.NOT. LAERO_FT) .AND. (.NOT. LORILAM) &
.AND. (.NOT. LSALT) .AND. (.NOT. LDUST)) THEN
WRITE(UNIT=ILUOUT,FMT=9001) KMI
WRITE(UNIT=ILUOUT,FMT=*) '!!! WARNING !!! FOR REPRODUCTIBILITY BETWEEN START and START+RESTART,'
WRITE(UNIT=ILUOUT,FMT=*) 'YOU MUST USE LAERO_FT=T WITH CAER=TEGE IF CCONF=RESTA IN ALL SEGMENTS'
WRITE(UNIT=ILUOUT,FMT=*) 'TO UPDATE THE OZONE AND AEROSOLS CLIMATOLOGY USED BY THE RADIATION CODE;'
END IF
END IF END IF
! !
! 3.6 check the initialization of the deep convection scheme ! 3.6 check the initialization of the deep convection scheme
......
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