My Project
 All Classes Files Functions Variables Macros
gr_v_scal_loc.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE gr_v_scal_loc(nx,x_v,x_scal)
5 c%W% %G%
6 c=======================================================================
7 c
8 c Author: Frederic Hourdin original: 11/11/92
9 c -------
10 c
11 c Subject:
12 c ------
13 c
14 c Method:
15 c --------
16 c
17 c Interface:
18 c ----------
19 c
20 c Input:
21 c ------
22 c
23 c Output:
24 c -------
25 c
26 c=======================================================================
27  USE parallel
28  IMPLICIT NONE
29 c-----------------------------------------------------------------------
30 c Declararations:
31 c ---------------
32 
33 #include "dimensions.h"
34 #include "paramet.h"
35 #include "comgeom.h"
36 
37 c Arguments:
38 c ----------
39 
40  INTEGER nx
41  REAL x_v(ijb_v:ije_v,nx),x_scal(ijb_v:ije_v,nx)
42 
43 c Local:
44 c ------
45 
46  INTEGER l,ij
47  INTEGER :: ijb,ije
48 c-----------------------------------------------------------------------
49  ijb=ij_begin
50  ije=ij_end
51  if (pole_nord) ijb=ij_begin+iip1
52  if (pole_sud) ije=ij_end-iip1
53 
54 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
55  DO l=1,nx
56  DO ij=ijb,ije
57  x_scal(ij,l)=
58  s(airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l))
59  s /(airev(ij-iip1)+airev(ij))
60  ENDDO
61  ENDDO
62 !$OMP ENDDO NOWAIT
63 
64  if (pole_nord) then
65 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
66  DO l=1,nx
67  DO ij=1,iip1
68  x_scal(ij,l)=0.
69  ENDDO
70  ENDDO
71 !$OMP ENDDO NOWAIT
72  endif
73 
74  if (pole_sud) then
75 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
76  DO l=1,nx
77  DO ij=ip1jm+1,ip1jmp1
78  x_scal(ij,l)=0.
79  ENDDO
80  ENDDO
81 !$OMP ENDDO NOWAIT
82  endif
83 
84  RETURN
85  END