/[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 190 by guez, Thu Apr 14 15:15:56 2016 UTC revision 208 by guez, Wed Dec 7 16:44:53 2016 UTC
# Line 8  contains Line 8  contains
8    
9      ! Author: F. Hourdin      ! Author: F. Hourdin
10    
11        USE comgeom, ONLY: aire
12        USE dimens_m, ONLY: jjm, llm
13        use groupe_m, only: groupe
14      USE histwrite_m, ONLY: histwrite      USE histwrite_m, ONLY: histwrite
15      use initfluxsto_m, only: initfluxsto      use initfluxsto_m, only: initfluxsto
     USE dimens_m, ONLY: jjm, llm  
16      USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1      USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1
     USE comgeom, ONLY: aire  
17      USE tracstoke, ONLY: istdyn, istphy      USE tracstoke, ONLY: istdyn, istphy
18    
19      REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)      REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
# Line 75  contains Line 76  contains
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)/real(istdyn)
80               tetac(ij, l) = tetac(ij, l)/float(istdyn)               tetac(ij, l) = tetac(ij, l)/real(istdyn)
81               phic(ij, l) = phic(ij, l)/float(istdyn)               phic(ij, l) = phic(ij, l)/real(istdyn)
82            END DO            END DO
83            DO ij = 1, ip1jm            DO ij = 1, ip1jm
84               pbarvc(ij, l) = pbarvc(ij, l)/float(istdyn)               pbarvc(ij, l) = pbarvc(ij, l)/real(istdyn)
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 pr\`es du p\^ole.         ! 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.190  
changed lines
  Added in v.208
  ViewVC Help
Powered by ViewVC 1.1.21