Changeset 4653 for LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90

Timestamp:

Aug 30, 2023, 9:15:43 AM (15 months ago)

Author:

evignon

Message:

prise en compte de l'humidite pour le calcul du flux de flottabilite dans atke
+ petits ajustements

File:

• LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90

@@ -13 +13 @@
   integer :: lunout,prt_level
   !$OMP THREADPRIVATE(lunout,prt_level)
-  real :: rg, rd, rpi, rcpd
-  !$OMP THREADPRIVATE(rg, rd, rpi, rcpd)
+  real :: rg, rd, rpi, rcpd, rv
+  !$OMP THREADPRIVATE(rg, rd, rpi, rcpd, rv)
 
   real :: viscom, viscoh
@@ -22 +22 @@
   !$OMP THREADPRIVATE(lmin)
 
-  logical :: ok_vdiff_tke
-  !$OMP THREADPRIVATE(ok_vdiff_tke)
+  logical :: atke_ok_vdiff, atke_ok_virtual
+  !$OMP THREADPRIVATE(atke_ok_vdiff,atke_ok_virtual)
 
 CONTAINS
 
-SUBROUTINE atke_ini(prt_level_in, lunout_in, rg_in, rd_in, rpi_in, rcpd_in)
+SUBROUTINE atke_ini(prt_level_in, lunout_in, rg_in, rd_in, rpi_in, rcpd_in, rv_in)
 
    USE ioipsl_getin_p_mod, ONLY : getin_p
 
   integer, intent(in) :: lunout_in,prt_level_in
-  real, intent(in) :: rg_in, rd_in, rpi_in, rcpd_in
+  real, intent(in) :: rg_in, rd_in, rpi_in, rcpd_in, rv_in
 
 
@@ -41 +41 @@
   rpi=rpi_in
   rcpd=rcpd_in
+  rv=rv_in
 
   viscom=1.46E-5
@@ -59 +60 @@
 
   ! activate vertical diffusion of TKE or not
-  ok_vdiff_tke=.false.
-  CALL getin_p('atke_ok_vdiff_tke',ok_vdiff_tke)
+  atke_ok_vdiff=.false.
+  CALL getin_p('atke_ok_vdiff',atke_ok_vdiff)
+
+
+  ! account for vapor for flottability
+  atke_ok_virtual=.true.
+  CALL getin_p('atke_ok_virtual',atke_ok_virtual)
+
 
   ! flag that controls the numerical treatment of diffusion coeffiient calculation
@@ -82 +89 @@
 
   ! constant for tke dissipation calculation
-  cepsilon=16.6/2./sqrt(2.) ! default value as in yamada4
+  cepsilon=5.87 ! default value as in yamada4
   CALL getin_p('atke_cepsilon',cepsilon)
+
+
+ ! coefficient for surface TKE 
+ ! following Lenderink & Holtslag 2004, ctkes=sqrt(cepsilon)
+ ! (provided by limit condition in neutral conditions)
+  ctkes=sqrt(cepsilon)
+
 
   ! slope of Pr=f(Ri) for stable conditions
@@ -111 +125 @@
   CALL getin_p('atke_smmin',smmin)
 
- ! coefficient for surface TKE (default value from Arpege, see E. Bazile note)
-  ctkes=3.75
-  CALL getin_p('atke_ctkes',ctkes)
 
   ! ratio between the eddy diffusivity coeff for tke wrt that for momentum