Changeset 3838
- Timestamp:
- Feb 11, 2021, 10:39:04 AM (4 years ago)
- Location:
- LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/pbl_surface_mod.F90
r3825 r3838 1839 1839 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1840 1840 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1841 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 1841 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 1842 yn2mout(:, nsrf, :)) 1842 1843 ELSE 1843 1844 CALL stdlevvar(klon, knon, is_ter, zxli, & … … 2638 2639 uzon, vmer, tair1, qair1, zgeo1, & 2639 2640 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2640 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 2641 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 2642 yn2mout(:, nsrf, :)) 2641 2643 ELSE 2642 2644 CALL stdlevvar(klon, knon, nsrf, zxli, & … … 2650 2652 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & 2651 2653 tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2652 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :)) 2654 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, & 2655 yn2mout_x(:, nsrf, :)) 2653 2656 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2654 2657 uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, & 2655 2658 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2656 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :)) 2659 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, & 2660 yn2mout_w(:, nsrf, :)) 2657 2661 ELSE 2658 2662 CALL stdlevvar(klon, knon, nsrf, zxli, & -
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/stdlevvar_mod.F90
r3832 r3838 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.