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

Last change on this file since 5118 was 5099, checked in by abarral, 4 months ago

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

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.