LMDZ
covcont_loc.F
Go to the documentation of this file.
1  SUBROUTINE covcont_loc (klevel,ucov, vcov, ucont, vcont )
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( ijb_u:ije_u,klevel ), vcov( ijb_v:ije_v,klevel )
25  REAL ucont( ijb_u:ije_u,klevel ), vcont( ijb_v:ije_v,klevel )
26  INTEGER l,ij
27  INTEGER ijbu,ijbv,ijeu,ijev
28 
29 
30  ijbu=ij_begin-iip1
31  ijbv=ij_begin-iip1
32  ijeu=ij_end+iip1
33  ijev=ij_end+iip1
34 
35  if (pole_nord) then
36  ijbu=ij_begin+iip1
37  ijbv=ij_begin
38  endif
39 
40  if (pole_sud) then
41  ijeu=ij_end-iip1
42  ijev=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
integer, save ij_end
logical, save pole_sud
!$Header!CDK comgeom COMMON comgeom unscu2
Definition: comgeom.h:25
integer, save ijb_v
!$Header!CDK comgeom COMMON comgeom unscv2
Definition: comgeom.h:25
logical, save pole_nord
integer, save ij_begin
integer, save ije_v
subroutine covcont_loc(klevel, ucov, vcov, ucont, vcont)
Definition: covcont_loc.F:2
integer, save ije_u
integer, save ijb_u