248 |
|
|
249 |
REAL flux_q(klon, nbsrf) ! flux turbulent d'humidite à la surface |
REAL flux_q(klon, nbsrf) ! flux turbulent d'humidite à la surface |
250 |
REAL flux_t(klon, nbsrf) ! flux turbulent de chaleur à la surface |
REAL flux_t(klon, nbsrf) ! flux turbulent de chaleur à la surface |
251 |
REAL flux_u(klon, nbsrf) ! flux turbulent de vitesse u à la surface |
|
252 |
REAL flux_v(klon, nbsrf) ! flux turbulent de vitesse v à la surface |
REAL flux_u(klon, nbsrf), flux_v(klon, nbsrf) |
253 |
|
! tension du vent (flux turbulent de vent) à la surface, en Pa |
254 |
|
|
255 |
! Le rayonnement n'est pas calcul\'e tous les pas, il faut donc que |
! Le rayonnement n'est pas calcul\'e tous les pas, il faut donc que |
256 |
! les variables soient r\'emanentes. |
! les variables soient r\'emanentes. |
370 |
|
|
371 |
REAL zustrdr(klon), zvstrdr(klon) |
REAL zustrdr(klon), zvstrdr(klon) |
372 |
REAL zustrli(klon), zvstrli(klon) |
REAL zustrli(klon), zvstrli(klon) |
|
REAL zustrph(klon), zvstrph(klon) |
|
373 |
REAL aam, torsfc |
REAL aam, torsfc |
374 |
|
|
375 |
REAL ve_lay(klon, llm) ! transport meri. de l'energie a chaque niveau vert. |
REAL ve_lay(klon, llm) ! transport meri. de l'energie a chaque niveau vert. |
945 |
ENDDO |
ENDDO |
946 |
ENDIF |
ENDIF |
947 |
|
|
948 |
! Stress n\'ecessaires : toute la physique |
CALL aaam_bud(rg, romega, pphis, zustrdr, zustrli, & |
949 |
|
sum((u_seri - u) / dtphys * zmasse, dim = 2), zvstrdr, & |
950 |
DO i = 1, klon |
zvstrli, sum((v_seri - v) / dtphys * zmasse, dim = 2), paprs, u, v, & |
951 |
zustrph(i) = 0. |
aam, torsfc) |
|
zvstrph(i) = 0. |
|
|
ENDDO |
|
|
DO k = 1, llm |
|
|
DO i = 1, klon |
|
|
zustrph(i) = zustrph(i) + (u_seri(i, k) - u(i, k)) / dtphys & |
|
|
* zmasse(i, k) |
|
|
zvstrph(i) = zvstrph(i) + (v_seri(i, k) - v(i, k)) / dtphys & |
|
|
* zmasse(i, k) |
|
|
ENDDO |
|
|
ENDDO |
|
|
|
|
|
CALL aaam_bud(rg, romega, pphis, zustrdr, zustrli, zustrph, zvstrdr, & |
|
|
zvstrli, zvstrph, paprs, u, v, aam, torsfc) |
|
952 |
|
|
953 |
! Calcul des tendances traceurs |
! Calcul des tendances traceurs |
954 |
call phytrac(julien, time, firstcal, lafin, dtphys, t, paprs, play, mfu, & |
call phytrac(julien, time, firstcal, lafin, dtphys, t, paprs, play, mfu, & |