26 |
use dynredem1_m, only: dynredem1 |
use dynredem1_m, only: dynredem1 |
27 |
USE exner_hyb_m, ONLY: exner_hyb |
USE exner_hyb_m, ONLY: exner_hyb |
28 |
use filtreg_m, only: filtreg |
use filtreg_m, only: filtreg |
29 |
|
use fluxstokenc_m, only: fluxstokenc |
30 |
use geopot_m, only: geopot |
use geopot_m, only: geopot |
31 |
USE guide_m, ONLY: guide |
USE guide_m, ONLY: guide |
32 |
use inidissip_m, only: idissip |
use inidissip_m, only: idissip |
35 |
USE pressure_var, ONLY: p3d |
USE pressure_var, ONLY: p3d |
36 |
USE temps, ONLY: itau_dyn |
USE temps, ONLY: itau_dyn |
37 |
use writedynav_m, only: writedynav |
use writedynav_m, only: writedynav |
38 |
|
use writehist_m, only: writehist |
39 |
|
|
40 |
! Variables dynamiques: |
! Variables dynamiques: |
41 |
REAL, intent(inout):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) vent covariant |
REAL, intent(inout):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) vent covariant |
46 |
|
|
47 |
REAL, intent(inout):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol, en Pa |
REAL, intent(inout):: ps(:, :) ! (iim + 1, jjm + 1) pression au sol, en Pa |
48 |
REAL masse(:, :, :) ! (iim + 1, jjm + 1, llm) masse d'air |
REAL masse(:, :, :) ! (iim + 1, jjm + 1, llm) masse d'air |
49 |
REAL phis(:, :) ! (iim + 1, jjm + 1) geopotentiel au sol |
REAL, intent(in):: phis(:, :) ! (iim + 1, jjm + 1) surface geopotential |
50 |
|
|
51 |
REAL, intent(inout):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx) |
REAL, intent(inout):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx) |
52 |
! mass fractions of advected fields |
! mass fractions of advected fields |
170 |
end if |
end if |
171 |
|
|
172 |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
173 |
! calcul des tendances physiques: |
! Calcul des tendances physiques: |
174 |
|
|
175 |
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
forall (l = 1: llm + 1) p3d(:, :, l) = ap(l) + bp(l) * ps |
176 |
CALL exner_hyb(ps, p3d, pks, pk, pkf) |
CALL exner_hyb(ps, p3d, pks, pk, pkf) |
184 |
phis, phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, & |
phis, phi, dudyn, dv, dq, w, dufi, dvfi, dtetafi, dqfi, dpfi, & |
185 |
lafin = itau + 1 == itaufin) |
lafin = itau + 1 == itaufin) |
186 |
|
|
187 |
! ajout des tendances physiques: |
! Ajout des tendances physiques: |
188 |
CALL addfi(nqmx, dtphys, ucov, vcov, teta, q, ps, dufi, dvfi, & |
CALL addfi(nqmx, dtphys, ucov, vcov, teta, q, ps, dufi, dvfi, & |
189 |
dtetafi, dqfi, dpfi) |
dtetafi, dqfi, dpfi) |
190 |
ENDIF |
ENDIF |
234 |
|
|
235 |
IF (MOD(itau + 1, iecri * day_step) == 0) THEN |
IF (MOD(itau + 1, iecri * day_step) == 0) THEN |
236 |
CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi) |
CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi) |
237 |
CALL writehist(itau, vcov, ucov, teta, phi, q, masse, ps, phis) |
CALL writehist(itau, vcov, ucov, teta, phi, q, masse, ps) |
238 |
END IF |
END IF |
239 |
end do time_integration |
end do time_integration |
240 |
|
|