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

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

Modifications sur l'albedo JG
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 70.9 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  temperature tendency due to moist convective processes
12       DO l=1, klev
13       DO i=1, klon
14        zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
15       ENDDO !i
16       ENDDO !l
17c
18
19cIM: 101003 : K/30min ==> K/s
20      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
21cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
22      CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d)
23c
24cIM: 101003 : K/30min ==> K/s
25      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
26cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
27      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
28c
29c
30      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
31c
32      IF(lev_histmth.GE.1) THEN
33c
34c Champs 2D:
35c
36cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
37      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis)
38C
39cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
40cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
41      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
42c
43      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
44cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
45      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
46     $                   zx_tmp_fi2d)
47c
48      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
49cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
50      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
51     $                   zx_tmp_fi2d)
52c
53cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
54      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
55
56cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
57      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
58c
59cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
60      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
61c
62c ENSEMBLES BEG
63cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
64c     CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
65c
66cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
67c     CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
68c
69c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
70c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
71c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
72c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
73c    .               iim*jjmp1,ndex2d)
74c
75      DO i=1, klon
76       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
77      ENDDO
78cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
79      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
80c
81      DO i = 1, klon
82         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
83      ENDDO
84cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
85      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
86c
87c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
88c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
89c    .                                   iim*jjmp1*klev,ndex3d)
90c
91c ENSEMBLES END
92c
93cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
94      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
95c
96cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
97      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
98C
99cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
100      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
101C
102      DO i = 1, klon
103         zx_tmp_fi2d(i) = paprs(i,1)
104      ENDDO
105cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
106      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
107c
108cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
109      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
110c
111      if (.not. ok_veget) then
112cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
113          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
114      endif
115c
116      DO i = 1, klon
117         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
118      ENDDO
119cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
120      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
121c
122cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
123c     CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
124c
125      DO i = 1, klon
126         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
127      ENDDO
128cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
129      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
130c
131      DO i = 1, klon
132         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
133      ENDDO
134cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
135      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
136c
137cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
138      CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall)
139c
140c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
141c     CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
142c    .               ndex2d)
143c
144cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
145      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
146c
147cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
148      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
149c
150cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
151      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
152c
153cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
154      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
155c
156cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
157      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
158c
159      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
160cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
161      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
162c
163c
164      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
165cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
166      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
167c
168      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
169cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
170      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
171c
172      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
173cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
174      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
175c
176cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
177      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
178c
179cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
180      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
181c
182cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
183      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
184c
185cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
186      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
187c
188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
189      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
190c
191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
192      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
193c
194      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
195cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
196      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
197c
198      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
199cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
200      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
201c
202c200 END
203cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
204      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
205c
206cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
207      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
208c
209cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
210      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
211c
212cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
213      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
214c
215cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
216      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
217c
218c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
219      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
220cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
221      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
222c
223c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
224      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
225cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
226      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
227c
228c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
229      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
230cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
231      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
232c
233c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
234      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
235cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
236      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
237c
238cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
239cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
240      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
241cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
242      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
243c
244      DO i=1, klon
245       sollwdownclr(i)=-1.*lwdn0(i,1)
246      ENDDO !i=1, klon
247
248cIM 080904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
249cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
250      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
251cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
252      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
253c
254cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
255      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
256c
257c     DO i=1, klon
258cIM 080904     sollwdownclr(i)=lwdn0(i,1)
259c      sollwdownclr(i)=-1.*lwdn0(i,1)
260c     ENDDO !i=1, klon
261c
262cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
263      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
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.GT.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) = falb1( 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, albsol1,zx_tmp_2d)
691      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol1)
692c
693cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol2,zx_tmp_2d)
694      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsol2)
695c
696      ENDIF !lev_histmth.GE.3
697c
698c
699      IF(lev_histmth.GE.4) THEN
700c
701c  FH Sorties pour la couche limite
702      CALL histwrite_phy(nid_mth,"kz",itau_w,ycoefh)
703      CALL histwrite_phy(nid_mth,"kz_max",itau_w,ycoefh)
704
705      if(iflag_pbl>1) then
706      zx_tmp_fi3d=0.
707      do nsrf=1,nbsrf
708         do k=1,klev
709          zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
710     ,    +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
711         enddo
712      enddo
713      CALL histwrite_phy(nid_mth,"tke",itau_w,zx_tmp_fi3d)
714      CALL histwrite_phy(nid_mth,"tke_max",itau_w,zx_tmp_fi3d)
715      endif
716
717
718
719cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
720      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
721
722      IF (iflag_con.GE.3) THEN
723c
724cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
725      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
726c
727cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
728      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
729c
730cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
731      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
732c
733cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
734      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
735c
736      ENDIF !iflag_con.GE.3
737c
738cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
739      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
740c
741cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
742      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
743c
744cIM: 101003 : K/30min ==> K/s
745      zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
746cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
747      CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
748c
749      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
750cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
751      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
752c
753cIM: 101003 : K/30min ==> K/s
754      zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
755cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
756      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
757cIM: 071003
758cIM: 101003 : K/30min ==> K/s
759      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
760     .                             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,"dtlschr",itau_w,zx_tmp_fi3d)
763c
764      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(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,"dqlsc",itau_w,zx_tmp_fi3d)
767c
768cIM: 101003 : K/30min ==> K/s
769      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
770cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
771      CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
772c
773      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
774cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
775      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
776c
777cIM: 101003 : K/30min ==> K/s
778      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
779cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
780      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
781c
782      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
783cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
784      CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
785c
786      zpt_conv = 0.
787      where (ptconv) zpt_conv = 1.
788cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
789      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
790c
791cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
792      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
793c
794cIM: 101003 : K/30min ==> K/s
795      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
796cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
797      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
798c
799      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
800cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
801      CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
802c
803cIM: 101003 : K/30min ==> K/s
804      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
805cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
806      CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d)
807c
808      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
809cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
810      CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d)
811c
812cIM: 101003 : K/day ==> K/s
813      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
814cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
815      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
816c
817cIM: 101003 : K/day ==> K/s     
818      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
819cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
820      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
821c
822cIM: 101003 : K/day ==> K/s     
823      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
824cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
825      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
826c
827cIM: 101003 : K/day ==> K/s     
828      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
829cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
830      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
831c
832cIM: 101003 : deja en K/s     
833      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
834cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
835      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
836c
837      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
838cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
839      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
840c
841      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
842cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
843      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
844c
845      IF (ok_orodr) THEN
846      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
847cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
848      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
849c
850      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
851cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
852      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
853c
854      ENDIF
855C
856      IF (ok_orolf) THEN
857      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
858cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
859      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
860c
861      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
862cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
863      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
864      ENDIF
865C
866#if defined(INCA_AER) && defined(CPP_COUPLE)
867      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
868      CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
869      CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
870      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
871      CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
872      CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
873      CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
874      CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
875      CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
876      CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
877      CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
878      CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
879      CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
880      CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
881#endif
882#ifndef INCA
883      IF (nqmax.GE.3) THEN
884      DO iq=3,nqmax
885cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
886         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
887      ENDDO
888      ENDIF
889#endif
890c
891      ENDIF !lev_histmth.GE.4
892c
893      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
894c
895      IF(lev_histmth.GE.1) THEN
896c
897c Champs 2D:
898c
899cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
900      CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG)
901C
902cIM 310804   CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
903cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
904      CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
905c
906      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
907cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
908      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
909     $      zx_tmp_fi2d)
910c
911      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
912cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
913      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
914     $      zx_tmp_fi2d)
915c
916      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100.
917cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
918      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
919     $      zx_tmp_fi2d)
920c
921      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)
922cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
923      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
924     $      zx_tmp_fi2d)
925c
926      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100.
927cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
928      CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
929     $      zx_tmp_fi2d)
930c
931      zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)
932cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
933      CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
934     $      zx_tmp_fi2d)
935c
936cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
937      CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat)
938c
939      zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon)
940cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
941      CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d)
942c
943c     IF(itap.EQ.ecrit_mth) THEN
944c       DO i=1, klon
945c         IF(nday_clr(i).GT.0.) THEN
946c          tops01(i)=tops01(i)/nday_clr(i)
947c         ELSE
948c          tops01(i)=1.e+20
949c         ENDIF
950c       ENDDO
951c       CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d)
952c       CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d,
953c    .       iim*jjmp1,ndex2d)
954c     ENDIF !itap.EQ.ecrit_mth
955c
956c220404 BEG
957c
958c temperature tendency due to total diabatic heating
959       DO l=1, klev
960       DO i=1, klon
961       d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime
962       ENDDO !i
963       ENDDO !l
964c
965cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
966      CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t)
967c
968c temperature tendency due to SW radiation
969       DO l=1, klev
970       DO i=1, klon
971       zx_tmp_fi3d(i,l)=heat(i,l)/RDAY
972       ENDDO !i
973       ENDDO !l
974c
975cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
976      CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d)
977c
978c temperature tendency due to LW radiation
979       DO l=1, klev
980       DO i=1, klon
981        zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY
982       ENDDO !i
983       ENDDO !l
984c
985cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
986      CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
987c
988
989cIM: 101003 : K/30min ==> K/s
990      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
991cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
992      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
993
994c
995      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
996cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
997      CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d)
998c
999
1000
1001
1002
1003c
1004c  temperature tendency due to large scale precipitation
1005       DO l=1, klev
1006       DO i=1, klon
1007       zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys
1008       ENDDO !i
1009       ENDDO !l
1010c
1011cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
1012      CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
1013c
1014c  total moisture tendency due to diabatic processes
1015       DO l=1, klev
1016       DO i=1, klon
1017        d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime
1018        zx_tmp_fi3d(i,l)=d_qx(i,l,ivap)
1019       ENDDO !i
1020       ENDDO !l
1021c
1022cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
1023      CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d)
1024c
1025cBEG 210404 BEG
1026c
1027      IF(ok_orodr.AND.ok_orolf) THEN
1028      DO k = 1, klev
1029      DO i = 1, klon
1030        d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys
1031        d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys
1032      ENDDO
1033      ENDDO
1034      ENDIF !(ok_orodr.AND.ok_orolf) THEN
1035c
1036cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
1037      CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli)
1038c
1039cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
1040      CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli)
1041c
1042      zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1043cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1044      CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d)
1045c
1046      zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1047cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1048      CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d)
1049c
1050cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
1051      CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u)
1052c
1053cym      CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
1054      CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v)
1055
1056cEND 210404 END
1057c
1058c cldtrue
1059c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d)
1060c     CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
1061c    $       iim*jjmp1*klev,ndex3d)
1062c
1063c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d)
1064c     CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
1065c    $       iim*jjmp1*klev,ndex3d)
1066c
1067c     CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d)
1068c     CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
1069c    $       iim*jjmp1*klev,ndex3d)
1070c
1071      DO i=1, klon
1072cIM 151004    zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1)
1073       zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i)
1074      ENDDO
1075cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1076      CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d)
1077c
1078      DO i=1, klon
1079cIM 151004    zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1)
1080       zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i)
1081      ENDDO
1082cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1083      CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d)
1084c
1085c220404 END
1086cAMIP2 END
1087c
1088cym     CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
1089      CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
1090
1091cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
1092      CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
1093c
1094cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
1095      CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
1096c
1097c ENSEMBLES BEG
1098cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
1099c     CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
1100c
1101cym     CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
1102c     CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
1103c
1104c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
1105c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
1106c     CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
1107c     CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
1108c    .               iim*jjmp1,ndex2d)
1109c
1110      DO i=1, klon
1111       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
1112      ENDDO
1113cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1114      CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
1115c
1116      DO i = 1, klon
1117         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
1118      ENDDO
1119cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1120      CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
1121c
1122c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
1123c     CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
1124c    .                                   iim*jjmp1*klev,ndex3d)
1125c
1126c ENSEMBLES END
1127c
1128cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
1129      CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
1130c
1131cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
1132      CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
1133C
1134cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
1135      CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
1136C
1137      DO i = 1, klon
1138         zx_tmp_fi2d(i) = paprs(i,1)
1139      ENDDO
1140cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1141      CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
1142c
1143cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
1144      CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
1145c
1146      if (.not. ok_veget) then
1147cym          CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
1148          CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
1149      endif
1150c
1151      DO i = 1, klon
1152         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
1153      ENDDO
1154cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1155      CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
1156c
1157cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
1158c     CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
1159c
1160      DO i = 1, klon
1161         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
1162      ENDDO
1163cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1164      CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
1165c
1166      DO i = 1, klon
1167         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
1168      ENDDO
1169cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1170      CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
1171c
1172cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
1173      CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall)
1174c
1175c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
1176c     CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
1177c    .               ndex2d)
1178c
1179cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
1180      CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
1181c
1182cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
1183      CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
1184c
1185cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
1186      CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
1187c
1188cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
1189      CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
1190c
1191cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
1192      CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
1193c
1194      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
1195cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1196      CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
1197c
1198c
1199      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
1200cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1201      CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
1202c
1203      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
1204cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1205      CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
1206c
1207      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
1208cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1209      CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
1210c
1211cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
1212      CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
1213c
1214cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
1215      CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
1216c
1217cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
1218      CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
1219c
1220cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
1221      CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
1222c
1223cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
1224      CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
1225c
1226cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
1227      CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
1228c
1229      zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
1230cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1231      CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
1232c
1233      zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
1234cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1235      CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
1236c
1237c200 END
1238cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
1239      CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
1240c
1241cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
1242      CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
1243c
1244cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
1245      CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
1246c
1247cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
1248      CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
1249c
1250cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
1251      CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
1252c
1253c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
1254      zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
1255cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1256      CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
1257c
1258c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
1259      zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
1260cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1261      CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
1262c
1263c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
1264      zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
1265cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1266      CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
1267c
1268c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
1269      zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
1270cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1271      CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
1272c
1273cIM 080904   zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
1274cIM 130904   zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
1275      zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
1276cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
1277      CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
1278c
1279      DO i=1, klon
1280       sollwdownclr(i)=-1.*lwdn0(i,1)
1281      ENDDO !i=1, klon
1282c
1283cIM 080904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
1284cIM 130904   zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
1285      zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
1286cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
1287      CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
1288c
1289cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
1290      CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
1291c
1292c     DO i=1, klon
1293cIM 08004    sollwdownclr(i)=lwdn0(i,1)
1294c      sollwdownclr(i)=-1.*lwdn0(i,1)
1295c     ENDDO !i=1, klon
1296c
1297cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
1298      CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
1299c
1300cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
1301      CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
1302c
1303      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
1304c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
1305cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1306      CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
1307c
1308cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
1309      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
1310c
1311c
1312c      DO i = 1, klon
1313c         zx_tmp_fi2d(i) = fluxu(i,1)
1314c      ENDDO
1315c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1316c      CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1317c
1318c      DO i = 1, klon
1319c         zx_tmp_fi2d(i) = fluxv(i,1)
1320c      ENDDO
1321c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1322c      CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
1323c
1324c
1325cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
1326       CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
1327c
1328cym       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
1329       CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
1330c
1331       DO nsrf = 1, nbsrf
1332        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
1333cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1334        CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
1335     $      zx_tmp_fi2d) 
1336C     
1337        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
1338cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1339        CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
1340     $      zx_tmp_fi2d)
1341       ENDDO
1342c
1343        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
1344cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1345        CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
1346     $      zx_tmp_fi2d) 
1347c
1348        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
1349cym        CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
1350        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
1351     $      zx_tmp_fi2d) 
1352c
1353      DO nsrf = 1, nbsrf
1354CYYY
1355C
1356        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
1357cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1358        CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
1359     $      zx_tmp_fi2d)
1360C
1361        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
1362cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1363        CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
1364     $      zx_tmp_fi2d)
1365C
1366        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
1367cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1368        CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
1369     $      zx_tmp_fi2d) 
1370C
1371        zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
1372cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1373        CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
1374     $      zx_tmp_fi2d)
1375c
1376        zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
1377cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1378        CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
1379     $      zx_tmp_fi2d)
1380C
1381        zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
1382cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1383        CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
1384     $      zx_tmp_fi2d)
1385C
1386        zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
1387cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1388        CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
1389     $      zx_tmp_fi2d)
1390C
1391       IF(1.EQ.0) THEN
1392C
1393        zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
1394cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1395        CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
1396     $      zx_tmp_fi2d)
1397C
1398        zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
1399     $                        *pctsrf(1:klon,nsrf)
1400cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1401        CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
1402     $      zx_tmp_fi2d)
1403C
1404       ENDIF !(1.EQ.0) THEN
1405C
1406      END DO 
1407c
1408cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
1409      CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
1410c
1411cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
1412      CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
1413c
1414cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
1415      CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
1416c
1417cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
1418      CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
1419c
1420cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
1421      CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
1422c
1423cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
1424      CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
1425c
1426cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
1427      CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
1428c
1429      zx_tmp_fi2d(1:klon) = flwp(1:klon)
1430cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1431      CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
1432c
1433      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
1434cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1435      CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
1436c
1437cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
1438      CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
1439c
1440cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
1441      CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
1442c
1443cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
1444      CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
1445c
1446cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
1447      CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
1448cKE43
1449      IF (iflag_con.GE.3) THEN ! sb
1450c
1451cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
1452      CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
1453c
1454cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
1455      CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
1456c
1457cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
1458      CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
1459c
1460cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
1461      CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
1462c
1463c34EK
1464c
1465cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
1466      CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
1467c
1468      ENDIF !iflag_con.GE.3
1469c
1470cIM cf. AM 081204 BEG
1471cHBTM2
1472c
1473c
1474cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
1475      CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
1476c
1477cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
1478      CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
1479c
1480cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
1481      CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
1482c
1483cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
1484      CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
1485c
1486cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
1487      CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
1488c
1489cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
1490      CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
1491c
1492cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
1493      CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
1494c
1495cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
1496      CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
1497c
1498cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
1499      CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
1500c
1501cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
1502      CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
1503c
1504cIM cf. AM 081204 END
1505c
1506c  Champs interpolles sur des niveaux de pression
1507c     DO k=1, nlevENS
1508      DO k=1, nlevSTD
1509c
1510c        bb=clevSTD(k)
1511c
1512c        IF(k.GE.2) THEN
1513         IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
1514         IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
1515c         aa=clevSTD(k)
1516c         bb=aa(1:lnblnk1(aa))
1517c        ENDIF
1518c
1519       IF(bb2.EQ."850".OR.bb2.EQ."700"
1520     $   .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
1521c
1522cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
1523        CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
1524c
1525cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
1526        CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
1527c
1528cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
1529      CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
1530c
1531cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
1532      CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
1533c
1534      ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
1535      ENDDO !nlevSTD
1536c
1537      DO i=1, klon
1538       IF (pctsrf(i,is_oce).GT.epsfra) THEN
1539        zx_tmp_fi2d(i) = fluxo(i)
1540       ELSE
1541        zx_tmp_fi2d(i) = 0.
1542       ENDIF
1543      ENDDO
1544c
1545cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1546      CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
1547c
1548      DO i=1, klon
1549       IF (pctsrf(i,is_sic).GT.epsfra) THEN
1550        zx_tmp_fi2d(i) = fluxg(i)
1551       ELSE
1552        zx_tmp_fi2d(i) = 0.
1553       ENDIF
1554      ENDDO
1555c
1556cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1557      CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
1558c
1559      DO i=1, klon
1560       IF (pctsrf(i,is_oce).GT.epsfra.OR.
1561     .     pctsrf(i,is_sic).GT.epsfra) THEN
1562        zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ 
1563     .                   ftsol(i, is_sic) * pctsrf(i,is_sic))/
1564     .                   (pctsrf(i,is_oce)+pctsrf(i,is_sic))
1565       ELSE
1566        zx_tmp_fi2d(i) = 273.15
1567       ENDIF
1568      ENDDO
1569c
1570cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1571      CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
1572c
1573      IF (OCEAN.EQ.'force ') THEN
1574c
1575      DO i=1, klon
1576      IF((pctsrf(i,is_oce).GT.epsfra).OR.
1577     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
1578       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
1579     .                  fluxg(i)*pctsrf(i,is_sic)
1580      ELSE
1581       zx_tmp_fi2d(i) = 1.E+20
1582      ENDIF
1583      ENDDO
1584cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1585      CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
1586     .     zx_tmp_fi2d)
1587c
1588      ELSE IF (OCEAN.EQ.'slab  ') THEN
1589c
1590c
1591      DO i=1, klon
1592      IF((pctsrf(i,is_oce).GT.epsfra).OR.
1593     .   (pctsrf(i,is_sic).GT.epsfra)) THEN
1594       zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
1595     .                  fluxg(i)*pctsrf(i,is_sic)
1596      ELSE
1597       zx_tmp_fi2d(i) = 1.E+20
1598      ENDIF
1599      ENDDO
1600cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1601      CALL histwrite_phy(nid_mth,"slab_bils",itau_w,
1602     .     zx_tmp_fi2d)
1603c
1604      DO i=1, klon
1605       IF(pctsrf(i,is_oce).GT.epsfra.OR.
1606     $    pctsrf(i,is_sic).GT.epsfra) THEN
1607        zx_tmp_fi2d(i)=tslab(i)
1608       ELSE
1609        zx_tmp_fi2d(i) = 1.E+20
1610       ENDIF
1611      ENDDO !i=1, klon
1612cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
1613      CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
1614c
1615cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
1616      CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
1617c
1618cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
1619      CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
1620c
1621      ENDIF !OCEAN.EQ.force/slab
1622c
1623      ENDIF !lev_histmth.GE.1
1624c
1625      IF(lev_histmth.GE.2) THEN
1626c
1627c Champs 3D:
1628C
1629c
1630cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
1631      CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
1632c
1633cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
1634      CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
1635c
1636cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
1637      CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
1638c
1639c#ifdef histmthNMC
1640c#include "write_histmthNMC.h"
1641c#endif
1642c
1643cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
1644      CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
1645c
1646c     zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
1647c    .                         (1-qx(1:klon,1:klev,ivap))
1648c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1649c     CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
1650c    .                                   iim*jjmp1*klev,ndex3d)
1651c
1652cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
1653      CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
1654c
1655cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
1656      CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
1657c
1658cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
1659      CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
1660c
1661cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
1662      CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
1663c
1664cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
1665      CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
1666c
1667cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
1668      CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
1669c
1670cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
1671      CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
1672c
1673cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
1674      CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
1675c
1676      DO k=1, klev
1677       DO i=1, klon
1678cIM 151004        zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
1679c    $                    /(paprs(i,k)-paprs(i,k+1))
1680c    $                    *(paprs(i,1)/101325.0)
1681c        zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968
1682         zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
1683     $                    /(paprs(i,k)-paprs(i,k+1))
1684     $                    *(paprs(i,1)/101325.0)
1685       ENDDO !i
1686      ENDDO !k
1687c
1688c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
1689c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
1690cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1691      CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
1692c
1693cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
1694      CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
1695c
1696cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
1697      CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
1698c
1699      ENDIF !lev_histmth.GE.2
1700c
1701      IF(lev_histmth.GE.3) THEN
1702c
1703       DO nsrf=1, nbsrf
1704c
1705        zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1706cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1707        CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
1708     $      zx_tmp_fi2d)
1709C
1710        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1711cym        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
1712        CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
1713     $      zx_tmp_fi2d)
1714c
1715      zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1716cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
1717      CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
1718     $    ,agesno)
1719c
1720      ENDDO !nsrf=1, nbsrf
1721c
1722cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol1,zx_tmp_2d)
1723      CALL histwrite_phy(nid_mth,"albs",itau_w,albsol1)
1724c
1725cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol2,zx_tmp_2d)
1726      CALL histwrite_phy(nid_mth,"albslw",itau_w,albsol2)
1727c
1728      ENDIF !lev_histmth.GE.3
1729c
1730c
1731      IF(lev_histmth.GE.4) THEN
1732c
1733c  FH Sorties pour la couche limite
1734      CALL histwrite_phy(nid_mth,"kz",itau_w,ycoefh)
1735      CALL histwrite_phy(nid_mth,"kz_max",itau_w,ycoefh)
1736
1737      if(iflag_pbl>1) then
1738      zx_tmp_fi3d=0.
1739      do nsrf=1,nbsrf
1740         do k=1,klev
1741          zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
1742     ,    +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1743         enddo
1744      enddo
1745      CALL histwrite_phy(nid_mth,"tke",itau_w,zx_tmp_fi3d)
1746      CALL histwrite_phy(nid_mth,"tke_max",itau_w,zx_tmp_fi3d)
1747      endif
1748
1749
1750
1751cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
1752      CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
1753c
1754      IF (iflag_con.GE.3) THEN
1755c
1756cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
1757      CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
1758c
1759cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
1760      CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
1761c
1762cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
1763      CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
1764c
1765cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
1766      CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
1767c
1768      ENDIF !iflag_con.GE.3
1769c
1770cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
1771      CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
1772c
1773cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
1774      CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
1775c
1776cIM: 101003 : K/30min ==> K/s
1777c
1778      zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1779cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1780      CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
1781c
1782cIM: 101003 : K/30min ==> K/s
1783      zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1784     .                             d_t_eva(1:klon,1:klev))/pdtphys
1785cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
1786      CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
1787c
1788      zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1789cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1790      CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
1791c
1792cIM: 101003 : K/30min ==> K/s
1793      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(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,"dtvdf",itau_w,zx_tmp_fi3d)
1796c
1797      zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1798cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1799      CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
1800c
1801cIM: 101003 : K/30min ==> K/s
1802      zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1803cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1804      CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
1805c
1806      zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(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,"dqeva",itau_w,zx_tmp_fi3d)
1809c
1810      zpt_conv = 0.
1811      where (ptconv) zpt_conv = 1.
1812cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
1813      CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
1814c
1815cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
1816      CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
1817cIM: 101003 : K/30min ==> K/s
1818      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1819cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1820      CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
1821c
1822      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(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,"dqajs",itau_w,zx_tmp_fi3d)
1825c
1826cIM: 101003 : K/30min ==> K/s
1827c
1828      zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
1829cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1830      CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d)
1831c
1832c
1833cIM: 101003 : K/day ==> K/s
1834      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1835cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1836      CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
1837c
1838cIM: 101003 : K/day ==> K/s     
1839      zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1840cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1841      CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
1842c
1843cIM: 101003 : K/day ==> K/s     
1844      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1845cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1846      CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
1847c
1848cIM: 101003 : K/day ==> K/s     
1849      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1850cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1851      CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
1852c
1853cIM: 101003 : deja en K/s     
1854      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1855cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1856      CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
1857c
1858      zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1859cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1860      CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
1861c
1862      zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1863cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1864      CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
1865c
1866      IF (ok_orodr) THEN
1867      zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1868cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1869      CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
1870c
1871      zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1872cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1873      CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
1874c
1875      ENDIF
1876C
1877      IF (ok_orolf) THEN
1878      zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1879cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1880      CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
1881c
1882      zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1883cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
1884      CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
1885      ENDIF
1886C
1887#if defined(INCA_AER) && defined(CPP_COUPLE)
1888      CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
1889      CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
1890      CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
1891      CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
1892      CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
1893      CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
1894      CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
1895      CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
1896      CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
1897      CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
1898      CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
1899      CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
1900      CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
1901      CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
1902#endif
1903#ifndef INCA
1904      IF (nqmax.GE.3) THEN
1905      DO iq=3,nqmax
1906cym         CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
1907         CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
1908      ENDDO
1909      ENDIF
1910#endif
1911c
1912      ENDIF !lev_histmth.GE.4
1913c
1914      ENDIF !(type_run.EQ."AMIP")
1915c
1916      if (ok_sync) then
1917c$OMP MASTER
1918        call histsync(nid_mth)
1919c$OMP END MASTER
1920      endif
1921c
1922      ENDIF !ok_mensuel
Note: See TracBrowser for help on using the repository browser.