6 |
|
|
7 |
contains |
contains |
8 |
|
|
9 |
SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, clesphy0, & |
SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, time_0) |
|
time_0) |
|
10 |
|
|
11 |
! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34 |
! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34 |
12 |
|
|
46 |
use pression_m, only: pression |
use pression_m, only: pression |
47 |
use pressure_var, only: p3d |
use pressure_var, only: p3d |
48 |
|
|
49 |
integer nq |
integer, intent(in):: nq |
|
REAL, intent(in):: clesphy0(:) |
|
50 |
|
|
51 |
! Variables dynamiques: |
! Variables dynamiques: |
52 |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants |
94 |
|
|
95 |
INTEGER itau, itaufinp1 |
INTEGER itau, itaufinp1 |
96 |
INTEGER iday ! jour julien |
INTEGER iday ! jour julien |
97 |
REAL time ! Heure de la journee en fraction d'1 jour |
REAL time ! time of day, as a fraction of day length |
98 |
|
|
99 |
REAL SSUM |
REAL SSUM |
100 |
real finvmaold(ip1jmp1, llm) |
real finvmaold(ip1jmp1, llm) |
218 |
|
|
219 |
CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, & |
CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, & |
220 |
masse, ps, pk, phis, phi, du, dv, dteta, dq, w, & |
masse, ps, pk, phis, phi, du, dv, dteta, dq, w, & |
221 |
clesphy0, dufi, dvfi, dtetafi, dqfi, dpfi) |
dufi, dvfi, dtetafi, dqfi, dpfi) |
222 |
|
|
223 |
! ajout des tendances physiques: |
! ajout des tendances physiques: |
224 |
CALL addfi(nqmx, dtphys, & |
CALL addfi(nqmx, dtphys, & |
335 |
dt = 2. * dtvr |
dt = 2. * dtvr |
336 |
END IF |
END IF |
337 |
ELSE |
ELSE |
338 |
! ...... pas leapfrog ..... |
! pas leapfrog |
339 |
leapf = .TRUE. |
leapf = .TRUE. |
340 |
dt = 2. * dtvr |
dt = 2. * dtvr |
341 |
END IF |
END IF |