18 |
USE conf_gcm_m, ONLY: day_step, iperiod |
USE conf_gcm_m, ONLY: day_step, iperiod |
19 |
use conf_guide_m, only: conf_guide, guide_u, guide_v, guide_t, guide_q, & |
use conf_guide_m, only: conf_guide, guide_u, guide_v, guide_t, guide_q, & |
20 |
ncep, ini_anal, tau_min_u, tau_max_u, tau_min_v, tau_max_v, & |
ncep, ini_anal, tau_min_u, tau_max_u, tau_min_v, tau_max_v, & |
21 |
tau_min_t, tau_max_t, tau_min_q, tau_max_q, tau_min_p, tau_max_p, & |
tau_min_t, tau_max_t, tau_min_q, tau_max_q, online |
|
online |
|
22 |
USE dimens_m, ONLY: iim, jjm, llm |
USE dimens_m, ONLY: iim, jjm, llm |
23 |
USE disvert_m, ONLY: ap, bp, preff, presnivs |
USE disvert_m, ONLY: ap, bp, preff, presnivs |
24 |
use dump2d_m, only: dump2d |
use dump2d_m, only: dump2d |
25 |
USE exner_hyb_m, ONLY: exner_hyb |
USE exner_hyb_m, ONLY: exner_hyb |
26 |
USE inigrads_m, ONLY: inigrads |
USE inigrads_m, ONLY: inigrads |
|
use massdair_m, only: massdair |
|
27 |
use netcdf, only: nf90_nowrite, nf90_close, nf90_inq_dimid |
use netcdf, only: nf90_nowrite, nf90_close, nf90_inq_dimid |
28 |
use netcdf95, only: nf95_inquire_dimension, nf95_open |
use netcdf95, only: nf95_inquire_dimension, nf95_open |
29 |
use nr_util, only: pi |
use nr_util, only: pi |
30 |
USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1, jjp1, llmp1 |
USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1 |
31 |
USE q_sat_m, ONLY: q_sat |
USE q_sat_m, ONLY: q_sat |
32 |
use read_reanalyse_m, only: read_reanalyse |
use read_reanalyse_m, only: read_reanalyse |
33 |
USE serre, ONLY: clat, clon |
USE serre, ONLY: clat, clon |
65 |
! alpha=1 signifie pas d'injection |
! alpha=1 signifie pas d'injection |
66 |
! alpha=0 signifie injection totale |
! alpha=0 signifie injection totale |
67 |
REAL, save:: alpha_q(iim + 1, jjm + 1) |
REAL, save:: alpha_q(iim + 1, jjm + 1) |
68 |
REAL, save:: alpha_t(iim + 1, jjm + 1), alpha_p(ip1jmp1) |
REAL, save:: alpha_t(iim + 1, jjm + 1) |
69 |
REAL, save:: alpha_u(iim + 1, jjm + 1), alpha_v(iim + 1, jjm) |
REAL, save:: alpha_u(iim + 1, jjm + 1), alpha_v(iim + 1, jjm) |
70 |
|
|
71 |
INTEGER, save:: step_rea, count_no_rea |
INTEGER, save:: step_rea, count_no_rea |
113 |
|
|
114 |
factt = dtvr * iperiod / daysec |
factt = dtvr * iperiod / daysec |
115 |
|
|
116 |
CALL tau2alpha(3, iip1, jjm, factt, tau_min_v, tau_max_v, alpha_v) |
CALL tau2alpha(3, factt, tau_min_v, tau_max_v, alpha_v) |
117 |
CALL tau2alpha(2, iip1, jjp1, factt, tau_min_u, tau_max_u, alpha_u) |
CALL tau2alpha(2, factt, tau_min_u, tau_max_u, alpha_u) |
118 |
CALL tau2alpha(1, iip1, jjp1, factt, tau_min_t, tau_max_t, alpha_t) |
CALL tau2alpha(1, factt, tau_min_t, tau_max_t, alpha_t) |
119 |
CALL tau2alpha(1, iip1, jjp1, factt, tau_min_p, tau_max_p, alpha_p) |
CALL tau2alpha(1, factt, tau_min_q, tau_max_q, alpha_q) |
|
CALL tau2alpha(1, iip1, jjp1, factt, tau_min_q, tau_max_q, alpha_q) |
|
120 |
|
|
121 |
CALL dump2d(iip1, jjp1, aire, 'AIRE MAILLe ') |
CALL dump2d(iip1, jjp1, aire, 'AIRE MAILLe ') |
122 |
CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ') |
CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ') |
126 |
alpha_t = 0. |
alpha_t = 0. |
127 |
alpha_u = 0. |
alpha_u = 0. |
128 |
alpha_v = 0. |
alpha_v = 0. |
129 |
alpha_p = 0. |
alpha_q = 0. |
130 |
END IF |
END IF |
131 |
|
|
132 |
step_rea = 1 |
step_rea = 1 |