source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmth.h @ 506

Last change on this file since 506 was 506, checked in by lmdzadmin, 20 years ago

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

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