48 |
|
|
49 |
knon = size(ypaprs, 1) |
knon = size(ypaprs, 1) |
50 |
|
|
51 |
! Debut de la partie qui doit etre incluse a terme dans clmain. |
! Debut de la partie qui doit etre incluse a terme dans pbl_surface. |
52 |
|
|
53 |
do i = 1, knon |
do i = 1, knon |
54 |
zlay(i, 1) = RD * yt(i, 1) * 2 / (ypaprs(i, 1) + ypplay(i, 1)) & |
zlay(i, 1) = RD * yt(i, 1) * 2 / (ypaprs(i, 1) + ypplay(i, 1)) & |
74 |
forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1)) |
forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1)) |
75 |
ustar = SQRT(ycoefm * (yu(:, 1)**2 + yv(:, 1)**2)) |
ustar = SQRT(ycoefm * (yu(:, 1)**2 + yv(:, 1)**2)) |
76 |
|
|
77 |
! Fin de la partie qui doit être incluse à terme dans clmain |
! Fin de la partie qui doit être incluse à terme dans pbl_surface |
78 |
|
|
79 |
! Cette routine est ecrite pour avoir en entree ustar, teta et zlev |
! Cette routine est ecrite pour avoir en entree ustar, teta et zlev |
80 |
! Ici, on a inclus le calcul de ces trois variables dans la routine |
! Ici, on a inclus le calcul de ces trois variables dans la routine |