My Project
 All Classes Files Functions Variables Macros
gr_u_scal_p.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE gr_u_scal_p(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(ip1jmp1,nx),x_scal(ip1jmp1,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  DO l=1,nx
54  DO ij=ijb+1,ije
55  x_scal(ij,l)=
56  s(aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l))
57  s /(aireu(ij)+aireu(ij-1))
58  ENDDO
59  ENDDO
60 
61 cym CALL SCOPY(nx*jjp1,x_scal(iip1,1),iip1,x_scal(1,1),iip1)
62  ijb=ij_begin
63  ije=ij_end
64 
65  DO l=1,nx
66  DO ij=ijb,ije-iip1+1,iip1
67  x_scal(ij,l)=x_scal(ij+iip1-1,l)
68  ENDDO
69  ENDDO
70  RETURN
71 
72  END