Context Navigation

source: LMDZ5/trunk/libf/phylmd/phys_output_write.h @ 1518

Last change on this file since 1518 was 1518, checked in by idelkadi, 13 years ago

Modifications des routines de convection :

correction de bug : le champ sigd n'etait pas decompresse avant d'etre sorti de la convection.
sortir de nouveaux champs de convection (plcl, plfc, wb)

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 54.5 KB

Line
1	itau_w = itau_phy + itap
2
3	DO iff=1,nfiles
4
5	IF (clef_files(iff)) THEN
6	ndex2d = 0
7	ndex3d = 0
8
9	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
10	IF (o_phis%flag(iff)<=lev_files(iff)) THEN
11	CALL histwrite_phy(nid_files(iff),
12	$ o_phis%name,itau_w,pphis)
13	ENDIF
14
15	IF (o_aire%flag(iff)<=lev_files(iff)) THEN
16	CALL histwrite_phy(nid_files(iff),o_aire%name,itau_w,airephy)
17	ENDIF
18
19	IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN
20	DO i=1, klon
21	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
22	ENDDO
23	CALL histwrite_phy(nid_files(iff),
24	$ o_contfracATM%name,itau_w,zx_tmp_fi2d)
25	ENDIF
26
27	IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN
28	CALL histwrite_phy(nid_files(iff),o_contfracOR%name,itau_w,
29	$ pctsrf(:,is_ter))
30	ENDIF
31
32	IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN
33	CALL histwrite_phy(nid_files(iff),
34	$ o_aireTER%name,itau_w,paire_ter)
35	ENDIF
36
37	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
38
39	IF (o_flat%flag(iff)<=lev_files(iff)) THEN
40	CALL histwrite_phy(nid_files(iff),o_flat%name,itau_w,zxfluxlat)
41	ENDIF
42
43	IF (o_slp%flag(iff)<=lev_files(iff)) THEN
44	CALL histwrite_phy(nid_files(iff),o_slp%name,itau_w,slp)
45	ENDIF
46
47	IF (o_tsol%flag(iff)<=lev_files(iff)) THEN
48	CALL histwrite_phy(nid_files(iff),o_tsol%name,itau_w,zxtsol)
49	ENDIF
50
51	IF (o_t2m%flag(iff)<=lev_files(iff)) THEN
52	CALL histwrite_phy(nid_files(iff),o_t2m%name,itau_w,zt2m)
53	ENDIF
54
55	IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN
56	CALL histwrite_phy(nid_files(iff),o_t2m_min%name,itau_w,zt2m)
57	ENDIF
58
59	IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN
60	CALL histwrite_phy(nid_files(iff),o_t2m_max%name,itau_w,zt2m)
61	ENDIF
62
63	IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN
64	DO i=1, klon
65	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
66	ENDDO
67	CALL histwrite_phy(nid_files(iff),
68	s o_wind10m%name,itau_w,zx_tmp_fi2d)
69	ENDIF
70
71	IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN
72	DO i=1, klon
73	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
74	ENDDO
75	CALL histwrite_phy(nid_files(iff),o_wind10max%name,
76	$ itau_w,zx_tmp_fi2d)
77	ENDIF
78
79	IF (o_sicf%flag(iff)<=lev_files(iff)) THEN
80	DO i = 1, klon
81	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
82	ENDDO
83	CALL histwrite_phy(nid_files(iff),
84	$ o_sicf%name,itau_w,zx_tmp_fi2d)
85	ENDIF
86
87	IF (o_q2m%flag(iff)<=lev_files(iff)) THEN
88	CALL histwrite_phy(nid_files(iff),o_q2m%name,itau_w,zq2m)
89	ENDIF
90
91	IF (o_u10m%flag(iff)<=lev_files(iff)) THEN
92	CALL histwrite_phy(nid_files(iff),o_u10m%name,itau_w,zu10m)
93	ENDIF
94
95	IF (o_v10m%flag(iff)<=lev_files(iff)) THEN
96	CALL histwrite_phy(nid_files(iff),o_v10m%name,itau_w,zv10m)
97	ENDIF
98
99	IF (o_psol%flag(iff)<=lev_files(iff)) THEN
100	DO i = 1, klon
101	zx_tmp_fi2d(i) = paprs(i,1)
102	ENDDO
103	CALL histwrite_phy(nid_files(iff),
104	s o_psol%name,itau_w,zx_tmp_fi2d)
105	ENDIF
106
107	IF (o_mass%flag(iff)<=lev_files(iff)) THEN
108	CALL histwrite_phy(nid_files(iff),o_mass%name,itau_w,zmasse)
109	ENDIF
110
111
112	IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN
113	CALL histwrite_phy(nid_files(iff),o_qsurf%name,itau_w,zxqsurf)
114	ENDIF
115
116	if (.not. ok_veget) then
117	IF (o_qsol%flag(iff)<=lev_files(iff)) THEN
118	CALL histwrite_phy(nid_files(iff),o_qsol%name,itau_w,qsol)
119	ENDIF
120	endif
121
122	IF (o_precip%flag(iff)<=lev_files(iff)) THEN
123	DO i = 1, klon
124	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
125	ENDDO
126	CALL histwrite_phy(nid_files(iff),o_precip%name,
127	s itau_w,zx_tmp_fi2d)
128	ENDIF
129
130	IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN
131	CALL histwrite_phy(nid_files(iff),o_ndayrain%name,
132	s itau_w,nday_rain)
133	ENDIF
134
135	IF (o_plul%flag(iff)<=lev_files(iff)) THEN
136	DO i = 1, klon
137	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
138	ENDDO
139	CALL histwrite_phy(nid_files(iff),o_plul%name,itau_w,zx_tmp_fi2d)
140	ENDIF
141
142	IF (o_pluc%flag(iff)<=lev_files(iff)) THEN
143	DO i = 1, klon
144	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
145	ENDDO
146	CALL histwrite_phy(nid_files(iff),o_pluc%name,itau_w,zx_tmp_fi2d)
147	ENDIF
148
149	IF (o_snow%flag(iff)<=lev_files(iff)) THEN
150	CALL histwrite_phy(nid_files(iff),o_snow%name,itau_w,snow_fall)
151	ENDIF
152
153	IF (o_msnow%flag(iff)<=lev_files(iff)) THEN
154	CALL histwrite_phy(nid_files(iff),o_msnow%name,itau_w,snow_o)
155	ENDIF
156
157	IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN
158	CALL histwrite_phy(nid_files(iff),o_fsnow%name,itau_w,zfra_o)
159	ENDIF
160
161	IF (o_evap%flag(iff)<=lev_files(iff)) THEN
162	CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap)
163	ENDIF
164
165	IF (o_tops%flag(iff)<=lev_files(iff)) THEN
166	CALL histwrite_phy(nid_files(iff),o_tops%name,itau_w,topsw)
167	ENDIF
168
169	IF (o_tops0%flag(iff)<=lev_files(iff)) THEN
170	CALL histwrite_phy(nid_files(iff),o_tops0%name,itau_w,topsw0)
171	ENDIF
172
173	IF (o_topl%flag(iff)<=lev_files(iff)) THEN
174	CALL histwrite_phy(nid_files(iff),o_topl%name,itau_w,toplw)
175	ENDIF
176
177	IF (o_topl0%flag(iff)<=lev_files(iff)) THEN
178	CALL histwrite_phy(nid_files(iff),o_topl0%name,itau_w,toplw0)
179	ENDIF
180
181	IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN
182	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
183	CALL histwrite_phy(nid_files(iff),o_SWupTOA%name,
184	s itau_w,zx_tmp_fi2d)
185	ENDIF
186
187	IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN
188	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
189	CALL histwrite_phy(nid_files(iff),
190	$ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d)
191	ENDIF
192
193	IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN
194	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
195	CALL histwrite_phy(nid_files(iff),
196	s o_SWdnTOA%name,itau_w,zx_tmp_fi2d)
197	ENDIF
198
199	IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN
200	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
201	CALL histwrite_phy(nid_files(iff),
202	$ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d)
203	ENDIF
204
205	IF (o_nettop%flag(iff)<=lev_files(iff)) THEN
206	zx_tmp_fi2d(:) = topsw(:)-toplw(:)
207	CALL histwrite_phy(nid_files(iff),
208	$ o_nettop%name,itau_w,zx_tmp_fi2d)
209	ENDIF
210
211	IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
212	CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200)
213	ENDIF
214
215	IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN
216	CALL histwrite_phy(nid_files(iff),
217	s o_SWup200clr%name,itau_w,SWup200clr)
218	ENDIF
219
220	IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN
221	CALL histwrite_phy(nid_files(iff),o_SWdn200%name,itau_w,SWdn200)
222	ENDIF
223
224	IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN
225	CALL histwrite_phy(nid_files(iff),
226	s o_SWdn200clr%name,itau_w,SWdn200clr)
227	ENDIF
228
229	IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN
230	CALL histwrite_phy(nid_files(iff),o_LWup200%name,itau_w,LWup200)
231	ENDIF
232
233	IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN
234	CALL histwrite_phy(nid_files(iff),
235	s o_LWup200clr%name,itau_w,LWup200clr)
236	ENDIF
237
238	IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN
239	CALL histwrite_phy(nid_files(iff),
240	s o_LWdn200%name,itau_w,LWdn200)
241	ENDIF
242
243	IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN
244	CALL histwrite_phy(nid_files(iff),
245	s o_LWdn200clr%name,itau_w,LWdn200clr)
246	ENDIF
247
248	IF (o_sols%flag(iff)<=lev_files(iff)) THEN
249	CALL histwrite_phy(nid_files(iff),o_sols%name,itau_w,solsw)
250	ENDIF
251
252	IF (o_sols0%flag(iff)<=lev_files(iff)) THEN
253	CALL histwrite_phy(nid_files(iff),o_sols0%name,itau_w,solsw0)
254	ENDIF
255
256	IF (o_soll%flag(iff)<=lev_files(iff)) THEN
257	CALL histwrite_phy(nid_files(iff),o_soll%name,itau_w,sollw)
258	ENDIF
259
260	IF (o_radsol%flag(iff)<=lev_files(iff)) THEN
261	CALL histwrite_phy(nid_files(iff),o_radsol%name,itau_w,radsol)
262	ENDIF
263
264	IF (o_soll0%flag(iff)<=lev_files(iff)) THEN
265	CALL histwrite_phy(nid_files(iff),o_soll0%name,itau_w,sollw0)
266	ENDIF
267
268	IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN
269	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
270	CALL histwrite_phy(nid_files(iff),
271	s o_SWupSFC%name,itau_w,zx_tmp_fi2d)
272	ENDIF
273
274	IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN
275	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
276	CALL histwrite_phy(nid_files(iff),
277	$ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d)
278	ENDIF
279
280	IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN
281	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
282	CALL histwrite_phy(nid_files(iff),
283	$ o_SWdnSFC%name,itau_w,zx_tmp_fi2d)
284	ENDIF
285
286	IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
287	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
288	CALL histwrite_phy(nid_files(iff),
289	$ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d)
290	ENDIF
291
292	IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN
293	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
294	CALL histwrite_phy(nid_files(iff),
295	$ o_LWupSFC%name,itau_w,zx_tmp_fi2d)
296	ENDIF
297
298	IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN
299	CALL histwrite_phy(nid_files(iff),
300	$ o_LWdnSFC%name,itau_w,sollwdown)
301	ENDIF
302
303	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
304	IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN
305	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
306	CALL histwrite_phy(nid_files(iff),
307	$ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d)
308	ENDIF
309
310	IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
311	CALL histwrite_phy(nid_files(iff),
312	$ o_LWdnSFCclr%name,itau_w,sollwdownclr)
313	ENDIF
314
315	IF (o_bils%flag(iff)<=lev_files(iff)) THEN
316	CALL histwrite_phy(nid_files(iff),o_bils%name,itau_w,bils)
317	ENDIF
318
319	IF (o_sens%flag(iff)<=lev_files(iff)) THEN
320	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
321	CALL histwrite_phy(nid_files(iff),o_sens%name,itau_w,zx_tmp_fi2d)
322	ENDIF
323
324	IF (o_fder%flag(iff)<=lev_files(iff)) THEN
325	CALL histwrite_phy(nid_files(iff),o_fder%name,itau_w,fder)
326	ENDIF
327
328	IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN
329	CALL histwrite_phy(nid_files(iff),o_ffonte%name,itau_w,zxffonte)
330	ENDIF
331
332	IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN
333	CALL histwrite_phy(nid_files(iff),
334	$ o_fqcalving%name,itau_w,zxfqcalving)
335	ENDIF
336
337	IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN
338	CALL histwrite_phy(nid_files(iff),
339	$ o_fqfonte%name,itau_w,zxfqfonte)
340	ENDIF
341
342	IF (o_taux%flag(iff)<=lev_files(iff)) THEN
343	zx_tmp_fi2d=0.
344	do nsrf=1,nbsrf
345	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
346	enddo
347	CALL histwrite_phy(nid_files(iff),
348	$ o_taux%name,itau_w,zx_tmp_fi2d)
349	ENDIF
350
351	IF (o_tauy%flag(iff)<=lev_files(iff)) THEN
352	zx_tmp_fi2d=0.
353	do nsrf=1,nbsrf
354	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
355	enddo
356	CALL histwrite_phy(nid_files(iff),
357	$ o_tauy%name,itau_w,zx_tmp_fi2d)
358	ENDIF
359
360	DO nsrf = 1, nbsrf
361	! IF(nsrf.GE.2) THEN
362	IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
363	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
364	CALL histwrite_phy(nid_files(iff),
365	$ o_pourc_srf(nsrf)%name,itau_w,
366	$ zx_tmp_fi2d)
367	ENDIF
368
369	IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
370	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
371	CALL histwrite_phy(nid_files(iff),
372	$ o_fract_srf(nsrf)%name,itau_w,
373	$ zx_tmp_fi2d)
374	ENDIF
375	! ENDIF !nsrf.GT.2
376
377	IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
378	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
379	CALL histwrite_phy(nid_files(iff),
380	$ o_taux_srf(nsrf)%name,itau_w,
381	$ zx_tmp_fi2d)
382	ENDIF
383
384	IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
385	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
386	CALL histwrite_phy(nid_files(iff),
387	$ o_tauy_srf(nsrf)%name,itau_w,
388	$ zx_tmp_fi2d)
389	ENDIF
390
391	IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
392	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
393	CALL histwrite_phy(nid_files(iff),
394	$ o_tsol_srf(nsrf)%name,itau_w,
395	$ zx_tmp_fi2d)
396	ENDIF
397
398	IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
399	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
400	CALL histwrite_phy(nid_files(iff),o_u10m_srf(nsrf)%name,
401	$ itau_w,zx_tmp_fi2d)
402	ENDIF
403
404	IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
405	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
406	CALL histwrite_phy(nid_files(iff),o_v10m_srf(nsrf)%name,
407	$ itau_w,zx_tmp_fi2d)
408	ENDIF
409
410	IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
411	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
412	CALL histwrite_phy(nid_files(iff),o_t2m_srf(nsrf)%name,
413	$ itau_w,zx_tmp_fi2d)
414	ENDIF
415
416	IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
417	zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
418	CALL histwrite_phy(nid_files(iff),o_evap_srf(nsrf)%name,
419	$ itau_w,zx_tmp_fi2d)
420	ENDIF
421
422	IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
423	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
424	CALL histwrite_phy(nid_files(iff),
425	$ o_sens_srf(nsrf)%name,itau_w,
426	$ zx_tmp_fi2d)
427	ENDIF
428
429	IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
430	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
431	CALL histwrite_phy(nid_files(iff),
432	$ o_lat_srf(nsrf)%name,itau_w,
433	$ zx_tmp_fi2d)
434	ENDIF
435
436	IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
437	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
438	CALL histwrite_phy(nid_files(iff),
439	$ o_flw_srf(nsrf)%name,itau_w,
440	$ zx_tmp_fi2d)
441	ENDIF
442
443	IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
444	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
445	CALL histwrite_phy(nid_files(iff),
446	$ o_fsw_srf(nsrf)%name,itau_w,
447	$ zx_tmp_fi2d)
448	ENDIF
449
450	IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
451	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
452	CALL histwrite_phy(nid_files(iff),
453	$ o_wbils_srf(nsrf)%name,itau_w,
454	$ zx_tmp_fi2d)
455	ENDIF
456
457	IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
458	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
459	CALL histwrite_phy(nid_files(iff),
460	$ o_wbilo_srf(nsrf)%name,itau_w,
461	$ zx_tmp_fi2d)
462	ENDIF
463
464	if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
465	IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
466	CALL histwrite_phy(nid_files(iff),
467	$ o_tke_srf(nsrf)%name,itau_w,
468	$ pbl_tke(:,1:klev,nsrf))
469	ENDIF
470
471	IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
472	CALL histwrite_phy(nid_files(iff),
473	$ o_tke_max_srf(nsrf)%name,itau_w,
474	$ pbl_tke(:,1:klev,nsrf))
475	ENDIF
476	endif
477	ENDDO
478
479	IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN
480	CALL histwrite_phy(nid_files(iff),o_cdrm%name,itau_w,cdragm)
481	ENDIF
482
483	IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN
484	CALL histwrite_phy(nid_files(iff),o_cdrh%name,itau_w,cdragh)
485	ENDIF
486
487	IF (o_cldl%flag(iff)<=lev_files(iff)) THEN
488	CALL histwrite_phy(nid_files(iff),o_cldl%name,itau_w,cldl)
489	ENDIF
490
491	IF (o_cldm%flag(iff)<=lev_files(iff)) THEN
492	CALL histwrite_phy(nid_files(iff),o_cldm%name,itau_w,cldm)
493	ENDIF
494
495	IF (o_cldh%flag(iff)<=lev_files(iff)) THEN
496	CALL histwrite_phy(nid_files(iff),o_cldh%name,itau_w,cldh)
497	ENDIF
498
499	IF (o_cldt%flag(iff)<=lev_files(iff)) THEN
500	CALL histwrite_phy(nid_files(iff),o_cldt%name,
501	& itau_w,cldt)
502	ENDIF
503
504	IF (o_cldq%flag(iff)<=lev_files(iff)) THEN
505	CALL histwrite_phy(nid_files(iff),o_cldq%name,itau_w,cldq)
506	ENDIF
507
508	IF (o_lwp%flag(iff)<=lev_files(iff)) THEN
509	zx_tmp_fi2d(1:klon) = flwp(1:klon)
510	CALL histwrite_phy(nid_files(iff),
511	s o_lwp%name,itau_w,zx_tmp_fi2d)
512	ENDIF
513
514	IF (o_iwp%flag(iff)<=lev_files(iff)) THEN
515	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
516	CALL histwrite_phy(nid_files(iff),
517	s o_iwp%name,itau_w,zx_tmp_fi2d)
518	ENDIF
519
520	IF (o_ue%flag(iff)<=lev_files(iff)) THEN
521	CALL histwrite_phy(nid_files(iff),o_ue%name,itau_w,ue)
522	ENDIF
523
524	IF (o_ve%flag(iff)<=lev_files(iff)) THEN
525	CALL histwrite_phy(nid_files(iff),o_ve%name,itau_w,ve)
526	ENDIF
527
528	IF (o_uq%flag(iff)<=lev_files(iff)) THEN
529	CALL histwrite_phy(nid_files(iff),o_uq%name,itau_w,uq)
530	ENDIF
531
532	IF (o_vq%flag(iff)<=lev_files(iff)) THEN
533	CALL histwrite_phy(nid_files(iff),o_vq%name,itau_w,vq)
534	ENDIF
535
536	IF(iflag_con.GE.3) THEN ! sb
537	IF (o_cape%flag(iff)<=lev_files(iff)) THEN
538	CALL histwrite_phy(nid_files(iff),o_cape%name,itau_w,cape)
539	ENDIF
540
541	IF (o_pbase%flag(iff)<=lev_files(iff)) THEN
542	CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w,ema_pcb)
543	ENDIF
544
545	IF (o_ptop%flag(iff)<=lev_files(iff)) THEN
546	CALL histwrite_phy(nid_files(iff),o_ptop%name,itau_w,ema_pct)
547	ENDIF
548
549	IF (o_fbase%flag(iff)<=lev_files(iff)) THEN
550	CALL histwrite_phy(nid_files(iff),o_fbase%name,itau_w,ema_cbmf)
551	ENDIF
552
553	IF (o_plcl%flag(iff)<=lev_files(iff)) THEN
554	CALL histwrite_phy(nid_files(iff),o_plcl%name,itau_w,plcl)
555	ENDIF
556
557	IF (o_plfc%flag(iff)<=lev_files(iff)) THEN
558	CALL histwrite_phy(nid_files(iff),o_plfc%name,itau_w,plfc)
559	ENDIF
560
561	IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN
562	CALL histwrite_phy(nid_files(iff),o_wbeff%name,itau_w,wbeff)
563	ENDIF
564
565
566	IF (o_prw%flag(iff)<=lev_files(iff)) THEN
567	CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw)
568	ENDIF
569
570	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
571	CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)
572	ENDIF
573
574	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
575	CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)
576	ENDIF
577
578	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
579	CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)
580	ENDIF
581
582	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
583	CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)
584	ENDIF
585
586	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
587	CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)
588	ENDIF
589
590	IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN
591	zx_tmp_fi2d=float(itau_con)/float(itap)
592	CALL histwrite_phy(nid_files(iff),o_ftime_con%name,
593	s itau_w,zx_tmp_fi2d)
594	ENDIF
595
596	IF (o_mc%flag(iff)<=lev_files(iff)) THEN
597	if(iflag_thermals.gt.1)then
598	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
599	else
600	zx_tmp_fi3d=dnwd+dnwd0+upwd
601	endif
602	CALL histwrite_phy(nid_files(iff),o_mc%name,itau_w,zx_tmp_fi3d)
603	ENDIF
604
605	ENDIF !iflag_con .GE. 3
606
607	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
608	CALL histwrite_phy(nid_files(iff),o_s_pblh%name,itau_w,s_pblh)
609	ENDIF
610
611	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
612	CALL histwrite_phy(nid_files(iff),o_s_pblt%name,itau_w,s_pblt)
613	ENDIF
614
615	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
616	CALL histwrite_phy(nid_files(iff),o_s_lcl%name,itau_w,s_lcl)
617	ENDIF
618
619	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
620	CALL histwrite_phy(nid_files(iff),o_s_therm%name,itau_w,s_therm)
621	ENDIF
622
623	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
624	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
625	! CALL histwrite_phy(nid_files(iff),o_s_capCL%name,itau_w,s_capCL)
626	! ENDIF
627
628	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
629	! CALL histwrite_phy(nid_files(iff),o_s_oliqCL%name,itau_w,s_oliqCL)
630	! ENDIF
631
632	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
633	! CALL histwrite_phy(nid_files(iff),o_s_cteiCL%name,itau_w,s_cteiCL)
634	! ENDIF
635
636	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
637	! CALL histwrite_phy(nid_files(iff),o_s_trmb1%name,itau_w,s_trmb1)
638	! ENDIF
639
640	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
641	! CALL histwrite_phy(nid_files(iff),o_s_trmb2%name,itau_w,s_trmb2)
642	! ENDIF
643
644	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
645	! CALL histwrite_phy(nid_files(iff),o_s_trmb3%name,itau_w,s_trmb3)
646	! ENDIF
647
648	! Champs interpolles sur des niveaux de pression
649
650	ll=0
651	DO k=1, nlevSTD
652	! IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
653	! IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
654	bb2=clevSTD(k)
655	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
656	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
657	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
658
659	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
660	ll=ll+1
661	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
662	CALL histwrite_phy(nid_files(iff),o_uSTDlevs(ll)%name,
663	& itau_w,uwriteSTD(:,k,iff))
664	ENDIF
665
666	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
667	CALL histwrite_phy(nid_files(iff),o_vSTDlevs(ll)%name,
668	& itau_w,vwriteSTD(:,k,iff))
669	ENDIF
670
671	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
672	CALL histwrite_phy(nid_files(iff),o_wSTDlevs(ll)%name,
673	& itau_w,wwriteSTD(:,k,iff))
674	ENDIF
675
676	IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
677	CALL histwrite_phy(nid_files(iff),o_zSTDlevs(ll)%name,
678	& itau_w,phiwriteSTD(:,k,iff))
679	ENDIF
680
681	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
682	CALL histwrite_phy(nid_files(iff),o_qSTDlevs(ll)%name,
683	& itau_w, qwriteSTD(:,k,iff))
684	ENDIF
685
686	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
687	CALL histwrite_phy(nid_files(iff),o_tSTDlevs(ll)%name,
688	& itau_w, twriteSTD(:,k,iff))
689	ENDIF
690
691	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
692	ENDDO
693
694	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
695	DO i=1, klon
696	IF (pctsrf(i,is_oce).GT.epsfra.OR.
697	$ pctsrf(i,is_sic).GT.epsfra) THEN
698	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
699	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
700	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
701	ELSE
702	zx_tmp_fi2d(i) = 273.15
703	ENDIF
704	ENDDO
705	CALL histwrite_phy(nid_files(iff),
706	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
707	ENDIF
708
709	! Couplage convection-couche limite
710	IF (iflag_con.GE.3) THEN
711	IF (iflag_coupl>=1) THEN
712	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
713	CALL histwrite_phy(nid_files(iff),o_ale_bl%name,itau_w,ale_bl)
714	ENDIF
715	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
716	CALL histwrite_phy(nid_files(iff),o_alp_bl%name,itau_w,alp_bl)
717	ENDIF
718	ENDIF !iflag_coupl>=1
719	ENDIF !(iflag_con.GE.3)
720
721	! Wakes
722	IF (iflag_con.EQ.3) THEN
723	IF (iflag_wake>=1) THEN
724	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
725	CALL histwrite_phy(nid_files(iff),o_ale_wk%name,itau_w,ale_wake)
726	ENDIF
727	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
728	CALL histwrite_phy(nid_files(iff),o_alp_wk%name,itau_w,alp_wake)
729	ENDIF
730
731	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
732	CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)
733	ENDIF
734	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
735	CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)
736	ENDIF
737	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
738	CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)
739	ENDIF
740	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
741	CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe)
742	ENDIF
743	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
744	CALL histwrite_phy(nid_files(iff),o_wake_h%name,itau_w,wake_h)
745	ENDIF
746
747	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
748	CALL histwrite_phy(nid_files(iff),o_wake_s%name,itau_w,wake_s)
749	ENDIF
750
751	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
752	CALL histwrite_phy(nid_files(iff),o_wake_deltat%name,
753	$ itau_w,wake_deltat)
754	ENDIF
755
756	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
757	CALL histwrite_phy(nid_files(iff),o_wake_deltaq%name,
758	$ itau_w,wake_deltaq)
759	ENDIF
760
761	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
762	CALL histwrite_phy(nid_files(iff),
763	s o_wake_omg%name,itau_w,wake_omg)
764	ENDIF
765
766	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
767	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
768	& /pdtphys
769	CALL histwrite_phy(nid_files(iff),
770	& o_dtwak%name,itau_w,zx_tmp_fi3d)
771	ENDIF
772
773	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
774	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
775	CALL histwrite_phy(nid_files(iff),
776	& o_dqwak%name,itau_w,zx_tmp_fi3d)
777	ENDIF
778	ENDIF ! iflag_wake>=1
779
780	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
781	CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip)
782	ENDIF
783
784	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
785	CALL histwrite_phy(nid_files(iff),o_ftd%name,itau_w,ftd)
786	ENDIF
787
788	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
789	CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd)
790	ENDIF
791	ENDIF !(iflag_con.EQ.3)
792
793	IF (type_ocean=='slab ') THEN
794	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
795	$ CALL histwrite_phy(
796	$ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils)
797
798	ENDIF !type_ocean == force/slab
799
800	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
801	CALL histwrite_phy(nid_files(iff),
802	s o_weakinv%name,itau_w,weak_inversion)
803	ENDIF
804
805	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
806	CALL histwrite_phy(nid_files(iff),o_dthmin%name,itau_w,dthmin)
807	ENDIF
808
809	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
810	CALL histwrite_phy(nid_files(iff),o_cldtau%name,itau_w,cldtau)
811	ENDIF
812
813	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
814	CALL histwrite_phy(nid_files(iff),o_cldemi%name,itau_w,cldemi)
815	ENDIF
816
817	IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
818	CALL histwrite_phy(nid_files(iff),
819	s o_pr_con_l%name,itau_w,pmflxr(:,1:klev))
820	ENDIF
821
822	IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
823	CALL histwrite_phy(nid_files(iff),
824	s o_pr_con_i%name,itau_w,pmflxs(:,1:klev))
825	ENDIF
826
827	IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
828	CALL histwrite_phy(nid_files(iff),
829	s o_pr_lsc_l%name,itau_w,prfl(:,1:klev))
830	ENDIF
831
832	IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
833	CALL histwrite_phy(nid_files(iff),
834	s o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
835	ENDIF
836
837	IF (o_re%flag(iff)<=lev_files(iff)) THEN
838	CALL histwrite_phy(nid_files(iff),o_re%name,itau_w,re)
839	ENDIF
840
841	IF (o_fl%flag(iff)<=lev_files(iff)) THEN
842	CALL histwrite_phy(nid_files(iff),o_fl%name,itau_w,fl)
843	ENDIF
844
845
846
847	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
848	DO i=1, klon
849	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
850	ENDDO
851	CALL histwrite_phy(nid_files(iff),o_rh2m%name,itau_w,zx_tmp_fi2d)
852	ENDIF
853
854	IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN
855	DO i=1, klon
856	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
857	ENDDO
858	CALL histwrite_phy(nid_files(iff),o_rh2m_min%name,
859	s itau_w,zx_tmp_fi2d)
860	ENDIF
861
862	IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN
863	DO i=1, klon
864	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
865	ENDDO
866	CALL histwrite_phy(nid_files(iff),o_rh2m_max%name,
867	s itau_w,zx_tmp_fi2d)
868	ENDIF
869
870
871	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
872	CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m)
873	ENDIF
874
875	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
876	CALL histwrite_phy(nid_files(iff),o_tpot%name,itau_w,tpot)
877	ENDIF
878
879	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
880	CALL histwrite_phy(nid_files(iff),o_tpote%name,itau_w,tpote)
881	ENDIF
882
883	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
884	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
885	CALL histwrite_phy(nid_files(iff),
886	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
887	ENDIF
888
889	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
890	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
891	CALL histwrite_phy(nid_files(iff),
892	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
893	ENDIF
894
895	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
896	CALL histwrite_phy(nid_files(iff),
897	s o_LWdownOR%name,itau_w,sollwdown)
898	ENDIF
899
900	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
901	CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc)
902	ENDIF
903
904	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
905	CALL histwrite_phy(nid_files(iff),
906	s o_solldown%name,itau_w,sollwdown)
907	ENDIF
908
909	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
910	CALL histwrite_phy(nid_files(iff),
911	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
912	ENDIF
913
914	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
915	CALL histwrite_phy(nid_files(iff),
916	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
917	ENDIF
918
919	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
920	CALL histwrite_phy(nid_files(iff),
921	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
922	ENDIF
923
924	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
925	CALL histwrite_phy(nid_files(iff),
926	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
927	ENDIF
928
929	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
930	CALL histwrite_phy(nid_files(iff),o_rugs%name,itau_w,zxrugs)
931	ENDIF
932
933	! OD550 per species
934	IF (new_aod .and. (.not. aerosol_couple)) THEN
935	IF (ok_ade.OR.ok_aie) THEN
936
937	IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN
938	CALL histwrite_phy(nid_files(iff),o_od550aer%name,itau_w,
939	$ od550aer)
940	ENDIF
941	IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN
942	CALL histwrite_phy(nid_files(iff),o_od865aer%name,itau_w,
943	$ od865aer)
944	ENDIF
945	IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN
946	CALL histwrite_phy(nid_files(iff),o_absvisaer%name,itau_w,
947	$ absvisaer)
948	ENDIF
949	IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN
950	CALL histwrite_phy(nid_files(iff),o_od550lt1aer%name,itau_w,
951	$ od550lt1aer)
952	ENDIF
953
954	IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN
955	CALL histwrite_phy(nid_files(iff),o_sconcso4%name,itau_w,
956	$ sconcso4)
957	ENDIF
958	IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN
959	CALL histwrite_phy(nid_files(iff),o_sconcoa%name,itau_w,
960	$ sconcoa)
961	ENDIF
962	IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN
963	CALL histwrite_phy(nid_files(iff),o_sconcbc%name,itau_w,
964	$ sconcbc)
965	ENDIF
966	IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN
967	CALL histwrite_phy(nid_files(iff),o_sconcss%name,itau_w,
968	$ sconcss)
969	ENDIF
970	IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN
971	CALL histwrite_phy(nid_files(iff),o_sconcdust%name,itau_w,
972	$ sconcdust)
973	ENDIF
974
975	IF (o_concso4%flag(iff)<=lev_files(iff)) THEN
976	CALL histwrite_phy(nid_files(iff),o_concso4%name,itau_w,
977	$ concso4)
978	ENDIF
979	IF (o_concoa%flag(iff)<=lev_files(iff)) THEN
980	CALL histwrite_phy(nid_files(iff),o_concoa%name,itau_w,
981	$ concoa)
982	ENDIF
983	IF (o_concbc%flag(iff)<=lev_files(iff)) THEN
984	CALL histwrite_phy(nid_files(iff),o_concbc%name,itau_w,
985	$ concbc)
986	ENDIF
987	IF (o_concss%flag(iff)<=lev_files(iff)) THEN
988	CALL histwrite_phy(nid_files(iff),o_concss%name,itau_w,
989	$ concss)
990	ENDIF
991	IF (o_concdust%flag(iff)<=lev_files(iff)) THEN
992	CALL histwrite_phy(nid_files(iff),o_concdust%name,itau_w,
993	$ concdust)
994	ENDIF
995
996	IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN
997	CALL histwrite_phy(nid_files(iff),o_loadso4%name,itau_w,
998	$ loadso4)
999	ENDIF
1000	IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN
1001	CALL histwrite_phy(nid_files(iff),o_loadoa%name,itau_w,
1002	$ loadoa)
1003	ENDIF
1004	IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN
1005	CALL histwrite_phy(nid_files(iff),o_loadbc%name,itau_w,
1006	$ loadbc)
1007	ENDIF
1008	IF (o_loadss%flag(iff)<=lev_files(iff)) THEN
1009	CALL histwrite_phy(nid_files(iff),o_loadss%name,itau_w,
1010	$ loadss)
1011	ENDIF
1012	IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN
1013	CALL histwrite_phy(nid_files(iff),o_loaddust%name,itau_w,
1014	$ loaddust)
1015	ENDIF
1016
1017	DO naero = 1, naero_spc
1018	IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
1019	CALL histwrite_phy(nid_files(iff),
1020	$ o_tausumaero(naero)%name,itau_w,
1021	$ tausum_aero(:,2,naero) )
1022	ENDIF
1023	END DO
1024	endif
1025	ENDIF
1026
1027	IF (ok_ade) THEN
1028	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
1029	CALL histwrite_phy(nid_files(iff),o_topswad%name,itau_w,
1030	$ topswad_aero)
1031	ENDIF
1032	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
1033	CALL histwrite_phy(nid_files(iff),o_solswad%name,itau_w,
1034	$ solswad_aero)
1035	ENDIF
1036
1037	!====MS forcing diagnostics
1038	if (new_aod) then
1039	IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
1040	CALL histwrite_phy(nid_files(iff),o_swtoaas_nat%name,itau_w,
1041	$ topsw_aero(:,1))
1042	ENDIF
1043
1044	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
1045	CALL histwrite_phy(nid_files(iff),o_swsrfas_nat%name,itau_w,
1046	$ solsw_aero(:,1))
1047	ENDIF
1048
1049	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
1050	CALL histwrite_phy(nid_files(iff),o_swtoacs_nat%name,itau_w,
1051	$ topsw0_aero(:,1))
1052	ENDIF
1053
1054	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
1055	CALL histwrite_phy(nid_files(iff),o_swsrfcs_nat%name,itau_w,
1056	$ solsw0_aero(:,1))
1057	ENDIF
1058
1059	!ant
1060	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
1061	CALL histwrite_phy(nid_files(iff),o_swtoaas_ant%name,itau_w,
1062	$ topsw_aero(:,2))
1063	ENDIF
1064
1065	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
1066	CALL histwrite_phy(nid_files(iff),o_swsrfas_ant%name,itau_w,
1067	$ solsw_aero(:,2))
1068	ENDIF
1069
1070	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
1071	CALL histwrite_phy(nid_files(iff),o_swtoacs_ant%name,itau_w,
1072	$ topsw0_aero(:,2))
1073	ENDIF
1074
1075	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
1076	CALL histwrite_phy(nid_files(iff),o_swsrfcs_ant%name,itau_w,
1077	$ solsw0_aero(:,2))
1078	ENDIF
1079
1080	!cf
1081
1082	if (.not. aerosol_couple) then
1083	IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
1084	CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
1085	$ topswcf_aero(:,1))
1086	ENDIF
1087
1088	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
1089	CALL histwrite_phy(nid_files(iff),o_swsrfcf_nat%name,itau_w,
1090	$ solswcf_aero(:,1))
1091	ENDIF
1092
1093	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
1094	CALL histwrite_phy(nid_files(iff),o_swtoacf_ant%name,itau_w,
1095	$ topswcf_aero(:,2))
1096	ENDIF
1097
1098	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
1099	CALL histwrite_phy(nid_files(iff),o_swsrfcf_ant%name,itau_w,
1100	$ solswcf_aero(:,2))
1101	ENDIF
1102
1103	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
1104	CALL histwrite_phy(nid_files(iff),o_swtoacf_zero%name,itau_w,
1105	$ topswcf_aero(:,3))
1106	ENDIF
1107
1108	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
1109	CALL histwrite_phy(nid_files(iff),o_swsrfcf_zero%name,itau_w,
1110	$ solswcf_aero(:,3))
1111	ENDIF
1112	endif
1113
1114	endif ! new_aod
1115	!====MS forcing diagnostics
1116
1117	ENDIF
1118
1119	IF (ok_aie) THEN
1120	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
1121	CALL histwrite_phy(nid_files(iff),o_topswai%name,itau_w,
1122	$ topswai_aero)
1123	ENDIF
1124	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
1125	CALL histwrite_phy(nid_files(iff),o_solswai%name,itau_w,
1126	$ solswai_aero)
1127	ENDIF
1128	IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN
1129	CALL histwrite_phy(nid_files(iff),o_scdnc%name,itau_w,
1130	$ scdnc)
1131	ENDIF
1132	IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN
1133	CALL histwrite_phy(nid_files(iff),o_cldncl%name,itau_w,
1134	$ cldncl)
1135	ENDIF
1136	IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN
1137	CALL histwrite_phy(nid_files(iff),o_reffclws%name,itau_w,
1138	$ reffclws)
1139	ENDIF
1140	IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN
1141	CALL histwrite_phy(nid_files(iff),o_reffclwc%name,itau_w,
1142	$ reffclwc)
1143	ENDIF
1144	IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN
1145	CALL histwrite_phy(nid_files(iff),o_cldnvi%name,itau_w,
1146	$ cldnvi)
1147	ENDIF
1148	IF (o_lcc%flag(iff)<=lev_files(iff)) THEN
1149	CALL histwrite_phy(nid_files(iff),o_lcc%name,itau_w,
1150	$ lcc)
1151	ENDIF
1152	IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN
1153	CALL histwrite_phy(nid_files(iff),o_lcc3d%name,itau_w,
1154	$ lcc3d)
1155	ENDIF
1156	IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN
1157	CALL histwrite_phy(nid_files(iff),o_lcc3dcon%name,itau_w,
1158	$ lcc3dcon)
1159	ENDIF
1160	IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN
1161	CALL histwrite_phy(nid_files(iff),o_lcc3dstra%name,itau_w,
1162	$ lcc3dstra)
1163	ENDIF
1164	IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN
1165	CALL histwrite_phy(nid_files(iff),o_reffclwtop%name,itau_w,
1166	$ reffclwtop)
1167	ENDIF
1168	ENDIF
1169
1170	! Champs 3D:
1171	IF (ok_ade .OR. ok_aie) then
1172	IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN
1173	CALL histwrite_phy(nid_files(iff),o_ec550aer%name,itau_w,
1174	& ec550aer)
1175	ENDIF
1176	ENDIF
1177
1178	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
1179	CALL histwrite_phy(nid_files(iff),o_lwcon%name,itau_w,flwc)
1180	ENDIF
1181
1182	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
1183	CALL histwrite_phy(nid_files(iff),o_iwcon%name,itau_w,fiwc)
1184	ENDIF
1185
1186	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
1187	CALL histwrite_phy(nid_files(iff),o_temp%name,itau_w,t_seri)
1188	ENDIF
1189
1190	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
1191	CALL histwrite_phy(nid_files(iff),o_theta%name,itau_w,theta)
1192	ENDIF
1193
1194	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
1195	CALL histwrite_phy(nid_files(iff),o_ovapinit%name,itau_w,
1196	$ qx(:,:,ivap))
1197	ENDIF
1198
1199	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
1200	CALL histwrite_phy(nid_files(iff),
1201	$ o_ovap%name,itau_w,q_seri)
1202	ENDIF
1203
1204	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
1205	CALL histwrite_phy(nid_files(iff),o_geop%name,itau_w,zphi)
1206	ENDIF
1207
1208	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
1209	CALL histwrite_phy(nid_files(iff),o_vitu%name,itau_w,u_seri)
1210	ENDIF
1211
1212	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
1213	CALL histwrite_phy(nid_files(iff),o_vitv%name,itau_w,v_seri)
1214	ENDIF
1215
1216	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
1217	CALL histwrite_phy(nid_files(iff),o_vitw%name,itau_w,omega)
1218	ENDIF
1219
1220	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
1221	CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay)
1222	ENDIF
1223
1224	IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
1225	CALL histwrite_phy(nid_files(iff),o_paprs%name,
1226	s itau_w,paprs(:,1:klev))
1227	ENDIF
1228
1229	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
1230	CALL histwrite_phy(nid_files(iff),o_rneb%name,itau_w,cldfra)
1231	ENDIF
1232
1233	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
1234	CALL histwrite_phy(nid_files(iff),o_rnebcon%name,itau_w,rnebcon)
1235	ENDIF
1236
1237	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
1238	CALL histwrite_phy(nid_files(iff),o_rhum%name,itau_w,zx_rh)
1239	ENDIF
1240
1241	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
1242	CALL histwrite_phy(nid_files(iff), o_ozone%name, itau_w,
1243	$ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1244	ENDIF
1245
1246	IF (o_ozone_light%flag(iff)<=lev_files(iff) .and.
1247	$ read_climoz == 2) THEN
1248	CALL histwrite_phy(nid_files(iff), o_ozone_light%name, itau_w,
1249	$ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1250	ENDIF
1251
1252	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
1253	CALL histwrite_phy(nid_files(iff),o_dtphy%name,itau_w,d_t)
1254	ENDIF
1255
1256	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
1257	CALL histwrite_phy(nid_files(iff),
1258	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
1259	ENDIF
1260
1261	DO nsrf=1, nbsrf
1262	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1263	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1264	CALL histwrite_phy(nid_files(iff),
1265	s o_albe_srf(nsrf)%name,itau_w,
1266	$ zx_tmp_fi2d)
1267	ENDIF
1268
1269	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1270	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1271	CALL histwrite_phy(nid_files(iff),
1272	s o_rugs_srf(nsrf)%name,itau_w,
1273	$ zx_tmp_fi2d)
1274	ENDIF
1275
1276	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1277	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1278	CALL histwrite_phy(nid_files(iff),
1279	s o_ages_srf(nsrf)%name,itau_w
1280	$ ,zx_tmp_fi2d)
1281	ENDIF
1282	ENDDO !nsrf=1, nbsrf
1283
1284	IF (o_alb1%flag(iff)<=lev_files(iff)) THEN
1285	CALL histwrite_phy(nid_files(iff),o_alb1%name,itau_w,albsol1)
1286	ENDIF
1287
1288	IF (o_alb2%flag(iff)<=lev_files(iff)) THEN
1289	CALL histwrite_phy(nid_files(iff),o_alb2%name,itau_w,albsol2)
1290	ENDIF
1291
1292	!FH Sorties pour la couche limite
1293	if (iflag_pbl>1) then
1294	zx_tmp_fi3d=0.
1295	do nsrf=1,nbsrf
1296	do k=1,klev
1297	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
1298	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1299	enddo
1300	enddo
1301	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
1302	CALL histwrite_phy(nid_files(iff),o_tke%name,itau_w,zx_tmp_fi3d)
1303	ENDIF
1304
1305	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
1306	CALL histwrite_phy(nid_files(iff),
1307	s o_tke_max%name,itau_w,zx_tmp_fi3d)
1308	ENDIF
1309	endif
1310
1311	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
1312	CALL histwrite_phy(nid_files(iff),o_kz%name,itau_w,coefh)
1313	ENDIF
1314
1315	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
1316	CALL histwrite_phy(nid_files(iff),o_kz_max%name,itau_w,coefh)
1317	ENDIF
1318
1319	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
1320	CALL histwrite_phy(nid_files(iff),o_clwcon%name,itau_w,clwcon0)
1321	ENDIF
1322
1323	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
1324	CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn)
1325	ENDIF
1326
1327	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
1328	CALL histwrite_phy(nid_files(iff),o_dqdyn%name,itau_w,d_q_dyn)
1329	ENDIF
1330
1331	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
1332	CALL histwrite_phy(nid_files(iff),o_dudyn%name,itau_w,d_u_dyn)
1333	ENDIF
1334
1335	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
1336	CALL histwrite_phy(nid_files(iff),o_dvdyn%name,itau_w,d_v_dyn)
1337	ENDIF
1338
1339	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
1340	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1341	CALL histwrite_phy(nid_files(iff),o_dtcon%name,itau_w,zx_tmp_fi3d)
1342	ENDIF
1343
1344	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
1345	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1346	CALL histwrite_phy(nid_files(iff),o_ducon%name,itau_w,zx_tmp_fi3d)
1347	ENDIF
1348
1349	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
1350	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1351	CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d)
1352	ENDIF
1353
1354	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
1355	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1356	CALL histwrite_phy(nid_files(iff),o_dtlsc%name,itau_w,zx_tmp_fi3d)
1357	ENDIF
1358
1359	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
1360	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1361	$ d_t_eva(1:klon,1:klev))/pdtphys
1362	CALL histwrite_phy(nid_files(iff),
1363	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
1364	ENDIF
1365
1366	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
1367	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1368	CALL histwrite_phy(nid_files(iff),o_dqlsc%name,itau_w,zx_tmp_fi3d)
1369	ENDIF
1370
1371	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1372	! Sorties specifiques a la separation thermiques/non thermiques
1373	if (iflag_thermals>1) then
1374
1375	IF (o_dtlscth%flag(iff)<=lev_files(iff)) THEN
1376	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1377	CALL
1378	s histwrite_phy(nid_files(iff),o_dtlscth%name,itau_w,zx_tmp_fi3d)
1379	ENDIF
1380
1381	IF (o_dtlscst%flag(iff)<=lev_files(iff)) THEN
1382	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1383	CALL
1384	s histwrite_phy(nid_files(iff),o_dtlscst%name,itau_w,zx_tmp_fi3d)
1385	ENDIF
1386
1387	IF (o_dqlscth%flag(iff)<=lev_files(iff)) THEN
1388	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1389	CALL
1390	s histwrite_phy(nid_files(iff),o_dqlscth%name,itau_w,zx_tmp_fi3d)
1391	ENDIF
1392
1393	IF (o_dqlscst%flag(iff)<=lev_files(iff)) THEN
1394	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1395	CALL
1396	s histwrite_phy(nid_files(iff),o_dqlscst%name,itau_w,zx_tmp_fi3d)
1397	ENDIF
1398
1399	IF (o_plulth%flag(iff)<=lev_files(iff)) THEN
1400	CALL histwrite_phy(nid_files(iff),o_plulth%name,itau_w,plul_th)
1401	ENDIF
1402
1403	IF (o_plulst%flag(iff)<=lev_files(iff)) THEN
1404	CALL histwrite_phy(nid_files(iff),o_plulst%name,itau_w,plul_st)
1405	ENDIF
1406
1407	do k=1,klev
1408	do i=1,klon
1409	if (ptconvth(i,k)) then
1410	zx_tmp_fi3d(i,k)=1.
1411	else
1412	zx_tmp_fi3d(i,k)=0.
1413	endif
1414	enddo
1415	enddo
1416	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1417	CALL
1418	s histwrite_phy(nid_files(iff),o_ptconvth%name,itau_w,zx_tmp_fi3d)
1419	ENDIF
1420
1421	do i=1,klon
1422	zx_tmp_fi2d(1:klon)=lmax_th(:)
1423	enddo
1424	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1425	CALL
1426	s histwrite_phy(nid_files(iff),o_lmaxth%name,itau_w,zx_tmp_fi2d)
1427	ENDIF
1428
1429	endif ! iflag_thermals>1
1430
1431	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1432
1433	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
1434	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1435	CALL histwrite_phy(nid_files(iff),o_dtvdf%name,itau_w,zx_tmp_fi3d)
1436	ENDIF
1437
1438	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
1439	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1440	CALL histwrite_phy(nid_files(iff),o_dqvdf%name,itau_w,zx_tmp_fi3d)
1441	ENDIF
1442
1443	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
1444	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1445	CALL histwrite_phy(nid_files(iff),o_dteva%name,itau_w,zx_tmp_fi3d)
1446	ENDIF
1447
1448	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
1449	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1450	CALL histwrite_phy(nid_files(iff),o_dqeva%name,itau_w,zx_tmp_fi3d)
1451	ENDIF
1452
1453	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
1454	zpt_conv = 0.
1455	where (ptconv) zpt_conv = 1.
1456	CALL histwrite_phy(nid_files(iff),o_ptconv%name,itau_w,zpt_conv)
1457	ENDIF
1458
1459	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
1460	CALL histwrite_phy(nid_files(iff),o_ratqs%name,itau_w,ratqs)
1461	ENDIF
1462
1463	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
1464	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
1465	CALL histwrite_phy(nid_files(iff),o_dtthe%name,itau_w,zx_tmp_fi3d)
1466	ENDIF
1467
1468	IF (iflag_thermals.gt.1) THEN
1469	IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN
1470	! Pour l instant 0 a y reflichir pour les thermiques
1471	zx_tmp_fi2d=0.
1472	CALL histwrite_phy(nid_files(iff),o_ftime_th%name,
1473	s itau_w,zx_tmp_fi2d)
1474	ENDIF
1475
1476	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
1477	CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm)
1478	ENDIF
1479
1480	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
1481	CALL histwrite_phy(nid_files(iff),o_e_th%name,itau_w,entr_therm)
1482	ENDIF
1483
1484	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
1485	CALL histwrite_phy(nid_files(iff),o_w_th%name,itau_w,zw2)
1486	ENDIF
1487
1488	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
1489	CALL histwrite_phy(nid_files(iff),o_q_th%name,itau_w,zqasc)
1490	ENDIF
1491
1492	IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN
1493	CALL histwrite_phy(nid_files(iff),
1494	s o_lambda_th%name,itau_w,lambda_th)
1495	ENDIF
1496
1497	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
1498	CALL histwrite_phy(nid_files(iff),o_a_th%name,itau_w,fraca)
1499	ENDIF
1500
1501	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
1502	CALL histwrite_phy(nid_files(iff),o_d_th%name,itau_w,detr_therm)
1503	ENDIF
1504
1505
1506	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1507	CALL histwrite_phy(nid_files(iff),o_f0_th%name,itau_w,f0)
1508	ENDIF
1509
1510	IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN
1511	CALL histwrite_phy(nid_files(iff),o_zmax_th%name,itau_w,zmax_th)
1512	ENDIF
1513
1514	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
1515	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
1516	CALL histwrite_phy(nid_files(iff),o_dqthe%name,itau_w,zx_tmp_fi3d)
1517	ENDIF
1518
1519	ENDIF !iflag_thermals
1520
1521	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
1522	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1523	CALL histwrite_phy(nid_files(iff),o_dtajs%name,itau_w,zx_tmp_fi3d)
1524	ENDIF
1525
1526	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
1527	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1528	CALL histwrite_phy(nid_files(iff),o_dqajs%name,itau_w,zx_tmp_fi3d)
1529	ENDIF
1530
1531	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
1532	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1533	CALL histwrite_phy(nid_files(iff),o_dtswr%name,itau_w,zx_tmp_fi3d)
1534	ENDIF
1535
1536	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
1537	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1538	CALL histwrite_phy(nid_files(iff),o_dtsw0%name,itau_w,zx_tmp_fi3d)
1539	ENDIF
1540
1541	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
1542	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1543	CALL histwrite_phy(nid_files(iff),o_dtlwr%name,itau_w,zx_tmp_fi3d)
1544	ENDIF
1545
1546	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
1547	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1548	CALL histwrite_phy(nid_files(iff),o_dtlw0%name,itau_w,zx_tmp_fi3d)
1549	ENDIF
1550
1551	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
1552	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1553	CALL histwrite_phy(nid_files(iff),o_dtec%name,itau_w,zx_tmp_fi3d)
1554	ENDIF
1555
1556	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
1557	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1558	CALL histwrite_phy(nid_files(iff),o_duvdf%name,itau_w,zx_tmp_fi3d)
1559	ENDIF
1560
1561	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
1562	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1563	CALL histwrite_phy(nid_files(iff),o_dvvdf%name,itau_w,zx_tmp_fi3d)
1564	ENDIF
1565
1566	IF (ok_orodr) THEN
1567	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
1568	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1569	CALL histwrite_phy(nid_files(iff),o_duoro%name,itau_w,zx_tmp_fi3d)
1570	ENDIF
1571
1572	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
1573	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1574	CALL histwrite_phy(nid_files(iff),o_dvoro%name,itau_w,zx_tmp_fi3d)
1575	ENDIF
1576	ENDIF
1577
1578	IF (ok_orolf) THEN
1579	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
1580	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1581	CALL histwrite_phy(nid_files(iff),o_dulif%name,itau_w,zx_tmp_fi3d)
1582	ENDIF
1583
1584	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
1585	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1586	CALL histwrite_phy(nid_files(iff),o_dvlif%name,itau_w,zx_tmp_fi3d)
1587	ENDIF
1588	ENDIF
1589
1590	if (nqtot.GE.3) THEN
1591	DO iq=3,nqtot
1592	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
1593	CALL histwrite_phy(nid_files(iff),
1594	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
1595	ENDIF
1596	ENDDO
1597	endif
1598
1599	if (ok_sync) then
1600	c$OMP MASTER
1601	call histsync(nid_files(iff))
1602	c$OMP END MASTER
1603	endif
1604
1605	ENDIF ! clef_files
1606
1607	ENDDO ! iff=1,nfiles

Note: See TracBrowser for help on using the repository browser.

Download in other formats: