9 |
c en faisant Qzon=Cv T + L * ... |
c en faisant Qzon=Cv T + L * ... |
10 |
c vQ..A=Cp T + L * ... |
c vQ..A=Cp T + L * ... |
11 |
|
|
12 |
USE IOIPSL |
USE histcom |
13 |
|
use calendar |
14 |
|
use histwrite_m |
15 |
use dimens_m |
use dimens_m |
16 |
use paramet_m |
use paramet_m |
17 |
use comconst |
use comconst |
142 |
c Initialisation du fichier contenant les moyennes zonales. |
c Initialisation du fichier contenant les moyennes zonales. |
143 |
c --------------------------------------------------------- |
c --------------------------------------------------------- |
144 |
|
|
|
character*10 infile |
|
|
|
|
145 |
integer fileid |
integer fileid |
146 |
integer thoriid, zvertiid |
integer thoriid, zvertiid |
147 |
save fileid |
save fileid |
150 |
|
|
151 |
C Variables locales |
C Variables locales |
152 |
C |
C |
|
integer tau0 |
|
153 |
real zjulian |
real zjulian |
154 |
character*3 str |
character*3 str |
155 |
character*10 ctrac |
character*10 ctrac |
209 |
c Initialisation du fichier contenant les moyennes zonales. |
c Initialisation du fichier contenant les moyennes zonales. |
210 |
c --------------------------------------------------------- |
c --------------------------------------------------------- |
211 |
|
|
|
infile='dynzon' |
|
|
|
|
212 |
zan = annee_ref |
zan = annee_ref |
213 |
dayref = day_ref |
dayref = day_ref |
214 |
CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) |
CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) |
|
tau0 = itau_dyn |
|
215 |
|
|
216 |
rlong=0. |
rlong=0. |
217 |
rlatg=rlatv*180./pi |
rlatg=rlatv*180./pi |
218 |
|
|
219 |
call histbeg_totreg(infile, rlong(:1), rlatg, |
call histbeg_totreg('dynzon', rlong(:1), rlatg, |
220 |
. 1, 1, 1, jjm, |
. 1, 1, 1, jjm, |
221 |
. tau0, zjulian, dt_cum, thoriid, fileid) |
. itau_dyn, zjulian, dt_cum, thoriid, fileid) |
222 |
|
|
223 |
C |
C |
224 |
C Appel a histvert pour la grille verticale |
C Appel a histvert pour la grille verticale |