1 |
|
2017 |
SUBROUTINE vitvert (convm, w) |
2 |
|
|
! |
3 |
|
|
!------------------------------------------------------------------------------- |
4 |
|
|
! Authors: P. Le Van , Fr. Hourdin. |
5 |
|
|
!------------------------------------------------------------------------------- |
6 |
|
|
! Purpose: Compute vertical speed at sigma levels. |
7 |
|
|
USE comvert_mod, ONLY: bp |
8 |
|
|
IMPLICIT NONE |
9 |
|
|
include "dimensions.h" |
10 |
|
|
include "paramet.h" |
11 |
|
|
!=============================================================================== |
12 |
|
|
! Arguments: |
13 |
|
|
REAL, INTENT(IN) :: convm(ip1jmp1,llm) |
14 |
|
|
REAL, INTENT(OUT) :: w (ip1jmp1,llm) |
15 |
|
|
!=============================================================================== |
16 |
|
|
! Notes: Vertical speed is oriented from bottom to top. |
17 |
|
|
! * At ground - level sigma(1): w(i,j,1) = 0. |
18 |
|
|
! * At top - level sigma(llm+1): w(i,j,l) = 0. (not stored in w) |
19 |
|
|
!=============================================================================== |
20 |
|
|
! Local variables: |
21 |
|
|
INTEGER :: l |
22 |
|
|
!=============================================================================== |
23 |
✓✓✓✓
|
83546157 |
DO l=1,llmm1; w(:,l+1)=convm(:,l+1)-bp(l+1)*convm(:,1); END DO |
24 |
✓✓ |
2198530 |
w(:,1)=0. |
25 |
|
|
|
26 |
|
2017 |
END SUBROUTINE vitvert |
27 |
|
|
|