106 |
! dflux_q derive du flux latent |
! dflux_q derive du flux latent |
107 |
! IM "slab" ocean |
! IM "slab" ocean |
108 |
|
|
109 |
REAL, intent(out):: ycoefh(:, :) ! (klon, klev) |
REAL, intent(out):: ycoefh(:, 2:) ! (klon, 2:klev) |
110 |
! Pour pouvoir extraire les coefficients d'\'echange, le champ |
! Pour pouvoir extraire les coefficients d'\'echange, le champ |
111 |
! "ycoefh" a \'et\'e cr\'e\'e. Nous avons moyenn\'e les valeurs de |
! "ycoefh" a \'et\'e cr\'e\'e. Nous avons moyenn\'e les valeurs de |
112 |
! ce champ sur les quatre sous-surfaces du mod\`ele. |
! ce champ sur les quatre sous-surfaces du mod\`ele. |
490 |
d_v(i, k) = d_v(i, k) + y_d_v(j, k) |
d_v(i, k) = d_v(i, k) + y_d_v(j, k) |
491 |
END DO |
END DO |
492 |
END DO |
END DO |
|
|
|
|
DO j = 1, knon |
|
|
i = ni(j) |
|
|
DO k = 2, klev |
|
|
ycoefh(i, k) = ycoefh(i, k) + coefh(j, k) |
|
|
END DO |
|
|
END DO |
|
493 |
|
|
494 |
DO j = 1, knon |
ycoefh(ni(:knon), :) = ycoefh(ni(:knon), :) + coefh(:knon, :) |
|
i = ni(j) |
|
|
ycoefh(i, 1) = ycoefh(i, 1) + ycdragh(j) |
|
|
END DO |
|
495 |
|
|
496 |
! diagnostic t, q a 2m et u, v a 10m |
! diagnostic t, q a 2m et u, v a 10m |
497 |
|
|