[2293] | 1 | SUBROUTINE 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] | 11 | USE paramet_mod_h |
---|
[5271] | 12 | IMPLICIT 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 | |
---|
| 27 | END SUBROUTINE pression |
---|
| 28 | |
---|
| 29 | |
---|