1 |
! |
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gr_v_scal.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $ |
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 |
use dimens_m |
28 |
use paramet_m |
29 |
use comgeom |
30 |
IMPLICIT NONE |
31 |
c----------------------------------------------------------------------- |
32 |
c Declararations: |
33 |
c --------------- |
34 |
|
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 |