136 |
dtvr, itau) |
dtvr, itau) |
137 |
|
|
138 |
! integrations dynamique et traceurs: |
! integrations dynamique et traceurs: |
139 |
CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, & |
CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, dp, & |
140 |
dp, vcov, ucov, teta, q, ps, masse, finvmaold, leapf, dt) |
vcov, ucov, teta, q(:, :, :2), ps, masse, finvmaold, dt, leapf) |
141 |
|
|
142 |
if (.not. leapf) then |
if (.not. leapf) then |
143 |
! Matsuno backward |
! Matsuno backward |
150 |
phi, .false., du, dv, dteta, dp, w, pbaru, pbarv, time_0) |
phi, .false., du, dv, dteta, dp, w, pbaru, pbarv, time_0) |
151 |
|
|
152 |
! integrations dynamique et traceurs: |
! integrations dynamique et traceurs: |
153 |
CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, & |
CALL integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, & |
154 |
dteta, dp, vcov, ucov, teta, q, ps, masse, finvmaold, .false., & |
dp, vcov, ucov, teta, q(:, :, :2), ps, masse, finvmaold, dtvr, & |
155 |
dtvr) |
leapf=.false.) |
156 |
end if |
end if |
157 |
|
|
158 |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |