Custom Query (96 matches)

cycle_diurne et iflag_cycle_diurne : la première est lue dans conf_phys et pas utilisée, la seconde en dur dans physiq_mod

je me suis donc aperçu cet été (2017) que la correction de la ALP des wakes, pour prendre en compte la probabilité de présence des fronts de rafales dans la maille considérée, était absente. Je ne sais absolument quand cette disparition s'est produite. Si ça se trouve ça s'est fait au moment de l'introduction des poches dans LMDZ.

Jean-Yves

Rapport de possible bug dans yamada4 identifié par Yann dans le prototype dynamico-aquaplanète:

Produit d'un svn diff dans corrections dans yamada4.F90:

@@ -246,7 +246,7 @@
    IF (iflag_pbl==6) THEN

      DO k = 2, klev
-      q2(:, k) = l(:, k)**2*zz(:, k)
+      q2(:ngrid, k) = l(:ngrid, k)**2*zz(:ngrid, k)
      END DO


@@ -339,12 +339,14 @@
      ! print*,'Schema mixte D'
      ! print*,'Longueur ',l(:,:)
      DO k = 2, klev - 1
-      l(:, k) = max(l(:,k), 1.)
-      km(:, k) = l(:, k)*sqrt(q2(:,k))*sm(:, k)
-      q2(:, k) = q2(:, k) + dt*km(:, k)*m2(:, k)*(1.-rif(:,k))
-      q2(:, k) = min(max(q2(:,k),1.E-10), 1.E4)
-      q2(:, k) = 1./(1./sqrt(q2(:,k))+dt/(2*l(:,k)*b1))
-      q2(:, k) = q2(:, k)*q2(:, k)
+      DO ig = 1, ngrid
+        l(ig, k) = max(l(ig,k), 1.)
+        km(ig, k) = l(ig, k)*sqrt(q2(ig,k))*sm(ig, k)
+        q2(ig, k) = q2(ig, k) + dt*km(ig, k)*m2(ig, k)*(1.-rif(ig,k))
+        q2(ig, k) = min(max(q2(ig,k),1.E-10), 1.E4)
+        q2(ig, k) = 1./(1./sqrt(q2(ig,k))+dt/(2*l(ig,k)*b1))
+        q2(ig, k) = q2(ig, k)*q2(ig, k)
+      END DO
      END DO


@@ -376,7 +378,7 @@
    ! Transport diffusif vertical de la TKE.
    IF (iflag_pbl>=12) THEN
      ! print*,'YAMADA VDIF'
-    q2(:, 1) = q2(:, 2)
+    q2(1:ngrid, 1) = q2(1:ngrid, 2)
      CALL vdif_q2(dt, g, rconst, ngrid, plev, temp, kq, q2)
    END IF