source: LMDZ4/trunk/libf/phylmd/write_histmth.h @ 602

Last change on this file since 602 was 602, checked in by Laurent Fairhead, 19 years ago

Sortie de prw du if sur iflag_con FH
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 30.7 KB
Line 
1!
2! $Header$
3!
4      IF (ok_mensuel) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8
9         zsto = dtime
10         zout = dtime * ecrit_mth
11         itau_w = itau_phy + itap
12
13      i = NINT(zout/zsto)
14c
15      IF(lev_histmth.GE.1) THEN
16c
17c Champs 2D:
18c
19      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
20      CALL histwrite(nid_mth,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
21C
22      i = NINT(zout/zsto)
23      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
24      CALL histwrite(nid_mth,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
25c
26      zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_ter)
27      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
28      CALL histwrite(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
29     $      zx_tmp_2d,iim*jjmp1,ndex2d)
30c
31      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
32      CALL histwrite(nid_mth,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
33
34      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
35      CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
36c
37      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
38      CALL histwrite(nid_mth,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
39c
40c ENSEMBLES BEG
41      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
42      CALL histwrite(nid_mth,"t2m_min",itau_w,zx_tmp_2d,
43     .               iim*jjmp1,ndex2d)
44c
45      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
46      CALL histwrite(nid_mth,"t2m_max",itau_w,zx_tmp_2d,
47     .               iim*jjmp1,ndex2d)
48c
49c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
50c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
51c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
52c     CALL histwrite(nid_mth,"tsoil",itau_w,zx_tmp_2d,
53c    .               iim*jjmp1,ndex2d)
54c
55      DO i=1, klon
56       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
57      ENDDO
58      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
59      CALL histwrite(nid_mth,"wind10m",itau_w,zx_tmp_2d,
60     $     iim*jjmp1,ndex2d)
61c
62      DO i = 1, klon
63         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
64      ENDDO
65      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
66      CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
67c
68c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
69c     CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
70c    .                                   iim*jjmp1*klev,ndex3d)
71c
72c ENSEMBLES END
73c
74      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
75      CALL histwrite(nid_mth,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
76c
77      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
78      CALL histwrite(nid_mth,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
79C
80      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
81      CALL histwrite(nid_mth,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
82C
83      DO i = 1, klon
84         zx_tmp_fi2d(i) = paprs(i,1)
85      ENDDO
86      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
87      CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
88c
89      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
90      CALL histwrite(nid_mth,"qsurf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
91c
92      if (.not. ok_veget) then
93          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
94          CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1
95     $        ,ndex2d)
96      endif
97c
98      DO i = 1, klon
99         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
100      ENDDO
101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
102      CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
103c
104      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
105      CALL histwrite(nid_mth,"ndayrain",itau_w,zx_tmp_2d,
106     $               iim*jjmp1,ndex2d)
107c
108      DO i = 1, klon
109         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
110      ENDDO
111      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
112      CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
113c
114      DO i = 1, klon
115         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
116      ENDDO
117      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
118      CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
119c
120      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
121      CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
122c
123      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
124      CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
125     .               ndex2d)
126c
127cIM: 071003
128c     zx_tmp_fi2d(1:klon)=evap(1:klon)*86400.
129c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
130      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
131      CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
132c
133      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
134      CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
135c
136      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
137      CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
138c
139      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
140      CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
141c
142      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
143      CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
144c
145c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
146cIM 080304   zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 2 )
147      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
148      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
149      CALL histwrite(nid_mth, "SWupTOA",itau_w,zx_tmp_2d,
150     .                               iim*jjmp1,ndex2d)
151c
152c
153c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
154cIM 080304   zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 2 )
155      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
156      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
157      CALL histwrite(nid_mth, "SWupTOAclr",itau_w,zx_tmp_2d,
158     .                               iim*jjmp1,ndex2d)
159c
160c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
161cIM 080304   zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 2 )
162      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
163      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
164      CALL histwrite(nid_mth, "SWdnTOA",itau_w,zx_tmp_2d,
165     .                               iim*jjmp1,ndex2d)
166c
167c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
168cIM 080304   zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 2 )
169      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
170      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
171      CALL histwrite(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_2d,
172     .                               iim*jjmp1,ndex2d)
173c
174      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
175      CALL histwrite(nid_mth,"SWup200",itau_w,zx_tmp_2d,
176     $     iim*jjmp1,ndex2d)
177c
178      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
179      CALL histwrite(nid_mth,"SWup200clr",itau_w,zx_tmp_2d,
180     $     iim*jjmp1,ndex2d)
181c
182      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
183      CALL histwrite(nid_mth,"SWdn200",itau_w,zx_tmp_2d,
184     $     iim*jjmp1,ndex2d)
185c
186      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
187      CALL histwrite(nid_mth,"SWdn200clr",itau_w,zx_tmp_2d,
188     $     iim*jjmp1,ndex2d)
189c
190      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
191      CALL histwrite(nid_mth,"LWup200",itau_w,zx_tmp_2d,
192     $     iim*jjmp1,ndex2d)
193c
194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
195      CALL histwrite(nid_mth,"LWup200clr",itau_w,zx_tmp_2d,
196     $     iim*jjmp1,ndex2d)
197c
198      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
199      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
200      CALL histwrite(nid_mth,"LWdn200",itau_w,zx_tmp_2d,
201     $     iim*jjmp1,ndex2d)
202c
203      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
204      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
205      CALL histwrite(nid_mth,"LWdn200clr",itau_w,zx_tmp_2d,
206     $     iim*jjmp1,ndex2d)
207c
208c200 END
209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
210      CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
211c
212      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
213      CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
214c
215      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
216      CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
217c
218      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
219      CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
220c
221c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
222      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
223      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
224      CALL histwrite(nid_mth, "SWupSFC",itau_w,zx_tmp_2d,
225     .                               iim*jjmp1,ndex2d)
226c
227c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
228      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
229      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
230      CALL histwrite(nid_mth, "SWupSFCclr",itau_w,zx_tmp_2d,
231     .                               iim*jjmp1,ndex2d)
232c
233c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
234      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
235      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
236      CALL histwrite(nid_mth, "SWdnSFC",itau_w,zx_tmp_2d,
237     .                               iim*jjmp1,ndex2d)
238c
239c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
240      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
241      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
242      CALL histwrite(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_2d,
243     .                               iim*jjmp1,ndex2d)
244c
245      zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
246      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
247      CALL histwrite(nid_mth,"LWupSFC",itau_w,zx_tmp_2d,iim*jjmp1,
248     .               ndex2d)
249c
250      zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
251      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
252      CALL histwrite(nid_mth,"LWupSFCclr",itau_w,zx_tmp_2d,iim*jjmp1,
253     .               ndex2d)
254c
255      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
256      CALL histwrite(nid_mth,"LWdnSFC",itau_w,zx_tmp_2d,iim*jjmp1,
257     .               ndex2d)
258c
259      DO i=1, klon
260       sollwdownclr(i)=lwdn0(i,1)
261      ENDDO !i=1, klon
262c
263      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
264      CALL histwrite(nid_mth,"LWdnSFCclr",itau_w,zx_tmp_2d,
265     $               iim*jjmp1,ndex2d)
266c
267c
268c effets des aerosols
269c
270c     IF (ok_ade.OR.ok_aie) THEN
271      zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
272c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswad,zx_tmp_2d)
273      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
274      CALL histwrite(nid_mth,"topsad",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
275c
276      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
277c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswad,zx_tmp_2d)
278      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
279      CALL histwrite(nid_mth,"solsad",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
280c
281      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
282c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topswai,zx_tmp_2d)
283      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
284      CALL histwrite(nid_mth,"topsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
285c
286      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
287c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solswai,zx_tmp_2d)
288      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d ,zx_tmp_2d)
289      CALL histwrite(nid_mth,"solsai",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
290c     endif
291c
292      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
293      CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
294c
295      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
296c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
297      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
298      CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
299c
300      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
301      CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
302c
303c
304c      DO i = 1, klon
305c         zx_tmp_fi2d(i) = fluxu(i,1)
306c      ENDDO
307c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
308c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
309c
310c      DO i = 1, klon
311c         zx_tmp_fi2d(i) = fluxv(i,1)
312c      ENDDO
313c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
314c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
315c
316c
317       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
318       CALL histwrite(nid_mth,"ffonte",itau_w,zx_tmp_2d,iim*jjmp1,
319     $                ndex2d)
320c
321       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
322       CALL histwrite(nid_mth,"fqcalving",itau_w,zx_tmp_2d,iim*jjmp1,
323     $                ndex2d)
324cIM: 171003
325       DO nsrf = 1, nbsrf
326        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
327        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
328        CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
329     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
330C     
331        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
332        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
333        CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
334     $      zx_tmp_2d,iim*jjmp1,ndex2d)
335       ENDDO
336cIM: 171003
337c
338cIM      if ( lev_histday.gt.1 ) then
339      DO nsrf = 1, nbsrf
340CYYY
341       IF(nsrf.GT.1) THEN
342        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
343        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
344        CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
345     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
346       ENDIF !nsrf.GT.1
347C
348        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
349        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
350        CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
351     $      zx_tmp_2d,iim*jjmp1,ndex2d)
352C
353        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
354        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
355        CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
356     $      zx_tmp_2d,iim*jjmp1,ndex2d)
357C
358        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
359        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
360        CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
361     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
362C
363        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
364        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
365        CALL histwrite(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
366     $      zx_tmp_2d,iim*jjmp1,ndex2d)
367c
368        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
369        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
370        CALL histwrite(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
371     $      zx_tmp_2d,iim*jjmp1,ndex2d)
372C
373        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
374        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
375        CALL histwrite(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
376     $      zx_tmp_2d,iim*jjmp1,ndex2d)
377C
378
379      END DO 
380cIM      endif !lev_histday
381
382c
383      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
384      CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
385c
386      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
387      CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
388c
389      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
390      CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
391c
392      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
393      CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
394c
395      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
396      CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
397c
398      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
399      CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
400c
401      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
402      CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
403c
404      zx_tmp_fi2d(1:klon) = flwp(1:klon)
405      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
406      CALL histwrite(nid_mth,"lwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
407c
408      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
409      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
410      CALL histwrite(nid_mth,"iwp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
411c
412      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
413      CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
414c
415      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
416      CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
417c
418      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
419      CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
420c
421      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
422      CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
423cKE43
424      IF (iflag_con.GE.3) THEN ! sb
425c
426      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
427      CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
428c
429      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
430      CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
431c
432      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
433      CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
434c
435      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
436      CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
437c
438c34EK
439c
440      ENDIF !iflag_con.GE.3
441c
442      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
443      CALL histwrite(nid_mth,"prw",itau_w,zx_tmp_2d,
444     .               iim*jjmp1,ndex2d)
445c
446c
447c
448c
449c  Champs interpolles sur des niveaux de pression
450      DO k=1, nlevENS
451c
452        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
453        CALL histwrite(nid_mth,"u"//clev(k),itau_w,zx_tmp_2d,
454     $       iim*jjmp1,ndex2d)
455c
456        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
457        CALL histwrite(nid_mth,"v"//clev(k),itau_w,zx_tmp_2d,
458     $       iim*jjmp1,ndex2d)
459c
460      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,k),zx_tmp_2d)
461      CALL histwrite(nid_mth,"w"//clev(k),itau_w,zx_tmp_2d,
462     $     iim*jjmp1,ndex2d)
463c
464      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
465      CALL histwrite(nid_mth,"phi"//clev(k),itau_w,zx_tmp_2d,
466     $     iim*jjmp1,ndex2d)
467c
468      ENDDO !nlevENS
469      ENDIF !lev_histmth.GE.1
470c
471      IF(lev_histmth.GE.2) THEN
472c
473c Champs 3D:
474C
475c
476      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
477      CALL histwrite(nid_mth,"lwcon",itau_w,zx_tmp_3d,
478     .               iim*jjmp1*klev,ndex3d)
479c
480      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
481      CALL histwrite(nid_mth,"iwcon",itau_w,zx_tmp_3d,
482     .               iim*jjmp1*klev,ndex3d)
483c
484      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
485      CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d,
486     .                                   iim*jjmp1*klev,ndex3d)
487c
488c#ifdef histmthNMC
489c#include "write_histmthNMC.h"
490c#endif
491c
492      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
493      CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d,
494     .                                   iim*jjmp1*klev,ndex3d)
495c
496cIM: 071003
497c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
498c    .                         (1-qx(1:klon,1:klev,ivap))
499c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
500c     CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
501c    .                                   iim*jjmp1*klev,ndex3d)
502c
503      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
504      CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d,
505     .                                   iim*jjmp1*klev,ndex3d)
506c
507      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
508      CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d,
509     .                                   iim*jjmp1*klev,ndex3d)
510c
511      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
512      CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d,
513     .                                   iim*jjmp1*klev,ndex3d)
514c
515      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
516      CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d,
517     .                                   iim*jjmp1*klev,ndex3d)
518c
519      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
520      CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d,
521     .                                   iim*jjmp1*klev,ndex3d)
522c
523      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
524      CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d,
525     .                                   iim*jjmp1*klev,ndex3d)
526c
527      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
528      CALL histwrite(nid_mth,"rnebcon",itau_w,zx_tmp_3d,
529     .                                   iim*jjmp1*klev,ndex3d)
530c
531      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
532      CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d,
533     .                                   iim*jjmp1*klev,ndex3d)
534c
535      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
536      CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d,
537     .                                   iim*jjmp1*klev,ndex3d)
538c
539cIM 22.03.04 BEG
540      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
541      CALL histwrite(nid_mth,"dtphy",itau_w,zx_tmp_3d,
542     .                                   iim*jjmp1*klev,ndex3d)
543c
544      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap), 
545     .                 zx_tmp_3d)
546      CALL histwrite(nid_mth,"dqphy",itau_w,zx_tmp_3d,
547     .                                   iim*jjmp1*klev,ndex3d)
548cIM 22.03.04 END
549c
550      ENDIF !lev_histmth.GE.2
551c
552      IF(lev_histmth.GE.3) THEN
553c
554       DO nsrf=1, nbsrf
555c
556        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
557        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
558        CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
559     $      zx_tmp_2d,iim*jjmp1,ndex2d)
560C
561        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
562        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
563        CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
564     $      zx_tmp_2d,iim*jjmp1,ndex2d)
565c
566      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
567      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
568      CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w
569     $    ,zx_tmp_2d,iim*jjmp1,ndex2d)
570c
571      ENDDO !nsrf=1, nbsrf
572c
573      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
574      CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
575c
576      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
577      CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
578c
579      ENDIF !lev_histmth.GE.3
580c
581c
582      IF(lev_histmth.GE.4) THEN
583c
584      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
585      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
586     .                                   iim*jjmp1*klev,ndex3d)
587c
588      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
589      CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d,
590     .                                   iim*jjmp1*klev,ndex3d)
591c
592      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
593      CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d,
594     .                                   iim*jjmp1*klev,ndex3d)
595c
596      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
597      CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d,
598     .                                   iim*jjmp1*klev,ndex3d)
599c
600      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
601      CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d,
602     .                                   iim*jjmp1*klev,ndex3d)
603c
604c
605      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
606      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
607     .                                   iim*jjmp1*klev,ndex3d)
608c
609      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
610      CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d,
611     .                                   iim*jjmp1*klev,ndex3d)
612c
613cIM: 101003 : K/30min ==> K/s
614      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
615      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
616      CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d,
617     .                                   iim*jjmp1*klev,ndex3d)
618c
619      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d)
620      CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d,
621     .                                   iim*jjmp1*klev,ndex3d)
622c
623cIM: 101003 : K/30min ==> K/s
624      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
625      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
626      CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d,
627     .                                   iim*jjmp1*klev,ndex3d)
628cIM: 071003
629cIM: 101003 : K/30min ==> K/s
630      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
631     .                             d_t_eva(1:klon,1:klev))/pdtphys
632      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
633      CALL histwrite(nid_mth,"dtlschr",itau_w,zx_tmp_3d,
634     .                                   iim*jjmp1*klev,ndex3d)
635c
636      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d)
637      CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d,
638     .                                   iim*jjmp1*klev,ndex3d)
639c
640cIM: 101003 : K/30min ==> K/s
641      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
642      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
643      CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d,
644     .                                   iim*jjmp1*klev,ndex3d)
645c
646      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
647      CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d,
648     .                                   iim*jjmp1*klev,ndex3d)
649c
650cIM: 101003 : K/30min ==> K/s
651      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
652      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
653      CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d,
654     .                                   iim*jjmp1*klev,ndex3d)
655c
656      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d)
657      CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d,
658     .                                   iim*jjmp1*klev,ndex3d)
659c
660      zpt_conv = 0.
661      where (ptconv) zpt_conv = 1.
662      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
663      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
664     .                                   iim*(jjmp1)*klev,ndex3d)
665c
666      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
667      CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d,
668     .                                   iim*(jjmp1)*klev,ndex3d)
669c
670cIM: 101003 : K/30min ==> K/s
671      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
672      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
673      CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d,
674     .                                   iim*jjmp1*klev,ndex3d)
675c
676      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d)
677      CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d,
678     .                                   iim*jjmp1*klev,ndex3d)
679c
680cIM: 101003 : K/day ==> K/s
681cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)*pdtphys/RDAY
682      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
683      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
684      CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d,
685     .                                   iim*jjmp1*klev,ndex3d)
686c
687cIM: 101003 : K/day ==> K/s     
688cIM: LMD_ARMIP3   zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)*pdtphys/RDAY
689      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
690      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
691      CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d,
692     .                                   iim*jjmp1*klev,ndex3d)
693c
694cIM: 101003 : K/day ==> K/s     
695cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)*pdtphys/RDAY
696      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
697      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
698      CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d,
699     .                                   iim*jjmp1*klev,ndex3d)
700c
701cIM: 101003 : K/day ==> K/s     
702cIM: LMD_ARMIP3     zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)*pdtphys/RDAY
703      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
704      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
705      CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d,
706     .                                   iim*jjmp1*klev,ndex3d)
707c
708cIM: 101003 : deja en K/s     
709      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
710      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
711      CALL histwrite(nid_mth,"dtec",itau_w,zx_tmp_3d,
712     .                                   iim*jjmp1*klev,ndex3d)
713c
714      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
715      CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d,
716     .                                   iim*jjmp1*klev,ndex3d)
717c
718      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
719      CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d,
720     .                                   iim*jjmp1*klev,ndex3d)
721c
722      IF (ok_orodr) THEN
723      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d)
724      CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d,
725     .                                   iim*jjmp1*klev,ndex3d)
726c
727      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d)
728      CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d,
729     .                                   iim*jjmp1*klev,ndex3d)
730c
731      ENDIF
732C
733      IF (ok_orolf) THEN
734      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d)
735      CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d,
736     .                                   iim*jjmp1*klev,ndex3d)
737c
738      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d)
739      CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d,
740     .                                   iim*jjmp1*klev,ndex3d)
741      ENDIF
742C
743c
744c effets des aerosols
745c
746c     IF (ok_ade.OR.ok_aie) THEN
747      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, re, zx_tmp_3d)
748      CALL histwrite(nid_mth,"re",itau_w,zx_tmp_3d,
749     .                                   iim*jjmp1*klev,ndex3d)
750c
751      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fl, zx_tmp_3d)
752      CALL histwrite(nid_mth,"redenom",itau_w,zx_tmp_3d,
753     .                                   iim*jjmp1*klev,ndex3d)
754c
755      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau, zx_tmp_3d)
756      CALL histwrite(nid_mth,"tau",itau_w,zx_tmp_3d,
757     .                                   iim*jjmp1*klev,ndex3d)
758c
759      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtaupi, zx_tmp_3d)
760      CALL histwrite(nid_mth,"taupi",itau_w,zx_tmp_3d,
761     .                                   iim*jjmp1*klev,ndex3d)
762c     endif
763c
764      IF (nqmax.GE.3) THEN
765      DO iq=1,nqmax-2
766      IF (iq.LE.99) THEN
767         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d)
768         WRITE(str2,'(i2.2)') iq
769         CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d,
770     .                                   iim*jjmp1*klev,ndex3d)
771      ELSE
772         PRINT*, "Trop de traceurs"
773         CALL abort
774      ENDIF
775      ENDDO
776      ENDIF
777c
778      ENDIF !lev_histmth.GE.4
779c
780      if (ok_sync) then
781        call histsync(nid_mth)
782      endif
783      ENDIF
Note: See TracBrowser for help on using the repository browser.