4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE bernoui(phi, ecin, bern) |
function bernoui(phi, ecin) |
8 |
|
|
9 |
! From LMDZ4/libf/dyn3d/bernoui.F, version 1.1.1.1 2004/05/19 12:53:06 |
! From LMDZ4/libf/dyn3d/bernoui.F, version 1.1.1.1 2004/05/19 12:53:06 |
10 |
|
|
11 |
! Author: P. Le Van |
! Author: P. Le Van |
12 |
|
|
13 |
! Objet: calcul de la fonction de Bernouilli aux niveaux s. |
! Objet : calcul de la fonction de Bernouilli aux niveaux s. |
14 |
! fonction de Bernouilli = bern = filtre de(geopotentiel + energ.cinet.) |
|
15 |
|
USE dimens_m, ONLY: iim, jjm, llm |
|
USE dimens_m, ONLY: jjm, llm |
|
16 |
USE filtreg_m, ONLY: filtreg |
USE filtreg_m, ONLY: filtreg |
17 |
|
|
18 |
REAL, INTENT(IN):: phi(:, :, :), ecin (:, :, :)! (iim + 1, jjm + 1, llm) |
REAL, INTENT(IN):: phi(:, :, :), ecin (:, :, :) ! (iim + 1, jjm + 1, llm) |
19 |
REAL, intent(out):: bern(:, :, :) ! (iim + 1, jjm + 1, llm) |
|
20 |
|
REAL bernoui(iim + 1, jjm + 1, llm) |
21 |
|
! fonction de Bernouilli = filtre de (géopotentiel + énergie cinétique) |
22 |
|
|
23 |
!----------------------------------------------------------------------- |
!----------------------------------------------------------------------- |
24 |
|
|
25 |
bern = phi + ecin |
bernoui = phi + ecin |
26 |
CALL filtreg(bern, direct = .true., intensive = .true.) |
CALL filtreg(bernoui, direct = .true., intensive = .true.) |
27 |
|
|
28 |
END SUBROUTINE bernoui |
END function bernoui |
29 |
|
|
30 |
end module bernoui_m |
end module bernoui_m |