1 |
guez |
3 |
|
2 |
guez |
81 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gr_u_scal.F,v 1.1.1.1 2004/05/19 |
3 |
|
|
! 12:53:06 lmdzadmin Exp $ |
4 |
guez |
3 |
|
5 |
guez |
81 |
SUBROUTINE gr_u_scal(nx, x_u, x_scal) |
6 |
|
|
! %W% %G% |
7 |
|
|
! ======================================================================= |
8 |
guez |
3 |
|
9 |
guez |
81 |
! Author: Frederic Hourdin original: 11/11/92 |
10 |
|
|
! ------- |
11 |
guez |
3 |
|
12 |
guez |
81 |
! Subject: |
13 |
|
|
! ------ |
14 |
guez |
3 |
|
15 |
guez |
81 |
! Method: |
16 |
|
|
! -------- |
17 |
guez |
3 |
|
18 |
guez |
81 |
! Interface: |
19 |
|
|
! ---------- |
20 |
guez |
3 |
|
21 |
guez |
81 |
! Input: |
22 |
|
|
! ------ |
23 |
guez |
3 |
|
24 |
guez |
81 |
! Output: |
25 |
|
|
! ------- |
26 |
guez |
3 |
|
27 |
guez |
81 |
! ======================================================================= |
28 |
|
|
USE dimens_m |
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_u(ip1jmp1, 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 = ip1jmp1, 2, -1 |
52 |
|
|
x_scal(ij, l) = (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l))/ & |
53 |
|
|
(aireu(ij)+aireu(ij-1)) |
54 |
|
|
END DO |
55 |
|
|
END DO |
56 |
|
|
|
57 |
|
|
CALL scopy(nx*jjp1, x_scal(iip1,1), iip1, x_scal(1,1), iip1) |
58 |
|
|
|
59 |
|
|
RETURN |
60 |
|
|
END SUBROUTINE gr_u_scal |