My Project
 All Classes Files Functions Variables Macros
dudv2.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE dudv2 ( teta, pkf, bern, du, dv )
5 
6  IMPLICIT NONE
7 c
8 c=======================================================================
9 c
10 c Auteur: P. Le Van
11 c -------
12 c
13 c Objet:
14 c ------
15 c
16 c *****************************************************************
17 c ..... calcul du terme de pression (gradient de p/densite ) et
18 c du terme de ( -gradient de la fonction de Bernouilli ) ...
19 c *****************************************************************
20 c Ces termes sont ajoutes a d(ucov)/dt et a d(vcov)/dt ..
21 c
22 c
23 c teta , pkf, bern sont des arguments d'entree pour le s-pg ....
24 c du et dv sont des arguments de sortie pour le s-pg ....
25 c
26 c=======================================================================
27 c
28 #include "dimensions.h"
29 #include "paramet.h"
30 #include "comvert.h"
31 
32  REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ),
33  * du( ip1jmp1,llm ), dv( ip1jm,llm )
34  INTEGER l,ij
35 c
36 c
37  DO 5 l = 1,llm
38 c
39  DO 2 ij = iip2, ip1jm - 1
40  du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
41  * ( pkf( ij,l ) - pkf(ij+1,l) ) + bern(ij,l) - bern(ij+1,l)
42  2 CONTINUE
43 c
44 c
45 c ..... correction pour du(iip1,j,l), j=2,jjm ......
46 c ... du(iip1,j,l) = du(1,j,l) ...
47 c
48 CDIR$ IVDEP
49  DO 3 ij = iip1+ iip1, ip1jm, iip1
50  du( ij,l ) = du( ij - iim,l )
51  3 CONTINUE
52 c
53 c
54  DO 4 ij = 1,ip1jm
55  dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
56  * ( pkf(ij+iip1,l) - pkf( ij,l ) )
57  * + bern( ij+iip1,l ) - bern( ij ,l )
58  4 CONTINUE
59 c
60  5 CONTINUE
61 c
62  RETURN
63  END