trunk/dyn3d/massbarxy.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/massbarxy.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
!
      SUBROUTINE massbarxy(  masse, massebxy )
c
c **********************************************************************
c
c  Calcule les moyennes en x et  y de la masse d'air dans chaque maille.
c **********************************************************************
c    Auteurs : P. Le Van , Fr. Hourdin  .
c   ..........
c
c  ..  masse          est  un  argum. d'entree  pour le s-pg ...
c  ..  massebxy       est  un  argum. de sortie pour le s-pg ...
c     
c
c     IMPLICIT NONE
c
      use dimens_m
      use paramet_m
      use comconst
      use comgeom
c
       REAL, intent(in)::  masse( ip1jmp1,llm )
       real massebxy( ip1jm,llm )
c

      DO   100    l = 1 , llm
c
      DO 5 ij = 1, ip1jm - 1
      massebxy( ij,l ) = masse(    ij  ,l ) * alpha2(   ij    )   +
     +                   masse(   ij+1 ,l ) * alpha3(  ij+1   )   +
     +                   masse( ij+iip1,l ) * alpha1( ij+iip1 )   +
     +                   masse( ij+iip2,l ) * alpha4( ij+iip2 )
   5  CONTINUE

c    ....  correction pour     massebxy( iip1,j )  ........

CDIR$ IVDEP

      DO 7 ij = iip1, ip1jm, iip1
      massebxy( ij,l ) = massebxy( ij - iim,l )
   7  CONTINUE

100   CONTINUE
c
      RETURN
      END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/massbarxy.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			SUBROUTINE massbarxy( masse, massebxy )
5			c
6			c **********************************************************************
7			c
8			c Calcule les moyennes en x et y de la masse d'air dans chaque maille.
9			c **********************************************************************
10			c Auteurs : P. Le Van , Fr. Hourdin .
11			c ..........
12			c
13			c .. masse est un argum. d'entree pour le s-pg ...
14			c .. massebxy est un argum. de sortie pour le s-pg ...
15			c
16			c
17			c IMPLICIT NONE
18			c
19			use dimens_m
20			use paramet_m
21			use comconst
22			use comgeom
23			c
24	guez	69	REAL, intent(in):: masse( ip1jmp1,llm )
25			real massebxy( ip1jm,llm )
26	guez	3	c
27
28			DO 100 l = 1 , llm
29			c
30			DO 5 ij = 1, ip1jm - 1
31			massebxy( ij,l ) = masse( ij ,l ) * alpha2( ij ) +
32			+ masse( ij+1 ,l ) * alpha3( ij+1 ) +
33			+ masse( ij+iip1,l ) * alpha1( ij+iip1 ) +
34			+ masse( ij+iip2,l ) * alpha4( ij+iip2 )
35			5 CONTINUE
36
37			c .... correction pour massebxy( iip1,j ) ........
38
39			CDIR$ IVDEP
40
41			DO 7 ij = iip1, ip1jm, iip1
42			massebxy( ij,l ) = massebxy( ij - iim,l )
43			7 CONTINUE
44
45			100 CONTINUE
46			c
47			RETURN
48			END