source: LMDZ5/trunk/libf/dyn3d_common/massbarxy.F @ 2100

Last change on this file since 2100 was 1945, checked in by lguez, 11 years ago

Duplicated files moved to dyn3d_common.
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE massbarxy(  masse, massebxy )
5c
6c **********************************************************************
7c
8c  Calcule les moyennes en x et  y de la masse d'air dans chaque maille.
9c **********************************************************************
10c    Auteurs : P. Le Van , Fr. Hourdin  .
11c   ..........
12c
13c  ..  masse          est  un  argum. d'entree  pour le s-pg ...
14c  ..  massebxy       est  un  argum. de sortie pour le s-pg ...
15c     
16c
17c     IMPLICIT NONE
18c
19#include "dimensions.h"
20#include "paramet.h"
21#include "comconst.h"
22#include "comgeom.h"
23c
24       REAL  masse( ip1jmp1,llm ), massebxy( ip1jm,llm )
25c
26
27      DO   100    l = 1 , llm
28c
29      DO 5 ij = 1, ip1jm - 1
30      massebxy( ij,l ) = masse(    ij  ,l ) * alpha2(   ij    )   +
31     +                   masse(   ij+1 ,l ) * alpha3(  ij+1   )   +
32     +                   masse( ij+iip1,l ) * alpha1( ij+iip1 )   +
33     +                   masse( ij+iip2,l ) * alpha4( ij+iip2 )
34   5  CONTINUE
35
36c    ....  correction pour     massebxy( iip1,j )  ........
37
38CDIR$ IVDEP
39
40      DO 7 ij = iip1, ip1jm, iip1
41      massebxy( ij,l ) = massebxy( ij - iim,l )
42   7  CONTINUE
43
44100   CONTINUE
45c
46      RETURN
47      END
Note: See TracBrowser for help on using the repository browser.