Changeset 4653 for LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_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_exchange_coeff_mod.F90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90

@@ -6 +6 @@
 
 subroutine atke_compute_km_kh(ngrid,nlay,dtime, &
-                        wind_u,wind_v,temp,play,pinterf,cdrag_uv, &
+                        wind_u,wind_v,temp,qvap,play,pinterf,cdrag_uv, &
                         tke,Km_out,Kh_out)
 
@@ -26 +26 @@
 
 
-USE atke_turbulence_ini_mod, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, ok_vdiff_tke
-USE atke_turbulence_ini_mod, ONLY : cepsilon, pr_slope, pr_asym, pr_neut, ctkes,rg, rd
+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, iflag_atke_lmix, lmin, smmin
 
@@ -43 +43 @@
 REAL, DIMENSION(ngrid,nlay), INTENT(IN)       :: wind_v   ! meridional velocity (m/s)
 REAL, DIMENSION(ngrid,nlay), INTENT(IN)       :: temp   ! temperature (K)
+REAL, DIMENSION(ngrid,nlay), INTENT(IN)       :: qvap   ! specific humidity (kg/kg)
 REAL, DIMENSION(ngrid,nlay), INTENT(IN)       :: play   ! pressure (Pa)
 REAL, DIMENSION(ngrid,nlay+1), INTENT(IN)     :: pinterf   ! pressure at interfaces(Pa)
@@ -81 +82 @@
 REAL    :: netsource  ! net source term of tke
 REAL    :: ustar      ! friction velocity estimation
+REAL    :: invtau     
+REAL    :: rvap
 
 ! Initializations:
@@ -90 +93 @@
 END DO
 
-! Calculation of potential temperature: (if vapor -> todo virtual potential temperature)
+! Calculation of potential temperature: (if vapor -> virtual potential temperature)
 !=====================================
 
 preff=100000.
-! The result should not depend on the choice of preff
+! results should not depend on the choice of preff
 DO ilay=1,nlay
      DO igrid = 1, ngrid
@@ -101 +104 @@
 END DO
 
+! account for water vapor mass for buoyancy calculation
+IF (atke_ok_virtual) THEN
+  DO ilay=1,nlay
+     DO igrid = 1, ngrid
+        rvap=max(0.,qvap(igrid,ilay)/(1.-qvap(igrid,ilay)))
+        theta(igrid,ilay)=theta(igrid,ilay)*(1.+rvap/(RD/RV))/(1.+rvap)
+     END DO
+  END DO
+ENDIF
 
 
@@ -287 +299 @@
 
 
+ELSE IF (iflag_atke == 5) THEN
+! semi implicit scheme from Arpege (V. Masson methodology with 
+! Taylor expansion of the dissipation term)
+! and implicit resolution when switch num (when we use the mixing length as a function of tke)
+    DO ilay=2,nlay
+        DO igrid=1,ngrid
+           qq=max(sqrt(2.*tke(igrid,ilay)),1.e-10)
+           if (switch_num(igrid,ilay) .and. N2(igrid,ilay)>0.) then
+             invtau=clmix*Sm(igrid,ilay)/sqrt(N2(igrid,ilay))*shear2(igrid,ilay)*(1.-Ri(igrid,ilay)/Prandtl(igrid,ilay))  & 
+                   -sqrt(N2(igrid,ilay))/(cepsilon*clmix)
+             !qq=qq/(1.-dtime*invtau)
+             !qq=qq*exp(dtime*invtau)
+             tke(igrid,ilay)=tke(igrid,ilay)/(1.-dtime*invtau)
+             tke(igrid,ilay)=max(0.,tke(igrid,ilay))
+           else
+             qq=(l_exchange(igrid,ilay)*Sm(igrid,ilay)/sqrt(2.)*shear2(igrid,ilay)*(1.-Ri(igrid,ilay)/Prandtl(igrid,ilay)) &
+             +qq*(1.+dtime*qq/(cepsilon*l_exchange(igrid,ilay)*2.*sqrt(2.)))) &
+             /(1.+2.*qq*dtime/(cepsilon*l_exchange(igrid,ilay)*2.*sqrt(2.)))
+             qq=max(0.,qq)
+             tke(igrid,ilay)=0.5*(qq**2)
+           endif
+        ENDDO
+    ENDDO
+
+
 ELSE
    call abort_physic("atke_compute_km_kh", &
@@ -312 +349 @@
 ! vertical diffusion of TKE 
 !==========================
-IF (ok_vdiff_tke) THEN
+IF (atke_ok_vdiff) THEN
    CALL atke_vdiff_tke(ngrid,nlay,dtime,z_lay,z_interf,temp,play,l_exchange,Sm,tke)
 ENDIF
@@ -338 +375 @@
 
 ! routine that computes the vertical diffusion of TKE by the turbulence
-! using an implicit resolution (See note by Dufresne and Ghattas
+! using an implicit resolution (See note by Dufresne and Ghattas (2009))
 ! E Vignon, July 2023