22 |
use filtreg_m, only: filtreg |
use filtreg_m, only: filtreg |
23 |
|
|
24 |
! Arguments: |
! Arguments: |
25 |
INTEGER, INTENT (IN) :: itau |
INTEGER, INTENT(IN) :: itau |
26 |
REAL :: ucov(ip1jmp1, llm), masse(ip1jmp1, llm) |
REAL :: ucov(ip1jmp1, llm), masse(ip1jmp1, llm) |
27 |
real, intent(in):: teta(ip1jmp1, llm) |
real, intent(in):: teta(ip1jmp1, llm) |
28 |
REAL :: ps(ip1jmp1), phis(ip1jmp1) |
REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1) |
29 |
REAL :: vorpot(ip1jm, llm) |
REAL :: vorpot(ip1jm, llm) |
30 |
REAL :: phi(ip1jmp1, llm), bern(ip1jmp1, llm) |
REAL, intent(in):: phi(ip1jmp1, llm) |
31 |
|
real bern(ip1jmp1, llm) |
32 |
REAL :: dp(ip1jmp1) |
REAL :: dp(ip1jmp1) |
33 |
REAL, INTENT (IN):: time_0 |
REAL, INTENT (IN):: time_0 |
34 |
REAL, INTENT (IN):: pk(ip1jmp1, llm) |
REAL, INTENT (IN):: pk(ip1jmp1, llm) |
45 |
|
|
46 |
!----------------------------------------------------------------------- |
!----------------------------------------------------------------------- |
47 |
|
|
48 |
|
print *, "Call sequence information: sortvarc" |
49 |
|
|
50 |
time = real(itau) / day_step + time_0 |
time = real(itau) / day_step + time_0 |
51 |
dtvrs1j = dtvr/daysec |
dtvrs1j = dtvr/daysec |
52 |
rjour = real(int(itau*dtvrs1j)) |
rjour = real(int(itau*dtvrs1j)) |
61 |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
62 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
63 |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
64 |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1) |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.) |
65 |
|
|
66 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
67 |
radsg = rad/g |
radsg = rad/g |