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

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

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

trunk/dyn3d/convmas.f revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/dyn3d/convmas.f90 revision 81 by guez, Wed Mar 5 14:38:41 2014 UTC
# Line 1  Line 1 
1  !  
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19
3  !  ! 12:53:07 lmdzadmin Exp $
4        SUBROUTINE convmas (pbaru, pbarv, convm )  
5  c  SUBROUTINE convmas(pbaru, pbarv, convm)
6        use dimens_m  
7        use paramet_m    USE dimens_m
8        use disvert_m    USE paramet_m
9        use conf_gcm_m    USE disvert_m
10        use filtreg_m, only: filtreg    USE conf_gcm_m
11        IMPLICIT NONE    USE filtreg_m, ONLY: filtreg
12      IMPLICIT NONE
13  c=======================================================================  
14  c    ! =======================================================================
15  c   Auteurs:  P. Le Van , F. Hourdin  .  
16  c   -------    ! Auteurs:  P. Le Van , F. Hourdin  .
17  c    ! -------
18  c   Objet:  
19  c   ------    ! Objet:
20  c    ! ------
21  c   ********************************************************************  
22  c   .... calcul de la convergence du flux de masse aux niveaux p ...    ! ********************************************************************
23  c   ********************************************************************    ! .... calcul de la convergence du flux de masse aux niveaux p ...
24  c    ! ********************************************************************
25  c  
26  c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....  
27  c      .....  convm      est  un argument de sortie pour le s-pg  ....    ! pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
28  c    ! .....  convm      est  un argument de sortie pour le s-pg  ....
29  c    le calcul se fait de haut en bas,  
30  c    la convergence de masse au niveau p(llm+1) est egale a 0. et    ! le calcul se fait de haut en bas,
31  c    n'est pas stockee dans le tableau convm .    ! la convergence de masse au niveau p(llm+1) est egale a 0. et
32  c    ! n'est pas stockee dans le tableau convm .
33  c  
34  c=======================================================================  
35  c    ! =======================================================================
36  c   Declarations:  
37  c   -------------    ! Declarations:
38      ! -------------
39    
40        REAL, intent(in):: pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm )  
41        real, intent(out):: convm(  ip1jmp1,llm )    REAL, INTENT (IN) :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
42        INTEGER   l,ij    REAL, INTENT (OUT) :: convm(ip1jmp1, llm)
43      INTEGER l, ij
44    
45  c-----------------------------------------------------------------------  
46  c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......    ! -----------------------------------------------------------------------
47      ! ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
48        CALL  convflu( pbaru, pbarv, llm, convm )  
49      CALL convflu(pbaru, pbarv, llm, convm)
50  c-----------------------------------------------------------------------  
51  c   filtrage:    ! -----------------------------------------------------------------------
52  c   ---------    ! filtrage:
53      ! ---------
54         CALL filtreg( convm, jjp1, llm, 2, 2, .true.)  
55      CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)
56  c    integration de la convergence de masse de haut  en bas ......  
57      ! integration de la convergence de masse de haut  en bas ......
58        DO      l      = llmm1, 1, -1  
59          DO    ij     = 1, ip1jmp1    DO l = llmm1, 1, -1
60           convm(ij,l) = convm(ij,l) + convm(ij,l+1)      DO ij = 1, ip1jmp1
61          ENDDO        convm(ij, l) = convm(ij, l) + convm(ij, l+1)
62        ENDDO      END DO
63  c    END DO
64        RETURN  
65        END    RETURN
66    END SUBROUTINE convmas
Legend:
Removed from v.76  
changed lines
  Added in v.81
  ViewVC Help
Powered by ViewVC 1.1.21