source: LMDZ6/trunk/libf/dyn3d_common/vitvert.f90

Last change on this file was 5272, checked in by abarral, 25 hours ago

Turn paramet.h into a module
  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1SUBROUTINE 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  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
9USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
10          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
11IMPLICIT NONE
12
13
14!===============================================================================
15! Arguments:
16  REAL, INTENT(IN)  :: convm(ip1jmp1,llm)
17  REAL, INTENT(OUT) :: w    (ip1jmp1,llm)
18!===============================================================================
19! Notes: Vertical speed is oriented from bottom to top.
20!   * At ground - level sigma(1):     w(i,j,1) = 0.
21!   * At top    - level sigma(llm+1): w(i,j,l) = 0. (not stored in w)
22!===============================================================================
23! Local variables:
24  INTEGER :: l
25!===============================================================================
26  DO l=1,llmm1; w(:,l+1)=convm(:,l+1)-bp(l+1)*convm(:,1); END DO
27  w(:,1)=0.
28
29END SUBROUTINE vitvert
30
Note: See TracBrowser for help on using the repository browser.