|
! |
|
|
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gr_v_scal.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $ |
|
|
! |
|
|
SUBROUTINE gr_v_scal(nx,x_v,x_scal) |
|
|
c%W% %G% |
|
|
c======================================================================= |
|
|
c |
|
|
c Author: Frederic Hourdin original: 11/11/92 |
|
|
c ------- |
|
|
c |
|
|
c Subject: |
|
|
c ------ |
|
|
c |
|
|
c Method: |
|
|
c -------- |
|
|
c |
|
|
c Interface: |
|
|
c ---------- |
|
|
c |
|
|
c Input: |
|
|
c ------ |
|
|
c |
|
|
c Output: |
|
|
c ------- |
|
|
c |
|
|
c======================================================================= |
|
|
use dimens_m |
|
|
use paramet_m |
|
|
use comgeom |
|
|
IMPLICIT NONE |
|
|
c----------------------------------------------------------------------- |
|
|
c Declararations: |
|
|
c --------------- |
|
|
|
|
|
|
|
|
c Arguments: |
|
|
c ---------- |
|
|
|
|
|
INTEGER nx |
|
|
REAL x_v(ip1jm,nx),x_scal(ip1jmp1,nx) |
|
|
|
|
|
c Local: |
|
|
c ------ |
|
|
|
|
|
INTEGER l,ij |
|
|
|
|
|
c----------------------------------------------------------------------- |
|
|
|
|
|
DO l=1,nx |
|
|
DO ij=iip2,ip1jm |
|
|
x_scal(ij,l)= |
|
|
s (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) |
|
|
s /(airev(ij-iip1)+airev(ij)) |
|
|
ENDDO |
|
|
DO ij=1,iip1 |
|
|
x_scal(ij,l)=0. |
|
|
ENDDO |
|
|
DO ij=ip1jm+1,ip1jmp1 |
|
|
x_scal(ij,l)=0. |
|
|
ENDDO |
|
|
ENDDO |
|
1 |
|
|
2 |
RETURN |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/gr_v_scal.F,v 1.1.1.1 2004/05/19 |
3 |
END |
! 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 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_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 |