12 |
|
|
13 |
USE comconst, ONLY: cpp, kappa |
USE comconst, ONLY: cpp, kappa |
14 |
USE comgeom, ONLY: aireu_2d, airev_2d, aire_2d |
USE comgeom, ONLY: aireu_2d, airev_2d, aire_2d |
15 |
USE dimens_m, ONLY: iim, jjm, llm |
USE dimens_m, ONLY: jjm, llm |
16 |
USE disvert_m, ONLY: ap, bp, preff |
USE disvert_m, ONLY: ap, bp, preff |
17 |
USE exner_hyb_m, ONLY: exner_hyb |
USE exner_hyb_m, ONLY: exner_hyb |
18 |
use massbar_m, only: massbar |
use massbar_m, only: massbar |
|
use massdair_m, only: massdair |
|
19 |
USE paramet_m, ONLY: iip1, jjp1, llmp1 |
USE paramet_m, ONLY: iip1, jjp1, llmp1 |
20 |
use pres2lev_m, only: pres2lev |
use pres2lev_m, only: pres2lev |
21 |
|
|
46 |
real p(iip1, jjp1, llmp1) |
real p(iip1, jjp1, llmp1) |
47 |
real pks(iip1, jjp1) |
real pks(iip1, jjp1) |
48 |
real pls(iip1, jjp1, llm) |
real pls(iip1, jjp1, llm) |
49 |
real prefkap, unskap |
real unskap |
50 |
|
|
51 |
integer i, j, l |
integer i, j, l |
52 |
|
|
58 |
|
|
59 |
! Calcul de pls, pression au milieu des couches, en Pascals |
! Calcul de pls, pression au milieu des couches, en Pascals |
60 |
unskap=1./kappa |
unskap=1./kappa |
|
prefkap = preff ** kappa |
|
61 |
DO l = 1, llm |
DO l = 1, llm |
62 |
DO j=1, jjp1 |
DO j=1, jjp1 |
63 |
DO i =1, iip1 |
DO i =1, iip1 |