source: LMDZ6/branches/SETHET_DECOUPLE/libf/dyn3dmem/vitvert_loc.F90 @ 5440

Last change on this file since 5440 was 2600, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn comvert.h into module comvert_mod.F90
EM
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.3 KB
Line 
1SUBROUTINE vitvert_loc(convm, w)
2!
3!-------------------------------------------------------------------------------
4! Authors: P. Le Van , Fr. Hourdin.
5!-------------------------------------------------------------------------------
6! Purpose: Compute vertical speed at sigma levels.
7  USE parallel_lmdz
8  USE comvert_mod, ONLY: bp
9 
10  IMPLICIT NONE
11  include "dimensions.h"
12  include "paramet.h"
13!===============================================================================
14! Arguments:
15  REAL, INTENT(IN)  :: convm(ijb_u:ije_u,llm)
16  REAL, INTENT(OUT) :: w    (ijb_u:ije_u,llm)
17!===============================================================================
18! Notes: Vertical speed is oriented from bottom to top.
19!   * At ground - level sigma(1):     w(i,j,1) = 0.
20!   * At top    - level sigma(llm+1): w(i,j,l) = 0. (not stored in w)
21!===============================================================================
22! Local variables:
23  INTEGER :: l, ijb, ije
24!===============================================================================
25  ijb=ij_begin
26  ije=ij_end+iip1
27  IF(pole_sud) ije=ij_end
28!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
29  DO l=1,llmm1
30    w(ijb:ije,l+1)=convm(ijb:ije,l+1)-bp(l+1)*convm(ijb:ije,1)
31  END DO
32!$OMP END DO
33!$OMP MASTER
34  w(ijb:ije,1)=0.
35!$OMP END MASTER
36!$OMP BARRIER
37
38END SUBROUTINE vitvert_loc
39
Note: See TracBrowser for help on using the repository browser.