Changeset 3200 for LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Feb 12, 2018, 10:01:04 AM (6 years ago)

Author:

Laurent Fairhead

Message:

Inclusion of r3198 from trunk
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'?\195?\169tait 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.

FH

File:

• LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/pbl_surface_mod.F90 (modified) (4 diffs)

LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/pbl_surface_mod.F90

@@ -1266 +1266 @@
        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)
@@ -1272 +1272 @@
 !!             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)
@@ -2398 +2403 @@
        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)
@@ -2411 +2416 @@
 
        ELSE  ! (iflag_split .eq.0)
-        DO k = 1, klev
+        DO k = 1, klev+1
           DO j = 1, knon
             i = ni(j)