19 |
USE comgeom, ONLY: aire_2d, apoln, apols |
USE comgeom, ONLY: aire_2d, apoln, apols |
20 |
USE disvert_m, ONLY: ap, bp |
USE disvert_m, ONLY: ap, bp |
21 |
USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, & |
USE conf_gcm_m, ONLY: day_step, iconser, iperiod, iphysiq, nday, offline, & |
22 |
iflag_phys, ok_guide |
iflag_phys, ok_guide, iecri |
23 |
USE dimens_m, ONLY: iim, jjm, llm, nqmx |
USE dimens_m, ONLY: iim, jjm, llm, nqmx |
24 |
use dissip_m, only: dissip |
use dissip_m, only: dissip |
25 |
USE dynetat0_m, ONLY: day_ini |
USE dynetat0_m, ONLY: day_ini |
229 |
call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, & |
call bilan_dyn(ps, masse, pk, pbaru, pbarv, teta, phi, ucov, vcov, & |
230 |
q(:, :, :, 1)) |
q(:, :, :, 1)) |
231 |
ENDIF |
ENDIF |
232 |
|
|
233 |
|
IF (MOD(itau + 1, iecri * day_step) == 0) THEN |
234 |
|
CALL geopot((iim + 1) * (jjm + 1), teta, pk, pks, phis, phi) |
235 |
|
CALL writehist(itau, vcov, ucov, teta, phi, q, masse, ps, phis) |
236 |
|
END IF |
237 |
end do time_integration |
end do time_integration |
238 |
|
|
239 |
CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps, & |
CALL dynredem1("restart.nc", vcov, ucov, teta, q, masse, ps, & |