9 |
USE comconst, ONLY : daysec, g, omeg, rad |
USE comconst, ONLY : daysec, g, omeg, rad |
10 |
USE comgeom, ONLY : aire, cu, rlatu |
USE comgeom, ONLY : aire, cu, rlatu |
11 |
USE ener, ONLY : ang0, etot0, ptot0, rmsdpdt, rmsv, stot0, ztot0 |
USE ener, ONLY : ang0, etot0, ptot0, rmsdpdt, rmsv, stot0, ztot0 |
12 |
|
use filtreg_m, only: filtreg |
13 |
|
|
14 |
IMPLICIT NONE |
IMPLICIT NONE |
15 |
|
|
46 |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
47 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1)) |
48 |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
49 |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1) |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.) |
50 |
|
|
51 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
52 |
|
|