My Project
 All Classes Files Functions Variables Macros
covcont_p.F
Go to the documentation of this file.
1  SUBROUTINE covcont_p (klevel,ucov, vcov, ucont, vcont )
2  USE parallel
3  IMPLICIT NONE
4 
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 des compos. contravariantes a partir des comp.covariantes
15 c ********************************************************************
16 c
17 c=======================================================================
18 
19 #include "dimensions.h"
20 #include "paramet.h"
21 #include "comgeom.h"
22 
23  INTEGER klevel
24  REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel )
25  REAL ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel )
26  INTEGER l,ij
27  INTEGER ijb_u,ijb_v,ije_u,ije_v
28 
29 
30  ijb_u=ij_begin-iip1
31  ijb_v=ij_begin-iip1
32  ije_u=ij_end+iip1
33  ije_v=ij_end+iip1
34 
35  if (pole_nord) then
36  ijb_u=ij_begin+iip1
37  ijb_v=ij_begin
38  endif
39 
40  if (pole_sud) then
41  ije_u=ij_end-iip1
42  ije_v=ij_end-iip1
43  endif
44 
45 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
46  DO 10 l = 1,klevel
47 
48  DO 2 ij = ijb_u,ije_u
49  ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
50  2 CONTINUE
51 
52  DO 4 ij = ijb_v,ije_v
53  vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
54  4 CONTINUE
55 
56  10 CONTINUE
57 c$OMP END DO NOWAIT
58  RETURN
59  END