Changeset 4688

Timestamp:

Sep 13, 2023, 8:36:25 PM (8 months ago)

Author:

evignon

Message:

correction formulation nombre de Prandlt dans ATKE pour atelier de demain

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

@@ -140 +140 @@
 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) / pr_slope
+Ri1 = -2./rpi * (pr_asym - pr_neut) 
 
 
@@ -157 +157 @@
         ELSE ! stable cases
             Sm(igrid,ilay) = max(smmin,cn*(1.-Ri(igrid,ilay)/Ric))
-            Prandtl(igrid,ilay) = pr_neut + Ri(igrid,ilay) * pr_slope
+            ! prandlt expression from venayagamoorthy and stretch 2010, Li et al 2019
+            Prandtl(igrid,ilay) = pr_neut*exp(-pr_slope/pr_neut*Ri(igrid,ilay)+Ri(igrid,ilay)/pr_neut) &
+                                + Ri(igrid,ilay) * pr_slope
             IF (Ri(igrid,ilay) .GE. Prandtl(igrid,ilay)) THEN
                call abort_physic("atke_compute_km_kh", &
@@ -182 +184 @@
              lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(N2(igrid,ilay))
              lstrat=max(lstrat,lmin)
-             !Inverse quadratic interpolation, Van de Wiel et al. 2010
+             !Inverse interpolation, Van de Wiel et al. 2010
              l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
           ENDIF 
@@ -197 +199 @@
                     clmixshear*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay)))
              lstrat=max(lstrat,lmin)
-             !Inverse quadratic interpolation, Van de Wiel et al. 2010   
+             !Inverse interpolation, Van de Wiel et al. 2010   
              l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
           ENDIF