Changeset 3198 for LMDZ6/trunk/libf/phylmd/tend_to_tke.F90

Timestamp:

Feb 12, 2018, 1:24:03 AM (7 years ago)

Author:

fhourdin

Message:

Retour vers l'insensibilite au decoupage en sous domaine.
Les routines gwd_rando incluait le calcul de niveaux de reference
sur la base d'un profile pris au milieu du domaine (en klon/2).
Rempace par un test en presnivs.

Une autre intercation entre routines concernant la tke a fait apparaitre
que la tke n'était pas passee correctement au niveau klev+1 au moment
du regroupement des mailles sous les sous surface.

Ces changements garantissent la convergence numerique si
addtkeoro=0
iflag_pbl<12
et
ok_gwd_rando=n
La convergence n'est pas garantie pour les dernieres versions des physiq.def
mais les differences devraient etre mineures.

File:

• LMDZ6/trunk/libf/phylmd/tend_to_tke.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/tend_to_tke.F90

@@ -32 +32 @@
 !**************************************************************************************
 
- SUBROUTINE tend_to_tke(dt,plev,exner,temp,windu,windv,dt_a,du_a,dv_a,tke)
+ SUBROUTINE tend_to_tke(dt,plev,exner,temp,windu,windv,dt_a,du_a,dv_a,pctsrf,tke)
 
  USE dimphy, ONLY: klon, klev
@@ -53 +53 @@
   REAL du_a(klon,klev)      ! Zonal wind speed tendency [m/s], grid-cell average or for a one subsurface
   REAL dv_a(klon,klev)      ! Meridional wind speed tendency [m/s], grid-cell average or for a one subsurface
+  REAL pctsrf(klon,nbsrf+1)       ! Turbulent Kinetic energy [m2/s2], grid-cell average or for a subsurface
 
 ! Inputs/Outputs
@@ -119 +120 @@
 
 
- DO isrf=1,nbsrf
+ DO isrf=1,nsrf
     DO k=1,klev
-       tke(:,k,isrf)= tke(:,k,isrf)+tendu(:,k)+tendv(:,k)+tendt(:,k)
-       tke(:,k,isrf)= max(tke(:,k,isrf),1.e-10)
+       DO i=1,klon
+          IF (pctsrf(i,isrf)>0.) THEN
+            tke(i,k,isrf)= tke(i,k,isrf)+tendu(i,k)+tendv(i,k)+tendt(i,k)
+            tke(i,k,isrf)= max(tke(i,k,isrf),1.e-10)
+          ENDIF
+       ENDDO
     ENDDO
  ENDDO
-
-! dtke_t(:,:)=tendt(:,:)
-! dtke_u(:,:)=tendu(:,:)
-! dtke_v(:,:)=tendv(:,:)