30 |
USE q_sat_m, ONLY: q_sat |
USE q_sat_m, ONLY: q_sat |
31 |
use read_reanalyse_m, only: read_reanalyse |
use read_reanalyse_m, only: read_reanalyse |
32 |
USE serre, ONLY: clat, clon |
USE serre, ONLY: clat, clon |
33 |
use tau2alpha_m, only: tau2alpha, dxdys |
use tau2alpha_m, only: tau2alpha |
34 |
use writefield_m, only: writefield |
use writefield_m, only: writefield |
35 |
|
|
36 |
INTEGER, INTENT(IN):: itau |
INTEGER, INTENT(IN):: itau |
71 |
|
|
72 |
INTEGER ilon, ilat |
INTEGER ilon, ilat |
73 |
REAL factt ! pas de temps entre deux appels au guidage, en fraction de jour |
REAL factt ! pas de temps entre deux appels au guidage, en fraction de jour |
|
real ztau(iim + 1, jjm + 1) |
|
74 |
|
|
75 |
INTEGER ij, l |
INTEGER ij, l |
76 |
INTEGER ncid, dimid |
INTEGER ncid, dimid |
138 |
CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, & |
CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, & |
139 |
masserea2, nlev) |
masserea2, nlev) |
140 |
qrea2 = max(qrea2, 0.1) |
qrea2 = max(qrea2, 0.1) |
141 |
|
|
142 |
|
CALL writefield("alpha_u", alpha_u) |
143 |
|
CALL writefield("alpha_t", alpha_t) |
144 |
END IF first_call |
END IF first_call |
145 |
|
|
146 |
! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS: |
! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS: |
159 |
masserea2, nlev) |
masserea2, nlev) |
160 |
qrea2 = max(qrea2, 0.1) |
qrea2 = max(qrea2, 0.1) |
161 |
factt = dtvr * iperiod / daysec |
factt = dtvr * iperiod / daysec |
|
ztau = factt / max(alpha_t, 1E-10) |
|
|
CALL writefield("aire", aire) |
|
|
CALL writefield("dxdys", dxdys) |
|
|
CALL writefield("alpha_u", alpha_u) |
|
|
CALL writefield("alpha_t", alpha_t) |
|
|
CALL writefield("ztau", ztau) |
|
162 |
CALL writefield("ucov", ucov) |
CALL writefield("ucov", ucov) |
163 |
CALL writefield("ucovrea2", ucovrea2) |
CALL writefield("ucovrea2", ucovrea2) |
164 |
CALL writefield("teta", teta) |
CALL writefield("teta", teta) |