Skip to content
Snippets Groups Projects
Commit da8fa56f authored by CHABOUREAU Jean-Pierre's avatar CHABOUREAU Jean-Pierre
Browse files

Jean-Pierre 24/03/2017: LIMA bug fixes

(cherry picked from commit 71693a7b)
parent 6160c71a
No related branches found
No related tags found
No related merge requests found
...@@ -128,6 +128,7 @@ END MODULE MODI_LIMA_ADJUST ...@@ -128,6 +128,7 @@ END MODULE MODI_LIMA_ADJUST
!! ------------- !! -------------
!! Original ??/??/13 !! Original ??/??/13
!! C. Barthe * LACy* jan. 2014 add budgets !! C. Barthe * LACy* jan. 2014 add budgets
!! JP Chaboureau *LA* March 2014 fix the calculation of icy cloud fraction
!! !!
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
...@@ -1105,7 +1106,7 @@ DEALLOCATE(ZCTMIN) ...@@ -1105,7 +1106,7 @@ DEALLOCATE(ZCTMIN)
!* 5.2 compute the cloud fraction PCLDFR (binary !!!!!!!) !* 5.2 compute the cloud fraction PCLDFR (binary !!!!!!!)
! !
IF ( .NOT. OSUBG_COND ) THEN IF ( .NOT. OSUBG_COND ) THEN
WHERE (PRCS(:,:,:) > 1.E-12 / ZDT) WHERE (PRCS(:,:,:) + PRIS(:,:,:) + PRSS(:,:,:) > 1.E-12 / ZDT)
ZW(:,:,:) = 1. ZW(:,:,:) = 1.
ELSEWHERE ELSEWHERE
ZW(:,:,:) = 0. ZW(:,:,:) = 0.
......
...@@ -284,6 +284,8 @@ END MODULE MODI_READ_EXSEG_n ...@@ -284,6 +284,8 @@ END MODULE MODI_READ_EXSEG_n
!! Modification 02/2016 (M.Leriche) treat gas and aq. chemicals separately !! Modification 02/2016 (M.Leriche) treat gas and aq. chemicals separately
!! Modification 10/2016 (C.LAC) Add OSPLIT_WENO + Add droplet !! Modification 10/2016 (C.LAC) Add OSPLIT_WENO + Add droplet
!! deposition + Add max values !! deposition + Add max values
!! Modification 03/2017 (JP Chaboureau) Fix the initialization of
!! LUSERx-type variables for LIMA
!!------------------------------------------------------------------------------ !!------------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -871,6 +873,20 @@ SELECT CASE ( CCLOUD ) ...@@ -871,6 +873,20 @@ SELECT CASE ( CCLOUD )
WRITE(UNIT=ILUOUT,FMT='("YOU WANT TO USE THE C3R5 MICROPHYS. SCHEME.",& WRITE(UNIT=ILUOUT,FMT='("YOU WANT TO USE THE C3R5 MICROPHYS. SCHEME.",&
&" YOU WILL HAVE VAPOR, CLOUD WATER/ICE, RAIN, SNOW AND GRAUPEL ",/, & &" YOU WILL HAVE VAPOR, CLOUD WATER/ICE, RAIN, SNOW AND GRAUPEL ",/, &
&"LUSERV, LUSERC, LUSERR, LUSERI, LUSERS, LUSERG ARE SET TO TRUE")' ) &"LUSERV, LUSERC, LUSERR, LUSERI, LUSERS, LUSERG ARE SET TO TRUE")' )
!
LUSERV=.TRUE. ; LUSERC=.TRUE. ; LUSERR=.TRUE.
LUSERI=.TRUE. ; LUSECI=.TRUE.
LUSERS=.TRUE. ; LUSERG=.TRUE.
LUSERH=.FALSE.
END IF
ELSE IF (CCLOUD == 'LIMA') THEN
IF (.NOT. ( LUSERV .AND. LUSERC .AND. LUSERR .AND. LUSERI .AND. &
LUSERS .AND. LUSERG .AND. (.NOT. LUSERH) &
) ) THEN
WRITE(UNIT=ILUOUT,FMT=9002) KMI
WRITE(UNIT=ILUOUT,FMT='("YOU WANT TO USE THE LIMA MICROPHYS. SCHEME.",&
&" YOU WILL HAVE VAPOR, CLOUD WATER/ICE, RAIN, SNOW AND GRAUPEL ",/, &
&"LUSERV, LUSERC, LUSERR, LUSERI, LUSERS, LUSERG ARE SET TO TRUE")' )
! !
LUSERV=.TRUE. ; LUSERC=.TRUE. ; LUSERR=.TRUE. LUSERV=.TRUE. ; LUSERC=.TRUE. ; LUSERR=.TRUE.
LUSERI=.TRUE. ; LUSECI=.TRUE. LUSERI=.TRUE. ; LUSECI=.TRUE.
......
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