source: LMDZ5/branches/LF-private/libf/dyn3dmem/geopot_loc.F @ 5018

Last change on this file since 5018 was 1823, checked in by Ehouarn Millour, 11 years ago

Remplacement de parallel.F90 (en conflit avec orchidée) par parallel_lmdz.F90.
UG
.........................................
Renaming parallel.F90 (conflicting with orchidée) into parallel_lmdz.F90.
UG

File size: 1.7 KB
Line 
1      SUBROUTINE geopot_loc ( ngrid, teta, pk, pks, phis, phi )
2      USE parallel_lmdz
3      IMPLICIT NONE
4     
5     
6c=======================================================================
7c
8c   Auteur:  P. Le Van
9c   -------
10c
11c   Objet:
12c   ------
13c
14c    *******************************************************************
15c    ....   calcul du geopotentiel aux milieux des couches    .....
16c    *******************************************************************
17c
18c     ....   l'integration se fait de bas en haut  ....
19c
20c     .. ngrid,teta,pk,pks,phis sont des argum. d'entree pour le s-pg ..
21c              phi               est un  argum. de sortie pour le s-pg .
22c
23c=======================================================================
24c-----------------------------------------------------------------------
25c   Declarations:
26c   -------------
27
28#include "dimensions.h"
29#include "paramet.h"
30#include "comvert.h"
31
32c   Arguments:
33c   ----------
34      INTEGER ngrid
35      REAL teta(ijb_u:ije_u,llm),pks(ijb_u:ije_u),phis(ijb_u:ije_u),
36     *     pk(ijb_u:ije_u,llm) , phi(ijb_u:ije_u,llm)
37
38
39c   Local:
40c   ------
41     
42      INTEGER  l, ij,ijb,ije
43
44
45c-----------------------------------------------------------------------
46c     calcul de phi au niveau 1 pres du sol  .....
47      ijb=ij_begin
48      ije=ij_end+iip1
49     
50      IF (pole_sud)  ije=ij_end
51
52      DO  ij  = ijb, ije
53      phi( ij,1 ) = phis( ij ) + teta(ij,1) * ( pks(ij) - pk(ij,1) )
54      ENDDO
55
56c     calcul de phi aux niveaux superieurs  .......
57
58      DO  l = 2,llm
59        DO  ij    = ijb,ije
60        phi(ij,l) = phi(ij,l-1) + 0.5 * ( teta(ij,l)  + teta(ij,l-1) )
61     *                              *   (  pk(ij,l-1) -  pk(ij,l)    )
62        ENDDO
63      ENDDO
64
65      RETURN
66      END
Note: See TracBrowser for help on using the repository browser.