/[lmdze]/trunk/Sources/dyn3d/vitvert.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/vitvert.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/dyn3d/vitvert.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/Sources/dyn3d/vitvert.f revision 252 by guez, Mon Jan 22 15:02:56 2018 UTC
# Line 1  Line 1 
1  SUBROUTINE vitvert(convm , w)  module vitvert_m
2    
3    ! From libf/dyn3d/vitvert.F, version 1.1.1.1 2004/05/19 12:53:05    IMPLICIT NONE
   ! Authors: P. Le Van , F. Hourdin  
4    
5    ! Objet : calcul de la vitesse verticale aux niveaux sigma  contains
6    
7    ! La vitesse verticale est orientee de haut en bas .    pure function vitvert(convm)
   ! au sol, au niveau sigma(1), w(i, j, 1) = 0.  
   ! au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi  
   ! egale a 0. et n'est pas stockee dans le tableau w .  
8    
9    USE dimens_m, ONLY : llm      ! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05
10    USE paramet_m, ONLY : ip1jmp1      ! Authors: P. Le Van, F. Hourdin
   USE disvert_m, ONLY : bp  
11    
12    IMPLICIT NONE      ! Purpose: Compute vertical speed at sigma levels.
13    
14        ! Vertical speed is oriented from bottom to top. At ground-level
15        ! sigma(1): vitvert(i, j, 1) = 0. At top-level sigma(llm + 1), vertical
16        ! speed is 0 too and is not stored in vitvert.
17        
18        USE dimens_m, ONLY : llm
19        USE disvert_m, ONLY : bp
20    
21        real, intent(in):: convm(:, :, :) ! (iim + 1, jjm + 1, llm)
22        REAL vitvert(size(convm, 1), size(convm, 2), size(convm, 3))
23    
24    real, intent(in):: convm(ip1jmp1, llm)      ! Local:
25    REAL, intent(out):: w(ip1jmp1, llm)      INTEGER l
26    
27    ! Local:      !------------------------------------------------------
   INTEGER l  
28    
29    !------------------------------------------------------      forall (l = 2: llm) &
30             vitvert(:, :, l) = convm(:, :, l) - bp(l) * convm(:, :, 1)
31        vitvert(:, :, 1) = 0.
32    
33    forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1)    END function vitvert
   w(:, 1) = 0.  
34    
35  END SUBROUTINE vitvert  end module vitvert_m
Legend:
Removed from v.76  
changed lines
  Added in v.252
  ViewVC Help
Powered by ViewVC 1.1.21