/[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/libf/dyn3d/vitvert.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/Sources/dyn3d/vitvert.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC
# Line 1  Line 1 
1  !  module vitvert_m
 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/vitvert.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  
 !  
       SUBROUTINE vitvert ( convm , w )  
 c  
   
 c=======================================================================  
 c  
 c   Auteurs:  P. Le Van , F. Hourdin .  
 c   -------  
 c  
 c   Objet:  
 c   ------  
 c  
 c    *******************************************************************  
 c  .... calcul de la vitesse verticale aux niveaux sigma  ....  
 c    *******************************************************************  
 c     convm   est un argument  d'entree pour le s-pg  ......  
 c       w     est un argument de sortie pour le s-pg  ......  
 c  
 c    la vitesse verticale est orientee de  haut en bas .  
 c    au sol, au niveau sigma(1),   w(i,j,1) = 0.  
 c    au sommet, au niveau sigma(llm+1) , la vit.verticale est aussi  
 c    egale a 0. et n'est pas stockee dans le tableau w  .  
 c  
 c  
 c=======================================================================  
   
       use dimens_m  
       use paramet_m  
       use comvert  
       IMPLICIT NONE  
   
       REAL w(ip1jmp1,llm),convm(ip1jmp1,llm)  
       INTEGER   l, ij  
2    
3      IMPLICIT NONE
4    
5    contains
6    
7        DO 2  l = 1,llmm1    SUBROUTINE vitvert(convm, w)
   
       DO 1 ij = 1,ip1jmp1  
       w( ij, l+1 ) = convm( ij, l+1 ) - bp(l+1) * convm( ij, 1 )  
    1  CONTINUE  
   
    2  CONTINUE  
   
       DO 5 ij  = 1,ip1jmp1  
       w(ij,1)  = 0.  
 5     CONTINUE  
8    
9        RETURN      ! From libf/dyn3d/vitvert.F, version 1.1.1.1, 2004/05/19 12:53:05
10        END      ! Authors: P. Le Van, F. Hourdin
11    
12        ! Objet : calcul de la vitesse verticale aux niveaux sigma
13    
14        ! La vitesse verticale est orientée de haut en bas. Au sol, au
15        ! niveau sigma(1), w(i, j, 1) = 0. Au sommet, au niveau
16        ! sigma(llm+1), la vitesse verticale est aussi égale à 0 et n'est
17        ! pas stockée dans le tableau w.
18    
19        USE dimens_m, ONLY : llm
20        USE disvert_m, ONLY : bp
21        USE paramet_m, ONLY : ip1jmp1
22    
23        real, intent(in):: convm(ip1jmp1, llm)
24        REAL, intent(out):: w(ip1jmp1, llm)
25    
26        ! Local:
27        INTEGER l
28    
29        !------------------------------------------------------
30    
31        forall (l = 2: llm) w(:, l) = convm(:, l) - bp(l) * convm(:, 1)
32        w(:, 1) = 0.
33    
34      END SUBROUTINE vitvert
35    
36    end module vitvert_m
Legend:
Removed from v.3  
changed lines
  Added in v.134
  ViewVC Help
Powered by ViewVC 1.1.21