source: trunk/LMDZ.GENERIC/libf/dyn3d/dudv2.F @ 1422

Last change on this file since 1422 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.7 KB
Line 
1      SUBROUTINE dudv2 ( teta, pkf, bern, du, dv  )
2
3      IMPLICIT NONE
4c
5c=======================================================================
6c
7c   Auteur:  P. Le Van
8c   -------
9c
10c   Objet:
11c   ------
12c
13c   *****************************************************************
14c   ..... calcul du terme de pression (gradient de p/densite )   et
15c          du terme de ( -gradient de la fonction de Bernouilli ) ...
16c   *****************************************************************
17c          Ces termes sont ajoutes a  d(ucov)/dt et a d(vcov)/dt  ..
18c
19c
20c    teta , pkf, bern  sont des arguments d'entree  pour le s-pg  ....
21c    du et dv          sont des arguments de sortie pour le s-pg  ....
22c
23c=======================================================================
24c
25#include "dimensions.h"
26#include "paramet.h"
27
28      REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ),
29     *         du( ip1jmp1,llm ),  dv( ip1jm,llm )
30      INTEGER  l,ij
31c
32c
33      DO 5 l = 1,llm
34c
35      DO 2  ij  = iip2, ip1jm - 1
36       du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
37     * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
38   2  CONTINUE
39c
40c
41c    .....  correction  pour du(iip1,j,l),  j=2,jjm   ......
42c    ...          du(iip1,j,l) = du(1,j,l)                 ...
43c
44CDIR$ IVDEP
45      DO 3 ij = iip1+ iip1, ip1jm, iip1
46      du( ij,l ) = du( ij - iim,l )
47   3  CONTINUE
48c
49c
50      DO 4 ij  = 1,ip1jm
51      dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
52     *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
53     *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
54   4  CONTINUE
55c
56   5  CONTINUE
57c
58      RETURN
59      END
Note: See TracBrowser for help on using the repository browser.