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