Changeset 2897 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

May 31, 2017, 12:34:09 AM (7 years ago)

Author:

fhourdin

Message:

Introduction d'une possible prise en compte de génération de TKE
par les ondes de relief.
Etienne Vignon et FH

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -69 +69 @@
        d_t_oli,d_u_oli,d_v_oli, &
        d_t_oro,d_u_oro,d_v_oro, &
+       d_t_oro_gw,d_u_oro_gw,d_v_oro_gw, &
        d_t_lif,d_u_lif,d_v_lif, &
        d_t_ec, &
@@ -434 +435 @@
     REAL d_qx(klon,klev,nqtot)
     REAL d_ps(klon)
+  ! variables pour tend_to_tke
+    REAL duadd(klon,klev)
+    REAL dvadd(klon,klev)
+    REAL dtadd(klon,klev)
+
     ! Variables pour le transport convectif
     real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
@@ -606 +612 @@
     REAL, SAVE :: alp_offset
     !$OMP THREADPRIVATE(alp_offset)
-
+  
     !
     !RR:fin declarations poches froides
@@ -640 +646 @@
     real therm_tke_max(klon,klev)   ! Profil de TKE dans les thermiques
     real env_tke_max(klon,klev)     ! Profil de TKE dans l'environnement
+
+    !-------Activer les tendances de TKE due a l'orograp??ie---------
+     INTEGER, SAVE :: addtkeoro
+    !$OMP THREADPRIVATE(addtkeoro) 
+     REAL, SAVE :: alphatkeoro
+    !$OMP THREADPRIVATE(alphatkeoro) 
+     LOGICAL, SAVE :: smallscales_tkeoro
+    !$OMP THREADPRIVATE(smallscales_tkeoro) 
+
 
 
@@ -780 +795 @@
     real zqsat(klon,klev)
     !
-    INTEGER i, k, iq, nsrf, l
+    INTEGER i, k, iq, j, nsrf, ll, l
     !
     REAL t_coup
@@ -2244 +2259 @@
 !>jyg
        ENDIF
+
+
+
 
 
@@ -3870 +3888 @@
        IF (ok_strato) THEN
 
-          CALL drag_noro_strato(klon,klev,dtime,paprs,pplay, &
+          CALL drag_noro_strato(0,klon,klev,dtime,paprs,pplay, &
                zmea,zstd, zsig, zgam, zthe,zpic,zval, &
                igwd,idx,itest, &
@@ -4045 +4063 @@
     !
     !
+
+!===============================================================
+!            Additional tendency of TKE due to orography
+!===============================================================
+!
+! Inititialization
+!------------------
+
+    
+
+       addtkeoro=0   
+       CALL getin_p('addtkeoro',addtkeoro) 
+      
+       IF (prt_level.ge.5) &
+            print*,'addtkeoro', addtkeoro
+            
+       alphatkeoro=1.   
+       CALL getin_p('alphatkeoro',alphatkeoro)
+       alphatkeoro=min(max(0.,alphatkeoro),1.)
+
+       smallscales_tkeoro=.false.   
+       CALL getin_p('smallscales_tkeoro',smallscales_tkeoro) 
+
+
+        dtadd(:,:)=0.
+        duadd(:,:)=0.
+        dvadd(:,:)=0.
+
+
+
+! Choices for addtkeoro:
+!      ** 0 no TKE tendency from orography    
+!      ** 1 we include a fraction alphatkeoro of the whole tendency duoro
+!      ** 2 we include a fraction alphatkeoro of the gravity wave part of duoro
+!
+
+       IF (addtkeoro .GT. 0 .AND. ok_orodr ) THEN
+!      -------------------------------------------
+
+
+       !  selection des points pour lesquels le schema est actif:
+
+
+
+  IF (addtkeoro .EQ. 1 ) THEN
+
+            duadd(:,:)=alphatkeoro*d_u_oro(:,:)
+            dvadd(:,:)=alphatkeoro*d_v_oro(:,:)
+
+  ELSE IF (addtkeoro .EQ. 2) THEN
+
+
+
+       IF (smallscales_tkeoro) THEN
+       igwd=0
+       DO i=1,klon
+          itest(i)=0
+! Etienne: ici je prends en compte plus de relief que la routine drag_noro_strato
+! car on peut s'attendre a ce que les petites echelles produisent aussi de la TKE
+! Mais attention, cela ne va pas dans le sens de la conservation de l'energie! 
+          IF (zstd(i).GT.1.0) THEN
+             itest(i)=1
+             igwd=igwd+1
+             idx(igwd)=i
+          ENDIF
+       ENDDO
+
+     ELSE 
+
+       igwd=0
+       DO i=1,klon
+          itest(i)=0
+        IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN
+             itest(i)=1
+             igwd=igwd+1
+             idx(igwd)=i
+          ENDIF
+       ENDDO
+
+       END IF
+
+
+
+
+       CALL drag_noro_strato(addtkeoro,klon,klev,dtime,paprs,pplay, &
+               zmea,zstd, zsig, zgam, zthe,zpic,zval, &
+               igwd,idx,itest, &
+               t_seri, u_seri, v_seri, &
+               zulow, zvlow, zustrdr, zvstrdr, &
+               d_t_oro_gw, d_u_oro_gw, d_v_oro_gw)
+
+            zustrdr(:)=0.
+            zvstrdr(:)=0.
+            zulow(:)=0.
+            zvlow(:)=0.
+
+            duadd(:,:)=alphatkeoro*d_u_oro_gw(:,:)
+            dvadd(:,:)=alphatkeoro*d_v_oro_gw(:,:)
+ END IF
+    
+
+
+   ! TKE update from subgrid temperature and wind tendencies
+   !----------------------------------------------------------
+    forall (k=1: nbp_lev) exner(:, k) = (pplay(:, k)/paprs(:,1))**RKAPPA
+
+
+    CALL tend_to_tke(pdtphys,paprs,exner,t_seri,u_seri,v_seri,dtadd,duadd,dvadd,pbl_tke)
+
+
+
+       ENDIF
+!      -----
+!===============================================================
+
+
+
     !====================================================================
     ! Interface Simulateur COSP (Calipso, ISCCP, MISR, ..)