19 |
use clesphys2, only: read_clesphys2 |
use clesphys2, only: read_clesphys2 |
20 |
use com_io_dyn, only: histid, histvid, histaveid |
use com_io_dyn, only: histid, histvid, histaveid |
21 |
use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize |
use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize |
|
use comdissnew, only: lstardis, nitergdiv, nitergrot, niterh, tetagdiv, & |
|
|
tetagrot, tetatemp |
|
22 |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
23 |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
24 |
nday, raz_date, periodav, conf_gcm |
nday, raz_date, periodav, conf_gcm |
146 |
! Initialisation du filtre : |
! Initialisation du filtre : |
147 |
CALL inifilr |
CALL inifilr |
148 |
|
|
149 |
! Initialisation de la dissipation : |
CALL inidissip |
|
CALL inidissip(lstardis, nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, & |
|
|
tetatemp) |
|
|
|
|
150 |
call init_dyn_phy |
call init_dyn_phy |
151 |
|
|
152 |
! Initialisation de la physique : |
! Initialisation de la physique : |
183 |
|
|
184 |
CALL dynredem0("restart.nc", day_end, phis) |
CALL dynredem0("restart.nc", day_end, phis) |
185 |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
186 |
infile="dyn_hist.nc", t_ops = iecri * daysec, t_wrt = iecri * daysec) |
t_ops = iecri * daysec, t_wrt = iecri * daysec) |
187 |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
188 |
infile='dyn_hist_ave.nc', t_ops = iperiod * zdtvr, & |
t_ops = iperiod * zdtvr, t_wrt = periodav * daysec) |
|
t_wrt = periodav * daysec) |
|
189 |
|
|
190 |
! Choix des fréquences de stockage pour le hors-ligne : |
! Choix des fréquences de stockage pour le hors-ligne : |
191 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |