Context Navigation

source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_write.h @ 1405

Last change on this file since 1405 was 1405, checked in by jghattas, 14 years ago

Added tracer "Age of stratospheric air", activated by putting Aga in tracer.def.
The logical rnpb is now controled during execution. rnpb is true only if both tracers RN and PB existe in tracer.def. RN and PB can now be removed from tracer.def
Modification in output. It is now possible to have 3 tracers. But still only tracer number 3 and tracer number 4(if existing) will be written.

F Lott, JG

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

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

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

Download in other formats: