54 |
|
|
55 |
!======================================================================= |
!======================================================================= |
56 |
|
|
57 |
! ... Dans inigeom , nouveaux calculs pour les elongations cu , cv |
! Dans inigeom , nouveaux calculs pour les elongations cu , cv |
58 |
! et possibilite d'appeler une fonction f(y) a derivee tangente |
! et possibilite d'appeler une fonction f(y) a derivee tangente |
59 |
! hyperbolique a la place de la fonction a derivee sinusoidale. |
! hyperbolique a la place de la fonction a derivee sinusoidale. |
60 |
|
|
61 |
! ... Possibilite de choisir le shema de Van-leer pour l'advection de |
! ... Possibilite de choisir le shema de Van-leer pour l'advection de |
62 |
! q , en faisant iadv = 10 dans traceur (29/04/97) . |
! q , en faisant iadv = 10 dans traceur (29/04/97) . |
249 |
PRINT *, 'nlev', nlev |
PRINT *, 'nlev', nlev |
250 |
rcod = nf90_close(ncidpl) |
rcod = nf90_close(ncidpl) |
251 |
! Lecture du premier etat des reanalyses. |
! Lecture du premier etat des reanalyses. |
252 |
CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, masserea2, & |
CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, & |
253 |
psrea2, 1, nlev) |
masserea2, psrea2, 1, nlev) |
254 |
qrea2(:, :) = max(qrea2(:, :), 0.1) |
qrea2(:, :) = max(qrea2(:, :), 0.1) |
255 |
|
|
256 |
|
|
286 |
count_no_rea, ' non lectures' |
count_no_rea, ' non lectures' |
287 |
step_rea = step_rea + 1 |
step_rea = step_rea + 1 |
288 |
itau_test = itau |
itau_test = itau |
289 |
CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, qrea2, & |
CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, & |
290 |
masserea2, psrea2, 1, nlev) |
qrea2, masserea2, psrea2, 1, nlev) |
291 |
qrea2(:, :) = max(qrea2(:, :), 0.1) |
qrea2(:, :) = max(qrea2(:, :), 0.1) |
292 |
factt = dtvr*iperiod/daysec |
factt = dtvr*iperiod/daysec |
293 |
ztau(:) = factt/max(alpha_t(:), 1.E-10) |
ztau(:) = factt/max(alpha_t(:), 1.E-10) |