1 |
guez |
91 |
module vitvert_m |
2 |
guez |
3 |
|
3 |
guez |
91 |
IMPLICIT NONE |
4 |
guez |
3 |
|
5 |
guez |
91 |
contains |
6 |
guez |
3 |
|
7 |
guez |
91 |
SUBROUTINE vitvert(convm, w) |
8 |
guez |
3 |
|
9 |
guez |
91 |
! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05 |
10 |
|
|
! Authors: P. Le Van, F. Hourdin |
11 |
guez |
3 |
|
12 |
guez |
91 |
! Objet : calcul de la vitesse verticale aux niveaux sigma |
13 |
guez |
3 |
|
14 |
guez |
91 |
! La vitesse verticale est orientée de haut en bas. Au sol, au |
15 |
|
|
! niveau sigma(1), w(i, j, 1) = 0. Au sommet, au niveau |
16 |
|
|
! sigma(llm+1), la vitesse verticale est aussi égale à 0 et n'est |
17 |
|
|
! pas stockée dans le tableau w. |
18 |
guez |
3 |
|
19 |
guez |
91 |
USE dimens_m, ONLY : llm |
20 |
|
|
USE disvert_m, ONLY : bp |
21 |
|
|
USE paramet_m, ONLY : ip1jmp1 |
22 |
guez |
3 |
|
23 |
guez |
91 |
real, intent(in):: convm(ip1jmp1, llm) |
24 |
|
|
REAL, intent(out):: w(ip1jmp1, llm) |
25 |
guez |
3 |
|
26 |
guez |
91 |
! Local: |
27 |
|
|
INTEGER l |
28 |
guez |
3 |
|
29 |
guez |
91 |
!------------------------------------------------------ |
30 |
|
|
|
31 |
|
|
forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1) |
32 |
|
|
w(:, 1) = 0. |
33 |
|
|
|
34 |
|
|
END SUBROUTINE vitvert |
35 |
|
|
|
36 |
|
|
end module vitvert_m |