Custom Query (96 matches)

C'est la version de LMDZ qui doit être incluse dans IPSLCM6.0.5

Définition:

• envisager d’introduire une résistance à l’évaporation (dû à sur-estimation de l’évaporation en été voir CR réunion LMDZ6 9 mai 2016 )
• s'attaquer au probleme de l'intermittence des precip
• régler le plantage de l'epluchage en B/w2.
• bascule XIOS2

Prévoir une nouvelle testing

Bug remonté des tests avec dynamico-aquaplanète de Yann. A vérifier, mais il semble qu'il manque dans physiq_mod.F90 les initialisations de mass_solu_aero() et mass_solu_aero_pi()...

Proposition de la modification:

Ligne 3301+ :

ELSE

tausum_aero(:,:,:) = 0.

+ mass_solu_aero(:,:)=0 + mass_solu_aero_pi(:,:)=0

IF (iflag_rrtm .EQ. 0) THEN !--old radiation

tau_aero(:,:,:,:) = 1.e-15 piz_aero(:,:,:,:) = 1. cg_aero(:,:,:,:) = 0.

ELSE

tau_aero_sw_rrtm(:,:,:,:) = 1.e-15 tau_aero_lw_rrtm(:,:,:,:) = 1.e-15 piz_aero_sw_rrtm(:,:,:,:) = 1.0 cg_aero_sw_rrtm(:,:,:,:) = 0.0

ENDIF

ENDIF

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