44 |
|
|
45 |
!----------------------------------------------------------------------- |
!----------------------------------------------------------------------- |
46 |
|
|
47 |
|
print *, "Call sequence information: sortvarc" |
48 |
|
|
49 |
time = real(itau) / day_step + time_0 |
time = real(itau) / day_step + time_0 |
50 |
dtvrs1j = dtvr/daysec |
dtvrs1j = dtvr/daysec |
51 |
rjour = real(int(itau*dtvrs1j)) |
rjour = real(int(itau*dtvrs1j)) |
60 |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1) |
61 |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1) |
62 |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
CALL scopy(ijp1llm, bern, 1, bernf, 1) |
63 |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1) |
CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.) |
64 |
|
|
65 |
! Calcul du moment angulaire |
! Calcul du moment angulaire |
66 |
radsg = rad/g |
radsg = rad/g |