1 |
|
|
! |
2 |
|
|
! $Header$ |
3 |
|
|
! |
4 |
|
|
SUBROUTINE gr_u_scal(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 |
|
|
IMPLICIT NONE |
28 |
|
|
c----------------------------------------------------------------------- |
29 |
|
|
c Declararations: |
30 |
|
|
c --------------- |
31 |
|
|
|
32 |
|
|
INCLUDE "dimensions.h" |
33 |
|
|
INCLUDE "paramet.h" |
34 |
|
|
INCLUDE "comgeom.h" |
35 |
|
|
|
36 |
|
|
c Arguments: |
37 |
|
|
c ---------- |
38 |
|
|
|
39 |
|
|
INTEGER nx |
40 |
|
|
REAL x_u(ip1jmp1,nx),x_scal(ip1jmp1,nx) |
41 |
|
|
|
42 |
|
|
c Local: |
43 |
|
|
c ------ |
44 |
|
|
|
45 |
|
|
INTEGER l,ij |
46 |
|
|
|
47 |
|
|
c----------------------------------------------------------------------- |
48 |
|
|
|
49 |
|
|
DO l=1,nx |
50 |
|
|
DO ij=ip1jmp1,2,-1 |
51 |
|
|
x_scal(ij,l)= |
52 |
|
|
s (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l)) |
53 |
|
|
s /(aireu(ij)+aireu(ij-1)) |
54 |
|
|
ENDDO |
55 |
|
|
ENDDO |
56 |
|
|
|
57 |
|
|
CALL SCOPY(nx*jjp1,x_scal(iip1,1),iip1,x_scal(1,1),iip1) |
58 |
|
|
|
59 |
|
|
RETURN |
60 |
|
|
END |