95 |
! flux de vapeur d'eau (kg / m2 / s) à la surface |
! flux de vapeur d'eau (kg / m2 / s) à la surface |
96 |
|
|
97 |
REAL, intent(out):: flux_u(klon, nbsrf), flux_v(klon, nbsrf) |
REAL, intent(out):: flux_u(klon, nbsrf), flux_v(klon, nbsrf) |
98 |
! tension du vent à la surface, en Pa |
! tension du vent (flux turbulent de vent) à la surface, en Pa |
99 |
|
|
100 |
REAL, INTENT(out):: cdragh(klon), cdragm(klon) |
REAL, INTENT(out):: cdragh(klon), cdragm(klon) |
101 |
real q2(klon, klev + 1, nbsrf) |
real q2(klon, klev + 1, nbsrf) |
369 |
END DO |
END DO |
370 |
|
|
371 |
ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), coefm(:knon, 1)) |
ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), coefm(:knon, 1)) |
|
|
|
|
! iflag_pbl peut \^etre utilis\'e comme longueur de m\'elange |
|
|
|
|
372 |
CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), & |
CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), & |
373 |
yu(:knon, :), yv(:knon, :), yteta(:knon, :), & |
yu(:knon, :), yv(:knon, :), yteta(:knon, :), & |
374 |
coefm(:knon, 1), yq2(:knon, :), ykmm(:knon, :), & |
coefm(:knon, 1), yq2(:knon, :), ykmm(:knon, :), & |
375 |
ykmn(:knon, :), ykmq(:knon, :), ustar(:knon), iflag_pbl) |
ykmn(:knon, :), ykmq(:knon, :), ustar(:knon)) |
|
|
|
376 |
coefm(:knon, 2:) = ykmm(:knon, 2:klev) |
coefm(:knon, 2:) = ykmm(:knon, 2:klev) |
377 |
coefh(:knon, 2:) = ykmn(:knon, 2:klev) |
coefh(:knon, 2:) = ykmn(:knon, 2:klev) |
378 |
END IF |
END IF |
379 |
|
|
380 |
! calculer la diffusion des vitesses "u" et "v" |
CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), & |
381 |
CALL clvent(knon, dtime, yu(:knon, 1), yv(:knon, 1), & |
yt(:knon, :), yu(:knon, :), ypaprs(:knon, :), & |
382 |
coefm(:knon, :), yt, yu, ypaprs, ypplay, ydelp, y_d_u, & |
ypplay(:knon, :), ydelp(:knon, :), y_d_u(:knon, :), & |
383 |
y_flux_u(:knon)) |
y_flux_u(:knon)) |
384 |
CALL clvent(knon, dtime, yu(:knon, 1), yv(:knon, 1), & |
CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), & |
385 |
coefm(:knon, :), yt, yv, ypaprs, ypplay, ydelp, y_d_v, & |
yt(:knon, :), yv(:knon, :), ypaprs(:knon, :), & |
386 |
|
ypplay(:knon, :), ydelp(:knon, :), y_d_v(:knon, :), & |
387 |
y_flux_v(:knon)) |
y_flux_v(:knon)) |
388 |
|
|
389 |
! calculer la diffusion de "q" et de "h" |
! calculer la diffusion de "q" et de "h" |