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

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

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

revision 150 by guez, Thu Jun 18 13:49:26 2015 UTC revision 207 by guez, Thu Sep 1 10:30:53 2016 UTC
# Line 8  contains Line 8  contains
8    
9      ! Author: F. Hourdin      ! Author: F. Hourdin
10    
     USE histwrite_m, ONLY: histwrite  
     USE dimens_m, ONLY: jjm, llm, nqmx  
     USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1, jjp1  
11      USE comgeom, ONLY: aire      USE comgeom, ONLY: aire
12        USE dimens_m, ONLY: jjm, llm
13        use groupe_m, only: groupe
14        USE histwrite_m, ONLY: histwrite
15        use initfluxsto_m, only: initfluxsto
16        USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1
17      USE tracstoke, ONLY: istdyn, istphy      USE tracstoke, ONLY: istdyn, istphy
18    
19      REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)      REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
20      REAL, intent(in):: masse(ip1jmp1, llm)      REAL, intent(in):: masse(ip1jmp1, llm)
21      real, intent(in):: phi(ip1jmp1, llm)      real, intent(in):: phi(ip1jmp1, llm)
22      real, intent(in):: teta(ip1jmp1, llm)      real, intent(in):: teta(ip1jmp1, llm)
# Line 22  contains Line 24  contains
24      REAL, intent(in):: time_step      REAL, intent(in):: time_step
25      INTEGER, INTENT (IN):: itau      INTEGER, INTENT (IN):: itau
26    
27      ! Variables local to the procedure:      ! Local:
28    
29      REAL, SAVE:: pbaruc(ip1jmp1, llm), pbarvc(ip1jm, llm)      REAL, SAVE:: pbaruc(ip1jmp1, llm), pbarvc(ip1jm, llm)
30      REAL, SAVE:: massem(ip1jmp1, llm)      REAL, SAVE:: massem(ip1jmp1, llm)
# Line 54  contains Line 56  contains
56         CALL initial0(ijmllm, pbarvc)         CALL initial0(ijmllm, pbarvc)
57      END IF      END IF
58    
59      !   accumulation des flux de masse horizontaux      ! Accumulation des flux de masse horizontaux
60      DO l = 1, llm      DO l = 1, llm
61         DO ij = 1, ip1jmp1         DO ij = 1, ip1jmp1
62            pbaruc(ij, l) = pbaruc(ij, l) + pbaru(ij, l)            pbaruc(ij, l) = pbaruc(ij, l) + pbaru(ij, l)
# Line 66  contains Line 68  contains
68         END DO         END DO
69      END DO      END DO
70    
71      !   selection de la masse instantannee des mailles avant le transport.      ! S\'election de la masse instantan\'ee des mailles avant le transport.
72      IF (itau == 0) massem = masse      IF (itau == 0) massem = masse
73    
74      IF (mod(itau + 1, istdyn) == 0) THEN      IF (mod(itau + 1, istdyn) == 0) THEN
75         ! on advecte a ce pas de temps         ! On advecte \`a ce pas de temps
76         !    normalisation         ! normalisation
77         DO l = 1, llm         DO l = 1, llm
78            DO ij = 1, ip1jmp1            DO ij = 1, ip1jmp1
79               pbaruc(ij, l) = pbaruc(ij, l)/float(istdyn)               pbaruc(ij, l) = pbaruc(ij, l)/float(istdyn)
# Line 83  contains Line 85  contains
85            END DO            END DO
86         END DO         END DO
87    
88         !   traitement des flux de masse avant advection.         ! Traitement des flux de masse avant advection.
89         !     1. calcul de w         ! 1. Calcul de w
90         !     2. groupement des mailles pres du pole.         ! 2. Groupement des mailles pr\`es du p\^ole.
   
91         CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg)         CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg)
92    
93         CALL histwrite(fluxid, 'masse', itau, massem)         CALL histwrite(fluxid, 'masse', itau, massem)
Legend:
Removed from v.150  
changed lines
  Added in v.207
  ViewVC Help
Powered by ViewVC 1.1.21