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, & |
58 |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
59 |
|
|
60 |
CALL conf_gcm |
CALL conf_gcm |
61 |
|
call read_comdissnew |
62 |
|
|
63 |
print *, "Enter namelist 'main_nml'." |
print *, "Enter namelist 'main_nml'." |
64 |
read (unit=*, nml=main_nml) |
read (unit=*, nml=main_nml) |
81 |
call init_dyn_phy |
call init_dyn_phy |
82 |
|
|
83 |
! Initialisation de la physique : |
! Initialisation de la physique : |
84 |
IF (iflag_phys == 1) THEN |
IF (iflag_phys) THEN |
85 |
airephy = pack(aire_2d, dyn_phy) |
airephy = pack(aire_2d, dyn_phy) |
86 |
CALL suphec |
CALL suphec |
87 |
call yoethf |
call yoethf |