35 |
use inidissip_m, only: idissip |
use inidissip_m, only: idissip |
36 |
use integrd_m, only: integrd |
use integrd_m, only: integrd |
37 |
use nr_util, only: assert |
use nr_util, only: assert |
|
USE pressure_var, ONLY: p3d |
|
38 |
USE temps, ONLY: itau_dyn |
USE temps, ONLY: itau_dyn |
39 |
use writedynav_m, only: writedynav |
use writedynav_m, only: writedynav |
40 |
use writehist_m, only: writehist |
use writehist_m, only: writehist |
97 |
logical leapf |
logical leapf |
98 |
real dt ! time step, in s |
real dt ! time step, in s |
99 |
|
|
100 |
|
REAL p3d(iim + 1, jjm + 1, llm+1) ! pressure at layer interfaces, in Pa |
101 |
|
! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", |
102 |
|
! for interface "l") |
103 |
|
|
104 |
!--------------------------------------------------- |
!--------------------------------------------------- |
105 |
|
|
106 |
print *, "Call sequence information: leapfrog" |
print *, "Call sequence information: leapfrog" |
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 |
CALL calfis(ucov, vcov, teta, q, pk, phis, phi, w, dufi, dvfi, & |
CALL calfis(ucov, vcov, teta, q, p3d, pk, phis, phi, w, dufi, dvfi, & |
174 |
dtetafi, dqfi, dayvrai = itau / day_step + day_ini, & |
dtetafi, dqfi, dayvrai = itau / day_step + day_ini, & |
175 |
time = REAL(mod(itau, day_step)) / day_step, & |
time = REAL(mod(itau, day_step)) / day_step, & |
176 |
lafin = itau + 1 == itaufin) |
lafin = itau + 1 == itaufin) |