9 |
! l'advection de "q", en modifiant "iadv" dans "traceur.def". |
! l'advection de "q", en modifiant "iadv" dans "traceur.def". |
10 |
|
|
11 |
use comconst, only: daysec, dtvr, iniconst |
use comconst, only: daysec, dtvr, iniconst |
12 |
|
use comdissnew, only: read_comdissnew |
13 |
use comgeom, only: aire_2d, inigeom |
use comgeom, only: aire_2d, inigeom |
14 |
use comgeomphy, only: airephy |
use comgeomphy, only: airephy |
15 |
use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, periodav, & |
use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, periodav, & |
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 |
57 |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
58 |
|
|
59 |
CALL conf_gcm |
CALL conf_gcm |
60 |
|
call read_comdissnew |
61 |
|
|
62 |
print *, "Enter namelist 'main_nml'." |
print *, "Enter namelist 'main_nml'." |
63 |
read (unit=*, nml=main_nml) |
read (unit=*, nml=main_nml) |
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 |
|
|