10 |
|
|
11 |
use clesphys2, only: read_clesphys2 |
use clesphys2, only: read_clesphys2 |
12 |
use com_io_dyn, only: histid, histvid, histaveid |
use com_io_dyn, only: histid, histvid, histaveid |
13 |
use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize |
use comconst, only: daysec, cpp, dtvr, g, rad, r |
14 |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv |
15 |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
16 |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
22 |
use grid_change, only: dyn_phy, init_dyn_phy |
use grid_change, only: dyn_phy, init_dyn_phy |
23 |
use iniadvtrac_m, only: iniadvtrac |
use iniadvtrac_m, only: iniadvtrac |
24 |
use inidissip_m, only: inidissip |
use inidissip_m, only: inidissip |
25 |
|
use inifilr_m, only: inifilr |
26 |
use inigeom_m, only: inigeom |
use inigeom_m, only: inigeom |
27 |
use initdynav_m, only: initdynav |
use initdynav_m, only: initdynav |
28 |
use inithist_m, only: inithist |
use inithist_m, only: inithist |
30 |
use histcom, only: histclo |
use histcom, only: histclo |
31 |
use leapfrog_m, only: leapfrog |
use leapfrog_m, only: leapfrog |
32 |
use logic, only: iflag_phys |
use logic, only: iflag_phys |
|
use paramet_m, only: ip1jm, ip1jmp1 |
|
33 |
use suphec_m, only: suphec |
use suphec_m, only: suphec |
34 |
use temps, only: day_ref, annee_ref, day_end, itau_dyn |
use temps, only: day_ref, annee_ref, day_end, itau_dyn |
35 |
use tracstoke, only: istdyn, istphy |
use tracstoke, only: istdyn, istphy |
40 |
REAL zdtvr ! time step for dynamics, in s |
REAL zdtvr ! time step for dynamics, in s |
41 |
|
|
42 |
! Variables dynamiques : |
! Variables dynamiques : |
43 |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants |
REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm) ! vent covariant |
44 |
REAL teta(ip1jmp1, llm) ! température potentielle |
REAL teta(iim + 1, jjm + 1, llm) ! température potentielle |
45 |
REAL q(ip1jmp1, llm, nqmx) ! champs advectés |
REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés |
46 |
REAL ps(ip1jmp1) ! pression au sol (Pa) |
REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa) |
47 |
|
REAL masse(iim + 1, jjm + 1, llm) ! masse d'air |
|
REAL masse(ip1jmp1, llm) ! masse d'air |
|
48 |
REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol |
REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol |
49 |
|
|
50 |
! Variables pour le fichier histoire : |
! Variables pour le fichier histoire : |
65 |
read (unit=*, nml=main_nml) |
read (unit=*, nml=main_nml) |
66 |
write(unit=*, nml=main_nml) |
write(unit=*, nml=main_nml) |
67 |
|
|
|
! Initialisations: |
|
|
call initialize |
|
|
|
|
68 |
! Choix du calendrier : |
! Choix du calendrier : |
69 |
if (true_calendar) then |
if (true_calendar) then |
70 |
call ioconf_calendar('gregorian') |
call ioconf_calendar('gregorian') |