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 |
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 |