source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/dudv2.F @ 420

Last change on this file since 420 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.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#include "comvert.h"
28
29      REAL teta( ip1jmp1,llm ),pkf( ip1jmp1,llm ) ,bern( ip1jmp1,llm ),
30     *         du( ip1jmp1,llm ),  dv( ip1jm,llm )
31      INTEGER  l,ij
32c
33c
34      DO 5 l = 1,llm
35c
36      DO 2  ij  = iip2, ip1jm - 1
37       du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
38     * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
39   2  CONTINUE
40c
41c
42c    .....  correction  pour du(iip1,j,l),  j=2,jjm   ......
43c    ...          du(iip1,j,l) = du(1,j,l)                 ...
44c
45CDIR$ IVDEP
46      DO 3 ij = iip1+ iip1, ip1jm, iip1
47      du( ij,l ) = du( ij - iim,l )
48   3  CONTINUE
49c
50c
51      DO 4 ij  = 1,ip1jm
52      dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
53     *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
54     *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
55   4  CONTINUE
56c
57   5  CONTINUE
58c
59      RETURN
60      END
Note: See TracBrowser for help on using the repository browser.