Changeset 4745 for LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90

Timestamp:

Nov 7, 2023, 8:27:07 AM (7 months ago)

Author:

evignon

Message:

nettoyage et corrections dans les routines atke pour utilisation en 3D (terre + mars)

File:

• LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90

@@ -3 +3 @@
 implicit none
 
+! declaration of constants and parameters
 save
 
@@ -15 +16 @@
   real :: rg, rd, rpi, rcpd, rv
   !$OMP THREADPRIVATE(rg, rd, rpi, rcpd, rv)
-
   real :: viscom, viscoh
   !$OMP THREADPRIVATE(viscom,viscoh)
-
   real :: lmin=0.01              ! minimum mixing length
   !$OMP THREADPRIVATE(lmin)
-
   logical :: atke_ok_vdiff, atke_ok_virtual
   !$OMP THREADPRIVATE(atke_ok_vdiff,atke_ok_virtual)
@@ -27 +25 @@
 CONTAINS
 
-SUBROUTINE atke_ini(prt_level_in, lunout_in, rg_in, rd_in, rpi_in, rcpd_in, rv_in)
+SUBROUTINE atke_ini(rg_in, rd_in, rpi_in, rcpd_in, rv_in, viscom_in, viscoh_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, rv_in
+  real, intent(in) :: rg_in, rd_in, rpi_in, rcpd_in, rv_in, viscom_in, viscoh_in
 
 
-  lunout=lunout_in
-  prt_level=prt_level_in
+  ! input arguments (universal constants for planet)
+  !-------------------------------------------------
+ 
+  ! gravity acceleration
+  rg=rg_in
+  ! dry gas constant (R/M, R=perfect gas constant and M is the molar mass of the fluid)
   rd=rd_in
-  rg=rg_in
+  ! Pi number
   rpi=rpi_in
+  ! cp per unit mass of the gas
   rcpd=rcpd_in
+  ! water vapor constant (for simulations in Earth's atmosphere)
   rv=rv_in
+  ! kinematic molecular viscosity for momentum
+  viscom=viscom_in
+  ! kinematic molecular viscosity for heat
+  viscoh=viscoh_in
 
-  viscom=1.46E-5
-  viscoh=2.06E-5
+
+  !viscom=1.46E-5
+  !viscoh=2.06E-5
+
+
+  ! Read flag values in .def files
+  !-------------------------------
+
 
   ! flag that controls options in atke_compute_km_kh
@@ -124 +137 @@
 
   ! minimum anisotropy coefficient (defined here as minsqrt(Ez/Ek)) at large Ri. 
-  ! From Zilitinkevich et al. 2013, it equals sqrt(0.03)~0.17
-  
+  ! From Zilitinkevich et al. 2013, it equals sqrt(0.03)~0.17  
   smmin=0.17
   CALL getin_p('atke_smmin',smmin)