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 |
20 |
use temps |
use temps |
21 |
use iniprint |
use iniprint |
22 |
use inigrads_m, only: inigrads |
use inigrads_m, only: inigrads |
23 |
|
use nr_util, only: pi |
24 |
|
|
25 |
IMPLICIT NONE |
IMPLICIT NONE |
26 |
|
|
143 |
c Initialisation du fichier contenant les moyennes zonales. |
c Initialisation du fichier contenant les moyennes zonales. |
144 |
c --------------------------------------------------------- |
c --------------------------------------------------------- |
145 |
|
|
|
character*10 infile |
|
|
|
|
146 |
integer fileid |
integer fileid |
147 |
integer thoriid, zvertiid |
integer thoriid, zvertiid |
148 |
save fileid |
save fileid |
151 |
|
|
152 |
C Variables locales |
C Variables locales |
153 |
C |
C |
|
integer tau0 |
|
154 |
real zjulian |
real zjulian |
155 |
character*3 str |
character*3 str |
156 |
character*10 ctrac |
character*10 ctrac |
186 |
|
|
187 |
if (i_sortie.eq.1) then |
if (i_sortie.eq.1) then |
188 |
file='dynzon' |
file='dynzon' |
189 |
call inigrads(ifile |
call inigrads(ifile ,(/0./),180./pi,0.,0.,rlatv,-90.,90., |
190 |
s ,(/0./),180./pi,0.,0.,rlatv,-90.,90.,180./pi |
$ 180./pi ,presnivs,1. ,dt_cum,file,'dyn_zon ') |
|
s ,presnivs,1. |
|
|
s ,dt_cum,file,'dyn_zon ') |
|
191 |
endif |
endif |
192 |
|
|
193 |
nom(itemp)='T' |
nom(itemp)='T' |
210 |
c Initialisation du fichier contenant les moyennes zonales. |
c Initialisation du fichier contenant les moyennes zonales. |
211 |
c --------------------------------------------------------- |
c --------------------------------------------------------- |
212 |
|
|
|
infile='dynzon' |
|
|
|
|
213 |
zan = annee_ref |
zan = annee_ref |
214 |
dayref = day_ref |
dayref = day_ref |
215 |
CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) |
CALL ymds2ju(zan, 1, dayref, 0.0, zjulian) |
|
tau0 = itau_dyn |
|
216 |
|
|
217 |
rlong=0. |
rlong=0. |
218 |
rlatg=rlatv*180./pi |
rlatg=rlatv*180./pi |
219 |
|
|
220 |
call histbeg_totreg(infile, rlong(:1), rlatg, |
call histbeg_totreg('dynzon', rlong(:1), rlatg, |
221 |
. 1, 1, 1, jjm, |
. 1, 1, 1, jjm, |
222 |
. tau0, zjulian, dt_cum, thoriid, fileid) |
. itau_dyn, zjulian, dt_cum, thoriid, fileid) |
223 |
|
|
224 |
C |
C |
225 |
C Appel a histvert pour la grille verticale |
C Appel a histvert pour la grille verticale |