98 |
logical leapf |
logical leapf |
99 |
real dt ! time step, in s |
real dt ! time step, in s |
100 |
|
|
101 |
REAL p3d(iim + 1, jjm + 1, llm+1) ! pressure at layer interfaces, in Pa |
REAL p3d(iim + 1, jjm + 1, llm + 1) ! pressure at layer interfaces, in Pa |
102 |
! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", |
! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", |
103 |
! for interface "l") |
! for interface "l") |
104 |
|
|
170 |
CALL filtreg_scal(pkf, direct = .true., intensive = .true.) |
CALL filtreg_scal(pkf, direct = .true., intensive = .true.) |
171 |
end if |
end if |
172 |
|
|
173 |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys /= 0) THEN |
IF (MOD(itau + 1, iphysiq) == 0 .AND. iflag_phys) THEN |
174 |
CALL calfis(ucov, vcov, teta, q, p3d, pk, phis, phi, w, dufi, dvfi, & |
CALL calfis(ucov, vcov, teta, q, p3d, pk, phis, phi, w, dufi, dvfi, & |
175 |
dtetafi, dqfi, dayvrai = itau / day_step + day_ini, & |
dtetafi, dqfi, dayvrai = itau / day_step + day_ini, & |
176 |
time = REAL(mod(itau, day_step)) / day_step, & |
time = REAL(mod(itau, day_step)) / day_step, & |
202 |
forall (l = 1: llm) |
forall (l = 1: llm) |
203 |
teta(:, 1, l) = SUM(aire_2d(:iim, 1) * teta(:iim, 1, l)) & |
teta(:, 1, l) = SUM(aire_2d(:iim, 1) * teta(:iim, 1, l)) & |
204 |
/ apoln |
/ apoln |
205 |
teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm+1) & |
teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm + 1) & |
206 |
* teta(:iim, jjm + 1, l)) / apols |
* teta(:iim, jjm + 1, l)) / apols |
207 |
END forall |
END forall |
208 |
END IF |
END IF |