source: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h @ 4116

Last change on this file since 4116 was 925, checked in by Laurent Fairhead, 17 years ago

Erreur lors de la reintegration des sorties aerosols SD
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 69.8 KB
Line 
1c
2c $Header$
3c
4      IF (ok_mensuel) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8
9         itau_w = itau_phy + itap
10
11c
12      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
13c
14      IF(lev_histmth.GE.1) THEN
15c
16c Champs 2D:
17c
18cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
19      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis)
20C
21cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
22cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
23      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
24c
25      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
26cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
27      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
28     $                   zx_tmp_fi2d)
29c
30      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
31cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
32      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
33     $                   zx_tmp_fi2d)
34c
35cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
36      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
37
38cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
39      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
40c
41cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
42      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
43c
44c ENSEMBLES BEG
45cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
46      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
47c
48cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
49      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
50c
51c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
52c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
53c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
54c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
55c    .               iim*jjmp1,ndex2d)
56c
57      DO i=1, klon
58       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
59      ENDDO
60cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
61      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
62c
63      DO i = 1, klon
64         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
65      ENDDO
66cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
67      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
68c
69c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
70c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
71c    .                                   iim*jjmp1*klev,ndex3d)
72c
73c ENSEMBLES END
74c
75cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
76      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
77c
78cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
79      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
80C
81cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
82      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
83C
84      DO i = 1, klon
85         zx_tmp_fi2d(i) = paprs(i,1)
86      ENDDO
87cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
88      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
89c
90cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
91      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
92c
93      if (.not. ok_veget) then
94cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
95          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
96      endif
97c
98      DO i = 1, klon
99         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
100      ENDDO
101cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
102      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
103c
104cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
105      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
106c
107      DO i = 1, klon
108         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
109      ENDDO
110cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
111      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
112c
113      DO i = 1, klon
114         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
115      ENDDO
116cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
117      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
118c
119cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
120      CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall)
121c
122c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
123c     CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
124c    .               ndex2d)
125c
126cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
127      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
128c
129cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
130      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
131c
132cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
133      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
134c
135cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
136      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
137c
138cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
139      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
140c
141      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
142cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
143      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
144c
145c
146      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
147cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
148      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
149c
150      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
151cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
152      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
153c
154      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
155cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
156      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
157c
158cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
159      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
160c
161cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
162      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
163c
164cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
165      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
166c
167cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
168      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
169c
170cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
171      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
172c
173cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
174      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
175c
176      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
177cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
178      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
179c
180      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
181cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
182      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
183c
184c200 END
185cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
186      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
187c
188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
189      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
190c
191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
192      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
193c
194cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
195      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
196c
197cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
198      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
199c
200c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
201      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
202cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
203      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
204c
205c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
206      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
207cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
208      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
209c
210c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
211      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
212cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
213      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
214c
215c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
216      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
217cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
218      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
219c
220cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
221cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
222      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
223cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
224      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
225c
226      DO i=1, klon
227       sollwdownclr(i)=-1.*lwdn0(i,1)
228      ENDDO !i=1, klon
229
230cIM 080904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
231cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
232      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
233cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
234      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
235c
236cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
237      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
238c
239c     DO i=1, klon
240cIM 080904     sollwdownclr(i)=lwdn0(i,1)
241c      sollwdownclr(i)=-1.*lwdn0(i,1)
242c     ENDDO !i=1, klon
243c
244cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
245      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
246c
247c
248c effets des aerosols
249c
250      IF (ok_ade.OR.ok_aie) THEN
251      zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
252      CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
253c
254      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
255      CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
256c
257      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
258      CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
259c
260      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
261      CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
262      endif
263
264c
265cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
266      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
267c
268      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
269c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
270cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
271      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
272c
273cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
274      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
275c
276c
277c      DO i = 1, klon
278c         zx_tmp_fi2d(i) = fluxu(i,1)
279c      ENDDO
280c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
281c      CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
282c
283c      DO i = 1, klon
284c         zx_tmp_fi2d(i) = fluxv(i,1)
285c      ENDDO
286c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
287c      CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
288c
289c
290cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
291       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
292c
293cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
294       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
295
296cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqfonte,zx_tmp_2d)
297       CALL histwrite_phy(nid_mth,"fqfonte",itau_w,zxfqfonte)
298
299       DO nsrf = 1, nbsrf
300        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
301cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
302        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
303     $      zx_tmp_fi2d) 
304C     
305        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
306cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
307        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
308     $      zx_tmp_fi2d)
309       ENDDO
310c
311      DO nsrf = 1, nbsrf
312CYYY
313       IF(nsrf.GE.2) THEN
314        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
315cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
316        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
317     &                     zx_tmp_fi2d) 
318c
319        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
320cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
321        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w,
322     $      zx_tmp_fi2d) 
323       ENDIF !nsrf.GE.2
324C
325        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
326cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
327        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
328     $      zx_tmp_fi2d)
329C
330        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
331cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
332        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
333     $      zx_tmp_fi2d)
334C
335        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
336cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
337        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
338     $      zx_tmp_fi2d) 
339C
340        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
341cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
342        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
343     $      zx_tmp_fi2d)
344c
345        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
346cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
347        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
348     $      zx_tmp_fi2d)
349C
350        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
351cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
352        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
353     $      zx_tmp_fi2d)
354C
355        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
356cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
357        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
358     $      zx_tmp_fi2d)
359C
360       IF(1.EQ.0) THEN
361C
362        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
363cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
364        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
365     $      zx_tmp_fi2d)
366C
367        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
368     $                        *pctsrf(1:klon,nsrf)
369cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
370        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
371     $      zx_tmp_fi2d)
372C
373       ENDIF !(1.EQ.0) THEN
374C
375      END DO 
376c
377cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
378      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
379c
380cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
381      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
382c
383cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
384      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
385c
386cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
387      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
388c
389cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
390      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
391c
392cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
393      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt)
394c
395cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
396      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
397c
398      zx_tmp_fi2d(1:klon) = flwp(1:klon)
399cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
400      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
401c
402      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
403cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
404      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
405c
406cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
407      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
408c
409cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
410      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
411c
412cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
413      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
414c
415cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
416      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
417cKE43
418      IF (iflag_con.GE.3) THEN ! sb
419c
420cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
421      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
422c
423cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
424      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
425c
426cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
427      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
428c
429cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
430      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
431c
432c34EK
433c
434cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
435      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
436c
437      ENDIF !iflag_con.GE.3
438c
439c
440cIM cf. AM 081204 BEG
441cHBTM2
442c
443c
444cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
445      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
446c
447cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
448      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
449c
450cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
451      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
452c
453cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
454      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
455c
456cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
457      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
458c
459cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
460      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
461c
462cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
463      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
464c
465cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
466      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
467c
468cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
469      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
470c
471cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
472      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
473c
474cIM cf. AM 081204 END
475c
476c  Champs interpolles sur des niveaux de pression
477      DO k=1, nlevSTD
478c
479c        bb=clevSTD(k)
480c
481c        IF(k.GE.2) THEN
482         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
483         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
484c         aa=clevSTD(k)
485c         bb=aa(1:lnblnk1(aa))
486c        ENDIF
487c
488      IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
489     $   bb2.EQ."500".OR.bb2.EQ."200") THEN
490c
491cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
492        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
493c
494cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
495        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
496c
497cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
498      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
499c
500cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
501      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
502c
503      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
504      ENDDO !nlevSTD
505c
506      DO i=1, klon
507       IF (pctsrf(i,is_oce).GT.epsfra) THEN
508        zx_tmp_fi2d(i) = fluxo(i)
509       ELSE
510        zx_tmp_fi2d(i) = 0.
511       ENDIF
512      ENDDO
513c
514cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
515      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
516c
517      DO i=1, klon
518       IF (pctsrf(i,is_sic).GT.epsfra) THEN
519        zx_tmp_fi2d(i) = fluxg(i)
520       ELSE
521        zx_tmp_fi2d(i) = 0.
522       ENDIF
523      ENDDO
524c
525cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
526      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
527c
528      DO i=1, klon
529       IF (pctsrf(i,is_oce).GT.epsfra.OR.
530     .     pctsrf(i,is_sic).GT.epsfra) THEN
531        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
532     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
533     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
534       ELSE
535        zx_tmp_fi2d(i) = 273.15
536       ENDIF
537      ENDDO
538c
539cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
540      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
541c
542      IF (OCEAN.EQ.'force ') THEN
543c
544      DO i=1, klon
545      IF((pctsrf(i,is_oce).GT.epsfra).OR.
546     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
547       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
548     .                  fluxg(i)*pctsrf(i,is_sic)
549      ELSE
550       zx_tmp_fi2d(i) = 1.E+20
551      ENDIF
552      ENDDO
553cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
554      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,zx_tmp_fi2d)
555c
556      ELSE IF (OCEAN.EQ.'slab  ') THEN
557c
558      DO i=1, klon
559      IF((pctsrf(i,is_oce).GT.epsfra).OR.
560     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
561       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
562     .                  fluxg(i)*pctsrf(i,is_sic)
563      ELSE
564       zx_tmp_fi2d(i) = 1.E+20
565      ENDIF
566      ENDDO
567cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
568      CALL histwrite_phy(nid_mth,"slab_bils",itau_w,zx_tmp_fi2d)
569c
570      DO i=1, klon
571       IF(pctsrf(i,is_oce).GT.epsfra.OR.
572     $    pctsrf(i,is_sic).GT.epsfra) THEN
573        zx_tmp_fi2d(i)=tslab(i)
574       ELSE
575        zx_tmp_fi2d(i) = 1.E+20
576       ENDIF
577      ENDDO !i=1, klon
578cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
579      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
580c
581cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
582      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
583c
584cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
585      CALL histwrite_phy(nid_mth,"siceh",itau_w, seaice/1000.)
586c
587      ENDIF !OCEAN.EQ.force/slab
588c
589      ENDIF !lev_histmth.GE.1
590c
591      IF(lev_histmth.GE.2) THEN
592c
593c Champs 3D:
594C
595c
596cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
597      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
598c
599cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
600      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
601c
602cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
603      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
604c
605c#ifdef histmthNMC
606c#include "write_histmthNMC.h"
607c#endif
608c
609cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
610      CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap))
611c
612c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
613c    .                         (1-qx(1:klon,1:klev,ivap))
614c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
615c     CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
616c    .                                   iim*jjmp1*klev,ndex3d)
617c
618cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
619      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
620c
621cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
622      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
623c
624cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
625      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
626c
627cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
628      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
629c
630cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
631      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
632c
633cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
634      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra)
635c
636cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
637      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
638c
639cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
640      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
641
642      DO k=1, klev
643       DO i=1, klon
644cIM 151004      zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
645c    $                    /(paprs(i,k)-paprs(i,k+1))
646c    $                    *(paprs(i,1)/101325.0)
647c        zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968
648         zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
649     $                    /(paprs(i,k)-paprs(i,k+1))
650     $                    *(paprs(i,1)/101325.0)
651       ENDDO !i
652      ENDDO !k
653c
654c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
655c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
656cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
657      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
658c
659cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
660      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
661c
662cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap), 
663cym     .                 zx_tmp_3d)
664      CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap))
665cIM 22.03.04 END
666c
667      ENDIF !lev_histmth.GE.2
668c
669      IF(lev_histmth.GE.3) THEN
670c
671       DO nsrf=1, nbsrf
672c
673        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
674cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
675        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
676     $      zx_tmp_fi2d)
677C
678        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
679cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
680        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
681     $      zx_tmp_fi2d)
682c
683      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
684cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
685      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
686     $    ,zx_tmp_fi2d)
687c
688      ENDDO !nsrf=1, nbsrf
689c
690cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
691      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
692c
693cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
694      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
695c
696      ENDIF !lev_histmth.GE.3
697c
698c
699      IF(lev_histmth.GE.4) THEN
700c
701cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
702      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
703
704      IF (iflag_con.GE.3) THEN
705c
706cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
707      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
708c
709cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
710      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
711c
712cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
713      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
714c
715cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
716      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
717c
718      ENDIF !iflag_con.GE.3
719c
720cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
721      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
722c
723cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
724      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
725c
726cIM: 101003 : K/30min ==> K/s
727      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
728cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
729      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
730c
731      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
732cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
733      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
734c
735cIM: 101003 : K/30min ==> K/s
736      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
737cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
738      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
739cIM: 071003
740cIM: 101003 : K/30min ==> K/s
741      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
742     .                             d_t_eva(1:klon,1:klev))/pdtphys
743cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
744      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
745c
746      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
747cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
748      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
749c
750cIM: 101003 : K/30min ==> K/s
751      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
752cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
753      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
754c
755      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
756cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
757      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
758c
759cIM: 101003 : K/30min ==> K/s
760      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
761cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
762      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
763c
764      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
765cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
766      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
767c
768      zpt_conv = 0.
769      where (ptconv) zpt_conv = 1.
770cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
771      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
772c
773cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
774      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
775c
776cIM: 101003 : K/30min ==> K/s
777      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
778cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
779      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
780c
781      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
782cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
783      CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
784c
785cIM: 101003 : K/day ==> K/s
786      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
787cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
788      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
789c
790cIM: 101003 : K/day ==> K/s     
791      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
792cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
793      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
794c
795cIM: 101003 : K/day ==> K/s     
796      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
797cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
798      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
799c
800cIM: 101003 : K/day ==> K/s     
801      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
802cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
803      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
804c
805cIM: 101003 : deja en K/s     
806      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
807cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
808      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
809c
810      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
811cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
812      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
813c
814      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
815cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
816      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
817c
818      IF (ok_orodr) THEN
819      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
820cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
821      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
822c
823      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
824cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
825      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
826c
827      ENDIF
828C
829      IF (ok_orolf) THEN
830      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
831cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
832      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
833c
834      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
835cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
836      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
837      ENDIF
838C
839
840#if defined(INCA_AER) && defined(CPP_COUPLE)
841      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
842      CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
843      CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
844      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
845      CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
846      CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
847      CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
848      CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
849      CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
850      CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
851      CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
852      CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
853      CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
854      CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
855#endif
856#ifndef INCA
857      IF (nqmax.GE.3) THEN
858      DO iq=3,nqmax
859cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
860         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
861      ENDDO
862      ENDIF
863#endif
864      ENDIF !lev_histmth.GE.4
865c
866      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
867c
868      IF(lev_histmth.GE.1) THEN
869c
870c Champs 2D:
871c
872cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
873      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG)
874C
875cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
876cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
877      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
878c
879      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
880cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
881      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
882     $      zx_tmp_fi2d)
883c
884      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
885cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
886      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
887     $      zx_tmp_fi2d)
888c
889      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100.
890cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
891      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
892     $      zx_tmp_fi2d)
893c
894      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)
895cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
896      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
897     $      zx_tmp_fi2d)
898c
899      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100.
900cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
901      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
902     $      zx_tmp_fi2d)
903c
904      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)
905cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
906      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
907     $      zx_tmp_fi2d)
908c
909cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
910      CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat)
911c
912      zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon)
913cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
914      CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d)
915c
916c     IF(itap.EQ.ecrit_mth) THEN
917c       DO i=1, klon
918c         IF(nday_clr(i).GT.0.) THEN
919c          tops01(i)=tops01(i)/nday_clr(i)
920c         ELSE
921c          tops01(i)=1.e+20
922c         ENDIF
923c       ENDDO
924c       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d)
925c       CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d,
926c    .       iim*jjmp1,ndex2d)
927c     ENDIF !itap.EQ.ecrit_mth
928c
929c220404 BEG
930c
931c temperature tendency due to total diabatic heating
932       DO l=1, klev
933       DO i=1, klon
934       d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime
935       ENDDO !i
936       ENDDO !l
937c
938cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
939      CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t)
940c
941c temperature tendency due to SW radiation
942       DO l=1, klev
943       DO i=1, klon
944       zx_tmp_fi3d(i,l)=heat(i,l)/RDAY
945       ENDDO !i
946       ENDDO !l
947c
948cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
949      CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d)
950c
951c temperature tendency due to LW radiation
952       DO l=1, klev
953       DO i=1, klon
954        zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY
955       ENDDO !i
956       ENDDO !l
957c
958cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
959      CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
960c
961c  temperature tendency due to moist convective processes
962       DO l=1, klev
963       DO i=1, klon
964        zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
965       ENDDO !i
966       ENDDO !l
967c
968cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
969      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
970c
971c temperature tendency due to dry convective processes
972       DO l=1, klev
973       DO i=1, klon
974        zx_tmp_fi3d(i,l)=d_t_ajs(i,l)/pdtphys
975       ENDDO !i
976       ENDDO !l
977c
978cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
979      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
980c
981c  temperature tendency due to large scale precipitation
982       DO l=1, klev
983       DO i=1, klon
984       zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys
985       ENDDO !i
986       ENDDO !l
987c
988cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
989      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
990c
991c  total moisture tendency due to diabatic processes
992       DO l=1, klev
993       DO i=1, klon
994        d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime
995        zx_tmp_fi3d(i,l)=d_qx(i,l,ivap)
996       ENDDO !i
997       ENDDO !l
998c
999cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
1000      CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d)
1001c
1002cBEG 210404 BEG
1003c
1004      IF(ok_orodr.AND.ok_orolf) THEN
1005      DO k = 1, klev
1006      DO i = 1, klon
1007        d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys
1008        d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys
1009      ENDDO
1010      ENDDO
1011      ENDIF !(ok_orodr.AND.ok_orolf) THEN
1012c
1013cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
1014      CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli)
1015c
1016cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
1017      CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli)
1018c
1019      zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1020cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1021      CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d)
1022c
1023      zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1024cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1025      CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d)
1026c
1027cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
1028      CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u)
1029c
1030cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
1031      CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v)
1032
1033cEND 210404 END
1034c
1035c cldtrue
1036c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d)
1037c     CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
1038c    $       iim*jjmp1*klev,ndex3d)
1039c
1040c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d)
1041c     CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
1042c    $       iim*jjmp1*klev,ndex3d)
1043c
1044c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d)
1045c     CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
1046c    $       iim*jjmp1*klev,ndex3d)
1047c
1048      DO i=1, klon
1049cIM 151004    zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1)
1050       zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i)
1051      ENDDO
1052cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1053      CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d)
1054c
1055      DO i=1, klon
1056cIM 151004    zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1)
1057       zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i)
1058      ENDDO
1059cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1060      CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d)
1061c
1062c220404 END
1063cAMIP2 END
1064c
1065cym     CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
1066      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
1067
1068cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
1069      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
1070c
1071cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
1072      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
1073c
1074c ENSEMBLES BEG
1075cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
1076      CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
1077c
1078cym     CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
1079      CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
1080c
1081c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
1082c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
1083c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
1084c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
1085c    .               iim*jjmp1,ndex2d)
1086c
1087      DO i=1, klon
1088       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
1089      ENDDO
1090cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1091      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
1092c
1093      DO i = 1, klon
1094         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
1095      ENDDO
1096cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1097      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
1098c
1099c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
1100c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
1101c    .                                   iim*jjmp1*klev,ndex3d)
1102c
1103c ENSEMBLES END
1104c
1105cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
1106      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
1107c
1108cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
1109      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
1110C
1111cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
1112      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
1113C
1114      DO i = 1, klon
1115         zx_tmp_fi2d(i) = paprs(i,1)
1116      ENDDO
1117cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1118      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
1119c
1120cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
1121      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
1122c
1123      if (.not. ok_veget) then
1124cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
1125          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
1126      endif
1127c
1128      DO i = 1, klon
1129         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
1130      ENDDO
1131cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1132      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
1133c
1134cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
1135      CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
1136c
1137      DO i = 1, klon
1138         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
1139      ENDDO
1140cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1141      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
1142c
1143      DO i = 1, klon
1144         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
1145      ENDDO
1146cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1147      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
1148c
1149cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
1150      CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall)
1151c
1152c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
1153c     CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
1154c    .               ndex2d)
1155c
1156cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
1157      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
1158c
1159cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
1160      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
1161c
1162cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
1163      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
1164c
1165cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
1166      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
1167c
1168cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
1169      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
1170c
1171      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
1172cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1173      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
1174c
1175c
1176      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
1177cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1178      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
1179c
1180      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
1181cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1182      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
1183c
1184      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
1185cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1186      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
1187c
1188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
1189      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
1190c
1191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
1192      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
1193c
1194cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
1195      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
1196c
1197cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
1198      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
1199c
1200cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
1201      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
1202c
1203cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
1204      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
1205c
1206      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
1207cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1208      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
1209c
1210      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
1211cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1212      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
1213c
1214c200 END
1215cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
1216      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
1217c
1218cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
1219      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
1220c
1221cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
1222      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
1223c
1224cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
1225      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
1226c
1227cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
1228      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
1229c
1230c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
1231      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
1232cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1233      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
1234c
1235c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
1236      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
1237cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1238      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
1239c
1240c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
1241      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
1242cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1243      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
1244c
1245c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
1246      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
1247cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1248      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
1249c
1250cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
1251cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
1252      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
1253cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
1254      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
1255c
1256      DO i=1, klon
1257       sollwdownclr(i)=-1.*lwdn0(i,1)
1258      ENDDO !i=1, klon
1259c
1260cIM 080904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
1261cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
1262      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
1263cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
1264      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
1265c
1266cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
1267      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
1268c
1269c     DO i=1, klon
1270cIM 08004    sollwdownclr(i)=lwdn0(i,1)
1271c      sollwdownclr(i)=-1.*lwdn0(i,1)
1272c     ENDDO !i=1, klon
1273c
1274cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
1275      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
1276c
1277cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
1278      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
1279c
1280      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
1281c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
1282cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1283      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
1284c
1285cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
1286      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
1287c
1288c Effets des aerosols
1289c
1290      IF (ok_ade.OR.ok_aie) THEN
1291        zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
1292        CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
1293c
1294        zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
1295        CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
1296c
1297        zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
1298        CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
1299c
1300        zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
1301        CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
1302      ENDIF
1303c
1304c
1305c      DO i = 1, klon
1306c         zx_tmp_fi2d(i) = fluxu(i,1)
1307c      ENDDO
1308c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1309c      CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1310c
1311c      DO i = 1, klon
1312c         zx_tmp_fi2d(i) = fluxv(i,1)
1313c      ENDDO
1314c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1315c      CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1316c
1317c
1318cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
1319       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
1320c
1321cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
1322       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
1323c
1324       DO nsrf = 1, nbsrf
1325        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
1326cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1327        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
1328     $      zx_tmp_fi2d) 
1329C     
1330        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
1331cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1332        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
1333     $      zx_tmp_fi2d)
1334       ENDDO
1335c
1336        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
1337cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1338        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
1339     $      zx_tmp_fi2d) 
1340c
1341        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
1342cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1343        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
1344     $      zx_tmp_fi2d) 
1345c
1346      DO nsrf = 1, nbsrf
1347CYYY
1348C
1349        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
1350cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1351        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
1352     $      zx_tmp_fi2d)
1353C
1354        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
1355cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1356        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
1357     $      zx_tmp_fi2d)
1358C
1359        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
1360cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1361        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
1362     $      zx_tmp_fi2d) 
1363C
1364        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
1365cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1366        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
1367     $      zx_tmp_fi2d)
1368c
1369        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
1370cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1371        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
1372     $      zx_tmp_fi2d)
1373C
1374        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
1375cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1376        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
1377     $      zx_tmp_fi2d)
1378C
1379        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
1380cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1381        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
1382     $      zx_tmp_fi2d)
1383C
1384       IF(1.EQ.0) THEN
1385C
1386        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
1387cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1388        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
1389     $      zx_tmp_fi2d)
1390C
1391        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
1392     $                        *pctsrf(1:klon,nsrf)
1393cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1394        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
1395     $      zx_tmp_fi2d)
1396C
1397       ENDIF !(1.EQ.0) THEN
1398C
1399      END DO 
1400c
1401cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
1402      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
1403c
1404cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
1405      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
1406c
1407cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
1408      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
1409c
1410cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
1411      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
1412c
1413cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
1414      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
1415c
1416cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
1417      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
1418c
1419cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
1420      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
1421c
1422      zx_tmp_fi2d(1:klon) = flwp(1:klon)
1423cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1424      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
1425c
1426      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
1427cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1428      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
1429c
1430cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
1431      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
1432c
1433cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
1434      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
1435c
1436cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
1437      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
1438c
1439cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
1440      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
1441cKE43
1442      IF (iflag_con.GE.3) THEN ! sb
1443c
1444cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
1445      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
1446c
1447cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
1448      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
1449c
1450cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
1451      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
1452c
1453cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
1454      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
1455c
1456c34EK
1457c
1458cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
1459      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
1460c
1461      ENDIF !iflag_con.GE.3
1462c
1463cIM cf. AM 081204 BEG
1464cHBTM2
1465c
1466c
1467cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
1468      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
1469c
1470cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
1471      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
1472c
1473cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
1474      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
1475c
1476cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
1477      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
1478c
1479cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
1480      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
1481c
1482cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
1483      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
1484c
1485cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
1486      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
1487c
1488cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
1489      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
1490c
1491cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
1492      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
1493c
1494cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
1495      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
1496c
1497cIM cf. AM 081204 END
1498c
1499c  Champs interpolles sur des niveaux de pression
1500c     DO k=1, nlevENS
1501      DO k=1, nlevSTD
1502c
1503c        bb=clevSTD(k)
1504c
1505c        IF(k.GE.2) THEN
1506         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1507         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1508c         aa=clevSTD(k)
1509c         bb=aa(1:lnblnk1(aa))
1510c        ENDIF
1511c
1512       IF(bb2.EQ."850".OR.bb2.EQ."700"
1513     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1514c
1515cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
1516        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
1517c
1518cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
1519        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
1520c
1521cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
1522      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
1523c
1524cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
1525      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
1526c
1527      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1528      ENDDO !nlevSTD
1529c
1530      DO i=1, klon
1531       IF (pctsrf(i,is_oce).GT.epsfra) THEN
1532        zx_tmp_fi2d(i) = fluxo(i)
1533       ELSE
1534        zx_tmp_fi2d(i) = 0.
1535       ENDIF
1536      ENDDO
1537c
1538cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1539      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
1540c
1541      DO i=1, klon
1542       IF (pctsrf(i,is_sic).GT.epsfra) THEN
1543        zx_tmp_fi2d(i) = fluxg(i)
1544       ELSE
1545        zx_tmp_fi2d(i) = 0.
1546       ENDIF
1547      ENDDO
1548c
1549cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1550      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
1551c
1552      DO i=1, klon
1553       IF (pctsrf(i,is_oce).GT.epsfra.OR.
1554     .     pctsrf(i,is_sic).GT.epsfra) THEN
1555        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
1556     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
1557     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
1558       ELSE
1559        zx_tmp_fi2d(i) = 273.15
1560       ENDIF
1561      ENDDO
1562c
1563cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1564      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
1565c
1566      IF (OCEAN.EQ.'force ') THEN
1567c
1568      DO i=1, klon
1569      IF((pctsrf(i,is_oce).GT.epsfra).OR.
1570     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
1571       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
1572     .                  fluxg(i)*pctsrf(i,is_sic)
1573      ELSE
1574       zx_tmp_fi2d(i) = 1.E+20
1575      ENDIF
1576      ENDDO
1577cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1578      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
1579     .     zx_tmp_fi2d)
1580c
1581      ELSE IF (OCEAN.EQ.'slab  ') THEN
1582c
1583c
1584      DO i=1, klon
1585      IF((pctsrf(i,is_oce).GT.epsfra).OR.
1586     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
1587       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
1588     .                  fluxg(i)*pctsrf(i,is_sic)
1589      ELSE
1590       zx_tmp_fi2d(i) = 1.E+20
1591      ENDIF
1592      ENDDO
1593cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1594      CALL histwrite_phy(nid_mth,"slab_bils",itau_w,
1595     .     zx_tmp_fi2d)
1596c
1597      DO i=1, klon
1598       IF(pctsrf(i,is_oce).GT.epsfra.OR.
1599     $    pctsrf(i,is_sic).GT.epsfra) THEN
1600        zx_tmp_fi2d(i)=tslab(i)
1601       ELSE
1602        zx_tmp_fi2d(i) = 1.E+20
1603       ENDIF
1604      ENDDO !i=1, klon
1605cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1606      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
1607c
1608cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
1609      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
1610c
1611cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
1612      CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
1613c
1614      ENDIF !OCEAN.EQ.force/slab
1615c
1616      ENDIF !lev_histmth.GE.1
1617c
1618      IF(lev_histmth.GE.2) THEN
1619c
1620c Champs 3D:
1621C
1622c
1623cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
1624      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
1625c
1626cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
1627      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
1628c
1629cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
1630      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
1631c
1632c#ifdef histmthNMC
1633c#include "write_histmthNMC.h"
1634c#endif
1635c
1636cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
1637      CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
1638c
1639c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
1640c    .                         (1-qx(1:klon,1:klev,ivap))
1641c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1642c     CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
1643c    .                                   iim*jjmp1*klev,ndex3d)
1644c
1645cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
1646      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
1647c
1648cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
1649      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
1650c
1651cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
1652      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
1653c
1654cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
1655      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
1656c
1657cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
1658      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
1659c
1660cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
1661      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
1662c
1663cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
1664      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
1665c
1666cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
1667      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
1668c
1669      DO k=1, klev
1670       DO i=1, klon
1671cIM 151004        zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
1672c    $                    /(paprs(i,k)-paprs(i,k+1))
1673c    $                    *(paprs(i,1)/101325.0)
1674c        zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968
1675         zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
1676     $                    /(paprs(i,k)-paprs(i,k+1))
1677     $                    *(paprs(i,1)/101325.0)
1678       ENDDO !i
1679      ENDDO !k
1680c
1681c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
1682c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
1683cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1684      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
1685c
1686cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
1687      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
1688c
1689cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
1690      CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
1691c
1692      ENDIF !lev_histmth.GE.2
1693c
1694      IF(lev_histmth.GE.3) THEN
1695c
1696       DO nsrf=1, nbsrf
1697c
1698        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
1699cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1700        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
1701     $      zx_tmp_fi2d)
1702C
1703        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1704cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1705        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
1706     $      zx_tmp_fi2d)
1707c
1708      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1709cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
1710      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
1711     $    ,agesno)
1712c
1713      ENDDO !nsrf=1, nbsrf
1714c
1715cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
1716      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
1717c
1718cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
1719      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
1720c
1721      ENDIF !lev_histmth.GE.3
1722c
1723c
1724      IF(lev_histmth.GE.4) THEN
1725c
1726cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
1727      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
1728c
1729      IF (iflag_con.GE.3) THEN
1730c
1731cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
1732      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
1733c
1734cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
1735      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
1736c
1737cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
1738      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
1739c
1740cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
1741      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
1742c
1743      ENDIF !iflag_con.GE.3
1744c
1745cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
1746      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
1747c
1748cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
1749      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
1750c
1751cIM: 101003 : K/30min ==> K/s
1752c
1753      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1754cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1755      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
1756c
1757cIM: 101003 : K/30min ==> K/s
1758      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1759     .                             d_t_eva(1:klon,1:klev))/pdtphys
1760cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1761      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
1762c
1763      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1764cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1765      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
1766c
1767cIM: 101003 : K/30min ==> K/s
1768      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1769cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1770      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
1771c
1772      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1773cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1774      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
1775c
1776cIM: 101003 : K/30min ==> K/s
1777      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1778cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1779      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
1780c
1781      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1782cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1783      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
1784c
1785      zpt_conv = 0.
1786      where (ptconv) zpt_conv = 1.
1787cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
1788      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
1789c
1790cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
1791      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
1792c
1793      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
1794cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1795      CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
1796c
1797cIM: 101003 : K/day ==> K/s
1798      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1799cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1800      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
1801c
1802cIM: 101003 : K/day ==> K/s     
1803      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1804cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1805      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
1806c
1807cIM: 101003 : K/day ==> K/s     
1808      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1809cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1810      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
1811c
1812cIM: 101003 : K/day ==> K/s     
1813      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1814cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1815      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
1816c
1817cIM: 101003 : deja en K/s     
1818      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1819cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1820      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
1821c
1822      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1823cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1824      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
1825c
1826      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1827cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1828      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
1829c
1830      IF (ok_orodr) THEN
1831      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1832cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1833      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
1834c
1835      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1836cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1837      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
1838c
1839      ENDIF
1840C
1841      IF (ok_orolf) THEN
1842      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1843cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1844      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
1845c
1846      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1847cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1848      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
1849      ENDIF
1850C
1851#if defined(INCA_AER) && defined(CPP_COUPLE)
1852      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
1853      CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
1854      CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
1855      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
1856      CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
1857      CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
1858      CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
1859      CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
1860      CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
1861      CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
1862      CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
1863      CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
1864      CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
1865      CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
1866#endif
1867#ifndef INCA
1868      IF (nqmax.GE.3) THEN
1869      DO iq=3,nqmax
1870cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
1871         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
1872      ENDDO
1873      ENDIF
1874c
1875#endif
1876      ENDIF !lev_histmth.GE.4
1877c
1878      ENDIF !(type_run.EQ."AMIP")
1879c
1880      if (ok_sync) then
1881c$OMP MASTER
1882        call histsync(nid_mth)
1883c$OMP END MASTER
1884      endif
1885c
1886      ENDIF !ok_mensuel
Note: See TracBrowser for help on using the repository browser.