184 |
REAL cdragh(klon) ! drag coefficient pour T and Q |
REAL cdragh(klon) ! drag coefficient pour T and Q |
185 |
REAL cdragm(klon) ! drag coefficient pour vent |
REAL cdragm(klon) ! drag coefficient pour vent |
186 |
|
|
187 |
REAL ycoefh(klon, llm) ! coef d'echange pour phytrac |
REAL coefh(klon, 2:llm) ! coef d'echange pour phytrac |
188 |
|
|
189 |
REAL, save:: ffonte(klon, nbsrf) |
REAL, save:: ffonte(klon, nbsrf) |
190 |
! flux thermique utilise pour fondre la neige |
! flux thermique utilise pour fondre la neige |
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. |
571 |
paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, & |
paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, & |
572 |
snow_fall, fsolsw, fsollw, frugs, agesno, rugoro, d_t_vdf, d_q_vdf, & |
snow_fall, fsolsw, fsollw, frugs, agesno, rugoro, d_t_vdf, d_q_vdf, & |
573 |
d_u_vdf, d_v_vdf, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, & |
d_u_vdf, d_v_vdf, d_ts, flux_t, flux_q, flux_u, flux_v, cdragh, & |
574 |
cdragm, q2, dsens, devap, ycoefh, t2m, q2m, u10m_srf, v10m_srf, & |
cdragm, q2, dsens, devap, coefh, t2m, q2m, u10m_srf, v10m_srf, & |
575 |
pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, trmb3, plcl, & |
pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, trmb3, plcl, & |
576 |
fqcalving, ffonte, run_off_lic_0) |
fqcalving, ffonte, run_off_lic_0) |
577 |
|
|
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, & |
955 |
mfd, pde_u, pen_d, ycoefh, fm_therm, entr_therm, u(:, 1), v(:, 1), & |
mfd, pde_u, pen_d, coefh, cdragh, fm_therm, entr_therm, u(:, 1), & |
956 |
ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, & |
v(:, 1), ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, & |
957 |
tr_seri, zmasse, ncid_startphy) |
dnwd, tr_seri, zmasse, ncid_startphy) |
958 |
|
|
959 |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
! Calculer le transport de l'eau et de l'energie (diagnostique) |
960 |
CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq) |
CALL transp(paprs, t_seri, q_seri, u_seri, v_seri, zphi, ve, vq, ue, uq) |