13 |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond |
14 |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & |
15 |
nday, raz_date, periodav, conf_gcm, iflag_phys |
nday, raz_date, periodav, conf_gcm, iflag_phys |
16 |
|
use conf_guide_m, only: conf_guide |
17 |
use dimens_m, only: iim, jjm, llm, nqmx |
use dimens_m, only: iim, jjm, llm, nqmx |
18 |
use dimphy, only: klon |
use dimphy, only: klon |
19 |
USE disvert_m, ONLY : disvert |
USE disvert_m, ONLY : disvert |
36 |
use tracstoke, only: istdyn, istphy |
use tracstoke, only: istdyn, istphy |
37 |
use unit_nml_m, only: unit_nml |
use unit_nml_m, only: unit_nml |
38 |
use yoethf_m, only: yoethf |
use yoethf_m, only: yoethf |
39 |
use write_field, only: NbField, FieldId |
use createnewfield_m, only: NbField, Ncid |
40 |
|
|
41 |
IMPLICIT NONE |
IMPLICIT NONE |
42 |
|
|
87 |
CALL iniconst |
CALL iniconst |
88 |
|
|
89 |
! Lecture du fichier "start.nc" : |
! Lecture du fichier "start.nc" : |
90 |
CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0) |
CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis) |
91 |
|
|
92 |
! On remet le calendrier à zéro si demandé : |
! On remet le calendrier à zéro si demandé : |
93 |
if (raz_date) then |
if (raz_date) then |
96 |
day_ref = dayref |
day_ref = dayref |
97 |
day_ini = dayref |
day_ini = dayref |
98 |
itau_dyn = 0 |
itau_dyn = 0 |
|
time_0 = 0. |
|
99 |
else |
else |
100 |
print *, 'On garde les dates du fichier "start".' |
print *, 'On garde les dates du fichier "start".' |
101 |
endif |
endif |
148 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 |
149 |
istphy = istdyn / iphysiq |
istphy = istdyn / iphysiq |
150 |
|
|
151 |
|
CALL conf_guide |
152 |
|
|
153 |
! Intégration temporelle du modèle : |
! Intégration temporelle du modèle : |
154 |
CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0) |
CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q) |
155 |
|
|
156 |
close(unit_nml) |
close(unit_nml) |
157 |
call histclo |
call histclo |
158 |
|
|
159 |
do i = 1, nbfield |
do i = 1, nbfield |
160 |
call nf95_close(FieldId(i)) |
call nf95_close(Ncid(i)) |
161 |
end do |
end do |
162 |
|
|
163 |
print *, 'Simulation finished' |
print *, 'Simulation finished' |