Changeset 3199 for LMDZ6/trunk/libf/phylmd

Timestamp:

Feb 12, 2018, 1:40:09 AM (7 years ago)

Author:

fhourdin

Message:

Pour la convergence independamment du decoupage de domaine.
Pour la convection, le niveau maximum k_upper_cv etait fixe
sur un profil d'altitude au milieu du domaine klon/2+1/klon
dependant donc du decoupage.
Remplace par un test en presnivs
-7*log(presnivs(k)/presnivs(1)) > 25.
independant du decoupage.
Garantit la convergence numerique pour l'ancienne physique.
Pas forcement pour la nouvelle si on a des valeurs
non nulles des variables convectives vers le sommet.

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -2635 +2635 @@
           ! Calculate the upmost level of deep convection loops: k_upper_cv
           !  (near 22 km)
-          izero = klon/2+1/klon
           k_upper_cv = klev
+          !izero = klon/2+1/klon
+          !DO k = klev,1,-1
+          !   IF (pphi(izero,k) > 22.e4) k_upper_cv = k
+          !ENDDO
+          ! FH : nouveau calcul base sur un profil global sans quoi
+          ! le modele etait sensible au decoupage de domaines
           DO k = klev,1,-1
-             IF (pphi(izero,k) > 22.e4) k_upper_cv = k
+             IF (-7*log(presnivs(k)/presnivs(1)) > 25.) k_upper_cv = k
           ENDDO
           IF (prt_level .ge. 5) THEN