Changeset 3839 for LMDZ6/trunk/libf/phylmd
- Timestamp:
- Feb 11, 2021, 10:48:55 AM (4 years ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
r3819 r3839 1876 1876 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1877 1877 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1878 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 1878 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 1879 yn2mout(:, nsrf, :)) 1879 1880 ELSE 1880 1881 CALL stdlevvar(klon, knon, is_ter, zxli, & … … 2712 2713 uzon, vmer, tair1, qair1, zgeo1, & 2713 2714 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2714 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 2715 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 2716 yn2mout(:, nsrf, :)) 2715 2717 ELSE 2716 2718 CALL stdlevvar(klon, knon, nsrf, zxli, & … … 2724 2726 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & 2725 2727 tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2726 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :)) 2728 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, & 2729 yn2mout_x(:, nsrf, :)) 2727 2730 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2728 2731 uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, & 2729 2732 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2730 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :)) 2733 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, & 2734 yn2mout_w(:, nsrf, :)) 2731 2735 ELSE 2732 2736 CALL stdlevvar(klon, knon, nsrf, zxli, & -
LMDZ6/trunk/libf/phylmd/stdlevvar_mod.F90
r3831 r3839 311 311 u1, v1, t1, q1, z1, & 312 312 ts1, qsurf, z0m, z0h, psol, pat1, & 313 t_2m, q_2m, t_10m, q_10m, u_10m, n2mout) 313 t_2m, q_2m, t_10m, q_10m, u_10m, ustar, & 314 n2mout) 314 315 ! 315 316 USE ioipsl_getin_p_mod, ONLY : getin_p … … 348 349 ! u_2m--output-R- vitesse du vent a 2m 349 350 ! u_10m--output-R- vitesse du vent a 10m 351 ! ustar--output-R- u* 350 352 !AM 351 353 ! t_10m--output-R- temperature de l'air a 10m … … 358 360 REAL, dimension(klon), intent(in) :: psol, pat1 359 361 ! 360 REAL, dimension(klon), intent(out) :: t_2m, q_2m 362 REAL, dimension(klon), intent(out) :: t_2m, q_2m, ustar 361 363 REAL, dimension(klon), intent(out) :: u_10m, t_10m, q_10m 362 364 INTEGER, dimension(klon, 6), intent(out) :: n2mout … … 449 451 tpot(i) = t1(i)* (psol(i)/pat1(i))**RKAPPA 450 452 zdu2 = MAX(CEPDUE*CEPDUE, speed(i)**2) 453 ustar(i) = sqrt(cdram(i) * zdu2) 451 454 ! 452 455 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.