9 |
! l'advection de "q", en modifiant "iadv" dans "traceur.def". |
! l'advection de "q", en modifiant "iadv" dans "traceur.def". |
10 |
|
|
11 |
use clesphys2, only: read_clesphys2 |
use clesphys2, only: read_clesphys2 |
|
use com_io_dyn, only: histid, histvid, histaveid |
|
12 |
use comconst, only: daysec, cpp, dtvr, g, rad, r |
use comconst, only: daysec, cpp, dtvr, g, rad, r |
13 |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
14 |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
149 |
print *, "day_end = ", day_end |
print *, "day_end = ", day_end |
150 |
|
|
151 |
CALL dynredem0("restart.nc", day_end, phis) |
CALL dynredem0("restart.nc", day_end, phis) |
152 |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, t_ops = iecri * daysec, & |
153 |
t_ops = iecri * daysec, t_wrt = iecri * daysec) |
t_wrt = iecri * daysec) |
154 |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, t_ops = iperiod * zdtvr, & |
155 |
t_ops = iperiod * zdtvr, t_wrt = periodav * daysec) |
t_wrt = periodav * daysec) |
156 |
|
|
157 |
! Choix des fréquences de stockage pour le hors-ligne : |
! Choix des fréquences de stockage pour le hors-ligne : |
158 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |