1 |
|
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gr_v_scal.F,v 1.1.1.1 2004/05/19 |
3 |
! 12:53:06 lmdzadmin Exp $ |
4 |
|
5 |
SUBROUTINE gr_v_scal(nx, x_v, x_scal) |
6 |
! %W% %G% |
7 |
! ======================================================================= |
8 |
|
9 |
! Author: Frederic Hourdin original: 11/11/92 |
10 |
! ------- |
11 |
|
12 |
! Subject: |
13 |
! ------ |
14 |
|
15 |
! Method: |
16 |
! -------- |
17 |
|
18 |
! Interface: |
19 |
! ---------- |
20 |
|
21 |
! Input: |
22 |
! ------ |
23 |
|
24 |
! Output: |
25 |
! ------- |
26 |
|
27 |
! ======================================================================= |
28 |
USE dimensions |
29 |
USE paramet_m |
30 |
USE comgeom |
31 |
IMPLICIT NONE |
32 |
! ----------------------------------------------------------------------- |
33 |
! Declararations: |
34 |
! --------------- |
35 |
|
36 |
|
37 |
! Arguments: |
38 |
! ---------- |
39 |
|
40 |
INTEGER nx |
41 |
REAL x_v(ip1jm, nx), x_scal(ip1jmp1, nx) |
42 |
|
43 |
! Local: |
44 |
! ------ |
45 |
|
46 |
INTEGER l, ij |
47 |
|
48 |
! ----------------------------------------------------------------------- |
49 |
|
50 |
DO l = 1, nx |
51 |
DO ij = iip2, ip1jm |
52 |
x_scal(ij, l) = (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l))/ & |
53 |
(airev(ij-iip1)+airev(ij)) |
54 |
END DO |
55 |
DO ij = 1, iip1 |
56 |
x_scal(ij, l) = 0. |
57 |
END DO |
58 |
DO ij = ip1jm + 1, ip1jmp1 |
59 |
x_scal(ij, l) = 0. |
60 |
END DO |
61 |
END DO |
62 |
|
63 |
RETURN |
64 |
END SUBROUTINE gr_v_scal |