310 |
END DO |
END DO |
311 |
END DO |
END DO |
312 |
|
|
|
! calculer Cdrag et les coefficients d'echange |
|
313 |
CALL coefkz(nsrf, ypaprs, ypplay, ksta, ksta_ter, yts(:knon), & |
CALL coefkz(nsrf, ypaprs, ypplay, ksta, ksta_ter, yts(:knon), & |
314 |
yrugos, yu, yv, yt, yq, yqsurf(:knon), coefm(:knon, :), & |
yrugos, yu, yv, yt, yq, yqsurf(:knon), coefm(:knon, :), & |
315 |
coefh(:knon, :), ycdragm(:knon), ycdragh(:knon)) |
coefh(:knon, :), ycdragm(:knon), ycdragh(:knon)) |
415 |
y_dflux_q(j) = y_dflux_q(j) * ypct(j) |
y_dflux_q(j) = y_dflux_q(j) * ypct(j) |
416 |
END DO |
END DO |
417 |
|
|
|
DO k = 2, klev |
|
|
DO j = 1, knon |
|
|
i = ni(j) |
|
|
coefh(j, k) = coefh(j, k) * ypct(j) |
|
|
coefm(j, k) = coefm(j, k) * ypct(j) |
|
|
END DO |
|
|
END DO |
|
|
DO j = 1, knon |
|
|
i = ni(j) |
|
|
ycdragh(j) = ycdragh(j) * ypct(j) |
|
|
ycdragm(j) = ycdragm(j) * ypct(j) |
|
|
END DO |
|
418 |
DO k = 1, klev |
DO k = 1, klev |
419 |
DO j = 1, knon |
DO j = 1, knon |
420 |
i = ni(j) |
i = ni(j) |
451 |
agesno(i, nsrf) = yagesno(j) |
agesno(i, nsrf) = yagesno(j) |
452 |
fqcalving(i, nsrf) = y_fqcalving(j) |
fqcalving(i, nsrf) = y_fqcalving(j) |
453 |
ffonte(i, nsrf) = y_ffonte(j) |
ffonte(i, nsrf) = y_ffonte(j) |
454 |
cdragh(i) = cdragh(i) + ycdragh(j) |
cdragh(i) = cdragh(i) + ycdragh(j) * ypct(j) |
455 |
cdragm(i) = cdragm(i) + ycdragm(j) |
cdragm(i) = cdragm(i) + ycdragm(j) * ypct(j) |
456 |
dflux_t(i) = dflux_t(i) + y_dflux_t(j) |
dflux_t(i) = dflux_t(i) + y_dflux_t(j) |
457 |
dflux_q(i) = dflux_q(i) + y_dflux_q(j) |
dflux_q(i) = dflux_q(i) + y_dflux_q(j) |
458 |
END DO |
END DO |
478 |
END DO |
END DO |
479 |
END DO |
END DO |
480 |
|
|
481 |
ycoefh(ni(:knon), :) = ycoefh(ni(:knon), :) + coefh(:knon, :) |
forall (k = 2:klev) ycoefh(ni(:knon), k) & |
482 |
|
= ycoefh(ni(:knon), k) + coefh(:knon, k) * ypct(:knon) |
483 |
|
|
484 |
! diagnostic t, q a 2m et u, v a 10m |
! diagnostic t, q a 2m et u, v a 10m |
485 |
|
|