32 |
use leapfrog_m, only: leapfrog |
use leapfrog_m, only: leapfrog |
33 |
use netcdf95, only: nf95_close |
use netcdf95, only: nf95_close |
34 |
use suphec_m, only: suphec |
use suphec_m, only: suphec |
|
use tracstoke, only: istdyn, istphy |
|
35 |
use unit_nml_m, only: unit_nml |
use unit_nml_m, only: unit_nml |
36 |
use yoethf_m, only: yoethf |
use yoethf_m, only: yoethf |
37 |
use createnewfield_m, only: NbField, Ncid |
use createnewfield_m, only: NbField, Ncid |
80 |
call init_dyn_phy |
call init_dyn_phy |
81 |
|
|
82 |
! Initialisation de la physique : |
! Initialisation de la physique : |
83 |
IF (iflag_phys == 1) THEN |
IF (iflag_phys) THEN |
84 |
airephy = pack(aire_2d, dyn_phy) |
airephy = pack(aire_2d, dyn_phy) |
85 |
CALL suphec |
CALL suphec |
86 |
call yoethf |
call yoethf |
92 |
CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec) |
CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec) |
93 |
call init_dynzon(dt_app = dtvr * iperiod) |
call init_dynzon(dt_app = dtvr * iperiod) |
94 |
|
|
|
! Choix des fr\'equences de stockage pour le hors-ligne : |
|
|
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
|
|
istphy = istdyn / iphysiq |
|
|
|
|
95 |
CALL conf_guide |
CALL conf_guide |
96 |
CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q) |
CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q) |
97 |
|
|