11 |
USE dimphy, ONLY: klev |
USE dimphy, ONLY: klev |
12 |
USE suphec_m, ONLY: rcpd, rd, rg, rkappa |
USE suphec_m, ONLY: rcpd, rd, rg, rkappa |
13 |
|
|
14 |
REAL, intent(out), dimension(:, :):: pkf, cq, dq, ch, dh ! (knon, klev) |
REAL, intent(in):: pkf(:, :) ! (knon, klev) |
15 |
|
REAL, intent(out), dimension(:, :):: cq, dq, ch, dh ! (knon, klev) |
16 |
|
|
17 |
REAL, intent(in):: paprs(:, :) ! (knon, klev + 1) |
REAL, intent(in):: paprs(:, :) ! (knon, klev + 1) |
18 |
! pression a inter-couche (Pa) |
! pression a inter-couche (Pa) |
47 |
|
|
48 |
!---------------------------------------------------------------- |
!---------------------------------------------------------------- |
49 |
|
|
|
forall (k = 1:klev) pkf(:, k) = (paprs(:, 1) / pplay(:, k))**RKAPPA |
|
|
! (La pression de r\'ef\'erence est celle au sol.) |
|
|
|
|
50 |
h = RCPD * t * pkf |
h = RCPD * t * pkf |
51 |
|
|
52 |
! Convertir les coefficients en variables convenables au calcul: |
! Convertir les coefficients en variables convenables au calcul: |