42 |
|
|
43 |
real p(iip1, jjp1, llmp1) |
real p(iip1, jjp1, llmp1) |
44 |
real pks(iip1, jjp1) |
real pks(iip1, jjp1) |
|
real pkf(iip1, jjp1, llm) |
|
45 |
real pls(iip1, jjp1, llm) |
real pls(iip1, jjp1, llm) |
46 |
real prefkap, unskap |
real prefkap, unskap |
47 |
|
|
52 |
! calcul de la pression au milieu des couches |
! calcul de la pression au milieu des couches |
53 |
forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi |
forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi |
54 |
call massdair(p, masse) |
call massdair(p, masse) |
55 |
CALL exner_hyb(psi, p, pks, pk, pkf) |
CALL exner_hyb(psi, p, pks, pk) |
56 |
|
|
57 |
! Calcul de pls, pression au milieu des couches, en Pascals |
! Calcul de pls, pression au milieu des couches, en Pascals |
58 |
unskap=1./kappa |
unskap=1./kappa |