Changeset 4441 for LMDZ6/trunk

Timestamp:

Feb 17, 2023, 5:11:37 PM (20 months ago)

Author:

fhourdin

Message:

Parametre controlant les downdrafts

Location:

LMDZ6/trunk/libf/phylmd

Files:

• thermcell_down.F90 (modified) (2 diffs)
• thermcell_ini_mod.F90 (modified) (4 diffs)
• thermcell_main.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/thermcell_down.F90

@@ -223 +223 @@
 
 
-   USE thermcell_ini_mod, ONLY : prt_level,RLvCp,RKAPPA,RETV
+   USE thermcell_ini_mod, ONLY : prt_level,RLvCp,RKAPPA,RETV,fact_thermals_down
    IMPLICIT NONE
 
@@ -271 +271 @@
       do ig=1,ngrid
          if (ilay.le.lmax(ig).and.lmax(ig)>1) then
-            edn(ig,ilay)=0.5*dup(ig,ilay)
-            ddn(ig,ilay)=0.5*eup(ig,ilay)
+            edn(ig,ilay)=fact_thermals_down*dup(ig,ilay)
+            ddn(ig,ilay)=fact_thermals_down*eup(ig,ilay)
             fdn(ig,ilay)=fdn(ig,ilay+1)+edn(ig,ilay)-ddn(ig,ilay)
             thetad(ig,ilay)=( fdn(ig,ilay+1)*thetad(ig,ilay+1) + edn(ig,ilay)*theta(ig,ilay) ) / (fdn(ig,ilay)+ddn(ig,ilay))

LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90

@@ -3 +3 @@
 
 save
+
    integer :: dvdq=1,dqimpl=-1,prt_level=0,lunout
    real RG,RD,RCPD,RKAPPA,RLVTT,RLvCp,RETV
@@ -8 +9 @@
    integer        :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure
    integer        :: iflag_thermals_down
-   real           :: entrdn_coef, detrdn_coef
+   real           :: fact_thermals_down
 
 !$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout)
@@ -15 +16 @@
 !$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure)
 !$OMP THREADPRIVATE(iflag_thermals_down)
-!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef)
+!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef,fact_thermals_down)
 
 
@@ -92 +93 @@
    CALL getin_p('iflag_thermals_down',iflag_thermals_down)
 
-
-   entrdn_coef = 0.5
-   CALL getin_p('entrdn_coef',entrdn_coef)
-
-   detrdn_coef = 0.5
-   CALL getin_p('detrdn_coef',detrdn_coef)
-
+   fact_thermals_down = 0.5
+   CALL getin_p('fact_thermals_down',fact_thermals_down)
 
      CALL getin_p('thermals_fact_epsilon',fact_epsilon)

LMDZ6/trunk/libf/phylmd/thermcell_main.F90

@@ -18 +18 @@
       USE thermcell_ini_mod, ONLY: thermcell_ini,dqimpl,dvdq,prt_level,lunout,prt_level
       USE thermcell_ini_mod, ONLY: iflag_thermals_closure,iflag_thermals_ed,tau_thermals,r_aspect_thermals
-      USE thermcell_ini_mod, ONLY: iflag_thermals_down
-      USE thermcell_ini_mod, ONLY: entrdn_coef, detrdn_coef
+      USE thermcell_ini_mod, ONLY: iflag_thermals_down, fact_thermals_down
       USE thermcell_ini_mod, ONLY: RD,RG
 
@@ -503 +502 @@
       IF (iflag_thermals_down .GT. 0) THEN
         if (debut) print*,'WARNING !!! routine thermcell_down en cours de developpement'
-        entrdn=entrdn_coef*detr0
-        detrdn=detrdn_coef*entr0 
+        entrdn=fact_thermals_down*detr0
+        detrdn=fact_thermals_down*entr0 
         ! we want to transport potential temperature, total water and momentum
         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zthl,zdthladj)