/[lmdze]/trunk/dyn3d/bilan_dyn.f
ViewVC logotype

Diff of /trunk/dyn3d/bilan_dyn.f

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

revision 54 by guez, Tue Dec 6 15:07:04 2011 UTC revision 55 by guez, Mon Dec 12 13:25:01 2011 UTC
# Line 10  contains Line 10  contains
10      ! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16      ! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16
11      ! 10:12:17 fairhead      ! 10:12:17 fairhead
12    
13      ! Sous-programme consacré à des diagnostics dynamiques de base      ! Sous-programme consacré à des diagnostics dynamiques de base.
14      ! De façon générale, les moyennes des scalaires Q sont pondérées par      ! De façon générale, les moyennes des scalaires Q sont pondérées
15      ! la masse. Les flux de masse sont, eux, simplement moyennés.      ! par la masse. Les flux de masse sont, eux, simplement moyennés.
16    
17      USE histcom, ONLY: histbeg_totreg, histdef, histend, histvert      USE histcom, ONLY: histbeg_totreg, histdef, histend, histvert
18      USE calendar, ONLY: ymds2ju      USE calendar, ONLY: ymds2ju
# Line 373  contains Line 373  contains
373         psi = psi*1.e-9         psi = psi*1.e-9
374         call histwrite(fileid, 'psi', itau, psi(:, :llm))         call histwrite(fileid, 'psi', itau, psi(:, :llm))
375    
376         ! Moyenne verticale         ! Intégrale verticale
377    
378         forall (iQ = 1: nQ, itr = 2: ntr) zavQ(:, itr, iQ) &         forall (iQ = 1: nQ, itr = 2: ntr) zavQ(:, itr, iQ) &
379              = sum(zvQ(:, :, itr, iQ) * zmasse, dim=2) / sum(zmasse, dim=2)              = sum(zvQ(:, :, itr, iQ) * zmasse, dim=2) / cv_2d(1, :)
380    
381         do iQ = 1, nQ         do iQ = 1, nQ
382            do itr = 2, ntr            do itr = 2, ntr
Legend:
Removed from v.54  
changed lines
  Added in v.55
  ViewVC Help
Powered by ViewVC 1.1.21