Changeset 3838 for LMDZ6

Timestamp:

Feb 11, 2021, 10:39:04 AM (3 years ago)

Author:

musat

Message:

Ajout ustar pour diagnostics pbl

Location:

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd

Files:

• pbl_surface_mod.F90 (modified) (3 diffs)
• stdlevvar_mod.F90 (modified) (4 diffs)

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/pbl_surface_mod.F90

@@ -1839 +1839 @@
                yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
                yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
-               yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
+               yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
+               yn2mout(:, nsrf, :))
           ELSE
           CALL stdlevvar(klon, knon, is_ter, zxli, &
@@ -2638 +2639 @@
             uzon, vmer, tair1, qair1, zgeo1, &
             tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
+            yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
+            yn2mout(:, nsrf, :))
         ELSE
         CALL stdlevvar(klon, knon, nsrf, zxli, &
@@ -2650 +2652 @@
             uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
             tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :))
+            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, &
+            yn2mout_x(:, nsrf, :))
          CALL stdlevvarn(klon, knon, nsrf, zxli, &
             uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, &
             tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :))
+            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, &
+            yn2mout_w(:, nsrf, :))
         ELSE
         CALL stdlevvar(klon, knon, nsrf, zxli, &

LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/stdlevvar_mod.F90

@@ -311 +311 @@
                            u1, v1, t1, q1, z1, &
                            ts1, qsurf, z0m, z0h, psol, pat1, &
-                           t_2m, q_2m, t_10m, q_10m, u_10m, n2mout)
+                           t_2m, q_2m, t_10m, q_10m, u_10m, ustar, &
+                           n2mout)
 !
       USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -348 +349 @@
 ! u_2m--output-R- vitesse du vent a 2m
 ! u_10m--output-R- vitesse du vent a 10m
+! ustar--output-R- u*
 !AM
 ! t_10m--output-R- temperature de l'air a 10m
@@ -358 +360 @@
       REAL, dimension(klon), intent(in) :: psol, pat1
 !
-      REAL, dimension(klon), intent(out) :: t_2m, q_2m
+      REAL, dimension(klon), intent(out) :: t_2m, q_2m, ustar
       REAL, dimension(klon), intent(out) :: u_10m, t_10m, q_10m
       INTEGER, dimension(klon, 6), intent(out) :: n2mout
@@ -449 +451 @@
         tpot(i) = t1(i)* (psol(i)/pat1(i))**RKAPPA
         zdu2 = MAX(CEPDUE*CEPDUE, speed(i)**2)
+        ustar(i) = sqrt(cdram(i) * zdu2)
 !
       ENDDO