1 |
SUBROUTINE vitvert(convm , w) |
module vitvert_m |
2 |
|
|
3 |
! From libf/dyn3d/vitvert.F, version 1.1.1.1 2004/05/19 12:53:05 |
IMPLICIT NONE |
|
! Authors: P. Le Van , F. Hourdin |
|
4 |
|
|
5 |
! Objet : calcul de la vitesse verticale aux niveaux sigma |
contains |
6 |
|
|
7 |
! La vitesse verticale est orientee de haut en bas . |
pure function vitvert(convm) |
|
! au sol, au niveau sigma(1), w(i, j, 1) = 0. |
|
|
! au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi |
|
|
! egale a 0. et n'est pas stockee dans le tableau w . |
|
8 |
|
|
9 |
USE dimens_m, ONLY : llm |
! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05 |
10 |
USE paramet_m, ONLY : ip1jmp1 |
! Authors: P. Le Van, F. Hourdin |
|
USE disvert_m, ONLY : bp |
|
11 |
|
|
12 |
IMPLICIT NONE |
! Purpose: Compute vertical speed at sigma levels. |
13 |
|
|
14 |
|
! Vertical speed is oriented from bottom to top. At ground-level |
15 |
|
! sigma(1): vitvert(i, j, 1) = 0. At top-level sigma(llm + 1), vertical |
16 |
|
! speed is 0 too and is not stored in vitvert. |
17 |
|
|
18 |
|
USE dimensions, ONLY : llm |
19 |
|
USE disvert_m, ONLY : bp |
20 |
|
|
21 |
|
real, intent(in):: convm(:, :, :) ! (iim + 1, jjm + 1, llm) |
22 |
|
REAL vitvert(size(convm, 1), size(convm, 2), size(convm, 3)) |
23 |
|
|
24 |
real, intent(in):: convm(ip1jmp1, llm) |
! Local: |
25 |
REAL, intent(out):: w(ip1jmp1, llm) |
INTEGER l |
26 |
|
|
27 |
! Local: |
!------------------------------------------------------ |
|
INTEGER l |
|
28 |
|
|
29 |
!------------------------------------------------------ |
forall (l = 2: llm) & |
30 |
|
vitvert(:, :, l) = convm(:, :, l) - bp(l) * convm(:, :, 1) |
31 |
|
vitvert(:, :, 1) = 0. |
32 |
|
|
33 |
forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1) |
END function vitvert |
|
w(:, 1) = 0. |
|
34 |
|
|
35 |
END SUBROUTINE vitvert |
end module vitvert_m |