11 |
|
|
12 |
! From dyn3d/sortvarc.F, version 1.1.1.1 2004/05/19 12:53:07 |
! From dyn3d/sortvarc.F, version 1.1.1.1 2004/05/19 12:53:07 |
13 |
! Author: P. Le Van |
! Author: P. Le Van |
14 |
! Objet : sortie des variables de contrôle |
! Objet : sortie des variables de contr\^ole |
15 |
|
|
16 |
USE comconst, ONLY: daysec, g, omeg, rad |
USE comconst, ONLY: daysec, g, omeg, rad |
17 |
USE comgeom, ONLY: aire_2d, cu_2d, rlatu |
USE comgeom, ONLY: aire_2d, cu_2d, rlatu |
18 |
USE dimens_m, ONLY: iim, jjm, llm |
USE dimens_m, ONLY: iim, jjm, llm |
19 |
USE ener, ONLY: ang0, etot0, ptot0, stot0, ztot0 |
USE ener, ONLY: ang0, etot0, ptot0, stot0, ztot0 |
20 |
use filtreg_m, only: filtreg |
use filtreg_scal_m, only: filtreg_scal |
21 |
use massbarxy_m, only: massbarxy |
use massbarxy_m, only: massbarxy |
22 |
USE paramet_m, ONLY: iip1, ip1jm, jjp1 |
USE paramet_m, ONLY: iip1, ip1jm, jjp1 |
23 |
|
|
53 |
rmsdpdt = sum(ge) - sum(ge(1, :)) |
rmsdpdt = sum(ge) - sum(ge(1, :)) |
54 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
55 |
bernf = bern |
bernf = bern |
56 |
CALL filtreg(bernf, direct = .false., intensive = .false.) |
CALL filtreg_scal(bernf, direct = .false., intensive = .false.) |
57 |
|
|
58 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
59 |
radsg = rad/g |
radsg = rad/g |