source: LMDZ.3.3/trunk/libf/dyn3d/massbarxy.F @ 507

Last change on this file since 507 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1      SUBROUTINE massbarxy(  masse, massebxy )
2c
3c **********************************************************************
4c
5c  Calcule les moyennes en x et  y de la masse d'air dans chaque maille.
6c **********************************************************************
7c    Auteurs : P. Le Van , Fr. Hourdin  .
8c   ..........
9c
10c  ..  masse          est  un  argum. d'entree  pour le s-pg ...
11c  ..  massebxy       est  un  argum. de sortie pour le s-pg ...
12c     
13c
14c     IMPLICIT NONE
15c
16#include "dimensions.h"
17#include "paramet.h"
18#include "comconst.h"
19#include "comgeom.h"
20c
21       REAL  masse( ip1jmp1,llm ), massebxy( ip1jm,llm )
22c
23
24      DO   100    l = 1 , llm
25c
26      DO 5 ij = 1, ip1jm - 1
27      massebxy( ij,l ) = masse(    ij  ,l ) * alpha2(   ij    )   +
28     +                   masse(   ij+1 ,l ) * alpha3(  ij+1   )   +
29     +                   masse( ij+iip1,l ) * alpha1( ij+iip1 )   +
30     +                   masse( ij+iip2,l ) * alpha4( ij+iip2 )
31   5  CONTINUE
32
33c    ....  correction pour     massebxy( iip1,j )  ........
34
35CDIR$ IVDEP
36
37      DO 7 ij = iip1, ip1jm, iip1
38      massebxy( ij,l ) = massebxy( ij - iim,l )
39   7  CONTINUE
40
41100   CONTINUE
42c
43      RETURN
44      END
Note: See TracBrowser for help on using the repository browser.