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

Last change on this file since 2306 was 2197, checked in by Ehouarn Millour, 10 years ago

Added 'implicit none' statements and proper variable definitions where they were missing.
EM

  • 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 )
5      IMPLICIT NONE
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 "comconst.h"
23#include "comgeom.h"
24c
25       REAL  masse( ip1jmp1,llm ), massebxy( ip1jm,llm )
26       INTEGER ij,l
27c
28
29      DO   100    l = 1 , llm
30c
31      DO 5 ij = 1, ip1jm - 1
32      massebxy( ij,l ) = masse(    ij  ,l ) * alpha2(   ij    )   +
33     +                   masse(   ij+1 ,l ) * alpha3(  ij+1   )   +
34     +                   masse( ij+iip1,l ) * alpha1( ij+iip1 )   +
35     +                   masse( ij+iip2,l ) * alpha4( ij+iip2 )
36   5  CONTINUE
37
38c    ....  correction pour     massebxy( iip1,j )  ........
39
40CDIR$ IVDEP
41
42      DO 7 ij = iip1, ip1jm, iip1
43      massebxy( ij,l ) = massebxy( ij - iim,l )
44   7  CONTINUE
45
46100   CONTINUE
47c
48      RETURN
49      END
Note: See TracBrowser for help on using the repository browser.