1 |
|
|
! |
2 |
|
|
! $Header$ |
3 |
|
|
! |
4 |
|
|
SUBROUTINE gr_v_scal(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 |
|
|
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_v(ip1jm,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=iip2,ip1jm |
51 |
|
|
x_scal(ij,l)= |
52 |
|
|
s (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) |
53 |
|
|
s /(airev(ij-iip1)+airev(ij)) |
54 |
|
|
ENDDO |
55 |
|
|
DO ij=1,iip1 |
56 |
|
|
x_scal(ij,l)=0. |
57 |
|
|
ENDDO |
58 |
|
|
DO ij=ip1jm+1,ip1jmp1 |
59 |
|
|
x_scal(ij,l)=0. |
60 |
|
|
ENDDO |
61 |
|
|
ENDDO |
62 |
|
|
|
63 |
|
|
RETURN |
64 |
|
|
END |