source: trunk/LMDZ.COMMON/libf/dyn3d_common/massbarxy.F @ 1453

Last change on this file since 1453 was 1422, checked in by milmd, 10 years ago

In GENERIC, MARS and COMMON models replace some include files by modules (usefull for decoupling physics with dynamics).

File size: 1.2 KB
RevLine 
[1]1!
2! $Header$
3!
4      SUBROUTINE massbarxy(  masse, massebxy )
[1391]5      IMPLICIT NONE
[1]6c
7c **********************************************************************
8c
9c  Calcule les moyennes en x et  y de la masse d'air dans chaque maille.
10c **********************************************************************
11c    Auteurs : P. Le Van , Fr. Hourdin  .
12c   ..........
13c
14c  ..  masse          est  un  argum. d'entree  pour le s-pg ...
15c  ..  massebxy       est  un  argum. de sortie pour le s-pg ...
16c     
17c
18c     IMPLICIT NONE
19c
20#include "dimensions.h"
21#include "paramet.h"
22#include "comgeom.h"
23c
24       REAL  masse( ip1jmp1,llm ), massebxy( ip1jm,llm )
[1391]25       INTEGER ij,l
[1]26c
27
28      DO   100    l = 1 , llm
29c
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
37c    ....  correction pour     massebxy( iip1,j )  ........
38
39CDIR$ IVDEP
40
41      DO 7 ij = iip1, ip1jm, iip1
42      massebxy( ij,l ) = massebxy( ij - iim,l )
43   7  CONTINUE
44
45100   CONTINUE
46c
47      RETURN
48      END
Note: See TracBrowser for help on using the repository browser.