18 |
USE exner_hyb_m, ONLY: exner_hyb |
USE exner_hyb_m, ONLY: exner_hyb |
19 |
use massbar_m, only: massbar |
use massbar_m, only: massbar |
20 |
use massdair_m, only: massdair |
use massdair_m, only: massdair |
21 |
|
use pres2lev_m, only: pres2lev |
22 |
|
|
23 |
integer nlevnc |
integer nlevnc |
24 |
real, intent(in):: psi(iip1, jjp1) |
real, intent(in):: psi(iip1, jjp1) |
43 |
|
|
44 |
real p(iip1, jjp1, llmp1) |
real p(iip1, jjp1, llmp1) |
45 |
real pks(iip1, jjp1) |
real pks(iip1, jjp1) |
|
real pkf(iip1, jjp1, llm) |
|
46 |
real pls(iip1, jjp1, llm) |
real pls(iip1, jjp1, llm) |
47 |
real prefkap, unskap |
real prefkap, unskap |
48 |
|
|
53 |
! calcul de la pression au milieu des couches |
! calcul de la pression au milieu des couches |
54 |
forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi |
forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi |
55 |
call massdair(p, masse) |
call massdair(p, masse) |
56 |
CALL exner_hyb(psi, p, pks, pk, pkf) |
CALL exner_hyb(psi, p, pks, pk) |
57 |
|
|
58 |
! Calcul de pls, pression au milieu des couches, en Pascals |
! Calcul de pls, pression au milieu des couches, en Pascals |
59 |
unskap=1./kappa |
unskap=1./kappa |