source: LMDZ5/trunk/libf/dyn3dmem/geopot_loc.F @ 1973

Last change on this file since 1973 was 1907, checked in by lguez, 12 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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
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.