source: LMDZ5/branches/IPSLCM6.0.10/libf/dyn3d_common/pression.F90 @ 4085

Last change on this file since 4085 was 2408, checked in by Laurent Fairhead, 9 years ago

Merged trunk changes r2298:2396 into testing branch
File size: 1.2 KB
Line 
1SUBROUTINE pression( ngrid, ap, bp, ps, p )
2!
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!-------------------------------------------------------------------------------
10  IMPLICIT NONE
11  include "dimensions.h"
12  include "paramet.h"
13!===============================================================================
14! Arguments:
15  INTEGER, INTENT(IN)  :: ngrid                 !--- NUMBER OF GRID POINTS
16  REAL,    INTENT(IN)  :: ap(llmp1), bp(llmp1)  !--- HYBRID COEFFICIENTS
17  REAL,    INTENT(IN)  :: ps(ngrid)             !--- SURFACE PRESSURE
18  REAL,    INTENT(OUT) :: p(ngrid,llmp1)        !--- 3D PRESSURE FIELD
19!===============================================================================
20! Local variables:
21  INTEGER :: l
22!===============================================================================
23  DO l=1,llmp1;  p(:,l) = ap(l) + bp(l) * ps(:);  END DO
24
25END SUBROUTINE pression
26
27
Note: See TracBrowser for help on using the repository browser.