From 6843cc1031a1ec688f042bdec5f1ba7d50cd552d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?S=C3=A9bastien=20Riette?= <sebastien.riette@meteo.fr> Date: Tue, 30 May 2023 12:21:27 +0200 Subject: [PATCH] STAT cloud scheme for LMDZ --- src/lmdz/ext/physiqex_mod.F90 | 17 +++++------------ 1 file changed, 5 insertions(+), 12 deletions(-) diff --git a/src/lmdz/ext/physiqex_mod.F90 b/src/lmdz/ext/physiqex_mod.F90 index 2cc767c31..545cab7c7 100644 --- a/src/lmdz/ext/physiqex_mod.F90 +++ b/src/lmdz/ext/physiqex_mod.F90 @@ -19,18 +19,8 @@ CONTAINS ! * Certains allocatable sont en klev, d'autres en klev+2. Il est possible de changer ceci mais il faut gérer ! les recopies pour que les params voient des tableaux en klev+2 (en effectuant des recopies des ! niveaux extrêmes comme fait pour le vent par exemple) - -!TODO à faire avant d'historiser dans LMDZ: -! * vérifier si rhodj ne devrait pas être calculé avec un rho humide ici (arome?) -! * brancher le sigma du schéma STAT -! * utiliser les variables ajustées - - - - - - - +! * L'eau tombée en surface (précipitations, sédimentation du nuage, terme de dépôt) se trouve dans +! les variables ZINPRC, ZINPRR, ZINPRS, ZINPRG @@ -518,6 +508,9 @@ ZTKES(:,:) = ZRTKES(:,:) / PRHODJ(:,:) ! Add tendencies of turb to total physics tendency d_u(:,1:klev) = d_u(:,1:klev) + ZRUS(:,2:klev+1)/PRHODJ(:,2:klev+1) d_v(:,1:klev) = d_v(:,1:klev) + ZRVS(:,2:klev+1)/PRHODJ(:,2:klev+1) +IF(PHYEX%PARAM_MFSHALLN%CMF_CLOUD=='STAT') THEN + PSIGS(:,:)=SQRT(PSIGS(:,:)**2 + PSIGMF(:,:)**2) +ENDIF !------------------------------------------------------------ ! Microphysics !------------------------------------------------------------ -- GitLab