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 |
32 |
use inigeom_m, only: inigeom |
use inigeom_m, only: inigeom |
33 |
use initdynav_m, only: initdynav |
use initdynav_m, only: initdynav |
34 |
use inithist_m, only: inithist |
use inithist_m, only: inithist |
35 |
USE IOIPSL, only: ioconf_calendar, histclo |
USE calendar, only: ioconf_calendar |
36 |
|
use histcom, only: histclo |
37 |
use leapfrog_m, only: leapfrog |
use leapfrog_m, only: leapfrog |
38 |
use logic, only: iflag_phys |
use logic, only: iflag_phys |
39 |
use paramet_m, only: ip1jm, ip1jmp1 |
use paramet_m, only: ip1jm, ip1jmp1 |
147 |
! Initialisation du filtre : |
! Initialisation du filtre : |
148 |
CALL inifilr |
CALL inifilr |
149 |
|
|
150 |
! Initialisation de la dissipation : |
CALL inidissip |
|
CALL inidissip(lstardis, nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, & |
|
|
tetatemp) |
|
|
|
|
151 |
call init_dyn_phy |
call init_dyn_phy |
152 |
|
|
153 |
! Initialisation de la physique : |
! Initialisation de la physique : |
184 |
|
|
185 |
CALL dynredem0("restart.nc", day_end, phis) |
CALL dynredem0("restart.nc", day_end, phis) |
186 |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
187 |
infile="dyn_hist.nc", t_ops = iecri * daysec, t_wrt = iecri * daysec) |
t_ops = iecri * daysec, t_wrt = iecri * daysec) |
188 |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
189 |
infile='dyn_hist_ave.nc', t_ops = iperiod * zdtvr, & |
t_ops = iperiod * zdtvr, t_wrt = periodav * daysec) |
|
t_wrt = periodav * daysec) |
|
190 |
|
|
191 |
! Choix des fréquences de stockage pour le hors-ligne : |
! Choix des fréquences de stockage pour le hors-ligne : |
192 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |