source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pression.F90 @ 5209

Last change on this file since 5209 was 5159, checked in by abarral, 3 months ago

Put dimensions.h and paramet.h into modules

File size: 1.3 KB
RevLine 
[2293]1SUBROUTINE pression( ngrid, ap, bp, ps, p )
[5099]2
[2293]3!-------------------------------------------------------------------------------
4! Authors: P. Le Van , Fr.Hourdin
5!-------------------------------------------------------------------------------
6! Purpose: Compute pressure p(l) at different levels from l = 1 (ground level)
7!          to l = llm +1. Those levels correspond to the llm layers interfaces,
8!          with p(ij,llm+1) = 0. and  p(ij,1) = ps(ij)  .   
9!-------------------------------------------------------------------------------
[5159]10USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
11  USE lmdz_paramet
[2293]12  IMPLICIT NONE
[5159]13
14
[2293]15!===============================================================================
16! Arguments:
17  INTEGER, INTENT(IN)  :: ngrid                 !--- NUMBER OF GRID POINTS
18  REAL,    INTENT(IN)  :: ap(llmp1), bp(llmp1)  !--- HYBRID COEFFICIENTS
19  REAL,    INTENT(IN)  :: ps(ngrid)             !--- SURFACE PRESSURE
20  REAL,    INTENT(OUT) :: p(ngrid,llmp1)        !--- 3D PRESSURE FIELD
21!===============================================================================
22! Local variables:
23  INTEGER :: l
24!===============================================================================
25  DO l=1,llmp1;  p(:,l) = ap(l) + bp(l) * ps(:);  END DO
26
27END SUBROUTINE pression
28
29
Note: See TracBrowser for help on using the repository browser.