trunk/dyn3d/gr_v_scal.f90


! $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)
  ! %W%    %G%
  ! =======================================================================

  ! Author:    Frederic Hourdin      original: 11/11/92
  ! -------

  ! Subject:
  ! ------

  ! Method:
  ! --------

  ! Interface:
  ! ----------

  ! Input:
  ! ------

  ! Output:
  ! -------

  ! =======================================================================
  USE dimens_m
  USE paramet_m
  USE comgeom
  IMPLICIT NONE
  ! -----------------------------------------------------------------------
  ! Declararations:
  ! ---------------


  ! Arguments:
  ! ----------

  INTEGER nx
  REAL x_v(ip1jm, nx), x_scal(ip1jmp1, nx)

  ! Local:
  ! ------

  INTEGER l, ij

  ! -----------------------------------------------------------------------

  DO l = 1, nx
    DO ij = iip2, ip1jm
      x_scal(ij, l) = (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l))/ &
        (airev(ij-iip1)+airev(ij))
    END DO
    DO ij = 1, iip1
      x_scal(ij, l) = 0.
    END DO
    DO ij = ip1jm + 1, ip1jmp1
      x_scal(ij, l) = 0.
    END DO
  END DO

  RETURN
END SUBROUTINE gr_v_scal
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 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