source: LMDZ6/trunk/libf/dyn3d_common/pression.f90 @ 5456

Last change on this file since 5456 was 5285, checked in by abarral, 2 months ago

As discussed internally, remove generic ONLY: ... for new _mod_h modules

File size: 1.3 KB
RevLine 
[2293]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!-------------------------------------------------------------------------------
[5271]10  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5285]11USE paramet_mod_h
[5271]12IMPLICIT NONE
13
[5272]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.