My Project
 All Classes Files Functions Variables Macros
gr_u_scal_loc.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE gr_u_scal_loc(nx,x_u,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_u(ijb_u:ije_u,nx),x_scal(ijb_u:ije_u,nx)
42 
43 c Local:
44 c ------
45 
46  INTEGER l,ij
47  INTEGER :: ijb,ije
48 
49 c-----------------------------------------------------------------------
50  ijb=ij_begin
51  ije=ij_end
52 
53 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
54  DO l=1,nx
55  DO ij=ijb+1,ije
56  x_scal(ij,l)=
57  s(aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l))
58  s /(aireu(ij)+aireu(ij-1))
59  ENDDO
60  ENDDO
61 !$OMP ENDDO NOWAIT
62 
63  ijb=ij_begin
64  ije=ij_end
65 
66 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
67  DO l=1,nx
68  DO ij=ijb,ije-iip1+1,iip1
69  x_scal(ij,l)=x_scal(ij+iip1-1,l)
70  ENDDO
71  ENDDO
72 !$OMP ENDDO NOWAIT
73  RETURN
74 
75  END