Changeset 2000 for LMDZ5/trunk/libf/phylmd

Timestamp:

Mar 26, 2014, 7:05:47 PM (11 years ago)

Author:

fhourdin

Message:

Introduction d'un parametre pour controler l'altitude a laquelle on va chercher l'air
pour l'entrainement pour obtenir des stratocumulus
dz = fact_thermals_ed_dz * z (these Arnaud Jam)

Location:

LMDZ5/trunk/libf/phylmd

Files:

• conf_phys_m.F90 (modified) (4 diffs)
• thermcell.h (modified) (2 diffs)
• thermcell_plume.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -115 +115 @@
 
     integer,SAVE :: iflag_thermals_ed_omp,iflag_thermals_optflux_omp,iflag_thermals_closure_omp
+    real, SAVE :: fact_thermals_ed_dz_omp
     integer,SAVE :: iflag_thermals_omp,nsplit_thermals_omp
     real,save :: tau_thermals_omp,alp_bl_k_omp
@@ -1141 +1142 @@
     iflag_thermals_omp = 0
     call getin('iflag_thermals',iflag_thermals_omp)
+    !
+    !Config Key  = iflag_thermals_ed
+    !Config Desc =
+    !Config Def  = 0
+    !Config Help =
+    !
+    fact_thermals_ed_dz_omp = 0.1
+
+    call getin('fact_thermals_ed_dz',fact_thermals_ed_dz_omp)
     !
     !
@@ -1823 +1833 @@
     iflag_thermals = iflag_thermals_omp
     iflag_thermals_ed = iflag_thermals_ed_omp
+    fact_thermals_ed_dz = fact_thermals_ed_dz_omp
     iflag_thermals_optflux = iflag_thermals_optflux_omp
     iflag_thermals_closure = iflag_thermals_closure_omp
@@ -2015 +2026 @@
     write(lunout,*)' iflag_thermals = ', iflag_thermals
     write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed
+    write(lunout,*)' fact_thermals_ed_dz = ', fact_thermals_ed_dz
     write(lunout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
     write(lunout,*)' iflag_thermals_closure = ', iflag_thermals_closure

LMDZ5/trunk/libf/phylmd/thermcell.h

@@ -9 +9 @@
       real,parameter     :: r_aspect_thermals=2.,l_mix_thermals=30.
       real               :: alp_bl_k
-      real               :: tau_thermals
+      real               :: tau_thermals,fact_thermals_ed_dz
       integer,parameter  :: w2di_thermals=0
       integer            :: isplit
@@ -17 +17 @@
 
       common/ctherm1/iflag_thermals,nsplit_thermals,iflag_thermals_closure
-      common/ctherm2/tau_thermals,alp_bl_k
+      common/ctherm2/tau_thermals,alp_bl_k,fact_thermals_ed_dz
       common/ctherm4/iflag_coupl,iflag_clos,iflag_wake
       common/ctherm5/iflag_thermals_ed,iflag_thermals_optflux

LMDZ5/trunk/libf/phylmd/thermcell_plume.F90

@@ -295 +295 @@
 !Modif AJAM
           
-        lmel=0.1*zlev(ig,l)
+        lmel=fact_thermals_ed_dz*zlev(ig,l)
 !        lmel=2.5*(zlev(ig,l)-zlev(ig,l-1))
         lt=l+1