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

Last change on this file since 5139 was 5134, checked in by abarral, 4 months ago

Replace academic.h, alpale.h, comdissip.h, comdissipn.h, comdissnew.h by modules
Remove unused clesph0.h
  • 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.