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

Timestamp:

Sep 11, 2023, 9:56:55 PM (13 months ago)

Author:

evignon

Message:

je renomme les routines atke suivant la convention decidee
pour la reecriture

File:

• LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90 (moved) (moved from LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90) (6 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

@@ -1 @@
-module atke_exchange_coeff_mod
+module lmdz_atke_exchange_coeff
 
 implicit none
@@ -26 +26 @@
 
 
-USE atke_turbulence_ini_mod, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual
-USE atke_turbulence_ini_mod, ONLY : cepsilon, pr_slope, pr_asym, pr_neut, ctkes,rg, rd, rv, atke_ok_vdiff
-USE atke_turbulence_ini_mod, ONLY : viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin
+USE lmdz_atke_turbulence_ini, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual
+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
 
 implicit none
@@ -182 +182 @@
              lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(N2(igrid,ilay))
              lstrat=max(lstrat,lmin)
-             !Monin-Obukhov consistent interpolation, Van de Wiel et al. 2010
-             l_exchange(igrid,ilay)=(1./(2.*l_exchange(igrid,ilay))+sqrt(1./(4.*l_exchange(igrid,ilay) &
-                                   *l_exchange(igrid,ilay))+1./(2.*lstrat*lstrat)))**(-1.0)
+             !Inverse quadratic interpolation, Van de Wiel et al. 2010
+             l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
           ENDIF 
       ENDDO
@@ -198 +197 @@
                     clmixshear*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay)))
              lstrat=max(lstrat,lmin)
-             !Monin-Obukhov consistent interpolation, Van de Wiel et al. 2010   
-             l_exchange(igrid,ilay)=(1./(2.*l_exchange(igrid,ilay))+sqrt(1./(4.*l_exchange(igrid,ilay) &
-                                   *l_exchange(igrid,ilay))+1./(2.*lstrat*lstrat)))**(-1.0)
-
+             !Inverse quadratic interpolation, Van de Wiel et al. 2010   
+             l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
           ENDIF
       ENDDO
@@ -356 +353 @@
 ! E Vignon, July 2023
 
-USE atke_turbulence_ini_mod, ONLY : rd, cke, viscom
+USE lmdz_atke_turbulence_ini, ONLY : rd, cke, viscom
 
 
@@ -449 +446 @@
 
 
-end module atke_exchange_coeff_mod
+end module lmdz_atke_exchange_coeff