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

Timestamp:

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

Author:

musat

Message:

Ajout ustar pour diagnostics pbl

File:

• LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/stdlevvar_mod.F90 (modified) (4 diffs)

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