Changeset 4029 for LMDZ6

Timestamp:

Nov 29, 2021, 12:56:23 PM (3 years ago)

Author:

musat

Message:

Ajout nouveaux flags wake par Jean-Yves Grandpeix

File:

• LMDZ6/trunk/DefLists/physiq.def_NPv6.3 (modified) (12 diffs)

LMDZ6/trunk/DefLists/physiq.def_NPv6.3

@@ -155 +155 @@
     # Controle de l'inertie des glaciers et banquises
     # inertie_sno : pour la neige (y compris sur les points "ter"
-    #      quand orchidee n'est pas actif utilisé)
+    #      quand orchidee n'est pas actif)
     # inertie_lic/sic : pour les "land ice" et "sea ice"
 inertie_sno=350
-    # iflag_sic=1 pour forcer à garder l'inertie des "sea ice"
+    # iflag_sic=1 pour forcer a garder l'inertie des "sea ice"
     # a inertie_sic meme si la surface est couverte de neige
 iflag_sic=1
@@ -164 +164 @@
 inertie_sic=2000
     # tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours)
-    #     # Utilise selon le forcage : tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2
-#tau_gl= _AUTO_ : DEFAULT = 5.
-tau_gl=5
+    # Utilise selon le forcage : tau_gl=5 pour Clim, tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2
+tau_gl=2
 
           # ???
@@ -179 +178 @@
 iflag_inertie=0
 
-          # Frequence appel rayonnement. Nombre appels par jour
+
+          # Frequence appel convection. Nombre appels par jour
+          # (nbapp_cv=0 ==> appel tous les pas de temps)
 nbapp_cv=48
 
@@ -191 +192 @@
 ###############################################################
 
-          #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
+          #seuils de la precipitation des nuages stratiformes (D: 2.6e-4 2.6e-4)
 cld_lc_lsc=0.00065
 cld_lc_con=0.00065
 
-          #constante de temps pour eleminer eau lsc et convective
+          #constante de temps pour eliminer eau lsc et convective
           # (D: 3600. 3600.)
 cld_tau_lsc=900
@@ -305 +306 @@
 iflag_con=3
 
+          # correct convection bug
+          #(default = y)
+keepbug_ice_frac=n
+
           #output level of energy conservation diagnostics
 if_ebil=0
@@ -381 +386 @@
           # Precipitation ejection from adiabatic ascents
           # cvflag_prec_eject: when TRUE ejection is activated. Default: False
+          # WARNING ! : cvflag_prec_eject=y is effective only if icvflag_Tpa=1
           # ejectliq : 0. =>liquid precipitation is carried up in the adiab updraughts
           #            1. =>liquid precipitation is ejected from adiab updraughts
@@ -392 +398 @@
 ejectice=0.
 
-          # Prevent some crashes
-          # Filter out bad wakes
-flag_wk_check_trgl=n
-iflag_wk_check_trgl=2
-iflag_alp_wk_cond=1
-
-
           # Check whether parcel level temperature and specific humidity
           # are reasonable and demand that LCL be high enough
@@ -429 +428 @@
 ##########################################################
 
-          # Appel separe des thermiques 
+          # Decimal code:  + first digit drives the splitting of thermals: 
+          #                                0 = thermals see the average profiles of T and q
+          #                                1 = thermals are only outside wakes
+          #                + second digit drives the splitting of diffuse PBL:
+          #                                0 = The diffuse PBL sees the average profiles of T and q
+          #                                1 = Two BLs, one in (w) and one in (x), with the same
+          #                                    surface temperature and humidity.
+          #                                2 = Two BLs, one in (w) and one in (x), with separate
+          #                                    surface temperature and humidity.
+          #                                3 = same as 1 over ocean; same as 0 everywhere else.
+          # Example: iflag_pbl_split=12 => thermals outside wakes, split BL with Tsurf and qsurf difference
+          # iflag_pbl_split=10 => thermals outside wakes
 iflag_pbl_split=10
 
@@ -511 +521 @@
 iflag_thermals_closure=2
 
+          # Flag permettant de garder le bug induisant des wake_deltaT non nuls
+          # a partir des tendances des thermiques, meme en l'absence de poches.
+          # Default : ok_bug_split_th=y ==> le bug est maintenu
+ok_bug_split_th=n
+
 
 ##########################################################
@@ -530 +545 @@
 ok_adjwk=y
 
+          # Prevent some crashes
+          # Filter out bad wakes
+flag_wk_check_trgl=n
+iflag_wk_check_trgl=2
+iflag_alp_wk_cond=1
+
+          #Wake population dynamic switch (D:0):
+          # iflag_wk_pop_dyn=0 ==> wake_dens is set to wdens_ref
+          # iflag_wk_pop_dyn=1 ==> wake_dens is a prognostic variable
+iflag_wk_pop_dyn=0
+
+          #Flag for the active wake ratio (D:1):
+          # iflag_wk_act=0: ==> active wake ratio = 0.
+          # iflag_wk_act=1: ==> active wake ratio = 1.
+          # iflag_wk_act=2: ==> active wake ratio = f(Wape)
+          # iflag_wk_act=3: ==> active wake ratio =
+          #     f(wake_dens, wake_radius, wape, Alp_wk, Alp_bl, birth rate fr cu)
+iflag_wk_act=1
+
 
 ####################################################################
@@ -535 +569 @@
 ####################################################################
           # The source terms of the TKE prognostic equation are diagnosed
-# from tendencies (du, dv, dT) 
+          # from tendencies (du, dv, dT) 
           # associated with subrgrid scaleamotions and treated as an 
           # additional heat source.
@@ -564 +598 @@
 ok_lic_cond=y
 ####################################################################
-
+          #min cloud droplet number concentration
 cdnc_min=10.
-          #min cloud droplet number concentration 
-ok_bug_split_th=n
-          # correct a bug in thermal splitting
-keepbug_ice_frac=n
-          # correct convection bug
-