Changeset 3198 for LMDZ6/trunk/libf/phylmd/pbl_surface_mod.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/pbl_surface_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1294 +1294 @@
        ENDDO
 !!! jyg le 07/02/2012 et le 10/04/2013
-        DO k = 1, klev
+        DO k = 1, klev+1
           DO j = 1, knon
              i = ni(j)
@@ -1300 +1300 @@
 !!             ytke(j,k)   = tke(i,k,nsrf)
              ytke(j,k)   = tke_x(i,k,nsrf)
+          ENDDO
+        ENDDO
 !>jyg
+        DO k = 1, klev
+          DO j = 1, knon
+             i = ni(j)
 !FC
              y_treedrg(j,k) =  treedrg(i,k,nsrf)
@@ -2408 +2413 @@
        IF (iflag_split .eq.0) THEN
         wake_dltke(:,:,nsrf) = 0.
-        DO k = 1, klev
+        DO k = 1, klev+1
            DO j = 1, knon
               i = ni(j)
@@ -2421 +2426 @@
 
        ELSE  ! (iflag_split .eq.0)
-        DO k = 1, klev
+        DO k = 1, klev+1
           DO j = 1, knon
             i = ni(j)