source: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vitvert_loc.F90 @ 5214

Last change on this file since 5214 was 5159, checked in by abarral, 3 months ago

Put dimensions.h and paramet.h into modules
  • 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 
10USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
11  USE lmdz_paramet
12  IMPLICIT NONE
13
14
15!===============================================================================
16! Arguments:
17  REAL, INTENT(IN)  :: convm(ijb_u:ije_u,llm)
18  REAL, INTENT(OUT) :: w    (ijb_u:ije_u,llm)
19!===============================================================================
20! Notes: Vertical speed is oriented from bottom to top.
21!   * At ground - level sigma(1):     w(i,j,1) = 0.
22!   * At top    - level sigma(llm+1): w(i,j,l) = 0. (not stored in w)
23!===============================================================================
24! Local variables:
25  INTEGER :: l, ijb, ije
26!===============================================================================
27  ijb=ij_begin
28  ije=ij_end+iip1
29  IF(pole_sud) ije=ij_end
30!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
31  DO l=1,llmm1
32    w(ijb:ije,l+1)=convm(ijb:ije,l+1)-bp(l+1)*convm(ijb:ije,1)
33  END DO
34!$OMP END DO
35!$OMP MASTER
36  w(ijb:ije,1)=0.
37!$OMP END MASTER
38!$OMP BARRIER
39
40END SUBROUTINE vitvert_loc
41
Note: See TracBrowser for help on using the repository browser.