source: LMDZ4/trunk/libf/dyn3dpar/massbarxy.F @ 1897

Last change on this file since 1897 was 774, checked in by Laurent Fairhead, 17 years ago

Suite du merge entre la version et la HEAD: quelques modifications de
Yann sur le

LF

  • 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.