572 |
paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, & |
paprs, play, fsnow, fqsurf, fevap, falbe, fluxlat, rain_fall, & |
573 |
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, & |
574 |
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, & |
575 |
cdragm, q2, dsens, devap, ycoefh, yu1, yv1, t2m, q2m, u10m_srf, & |
cdragm, q2, dsens, devap, ycoefh, t2m, q2m, u10m_srf, v10m_srf, & |
576 |
v10m_srf, pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, & |
pblh, capCL, oliqCL, cteiCL, pblT, therm, trmb1, trmb2, trmb3, plcl, & |
577 |
trmb3, plcl, fqcalving, ffonte, run_off_lic_0) |
fqcalving, ffonte, run_off_lic_0) |
578 |
|
|
579 |
! Incr\'ementation des flux |
! Incr\'ementation des flux |
580 |
|
|
582 |
evap = - sum(flux_q * pctsrf, dim = 2) |
evap = - sum(flux_q * pctsrf, dim = 2) |
583 |
fder = dlw + dsens + devap |
fder = dlw + dsens + devap |
584 |
|
|
585 |
|
yu1 = u_seri(:, 1) |
586 |
|
yv1 = v_seri(:, 1) |
587 |
|
|
588 |
DO k = 1, llm |
DO k = 1, llm |
589 |
DO i = 1, klon |
DO i = 1, klon |
590 |
t_seri(i, k) = t_seri(i, k) + d_t_vdf(i, k) |
t_seri(i, k) = t_seri(i, k) + d_t_vdf(i, k) |