8 |
|
|
9 |
! 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 |
10 |
! Authors: P. Le Van, L. Fairhead, F. Hourdin |
! Authors: P. Le Van, L. Fairhead, F. Hourdin |
11 |
! schema matsuno + leapfrog |
! Matsuno-leapfrog scheme. |
12 |
|
|
13 |
USE calfis_m, ONLY: calfis |
USE calfis_m, ONLY: calfis |
14 |
USE com_io_dyn, ONLY: histaveid |
USE com_io_dyn, ONLY: histaveid |
32 |
USE temps, ONLY: itau_dyn |
USE temps, ONLY: itau_dyn |
33 |
|
|
34 |
! Variables dynamiques: |
! Variables dynamiques: |
|
REAL, intent(inout):: vcov((iim + 1) * jjm, llm) ! vent covariant |
|
35 |
REAL, intent(inout):: ucov(ip1jmp1, llm) ! vent covariant |
REAL, intent(inout):: ucov(ip1jmp1, llm) ! vent covariant |
36 |
|
REAL, intent(inout):: vcov((iim + 1) * jjm, llm) ! vent covariant |
37 |
REAL, intent(inout):: teta(iim + 1, jjm + 1, llm) ! potential temperature |
REAL, intent(inout):: teta(iim + 1, jjm + 1, llm) ! potential temperature |
38 |
REAL ps(iim + 1, jjm + 1) ! pression au sol, en Pa |
REAL ps(iim + 1, jjm + 1) ! pression au sol, en Pa |
|
|
|
39 |
REAL masse(ip1jmp1, llm) ! masse d'air |
REAL masse(ip1jmp1, llm) ! masse d'air |
40 |
REAL phis(ip1jmp1) ! geopotentiel au sol |
REAL phis(ip1jmp1) ! geopotentiel au sol |
41 |
REAL q(ip1jmp1, llm, nqmx) ! mass fractions of advected fields |
REAL q(ip1jmp1, llm, nqmx) ! mass fractions of advected fields |
166 |
time = REAL(mod(itau, day_step)) / day_step + time_0 |
time = REAL(mod(itau, day_step)) / day_step + time_0 |
167 |
IF (time > 1.) time = time - 1. |
IF (time > 1.) time = time - 1. |
168 |
|
|
169 |
CALL calfis(nqmx, itau + 1 == itaufin, rdayvrai, time, ucov, vcov, & |
CALL calfis(itau + 1 == itaufin, rdayvrai, time, ucov, vcov, & |
170 |
teta, q, masse, ps, pk, phis, phi, du, dv, dteta, dq, w, dufi, & |
teta, q, masse, ps, pk, phis, phi, du, dv, dteta, dq, w, dufi, & |
171 |
dvfi, dtetafi, dqfi, dpfi) |
dvfi, dtetafi, dqfi, dpfi) |
172 |
|
|