39 |
use grid_change, only: dyn_phy, init_dyn_phy |
use grid_change, only: dyn_phy, init_dyn_phy |
40 |
use advtrac_m, only: iniadvtrac |
use advtrac_m, only: iniadvtrac |
41 |
use leapfrog_m, only: leapfrog |
use leapfrog_m, only: leapfrog |
42 |
|
use dynredem0_m, only: dynredem0 |
43 |
|
|
44 |
IMPLICIT NONE |
IMPLICIT NONE |
45 |
|
|
53 |
REAL ps(ip1jmp1) ! pression au sol (Pa) |
REAL ps(ip1jmp1) ! pression au sol (Pa) |
54 |
|
|
55 |
REAL masse(ip1jmp1, llm) ! masse d'air |
REAL masse(ip1jmp1, llm) ! masse d'air |
56 |
REAL phis(ip1jmp1) ! géopotentiel au sol |
REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol |
57 |
|
|
58 |
! Variables pour le fichier histoire : |
! Variables pour le fichier histoire : |
59 |
REAL time_0 |
REAL time_0 |
205 |
print *, "day_ini = ", day_ini |
print *, "day_ini = ", day_ini |
206 |
print *, "day_end = ", day_end |
print *, "day_end = ", day_end |
207 |
|
|
208 |
CALL dynredem0("restart.nc", day_end, phis, nqmx) |
CALL dynredem0("restart.nc", day_end, phis) |
209 |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, & |
210 |
infile="dyn_hist.nc", t_ops = iecri * daysec, t_wrt = iecri * daysec) |
infile="dyn_hist.nc", t_ops = iecri * daysec, t_wrt = iecri * daysec) |
211 |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |
CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, & |