1 |
! |
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/vitvert.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $ |
3 |
! |
4 |
SUBROUTINE vitvert ( convm , w ) |
5 |
c |
6 |
|
7 |
c======================================================================= |
8 |
c |
9 |
c Auteurs: P. Le Van , F. Hourdin . |
10 |
c ------- |
11 |
c |
12 |
c Objet: |
13 |
c ------ |
14 |
c |
15 |
c ******************************************************************* |
16 |
c .... calcul de la vitesse verticale aux niveaux sigma .... |
17 |
c ******************************************************************* |
18 |
c convm est un argument d'entree pour le s-pg ...... |
19 |
c w est un argument de sortie pour le s-pg ...... |
20 |
c |
21 |
c la vitesse verticale est orientee de haut en bas . |
22 |
c au sol, au niveau sigma(1), w(i,j,1) = 0. |
23 |
c au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi |
24 |
c egale a 0. et n'est pas stockee dans le tableau w . |
25 |
c |
26 |
c |
27 |
c======================================================================= |
28 |
|
29 |
use dimens_m |
30 |
use paramet_m |
31 |
use comvert |
32 |
IMPLICIT NONE |
33 |
|
34 |
REAL w(ip1jmp1,llm),convm(ip1jmp1,llm) |
35 |
INTEGER l, ij |
36 |
|
37 |
|
38 |
|
39 |
DO 2 l = 1,llmm1 |
40 |
|
41 |
DO 1 ij = 1,ip1jmp1 |
42 |
w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 ) |
43 |
1 CONTINUE |
44 |
|
45 |
2 CONTINUE |
46 |
|
47 |
DO 5 ij = 1,ip1jmp1 |
48 |
w(ij,1) = 0. |
49 |
5 CONTINUE |
50 |
|
51 |
RETURN |
52 |
END |