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

Last change on this file since 896 was 896, checked in by Laurent Fairhead, 16 years ago

Rajout de sorties concernant les aerosols qui avaient disparus lors d'un merge
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 69.2 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_2d)
253c
254      zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
255      CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_2d)
256c
257      zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
258      CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_2d)
259c
260      zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
261      CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_2d)
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.GT.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
1289c      DO i = 1, klon
1290c         zx_tmp_fi2d(i) = fluxu(i,1)
1291c      ENDDO
1292c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1293c      CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1294c
1295c      DO i = 1, klon
1296c         zx_tmp_fi2d(i) = fluxv(i,1)
1297c      ENDDO
1298c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1299c      CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1300c
1301c
1302cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
1303       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
1304c
1305cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
1306       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
1307c
1308       DO nsrf = 1, nbsrf
1309        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
1310cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1311        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
1312     $      zx_tmp_fi2d) 
1313C     
1314        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
1315cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1316        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
1317     $      zx_tmp_fi2d)
1318       ENDDO
1319c
1320        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
1321cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1322        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
1323     $      zx_tmp_fi2d) 
1324c
1325        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
1326cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1327        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
1328     $      zx_tmp_fi2d) 
1329c
1330      DO nsrf = 1, nbsrf
1331CYYY
1332C
1333        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
1334cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1335        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
1336     $      zx_tmp_fi2d)
1337C
1338        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
1339cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1340        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
1341     $      zx_tmp_fi2d)
1342C
1343        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
1344cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1345        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
1346     $      zx_tmp_fi2d) 
1347C
1348        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
1349cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1350        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
1351     $      zx_tmp_fi2d)
1352c
1353        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
1354cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1355        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
1356     $      zx_tmp_fi2d)
1357C
1358        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
1359cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1360        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
1361     $      zx_tmp_fi2d)
1362C
1363        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
1364cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1365        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
1366     $      zx_tmp_fi2d)
1367C
1368       IF(1.EQ.0) THEN
1369C
1370        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
1371cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1372        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
1373     $      zx_tmp_fi2d)
1374C
1375        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
1376     $                        *pctsrf(1:klon,nsrf)
1377cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1378        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
1379     $      zx_tmp_fi2d)
1380C
1381       ENDIF !(1.EQ.0) THEN
1382C
1383      END DO 
1384c
1385cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
1386      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
1387c
1388cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
1389      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
1390c
1391cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
1392      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
1393c
1394cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
1395      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
1396c
1397cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
1398      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
1399c
1400cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
1401      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
1402c
1403cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
1404      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
1405c
1406      zx_tmp_fi2d(1:klon) = flwp(1:klon)
1407cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1408      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
1409c
1410      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
1411cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1412      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
1413c
1414cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
1415      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
1416c
1417cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
1418      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
1419c
1420cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
1421      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
1422c
1423cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
1424      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
1425cKE43
1426      IF (iflag_con.GE.3) THEN ! sb
1427c
1428cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
1429      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
1430c
1431cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
1432      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
1433c
1434cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
1435      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
1436c
1437cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
1438      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
1439c
1440c34EK
1441c
1442cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
1443      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
1444c
1445      ENDIF !iflag_con.GE.3
1446c
1447cIM cf. AM 081204 BEG
1448cHBTM2
1449c
1450c
1451cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
1452      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
1453c
1454cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
1455      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
1456c
1457cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
1458      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
1459c
1460cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
1461      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
1462c
1463cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
1464      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
1465c
1466cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
1467      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
1468c
1469cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
1470      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
1471c
1472cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
1473      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
1474c
1475cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
1476      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
1477c
1478cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
1479      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
1480c
1481cIM cf. AM 081204 END
1482c
1483c  Champs interpolles sur des niveaux de pression
1484c     DO k=1, nlevENS
1485      DO k=1, nlevSTD
1486c
1487c        bb=clevSTD(k)
1488c
1489c        IF(k.GE.2) THEN
1490         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1491         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1492c         aa=clevSTD(k)
1493c         bb=aa(1:lnblnk1(aa))
1494c        ENDIF
1495c
1496       IF(bb2.EQ."850".OR.bb2.EQ."700"
1497     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1498c
1499cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
1500        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
1501c
1502cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
1503        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
1504c
1505cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
1506      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
1507c
1508cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
1509      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
1510c
1511      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1512      ENDDO !nlevSTD
1513c
1514      DO i=1, klon
1515       IF (pctsrf(i,is_oce).GT.epsfra) THEN
1516        zx_tmp_fi2d(i) = fluxo(i)
1517       ELSE
1518        zx_tmp_fi2d(i) = 0.
1519       ENDIF
1520      ENDDO
1521c
1522cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1523      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
1524c
1525      DO i=1, klon
1526       IF (pctsrf(i,is_sic).GT.epsfra) THEN
1527        zx_tmp_fi2d(i) = fluxg(i)
1528       ELSE
1529        zx_tmp_fi2d(i) = 0.
1530       ENDIF
1531      ENDDO
1532c
1533cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1534      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
1535c
1536      DO i=1, klon
1537       IF (pctsrf(i,is_oce).GT.epsfra.OR.
1538     .     pctsrf(i,is_sic).GT.epsfra) THEN
1539        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
1540     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
1541     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
1542       ELSE
1543        zx_tmp_fi2d(i) = 273.15
1544       ENDIF
1545      ENDDO
1546c
1547cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1548      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
1549c
1550      IF (OCEAN.EQ.'force ') THEN
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) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
1556     .                  fluxg(i)*pctsrf(i,is_sic)
1557      ELSE
1558       zx_tmp_fi2d(i) = 1.E+20
1559      ENDIF
1560      ENDDO
1561cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1562      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
1563     .     zx_tmp_fi2d)
1564c
1565      ELSE IF (OCEAN.EQ.'slab  ') THEN
1566c
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,"slab_bils",itau_w,
1579     .     zx_tmp_fi2d)
1580c
1581      DO i=1, klon
1582       IF(pctsrf(i,is_oce).GT.epsfra.OR.
1583     $    pctsrf(i,is_sic).GT.epsfra) THEN
1584        zx_tmp_fi2d(i)=tslab(i)
1585       ELSE
1586        zx_tmp_fi2d(i) = 1.E+20
1587       ENDIF
1588      ENDDO !i=1, klon
1589cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1590      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
1591c
1592cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
1593      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
1594c
1595cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
1596      CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
1597c
1598      ENDIF !OCEAN.EQ.force/slab
1599c
1600      ENDIF !lev_histmth.GE.1
1601c
1602      IF(lev_histmth.GE.2) THEN
1603c
1604c Champs 3D:
1605C
1606c
1607cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
1608      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
1609c
1610cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
1611      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
1612c
1613cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
1614      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
1615c
1616c#ifdef histmthNMC
1617c#include "write_histmthNMC.h"
1618c#endif
1619c
1620cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
1621      CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
1622c
1623c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
1624c    .                         (1-qx(1:klon,1:klev,ivap))
1625c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1626c     CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
1627c    .                                   iim*jjmp1*klev,ndex3d)
1628c
1629cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
1630      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
1631c
1632cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
1633      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
1634c
1635cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
1636      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
1637c
1638cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
1639      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
1640c
1641cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
1642      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
1643c
1644cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
1645      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
1646c
1647cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
1648      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
1649c
1650cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
1651      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
1652c
1653      DO k=1, klev
1654       DO i=1, klon
1655cIM 151004        zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
1656c    $                    /(paprs(i,k)-paprs(i,k+1))
1657c    $                    *(paprs(i,1)/101325.0)
1658c        zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968
1659         zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
1660     $                    /(paprs(i,k)-paprs(i,k+1))
1661     $                    *(paprs(i,1)/101325.0)
1662       ENDDO !i
1663      ENDDO !k
1664c
1665c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
1666c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
1667cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1668      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
1669c
1670cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
1671      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
1672c
1673cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
1674      CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
1675c
1676      ENDIF !lev_histmth.GE.2
1677c
1678      IF(lev_histmth.GE.3) THEN
1679c
1680       DO nsrf=1, nbsrf
1681c
1682        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
1683cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1684        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
1685     $      zx_tmp_fi2d)
1686C
1687        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1688cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1689        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
1690     $      zx_tmp_fi2d)
1691c
1692      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1693cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
1694      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
1695     $    ,agesno)
1696c
1697      ENDDO !nsrf=1, nbsrf
1698c
1699cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
1700      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
1701c
1702cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
1703      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
1704c
1705      ENDIF !lev_histmth.GE.3
1706c
1707c
1708      IF(lev_histmth.GE.4) THEN
1709c
1710cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
1711      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
1712c
1713      IF (iflag_con.GE.3) THEN
1714c
1715cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
1716      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
1717c
1718cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
1719      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
1720c
1721cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
1722      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
1723c
1724cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
1725      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
1726c
1727      ENDIF !iflag_con.GE.3
1728c
1729cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
1730      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
1731c
1732cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
1733      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
1734c
1735cIM: 101003 : K/30min ==> K/s
1736c
1737      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1738cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1739      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
1740c
1741cIM: 101003 : K/30min ==> K/s
1742      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1743     .                             d_t_eva(1:klon,1:klev))/pdtphys
1744cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1745      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
1746c
1747      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1748cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1749      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
1750c
1751cIM: 101003 : K/30min ==> K/s
1752      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1753cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1754      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
1755c
1756      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1757cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1758      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
1759c
1760cIM: 101003 : K/30min ==> K/s
1761      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1762cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1763      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
1764c
1765      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1766cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1767      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
1768c
1769      zpt_conv = 0.
1770      where (ptconv) zpt_conv = 1.
1771cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
1772      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
1773c
1774cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
1775      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
1776c
1777      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(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,"dqajs",itau_w,zx_tmp_fi3d)
1780c
1781cIM: 101003 : K/day ==> K/s
1782      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1783cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1784      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
1785c
1786cIM: 101003 : K/day ==> K/s     
1787      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1788cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1789      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
1790c
1791cIM: 101003 : K/day ==> K/s     
1792      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1793cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1794      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
1795c
1796cIM: 101003 : K/day ==> K/s     
1797      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1798cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1799      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
1800c
1801cIM: 101003 : deja en K/s     
1802      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1803cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1804      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
1805c
1806      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1807cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1808      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
1809c
1810      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1811cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1812      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
1813c
1814      IF (ok_orodr) THEN
1815      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1816cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1817      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
1818c
1819      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1820cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1821      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
1822c
1823      ENDIF
1824C
1825      IF (ok_orolf) THEN
1826      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(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,"dulif",itau_w,zx_tmp_fi3d)
1829c
1830      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1831cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1832      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
1833      ENDIF
1834C
1835#if defined(INCA_AER) && defined(CPP_COUPLE)
1836      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
1837      CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
1838      CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
1839      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
1840      CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
1841      CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
1842      CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
1843      CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
1844      CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
1845      CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
1846      CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
1847      CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
1848      CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
1849      CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
1850#endif
1851#ifndef INCA
1852      IF (nqmax.GE.3) THEN
1853      DO iq=3,nqmax
1854cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
1855         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
1856      ENDDO
1857      ENDIF
1858c
1859#endif
1860      ENDIF !lev_histmth.GE.4
1861c
1862      ENDIF !(type_run.EQ."AMIP")
1863c
1864      if (ok_sync) then
1865c$OMP MASTER
1866        call histsync(nid_mth)
1867c$OMP END MASTER
1868      endif
1869c
1870      ENDIF !ok_mensuel
Note: See TracBrowser for help on using the repository browser.