Changeset 4780 for LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

Timestamp:

Dec 21, 2023, 10:17:25 PM (9 months ago)

Author:

evignon

Message:

petites modifs dans les routines atke pour etre en accord avec les notations du papier

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

@@ -25 +25 @@
 
 
-USE lmdz_atke_turbulence_ini, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual
+USE lmdz_atke_turbulence_ini, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual, ri0, ri1
 USE lmdz_atke_turbulence_ini, ONLY : cepsilon, pr_slope, pr_asym, pr_neut, ctkes, rg, rd, rv, atke_ok_vdiff
-USE lmdz_atke_turbulence_ini, ONLY : viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin
+USE lmdz_atke_turbulence_ini, ONLY : viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin, cn
 
 implicit none
@@ -69 +69 @@
 
 INTEGER :: igrid,ilay ! horizontal,vertical index (flat grid)
-REAL    :: cn,Ri0,Ri1    ! parameter for Sm stability function and Prandlt
 REAL    :: preff      ! reference pressure for potential temperature calculations
 REAL    :: thetam     ! mean potential temperature at interface
@@ -132 +131 @@
 ! Computes the gradient Richardson's number and stability functions:
 !===================================================================
-
-! calculation of cn = Sm value at Ri=0
-! direct dependance on cepsilon to guarantee Fm=1 (first-order like stability function) at Ri=0
-cn=(1./sqrt(cepsilon))**(2/3)
-! calculation of Ri0 such that continuity in slope of Sm at Ri=0
-Ri0=2./rpi*(cinf - cn)*ric/cn
-! calculation of Ri1 to guarantee continuity in slope of Prandlt number at Ri=0
-Ri1 = -2./rpi * (pr_asym - pr_neut) 
-
 
 DO ilay=2,nlay
@@ -176 +166 @@
 
 IF (iflag_atke_lmix .EQ. 1 ) THEN
-
+! Blackadar formulation (~kappa l) + buoyancy length scale (Deardoff 1980) for very stable conditions
    DO ilay=2,nlay
       DO igrid=1,ngrid
@@ -223 +213 @@
 
 ELSE
-! default: neglect effect of local stratification and shear
+! default Blackadar formulation: neglect effect of local stratification and shear
 
    DO ilay=2,nlay+1
@@ -250 +240 @@
 
 ! full implicit scheme resolved with a second order polynomial equation
-
+! default solution which shows fair convergence properties
     DO ilay=2,nlay
         DO igrid=1,ngrid