Context Navigation

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

Last change on this file since 1633 was 1633, checked in by lguez, 12 years ago
"prw" can now be output even if "iflag_con < 3".
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 72.5 KB

Line
1	itau_w = itau_phy + itap + start_time * day_step / iphysiq
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),clef_stations(iff),
12	$ o_phis%name,itau_w,pphis)
13	ENDIF
14
15	IF (.NOT.clef_stations(iff)) THEN
16	IF (o_aire%flag(iff)<=lev_files(iff)) THEN
17	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
18	$o_aire%name,itau_w,airephy)
19	ENDIF
20
21	IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN
22	DO i=1, klon
23	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
24	ENDDO
25	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
26	$ o_contfracATM%name,itau_w,zx_tmp_fi2d)
27	ENDIF
28	ENDIF
29
30	IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN
31	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
32	$ o_contfracOR%name,itau_w,pctsrf(:,is_ter))
33	ENDIF
34
35	IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN
36	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
37	$ o_aireTER%name,itau_w,paire_ter)
38	ENDIF
39
40	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
41
42	IF (o_flat%flag(iff)<=lev_files(iff)) THEN
43	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
44	$ o_flat%name,itau_w,zxfluxlat)
45	ENDIF
46
47	IF (o_slp%flag(iff)<=lev_files(iff)) THEN
48	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
49	$ o_slp%name,itau_w,slp)
50	ENDIF
51
52	IF (o_tsol%flag(iff)<=lev_files(iff)) THEN
53	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
54	$ o_tsol%name,itau_w,zxtsol)
55	ENDIF
56
57	IF (o_t2m%flag(iff)<=lev_files(iff)) THEN
58	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
59	$ o_t2m%name,itau_w,zt2m)
60	ENDIF
61
62	IF (.NOT.clef_stations(iff)) THEN
63	IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN
64	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
65	$o_t2m_min%name,itau_w,zt2m)
66	ENDIF
67
68	IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN
69	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
70	$o_t2m_max%name,itau_w,zt2m)
71	ENDIF
72	ENDIF
73
74	IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN
75	DO i=1, klon
76	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
77	ENDDO
78	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
79	s o_wind10m%name,itau_w,zx_tmp_fi2d)
80	ENDIF
81
82	IF (.NOT.clef_stations(iff)) THEN
83	IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN
84	DO i=1, klon
85	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
86	ENDDO
87	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
88	$o_wind10max%name,itau_w,zx_tmp_fi2d)
89	ENDIF
90	ENDIF
91
92	IF (o_sicf%flag(iff)<=lev_files(iff)) THEN
93	DO i = 1, klon
94	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
95	ENDDO
96	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
97	$ o_sicf%name,itau_w,zx_tmp_fi2d)
98	ENDIF
99
100	IF (o_q2m%flag(iff)<=lev_files(iff)) THEN
101	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
102	$o_q2m%name,itau_w,zq2m)
103	ENDIF
104
105	IF (o_u10m%flag(iff)<=lev_files(iff)) THEN
106	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
107	$o_u10m%name,itau_w,zu10m)
108	ENDIF
109
110	IF (o_v10m%flag(iff)<=lev_files(iff)) THEN
111	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
112	$o_v10m%name,itau_w,zv10m)
113	ENDIF
114
115	IF (o_psol%flag(iff)<=lev_files(iff)) THEN
116	DO i = 1, klon
117	zx_tmp_fi2d(i) = paprs(i,1)
118	ENDDO
119	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
120	s o_psol%name,itau_w,zx_tmp_fi2d)
121	ENDIF
122
123	IF (o_mass%flag(iff)<=lev_files(iff)) THEN
124	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
125	s o_mass%name,itau_w,zmasse)
126	ENDIF
127
128
129	IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN
130	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
131	$o_qsurf%name,itau_w,zxqsurf)
132	ENDIF
133
134	if (.not. ok_veget) then
135	IF (o_qsol%flag(iff)<=lev_files(iff)) THEN
136	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
137	$o_qsol%name,itau_w,qsol)
138	ENDIF
139	endif
140
141	IF (o_precip%flag(iff)<=lev_files(iff)) THEN
142	DO i = 1, klon
143	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
144	ENDDO
145	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
146	$o_precip%name,itau_w,zx_tmp_fi2d)
147	ENDIF
148
149	IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN
150	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
151	$o_ndayrain%name,itau_w,nday_rain)
152	ENDIF
153
154	IF (o_plul%flag(iff)<=lev_files(iff)) THEN
155	DO i = 1, klon
156	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
157	ENDDO
158	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
159	$o_plul%name,itau_w,zx_tmp_fi2d)
160	ENDIF
161
162	IF (o_pluc%flag(iff)<=lev_files(iff)) THEN
163	DO i = 1, klon
164	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
165	ENDDO
166	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
167	$o_pluc%name,itau_w,zx_tmp_fi2d)
168	ENDIF
169
170	IF (o_snow%flag(iff)<=lev_files(iff)) THEN
171	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
172	$o_snow%name,itau_w,snow_fall)
173	ENDIF
174
175	IF (o_msnow%flag(iff)<=lev_files(iff)) THEN
176	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
177	$o_msnow%name,itau_w,snow_o)
178	ENDIF
179
180	IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN
181	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
182	$o_fsnow%name,itau_w,zfra_o)
183	ENDIF
184
185	IF (o_evap%flag(iff)<=lev_files(iff)) THEN
186	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
187	$o_evap%name,itau_w,evap)
188	ENDIF
189
190	IF (o_tops%flag(iff)<=lev_files(iff)) THEN
191	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
192	$o_tops%name,itau_w,topsw)
193	ENDIF
194
195	IF (o_tops0%flag(iff)<=lev_files(iff)) THEN
196	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
197	$o_tops0%name,itau_w,topsw0)
198	ENDIF
199
200	IF (o_topl%flag(iff)<=lev_files(iff)) THEN
201	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
202	$o_topl%name,itau_w,toplw)
203	ENDIF
204
205	IF (o_topl0%flag(iff)<=lev_files(iff)) THEN
206	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
207	$o_topl0%name,itau_w,toplw0)
208	ENDIF
209
210	IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN
211	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
212	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
213	$o_SWupTOA%name,itau_w,zx_tmp_fi2d)
214	ENDIF
215
216	IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN
217	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
218	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
219	$ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d)
220	ENDIF
221
222	IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN
223	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
224	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
225	s o_SWdnTOA%name,itau_w,zx_tmp_fi2d)
226	ENDIF
227
228	IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN
229	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
230	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
231	$ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d)
232	ENDIF
233
234	IF (o_nettop%flag(iff)<=lev_files(iff)) THEN
235	zx_tmp_fi2d(:) = topsw(:)-toplw(:)
236	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
237	$ o_nettop%name,itau_w,zx_tmp_fi2d)
238	ENDIF
239
240	IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
241	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
242	$o_SWup200%name,itau_w,SWup200)
243	ENDIF
244
245	IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN
246	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
247	s o_SWup200clr%name,itau_w,SWup200clr)
248	ENDIF
249
250	IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN
251	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
252	$o_SWdn200%name,itau_w,SWdn200)
253	ENDIF
254
255	IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN
256	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
257	s o_SWdn200clr%name,itau_w,SWdn200clr)
258	ENDIF
259
260	IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN
261	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
262	$o_LWup200%name,itau_w,LWup200)
263	ENDIF
264
265	IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN
266	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
267	s o_LWup200clr%name,itau_w,LWup200clr)
268	ENDIF
269
270	IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN
271	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
272	s o_LWdn200%name,itau_w,LWdn200)
273	ENDIF
274
275	IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN
276	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
277	s o_LWdn200clr%name,itau_w,LWdn200clr)
278	ENDIF
279
280	IF (o_sols%flag(iff)<=lev_files(iff)) THEN
281	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
282	$o_sols%name,itau_w,solsw)
283	ENDIF
284
285	IF (o_sols0%flag(iff)<=lev_files(iff)) THEN
286	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
287	$o_sols0%name,itau_w,solsw0)
288	ENDIF
289
290	IF (o_soll%flag(iff)<=lev_files(iff)) THEN
291	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
292	$o_soll%name,itau_w,sollw)
293	ENDIF
294
295	IF (o_radsol%flag(iff)<=lev_files(iff)) THEN
296	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
297	$o_radsol%name,itau_w,radsol)
298	ENDIF
299
300	IF (o_soll0%flag(iff)<=lev_files(iff)) THEN
301	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
302	$o_soll0%name,itau_w,sollw0)
303	ENDIF
304
305	IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN
306	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
307	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
308	s o_SWupSFC%name,itau_w,zx_tmp_fi2d)
309	ENDIF
310
311	IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN
312	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
313	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
314	$ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d)
315	ENDIF
316
317	IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN
318	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
319	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
320	$ o_SWdnSFC%name,itau_w,zx_tmp_fi2d)
321	ENDIF
322
323	IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
324	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
325	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
326	$ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d)
327	ENDIF
328
329	IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN
330	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
331	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
332	$ o_LWupSFC%name,itau_w,zx_tmp_fi2d)
333	ENDIF
334
335	IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN
336	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
337	$ o_LWdnSFC%name,itau_w,sollwdown)
338	ENDIF
339
340	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
341	IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN
342	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
343	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
344	$ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d)
345	ENDIF
346
347	IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
348	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
349	$ o_LWdnSFCclr%name,itau_w,sollwdownclr)
350	ENDIF
351
352	IF (o_bils%flag(iff)<=lev_files(iff)) THEN
353	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
354	$o_bils%name,itau_w,bils)
355	ENDIF
356
357	IF (o_sens%flag(iff)<=lev_files(iff)) THEN
358	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
359	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
360	$o_sens%name,itau_w,zx_tmp_fi2d)
361	ENDIF
362
363	IF (o_fder%flag(iff)<=lev_files(iff)) THEN
364	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
365	$o_fder%name,itau_w,fder)
366	ENDIF
367
368	IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN
369	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
370	$o_ffonte%name,itau_w,zxffonte)
371	ENDIF
372
373	IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN
374	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
375	$ o_fqcalving%name,itau_w,zxfqcalving)
376	ENDIF
377
378	IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN
379	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
380	$ o_fqfonte%name,itau_w,zxfqfonte)
381	ENDIF
382
383	IF (o_taux%flag(iff)<=lev_files(iff)) THEN
384	zx_tmp_fi2d=0.
385	do nsrf=1,nbsrf
386	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
387	enddo
388	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
389	$ o_taux%name,itau_w,zx_tmp_fi2d)
390	ENDIF
391
392	IF (o_tauy%flag(iff)<=lev_files(iff)) THEN
393	zx_tmp_fi2d=0.
394	do nsrf=1,nbsrf
395	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
396	enddo
397	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
398	$ o_tauy%name,itau_w,zx_tmp_fi2d)
399	ENDIF
400
401	DO nsrf = 1, nbsrf
402	! IF(nsrf.GE.2) THEN
403	IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
404	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
405	CALL histwrite_phy(nid_files(iff),
406	$ clef_stations(iff),
407	$ o_pourc_srf(nsrf)%name,itau_w,
408	$ zx_tmp_fi2d)
409	ENDIF
410
411	IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
412	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
413	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
414	$ o_fract_srf(nsrf)%name,itau_w,
415	$ zx_tmp_fi2d)
416	ENDIF
417	! ENDIF !nsrf.GT.2
418
419	IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
420	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
421	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
422	$ o_taux_srf(nsrf)%name,itau_w,
423	$ zx_tmp_fi2d)
424	ENDIF
425
426	IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
427	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
428	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
429	$ o_tauy_srf(nsrf)%name,itau_w,
430	$ zx_tmp_fi2d)
431	ENDIF
432
433	IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
434	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
435	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
436	$ o_tsol_srf(nsrf)%name,itau_w,
437	$ zx_tmp_fi2d)
438	ENDIF
439
440	IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
441	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
442	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
443	$o_u10m_srf(nsrf)%name,
444	$ itau_w,zx_tmp_fi2d)
445	ENDIF
446
447	IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
448	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
449	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
450	$o_v10m_srf(nsrf)%name,
451	$ itau_w,zx_tmp_fi2d)
452	ENDIF
453
454	IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
455	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
456	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
457	$o_t2m_srf(nsrf)%name,
458	$ itau_w,zx_tmp_fi2d)
459	ENDIF
460
461	IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
462	zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
463	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
464	$o_evap_srf(nsrf)%name,
465	$ itau_w,zx_tmp_fi2d)
466	ENDIF
467
468	IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
469	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
470	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
471	$ o_sens_srf(nsrf)%name,itau_w,
472	$ zx_tmp_fi2d)
473	ENDIF
474
475	IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
476	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
477	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
478	$ o_lat_srf(nsrf)%name,itau_w,
479	$ zx_tmp_fi2d)
480	ENDIF
481
482	IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
483	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
484	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
485	$ o_flw_srf(nsrf)%name,itau_w,
486	$ zx_tmp_fi2d)
487	ENDIF
488
489	IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
490	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
491	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
492	$ o_fsw_srf(nsrf)%name,itau_w,
493	$ zx_tmp_fi2d)
494	ENDIF
495
496	IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
497	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
498	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
499	$ o_wbils_srf(nsrf)%name,itau_w,
500	$ zx_tmp_fi2d)
501	ENDIF
502
503	IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
504	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
505	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
506	$ o_wbilo_srf(nsrf)%name,itau_w,
507	$ zx_tmp_fi2d)
508	ENDIF
509
510	if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
511	IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
512	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
513	$ o_tke_srf(nsrf)%name,itau_w,
514	$ pbl_tke(:,1:klev,nsrf))
515	ENDIF
516
517	IF (.NOT.clef_stations(iff)) THEN
518	IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
519	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
520	$ o_tke_max_srf(nsrf)%name,itau_w,
521	$ pbl_tke(:,1:klev,nsrf))
522	ENDIF
523	ENDIF
524	endif
525	ENDDO
526
527	IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN
528	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
529	$o_cdrm%name,itau_w,cdragm)
530	ENDIF
531
532	IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN
533	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
534	$o_cdrh%name,itau_w,cdragh)
535	ENDIF
536
537	IF (o_cldl%flag(iff)<=lev_files(iff)) THEN
538	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
539	$o_cldl%name,itau_w,cldl)
540	ENDIF
541
542	IF (o_cldm%flag(iff)<=lev_files(iff)) THEN
543	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
544	$o_cldm%name,itau_w,cldm)
545	ENDIF
546
547	IF (o_cldh%flag(iff)<=lev_files(iff)) THEN
548	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
549	$o_cldh%name,itau_w,cldh)
550	ENDIF
551
552	IF (o_cldt%flag(iff)<=lev_files(iff)) THEN
553	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
554	$o_cldt%name,
555	& itau_w,cldt)
556	ENDIF
557
558	IF (o_cldq%flag(iff)<=lev_files(iff)) THEN
559	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
560	$o_cldq%name,itau_w,cldq)
561	ENDIF
562
563	IF (o_lwp%flag(iff)<=lev_files(iff)) THEN
564	zx_tmp_fi2d(1:klon) = flwp(1:klon)
565	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
566	s o_lwp%name,itau_w,zx_tmp_fi2d)
567	ENDIF
568
569	IF (o_iwp%flag(iff)<=lev_files(iff)) THEN
570	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
571	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
572	s o_iwp%name,itau_w,zx_tmp_fi2d)
573	ENDIF
574
575	IF (o_ue%flag(iff)<=lev_files(iff)) THEN
576	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
577	$o_ue%name,itau_w,ue)
578	ENDIF
579
580	IF (o_ve%flag(iff)<=lev_files(iff)) THEN
581	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
582	$o_ve%name,itau_w,ve)
583	ENDIF
584
585	IF (o_uq%flag(iff)<=lev_files(iff)) THEN
586	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
587	$o_uq%name,itau_w,uq)
588	ENDIF
589
590	IF (o_vq%flag(iff)<=lev_files(iff)) THEN
591	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
592	$o_vq%name,itau_w,vq)
593	ENDIF
594
595	IF(iflag_con.GE.3) THEN ! sb
596	IF (o_cape%flag(iff)<=lev_files(iff)) THEN
597	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
598	$o_cape%name,itau_w,cape)
599	ENDIF
600
601	IF (o_pbase%flag(iff)<=lev_files(iff)) THEN
602	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
603	$o_pbase%name,itau_w,ema_pcb)
604	ENDIF
605
606	IF (o_ptop%flag(iff)<=lev_files(iff)) THEN
607	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
608	$o_ptop%name,itau_w,ema_pct)
609	ENDIF
610
611	IF (o_fbase%flag(iff)<=lev_files(iff)) THEN
612	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
613	$o_fbase%name,itau_w,ema_cbmf)
614	ENDIF
615
616	if (iflag_con /= 30) then
617	if (o_plcl%flag(iff)<=lev_files(iff)) THEN
618	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
619	$ o_plcl%name,itau_w,plcl)
620	ENDIF
621
622	IF (o_plfc%flag(iff)<=lev_files(iff)) THEN
623	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
624	$ o_plfc%name,itau_w,plfc)
625	ENDIF
626
627	IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN
628	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
629	$ o_wbeff%name,itau_w,wbeff)
630	ENDIF
631	end if
632
633	IF (.NOT.clef_stations(iff)) THEN
634	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
635	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
636	$o_cape_max%name,itau_w,cape)
637	ENDIF
638	ENDIF
639
640	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
641	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
642	$o_upwd%name,itau_w,upwd)
643	ENDIF
644
645	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
646	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
647	$o_Ma%name,itau_w,Ma)
648	ENDIF
649
650	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
651	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
652	$o_dnwd%name,itau_w,dnwd)
653	ENDIF
654
655	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
656	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
657	$o_dnwd0%name,itau_w,dnwd0)
658	ENDIF
659
660	IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN
661	zx_tmp_fi2d=float(itau_con)/float(itap)
662	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
663	$o_ftime_con%name,
664	s itau_w,zx_tmp_fi2d)
665	ENDIF
666
667	IF (o_mc%flag(iff)<=lev_files(iff)) THEN
668	if(iflag_thermals.gt.1)then
669	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
670	else
671	zx_tmp_fi3d=dnwd+dnwd0+upwd
672	endif
673	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
674	$o_mc%name,itau_w,zx_tmp_fi3d)
675	ENDIF
676
677	ENDIF !iflag_con .GE. 3
678
679	IF (o_prw%flag(iff) <= lev_files(iff))
680	$ CALL histwrite_phy(nid_files(iff), clef_stations(iff),
681	$ o_prw%name, itau_w, prw)
682
683	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
684	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
685	$o_s_pblh%name,itau_w,s_pblh)
686	ENDIF
687
688	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
689	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
690	$o_s_pblt%name,itau_w,s_pblt)
691	ENDIF
692
693	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
694	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
695	$o_s_lcl%name,itau_w,s_lcl)
696	ENDIF
697
698	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
699	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
700	$o_s_therm%name,itau_w,s_therm)
701	ENDIF
702
703	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
704	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
705	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
706	! $o_s_capCL%name,itau_w,s_capCL)
707	! ENDIF
708
709	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
710	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
711	! $o_s_oliqCL%name,itau_w,s_oliqCL)
712	! ENDIF
713
714	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
715	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
716	! $o_s_cteiCL%name,itau_w,s_cteiCL)
717	! ENDIF
718
719	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
720	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
721	! $o_s_trmb1%name,itau_w,s_trmb1)
722	! ENDIF
723
724	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
725	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
726	! $o_s_trmb2%name,itau_w,s_trmb2)
727	! ENDIF
728
729	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
730	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
731	! $o_s_trmb3%name,itau_w,s_trmb3)
732	! ENDIF
733
734	! Champs interpolles sur des niveaux de pression
735
736	ll=0
737	DO k=1, nlevSTD
738	bb2=clevSTD(k)
739	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
740	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
741	$ bb2.EQ."100".OR.
742	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
743
744	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
745	ll=ll+1
746	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
747	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
748	$o_uSTDlevs(ll)%name,
749	& itau_w,uwriteSTD(:,k,iff))
750	ENDIF
751
752	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
753	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
754	$o_vSTDlevs(ll)%name,
755	& itau_w,vwriteSTD(:,k,iff))
756	ENDIF
757
758	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
759	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
760	$o_wSTDlevs(ll)%name,
761	& itau_w,wwriteSTD(:,k,iff))
762	ENDIF
763
764	IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
765	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
766	$o_zSTDlevs(ll)%name,
767	& itau_w,phiwriteSTD(:,k,iff))
768	ENDIF
769
770	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
771	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
772	$o_qSTDlevs(ll)%name,
773	& itau_w, qwriteSTD(:,k,iff))
774	ENDIF
775
776	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
777	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
778	$o_tSTDlevs(ll)%name,
779	& itau_w, twriteSTD(:,k,iff))
780	ENDIF
781
782	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
783	ENDDO
784
785	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
786	DO i=1, klon
787	IF (pctsrf(i,is_oce).GT.epsfra.OR.
788	$ pctsrf(i,is_sic).GT.epsfra) THEN
789	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
790	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
791	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
792	ELSE
793	zx_tmp_fi2d(i) = 273.15
794	ENDIF
795	ENDDO
796	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
797	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
798	ENDIF
799
800	! Couplage convection-couche limite
801	IF (iflag_con.GE.3) THEN
802	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
803	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
804	$o_ale_bl%name,itau_w,ale_bl)
805	ENDIF
806	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
807	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
808	$o_alp_bl%name,itau_w,alp_bl)
809	ENDIF
810	ENDIF !(iflag_con.GE.3)
811
812	! Wakes
813	IF (iflag_con.EQ.3) THEN
814	IF (iflag_wake>=1) THEN
815	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
816	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
817	$o_ale_wk%name,itau_w,ale_wake)
818	ENDIF
819	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
820	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
821	$o_alp_wk%name,itau_w,alp_wake)
822	ENDIF
823
824	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
825	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
826	$o_WAPE%name,itau_w,wake_pe)
827	ENDIF
828	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
829	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
830	$o_wake_h%name,itau_w,wake_h)
831	ENDIF
832
833	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
834	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
835	$o_wake_s%name,itau_w,wake_s)
836	ENDIF
837
838	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
839	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
840	$o_wake_deltat%name,
841	$ itau_w,wake_deltat)
842	ENDIF
843
844	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
845	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
846	$o_wake_deltaq%name,
847	$ itau_w,wake_deltaq)
848	ENDIF
849
850	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
851	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
852	s o_wake_omg%name,itau_w,wake_omg)
853	ENDIF
854
855	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
856	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
857	& /pdtphys
858	CALL histwrite_phy(nid_files(iff),
859	$clef_stations(iff),
860	& o_dtwak%name,itau_w,zx_tmp_fi3d)
861	ENDIF
862
863	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
864	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
865	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
866	& o_dqwak%name,itau_w,zx_tmp_fi3d)
867	ENDIF
868	ENDIF ! iflag_wake>=1
869
870	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
871	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
872	$o_ale%name,itau_w,ale)
873	ENDIF
874	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
875	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
876	$o_alp%name,itau_w,alp)
877	ENDIF
878	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
879	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
880	$o_cin%name,itau_w,cin)
881	ENDIF
882	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
883	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
884	$o_Vprecip%name,itau_w,Vprecip)
885	ENDIF
886
887	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
888	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
889	$o_ftd%name,itau_w,ftd)
890	ENDIF
891
892	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
893	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
894	$o_fqd%name,itau_w,fqd)
895	ENDIF
896	ENDIF !(iflag_con.EQ.3)
897
898	IF (type_ocean=='slab ') THEN
899	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
900	$ CALL histwrite_phy(
901	$ nid_files(iff),clef_stations(iff),
902	$o_slab_bils%name,itau_w,slab_wfbils)
903
904	ENDIF !type_ocean == force/slab
905
906	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
907	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
908	s o_weakinv%name,itau_w,weak_inversion)
909	ENDIF
910
911	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
912	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
913	$o_dthmin%name,itau_w,dthmin)
914	ENDIF
915
916	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
917	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
918	$o_cldtau%name,itau_w,cldtau)
919	ENDIF
920
921	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
922	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
923	$o_cldemi%name,itau_w,cldemi)
924	ENDIF
925
926	IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
927	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
928	s o_pr_con_l%name,itau_w,pmflxr(:,1:klev))
929	ENDIF
930
931	IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
932	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
933	s o_pr_con_i%name,itau_w,pmflxs(:,1:klev))
934	ENDIF
935
936	IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
937	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
938	s o_pr_lsc_l%name,itau_w,prfl(:,1:klev))
939	ENDIF
940
941	IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
942	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
943	s o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
944	ENDIF
945
946	IF (o_re%flag(iff)<=lev_files(iff)) THEN
947	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
948	$o_re%name,itau_w,re)
949	ENDIF
950
951	IF (o_fl%flag(iff)<=lev_files(iff)) THEN
952	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
953	$o_fl%name,itau_w,fl)
954	ENDIF
955
956
957
958	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
959	DO i=1, klon
960	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
961	ENDDO
962	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
963	$o_rh2m%name,itau_w,zx_tmp_fi2d)
964	ENDIF
965
966	IF (.NOT.clef_stations(iff)) THEN
967	IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN
968	DO i=1, klon
969	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
970	ENDDO
971	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
972	$o_rh2m_min%name,
973	s itau_w,zx_tmp_fi2d)
974	ENDIF
975
976	IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN
977	DO i=1, klon
978	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
979	ENDDO
980	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
981	$o_rh2m_max%name,
982	s itau_w,zx_tmp_fi2d)
983	ENDIF
984	ENDIF
985
986
987	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
988	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
989	$o_qsat2m%name,itau_w,qsat2m)
990	ENDIF
991
992	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
993	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
994	$o_tpot%name,itau_w,tpot)
995	ENDIF
996
997	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
998	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
999	$o_tpote%name,itau_w,tpote)
1000	ENDIF
1001
1002	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
1003	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
1004	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1005	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
1006	ENDIF
1007
1008	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
1009	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
1010	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1011	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
1012	ENDIF
1013
1014	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
1015	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1016	s o_LWdownOR%name,itau_w,sollwdown)
1017	ENDIF
1018
1019	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
1020	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1021	$o_snowl%name,itau_w,snow_lsc)
1022	ENDIF
1023
1024	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
1025	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1026	s o_solldown%name,itau_w,sollwdown)
1027	ENDIF
1028
1029	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
1030	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1031	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
1032	ENDIF
1033
1034	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
1035	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1036	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
1037	ENDIF
1038
1039	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
1040	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1041	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
1042	ENDIF
1043
1044	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
1045	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1046	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
1047	ENDIF
1048
1049	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
1050	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1051	$o_rugs%name,itau_w,zxrugs)
1052	ENDIF
1053
1054	! OD550 per species
1055	IF (new_aod .and. (.not. aerosol_couple)) THEN
1056	IF (ok_ade.OR.ok_aie) THEN
1057
1058	IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN
1059	CALL histwrite_phy(nid_files(iff),
1060	$clef_stations(iff),
1061	$o_od550aer%name,itau_w,
1062	$ od550aer)
1063	ENDIF
1064	IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN
1065	CALL histwrite_phy(nid_files(iff),
1066	$clef_stations(iff),
1067	$o_od865aer%name,itau_w,
1068	$ od865aer)
1069	ENDIF
1070	IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN
1071	CALL histwrite_phy(nid_files(iff),
1072	$clef_stations(iff),
1073	$o_absvisaer%name,itau_w,
1074	$ absvisaer)
1075	ENDIF
1076	IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN
1077	CALL histwrite_phy(nid_files(iff),
1078	$clef_stations(iff),
1079	$o_od550lt1aer%name,itau_w,
1080	$ od550lt1aer)
1081	ENDIF
1082
1083	IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN
1084	CALL histwrite_phy(nid_files(iff),
1085	$clef_stations(iff),
1086	$o_sconcso4%name,itau_w,
1087	$ sconcso4)
1088	ENDIF
1089	IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN
1090	CALL histwrite_phy(nid_files(iff),
1091	$clef_stations(iff),
1092	$o_sconcoa%name,itau_w,
1093	$ sconcoa)
1094	ENDIF
1095	IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN
1096	CALL histwrite_phy(nid_files(iff),
1097	$clef_stations(iff),
1098	$o_sconcbc%name,itau_w,
1099	$ sconcbc)
1100	ENDIF
1101	IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN
1102	CALL histwrite_phy(nid_files(iff),
1103	$clef_stations(iff),
1104	$o_sconcss%name,itau_w,
1105	$ sconcss)
1106	ENDIF
1107	IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN
1108	CALL histwrite_phy(nid_files(iff),
1109	$clef_stations(iff),
1110	$o_sconcdust%name,itau_w,
1111	$ sconcdust)
1112	ENDIF
1113
1114	IF (o_concso4%flag(iff)<=lev_files(iff)) THEN
1115	CALL histwrite_phy(nid_files(iff),
1116	$clef_stations(iff),
1117	$o_concso4%name,itau_w,
1118	$ concso4)
1119	ENDIF
1120	IF (o_concoa%flag(iff)<=lev_files(iff)) THEN
1121	CALL histwrite_phy(nid_files(iff),
1122	$clef_stations(iff),
1123	$o_concoa%name,itau_w,
1124	$ concoa)
1125	ENDIF
1126	IF (o_concbc%flag(iff)<=lev_files(iff)) THEN
1127	CALL histwrite_phy(nid_files(iff),
1128	$clef_stations(iff),
1129	$o_concbc%name,itau_w,
1130	$ concbc)
1131	ENDIF
1132	IF (o_concss%flag(iff)<=lev_files(iff)) THEN
1133	CALL histwrite_phy(nid_files(iff),
1134	$clef_stations(iff),
1135	$o_concss%name,itau_w,
1136	$ concss)
1137	ENDIF
1138	IF (o_concdust%flag(iff)<=lev_files(iff)) THEN
1139	CALL histwrite_phy(nid_files(iff),
1140	$clef_stations(iff),
1141	$o_concdust%name,itau_w,
1142	$ concdust)
1143	ENDIF
1144
1145	IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN
1146	CALL histwrite_phy(nid_files(iff),
1147	$clef_stations(iff),
1148	$o_loadso4%name,itau_w,
1149	$ loadso4)
1150	ENDIF
1151	IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN
1152	CALL histwrite_phy(nid_files(iff),
1153	$clef_stations(iff),
1154	$o_loadoa%name,itau_w,
1155	$ loadoa)
1156	ENDIF
1157	IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN
1158	CALL histwrite_phy(nid_files(iff),
1159	$clef_stations(iff),
1160	$o_loadbc%name,itau_w,
1161	$ loadbc)
1162	ENDIF
1163	IF (o_loadss%flag(iff)<=lev_files(iff)) THEN
1164	CALL histwrite_phy(nid_files(iff),
1165	$clef_stations(iff),
1166	$o_loadss%name,itau_w,
1167	$ loadss)
1168	ENDIF
1169	IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN
1170	CALL histwrite_phy(nid_files(iff),
1171	$clef_stations(iff),
1172	$o_loaddust%name,itau_w,
1173	$ loaddust)
1174	ENDIF
1175
1176	DO naero = 1, naero_spc
1177	IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
1178	CALL histwrite_phy(nid_files(iff),
1179	$clef_stations(iff),
1180	$ o_tausumaero(naero)%name,itau_w,
1181	$ tausum_aero(:,2,naero) )
1182	ENDIF
1183	END DO
1184	endif
1185	ENDIF
1186
1187	IF (ok_ade) THEN
1188	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
1189	CALL histwrite_phy(nid_files(iff),
1190	$clef_stations(iff),
1191	$o_topswad%name,itau_w,
1192	$ topswad_aero)
1193	ENDIF
1194	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
1195	CALL histwrite_phy(nid_files(iff),
1196	$clef_stations(iff),
1197	$o_solswad%name,itau_w,
1198	$ solswad_aero)
1199	ENDIF
1200
1201	!====MS forcing diagnostics
1202	if (new_aod) then
1203	IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
1204	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1205	$o_swtoaas_nat%name,itau_w,
1206	$ topsw_aero(:,1))
1207	ENDIF
1208
1209	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
1210	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1211	$o_swsrfas_nat%name,itau_w,
1212	$ solsw_aero(:,1))
1213	ENDIF
1214
1215	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
1216	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1217	$o_swtoacs_nat%name,itau_w,
1218	$ topsw0_aero(:,1))
1219	ENDIF
1220
1221	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
1222	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1223	$o_swsrfcs_nat%name,itau_w,
1224	$ solsw0_aero(:,1))
1225	ENDIF
1226
1227	!ant
1228	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
1229	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1230	$o_swtoaas_ant%name,itau_w,
1231	$ topsw_aero(:,2))
1232	ENDIF
1233
1234	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
1235	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1236	$o_swsrfas_ant%name,itau_w,
1237	$ solsw_aero(:,2))
1238	ENDIF
1239
1240	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
1241	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1242	$o_swtoacs_ant%name,itau_w,
1243	$ topsw0_aero(:,2))
1244	ENDIF
1245
1246	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
1247	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1248	$o_swsrfcs_ant%name,itau_w,
1249	$ solsw0_aero(:,2))
1250	ENDIF
1251
1252	!cf
1253
1254	if (.not. aerosol_couple) then
1255	IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
1256	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1257	$o_swtoacf_nat%name,itau_w,
1258	$ topswcf_aero(:,1))
1259	ENDIF
1260
1261	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
1262	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1263	$o_swsrfcf_nat%name,itau_w,
1264	$ solswcf_aero(:,1))
1265	ENDIF
1266
1267	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
1268	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1269	$o_swtoacf_ant%name,itau_w,
1270	$ topswcf_aero(:,2))
1271	ENDIF
1272
1273	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
1274	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1275	$o_swsrfcf_ant%name,itau_w,
1276	$ solswcf_aero(:,2))
1277	ENDIF
1278
1279	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
1280	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1281	$o_swtoacf_zero%name,itau_w,
1282	$ topswcf_aero(:,3))
1283	ENDIF
1284
1285	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
1286	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1287	$o_swsrfcf_zero%name,itau_w,
1288	$ solswcf_aero(:,3))
1289	ENDIF
1290	endif
1291
1292	endif ! new_aod
1293	!====MS forcing diagnostics
1294
1295	ENDIF
1296
1297	IF (ok_aie) THEN
1298	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
1299	CALL histwrite_phy(nid_files(iff),
1300	$clef_stations(iff),
1301	$o_topswai%name,itau_w,
1302	$ topswai_aero)
1303	ENDIF
1304	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
1305	CALL histwrite_phy(nid_files(iff),
1306	$clef_stations(iff),
1307	$o_solswai%name,itau_w,
1308	$ solswai_aero)
1309	ENDIF
1310	IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN
1311	CALL histwrite_phy(nid_files(iff),
1312	$clef_stations(iff),
1313	$o_scdnc%name,itau_w,
1314	$ scdnc)
1315	ENDIF
1316	IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN
1317	CALL histwrite_phy(nid_files(iff),
1318	$clef_stations(iff),
1319	$o_cldncl%name,itau_w,
1320	$ cldncl)
1321	ENDIF
1322	IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN
1323	CALL histwrite_phy(nid_files(iff),
1324	$clef_stations(iff),
1325	$o_reffclws%name,itau_w,
1326	$ reffclws)
1327	ENDIF
1328	IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN
1329	CALL histwrite_phy(nid_files(iff),
1330	$clef_stations(iff),
1331	$o_reffclwc%name,itau_w,
1332	$ reffclwc)
1333	ENDIF
1334	IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN
1335	CALL histwrite_phy(nid_files(iff),
1336	$clef_stations(iff),
1337	$o_cldnvi%name,itau_w,
1338	$ cldnvi)
1339	ENDIF
1340	IF (o_lcc%flag(iff)<=lev_files(iff)) THEN
1341	CALL histwrite_phy(nid_files(iff),
1342	$clef_stations(iff),
1343	$o_lcc%name,itau_w,
1344	$ lcc)
1345	ENDIF
1346	IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN
1347	CALL histwrite_phy(nid_files(iff),
1348	$clef_stations(iff),
1349	$o_lcc3d%name,itau_w,
1350	$ lcc3d)
1351	ENDIF
1352	IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN
1353	CALL histwrite_phy(nid_files(iff),
1354	$clef_stations(iff),
1355	$o_lcc3dcon%name,itau_w,
1356	$ lcc3dcon)
1357	ENDIF
1358	IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN
1359	CALL histwrite_phy(nid_files(iff),
1360	$clef_stations(iff),
1361	$o_lcc3dstra%name,itau_w,
1362	$ lcc3dstra)
1363	ENDIF
1364	IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN
1365	CALL histwrite_phy(nid_files(iff),
1366	$clef_stations(iff),
1367	$o_reffclwtop%name,itau_w,
1368	$ reffclwtop)
1369	ENDIF
1370	ENDIF
1371
1372	! Champs 3D:
1373	IF (ok_ade .OR. ok_aie) then
1374	IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN
1375	CALL histwrite_phy(nid_files(iff),
1376	$clef_stations(iff),
1377	$o_ec550aer%name,itau_w,
1378	& ec550aer)
1379	ENDIF
1380	ENDIF
1381
1382	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
1383	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1384	$o_lwcon%name,itau_w,flwc)
1385	ENDIF
1386
1387	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
1388	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1389	$o_iwcon%name,itau_w,fiwc)
1390	ENDIF
1391
1392	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
1393	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1394	$o_temp%name,itau_w,t_seri)
1395	ENDIF
1396
1397	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
1398	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1399	$o_theta%name,itau_w,theta)
1400	ENDIF
1401
1402	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
1403	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1404	$o_ovapinit%name,itau_w,
1405	$ qx(:,:,ivap))
1406	ENDIF
1407
1408	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
1409	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1410	$ o_ovap%name,itau_w,q_seri)
1411	ENDIF
1412
1413	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
1414	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1415	$o_geop%name,itau_w,zphi)
1416	ENDIF
1417
1418	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
1419	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1420	$o_vitu%name,itau_w,u_seri)
1421	ENDIF
1422
1423	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
1424	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1425	$o_vitv%name,itau_w,v_seri)
1426	ENDIF
1427
1428	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
1429	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1430	$o_vitw%name,itau_w,omega)
1431	ENDIF
1432
1433	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
1434	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1435	$o_pres%name,itau_w,pplay)
1436	ENDIF
1437
1438	IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
1439	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1440	$o_paprs%name,
1441	s itau_w,paprs(:,1:klev))
1442	ENDIF
1443
1444	IF (o_zfull%flag(iff)<=lev_files(iff)) THEN
1445	DO i=1, klon
1446	zx_tmp_fi3d1(i,1)= pphis(i)/RG
1447	!020611 zx_tmp_fi3d(i,1)= pphis(i)/RG
1448	ENDDO
1449	DO k=1, klev
1450	!020611 DO k=1, klev-1
1451	DO i=1, klon
1452	!020611 zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (t_seri(i,k) RD
1453	zx_tmp_fi3d1(i,k+1)= zx_tmp_fi3d1(i,k) - (t_seri(i,k) RD
1454	$ (paprs(i,k+1) - paprs(i,k))) / ( pplay(i,k) * RG )
1455	ENDDO
1456	ENDDO
1457	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1458	$o_zfull%name,itau_w,zx_tmp_fi3d1(:,2:klevp1))
1459	!020611 $o_zfull%name,itau_w,zx_tmp_fi3d)
1460	ENDIF
1461
1462	IF (o_zhalf%flag(iff)<=lev_files(iff)) THEN
1463	DO i=1, klon
1464	zx_tmp_fi3d(i,1)= pphis(i)/RG - (
1465	$ (t_seri(i,1)+zxtsol(i))/2. RD
1466	$ (pplay(i,1) - paprs(i,1)))/( (paprs(i,1)+pplay(i,1))/2. * RG)
1467	ENDDO
1468	DO k=1, klev-1
1469	DO i=1, klon
1470	zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (
1471	$ (t_seri(i,k)+t_seri(i,k+1))/2. RD
1472	$ (pplay(i,k+1) - pplay(i,k))) / ( paprs(i,k) * RG )
1473	ENDDO
1474	ENDDO
1475	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1476	$o_zhalf%name,itau_w,zx_tmp_fi3d)
1477	ENDIF
1478
1479	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
1480	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1481	$o_rneb%name,itau_w,cldfra)
1482	ENDIF
1483
1484	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
1485	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1486	$o_rnebcon%name,itau_w,rnebcon)
1487	ENDIF
1488
1489	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
1490	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1491	$o_rhum%name,itau_w,zx_rh)
1492	ENDIF
1493
1494	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
1495	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1496	$o_ozone%name, itau_w,
1497	$ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1498	ENDIF
1499
1500	IF (o_ozone_light%flag(iff)<=lev_files(iff) .and.
1501	$ read_climoz == 2) THEN
1502	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1503	$o_ozone_light%name, itau_w,
1504	$ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1505	ENDIF
1506
1507	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
1508	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1509	$o_dtphy%name,itau_w,d_t)
1510	ENDIF
1511
1512	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
1513	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1514	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
1515	ENDIF
1516
1517	DO nsrf=1, nbsrf
1518	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1519	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1520	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1521	s o_albe_srf(nsrf)%name,itau_w,
1522	$ zx_tmp_fi2d)
1523	ENDIF
1524
1525	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1526	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1527	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1528	s o_rugs_srf(nsrf)%name,itau_w,
1529	$ zx_tmp_fi2d)
1530	ENDIF
1531
1532	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1533	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1534	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1535	s o_ages_srf(nsrf)%name,itau_w
1536	$ ,zx_tmp_fi2d)
1537	ENDIF
1538	ENDDO !nsrf=1, nbsrf
1539
1540	IF (o_alb1%flag(iff)<=lev_files(iff)) THEN
1541	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1542	$o_alb1%name,itau_w,albsol1)
1543	ENDIF
1544
1545	IF (o_alb2%flag(iff)<=lev_files(iff)) THEN
1546	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1547	$o_alb2%name,itau_w,albsol2)
1548	ENDIF
1549
1550	!FH Sorties pour la couche limite
1551	if (iflag_pbl>1) then
1552	zx_tmp_fi3d=0.
1553	do nsrf=1,nbsrf
1554	do k=1,klev
1555	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
1556	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1557	enddo
1558	enddo
1559	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
1560	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1561	$o_tke%name,itau_w,zx_tmp_fi3d)
1562	ENDIF
1563
1564	IF (.NOT.clef_stations(iff)) THEN
1565	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
1566	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1567	s o_tke_max%name,itau_w,zx_tmp_fi3d)
1568	ENDIF
1569	ENDIF
1570	endif
1571
1572	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
1573	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1574	$o_kz%name,itau_w,coefh)
1575	ENDIF
1576
1577	IF (.NOT.clef_stations(iff)) THEN
1578	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
1579	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1580	$o_kz_max%name,itau_w,coefh)
1581	ENDIF
1582	ENDIF
1583
1584	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
1585	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1586	$o_clwcon%name,itau_w,clwcon0)
1587	ENDIF
1588
1589	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
1590	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1591	$o_dtdyn%name,itau_w,d_t_dyn)
1592	ENDIF
1593
1594	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
1595	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1596	$o_dqdyn%name,itau_w,d_q_dyn)
1597	ENDIF
1598
1599	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
1600	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1601	$o_dudyn%name,itau_w,d_u_dyn)
1602	ENDIF
1603
1604	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
1605	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1606	$o_dvdyn%name,itau_w,d_v_dyn)
1607	ENDIF
1608
1609	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
1610	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1611	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1612	$o_dtcon%name,itau_w,zx_tmp_fi3d)
1613	ENDIF
1614
1615	if(iflag_thermals.eq.1)then
1616	IF (o_tntc%flag(iff)<=lev_files(iff)) THEN
1617	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys +
1618	$ d_t_ajsb(1:klon,1:klev)/pdtphys
1619	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1620	$o_tntc%name,itau_w,zx_tmp_fi3d)
1621	ENDIF
1622	else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then
1623	IF (o_tntc%flag(iff)<=lev_files(iff)) THEN
1624	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys +
1625	$ d_t_ajs(1:klon,1:klev)/pdtphys +
1626	$ d_t_wake(1:klon,1:klev)/pdtphys
1627	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1628	$o_tntc%name,itau_w,zx_tmp_fi3d)
1629	ENDIF
1630	endif
1631
1632	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
1633	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1634	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1635	$o_ducon%name,itau_w,zx_tmp_fi3d)
1636	ENDIF
1637
1638	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
1639	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1640	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1641	$o_dqcon%name,itau_w,zx_tmp_fi3d)
1642	ENDIF
1643
1644	if(iflag_thermals.eq.1)then
1645	IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN
1646	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1647	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1648	$o_tnhusc%name,itau_w,zx_tmp_fi3d)
1649	ENDIF
1650	else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then
1651	IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN
1652	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys +
1653	$ d_q_ajs(1:klon,1:klev)/pdtphys +
1654	$ d_q_wake(1:klon,1:klev)/pdtphys
1655	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1656	$o_tnhusc%name,itau_w,zx_tmp_fi3d)
1657	ENDIF
1658	endif
1659
1660	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
1661	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1662	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1663	$o_dtlsc%name,itau_w,zx_tmp_fi3d)
1664	ENDIF
1665
1666	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
1667	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1668	$ d_t_eva(1:klon,1:klev))/pdtphys
1669	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1670	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
1671	ENDIF
1672
1673	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
1674	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1675	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1676	$o_dqlsc%name,itau_w,zx_tmp_fi3d)
1677	ENDIF
1678
1679	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1680	! Sorties specifiques a la separation thermiques/non thermiques
1681	if (iflag_thermals>1) then
1682
1683	IF (o_dtlscth%flag(iff)<=lev_files(iff)) THEN
1684	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1685	CALL
1686	s histwrite_phy(nid_files(iff),clef_stations(iff),
1687	$o_dtlscth%name,itau_w,zx_tmp_fi3d)
1688	ENDIF
1689
1690	IF (o_dtlscst%flag(iff)<=lev_files(iff)) THEN
1691	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1692	CALL
1693	s histwrite_phy(nid_files(iff),clef_stations(iff),
1694	$o_dtlscst%name,itau_w,zx_tmp_fi3d)
1695	ENDIF
1696
1697	IF (o_dqlscth%flag(iff)<=lev_files(iff)) THEN
1698	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1699	CALL
1700	s histwrite_phy(nid_files(iff),clef_stations(iff),
1701	$o_dqlscth%name,itau_w,zx_tmp_fi3d)
1702	ENDIF
1703
1704	IF (o_dqlscst%flag(iff)<=lev_files(iff)) THEN
1705	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1706	CALL
1707	s histwrite_phy(nid_files(iff),clef_stations(iff),
1708	$o_dqlscst%name,itau_w,zx_tmp_fi3d)
1709	ENDIF
1710
1711	IF (o_plulth%flag(iff)<=lev_files(iff)) THEN
1712	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1713	$o_plulth%name,itau_w,plul_th)
1714	ENDIF
1715
1716	IF (o_plulst%flag(iff)<=lev_files(iff)) THEN
1717	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1718	$o_plulst%name,itau_w,plul_st)
1719	ENDIF
1720
1721	do k=1,klev
1722	do i=1,klon
1723	if (ptconvth(i,k)) then
1724	zx_tmp_fi3d(i,k)=1.
1725	else
1726	zx_tmp_fi3d(i,k)=0.
1727	endif
1728	enddo
1729	enddo
1730	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1731	CALL
1732	s histwrite_phy(nid_files(iff),clef_stations(iff),
1733	$o_ptconvth%name,itau_w,zx_tmp_fi3d)
1734	ENDIF
1735
1736	do i=1,klon
1737	zx_tmp_fi2d(1:klon)=lmax_th(:)
1738	enddo
1739	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1740	CALL
1741	s histwrite_phy(nid_files(iff),clef_stations(iff),
1742	$o_lmaxth%name,itau_w,zx_tmp_fi2d)
1743	ENDIF
1744
1745	endif ! iflag_thermals>1
1746
1747	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1748
1749	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
1750	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1751	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1752	$o_dtvdf%name,itau_w,zx_tmp_fi3d)
1753	ENDIF
1754
1755	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
1756	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1757	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1758	$o_dqvdf%name,itau_w,zx_tmp_fi3d)
1759	ENDIF
1760
1761	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
1762	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1763	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1764	$o_dteva%name,itau_w,zx_tmp_fi3d)
1765	ENDIF
1766
1767	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
1768	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1769	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1770	$o_dqeva%name,itau_w,zx_tmp_fi3d)
1771	ENDIF
1772
1773	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
1774	zpt_conv = 0.
1775	where (ptconv) zpt_conv = 1.
1776	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1777	$o_ptconv%name,itau_w,zpt_conv)
1778	ENDIF
1779
1780	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
1781	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1782	$o_ratqs%name,itau_w,ratqs)
1783	ENDIF
1784
1785	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
1786	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys -
1787	$ d_t_ajsb(1:klon,1:klev)/pdtphys
1788	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1789	$o_dtthe%name,itau_w,zx_tmp_fi3d)
1790	ENDIF
1791
1792	IF (iflag_thermals.gt.1) THEN
1793	IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN
1794	! Pour l instant 0 a y reflichir pour les thermiques
1795	zx_tmp_fi2d=0.
1796	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1797	$o_ftime_th%name,
1798	s itau_w,zx_tmp_fi2d)
1799	ENDIF
1800
1801	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
1802	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1803	$o_f_th%name,itau_w,fm_therm)
1804	ENDIF
1805
1806	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
1807	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1808	$o_e_th%name,itau_w,entr_therm)
1809	ENDIF
1810
1811	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
1812	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1813	$o_w_th%name,itau_w,zw2)
1814	ENDIF
1815
1816	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
1817	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1818	$o_q_th%name,itau_w,zqasc)
1819	ENDIF
1820
1821	IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN
1822	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1823	s o_lambda_th%name,itau_w,lambda_th)
1824	ENDIF
1825
1826	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
1827	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1828	$o_a_th%name,itau_w,fraca)
1829	ENDIF
1830
1831	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
1832	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1833	$o_d_th%name,itau_w,detr_therm)
1834	ENDIF
1835
1836	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1837	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1838	$o_f0_th%name,itau_w,f0)
1839	ENDIF
1840
1841	IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN
1842	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1843	$o_zmax_th%name,itau_w,zmax_th)
1844	ENDIF
1845
1846	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
1847	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys -
1848	$ d_q_ajsb(1:klon,1:klev)/pdtphys
1849	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1850	$o_dqthe%name,itau_w,zx_tmp_fi3d)
1851	ENDIF
1852
1853	ENDIF !iflag_thermals
1854
1855	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
1856	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1857	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1858	$o_dtajs%name,itau_w,zx_tmp_fi3d)
1859	ENDIF
1860
1861	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
1862	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1863	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1864	$o_dqajs%name,itau_w,zx_tmp_fi3d)
1865	ENDIF
1866
1867	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
1868	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1869	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1870	$o_dtswr%name,itau_w,zx_tmp_fi3d)
1871	ENDIF
1872
1873	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
1874	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1875	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1876	$o_dtsw0%name,itau_w,zx_tmp_fi3d)
1877	ENDIF
1878
1879	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
1880	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1881	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1882	$o_dtlwr%name,itau_w,zx_tmp_fi3d)
1883	ENDIF
1884
1885	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
1886	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1887	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1888	$o_dtlw0%name,itau_w,zx_tmp_fi3d)
1889	ENDIF
1890
1891	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
1892	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1893	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1894	$o_dtec%name,itau_w,zx_tmp_fi3d)
1895	ENDIF
1896
1897	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
1898	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1899	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1900	$o_duvdf%name,itau_w,zx_tmp_fi3d)
1901	ENDIF
1902
1903	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
1904	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1905	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1906	$o_dvvdf%name,itau_w,zx_tmp_fi3d)
1907	ENDIF
1908
1909	IF (ok_orodr) THEN
1910	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
1911	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1912	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1913	$o_duoro%name,itau_w,zx_tmp_fi3d)
1914	ENDIF
1915
1916	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
1917	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1918	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1919	$o_dvoro%name,itau_w,zx_tmp_fi3d)
1920	ENDIF
1921
1922	IF (o_dtoro%flag(iff)<=lev_files(iff)) THEN
1923	zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
1924	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1925	$o_dtoro%name,itau_w,zx_tmp_fi3d)
1926	ENDIF
1927	ENDIF
1928
1929	IF (ok_orolf) THEN
1930	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
1931	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1932	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1933	$o_dulif%name,itau_w,zx_tmp_fi3d)
1934	ENDIF
1935
1936	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
1937	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1938	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1939	$o_dvlif%name,itau_w,zx_tmp_fi3d)
1940	ENDIF
1941
1942	IF (o_dtlif%flag(iff)<=lev_files(iff)) THEN
1943	zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
1944	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1945	$o_dtlif%name,itau_w,zx_tmp_fi3d)
1946	ENDIF
1947	ENDIF
1948
1949	IF (ok_hines) THEN
1950	IF (o_duhin%flag(iff)<=lev_files(iff)) THEN
1951	zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys
1952	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1953	$o_duhin%name,itau_w,zx_tmp_fi3d)
1954	ENDIF
1955
1956
1957	IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN
1958	zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys
1959	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1960	$o_dvhin%name,itau_w,zx_tmp_fi3d)
1961	ENDIF
1962
1963	IF (o_dthin%flag(iff)<=lev_files(iff)) THEN
1964	zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
1965	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1966	$o_dthin%name,itau_w,zx_tmp_fi3d)
1967	ENDIF
1968	ENDIF
1969
1970	IF (o_rsu%flag(iff)<=lev_files(iff)) THEN
1971	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1972	$o_rsu%name,itau_w,swup)
1973	ENDIF
1974	IF (o_rsd%flag(iff)<=lev_files(iff)) THEN
1975	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1976	$o_rsd%name,itau_w,swdn)
1977	ENDIF
1978	IF (o_rlu%flag(iff)<=lev_files(iff)) THEN
1979	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1980	$o_rlu%name,itau_w,lwup)
1981	ENDIF
1982	IF (o_rld%flag(iff)<=lev_files(iff)) THEN
1983	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1984	$o_rld%name,itau_w,lwdn)
1985	ENDIF
1986
1987	IF (o_rsucs%flag(iff)<=lev_files(iff)) THEN
1988	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1989	$o_rsucs%name,itau_w,swup0)
1990	ENDIF
1991	IF (o_rsdcs%flag(iff)<=lev_files(iff)) THEN
1992	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1993	$o_rsdcs%name,itau_w,swdn0)
1994	ENDIF
1995	IF (o_rlucs%flag(iff)<=lev_files(iff)) THEN
1996	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1997	$o_rlucs%name,itau_w,lwup0)
1998	ENDIF
1999	IF (o_rldcs%flag(iff)<=lev_files(iff)) THEN
2000	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2001	$o_rldcs%name,itau_w,lwdn0)
2002	ENDIF
2003
2004	IF (o_tnt%flag(iff)<=lev_files(iff)) THEN
2005	zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+
2006	$d_t_dyn(1:klon,1:klev)
2007	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2008	$o_tnt%name,itau_w,zx_tmp_fi3d)
2009	ENDIF
2010
2011	IF (o_tntr%flag(iff)<=lev_files(iff)) THEN
2012	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY -
2013	$cool(1:klon,1:klev)/RDAY
2014	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2015	$o_tntr%name,itau_w,zx_tmp_fi3d)
2016	ENDIF
2017
2018	IF (o_tntscpbl%flag(iff)<=lev_files(iff)) THEN
2019	zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+
2020	$ d_t_eva(1:klon,1:klev)+
2021	$ d_t_vdf(1:klon,1:klev))/pdtphys
2022	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2023	$o_tntscpbl%name,itau_w,zx_tmp_fi3d)
2024	ENDIF
2025
2026	IF (o_tnhus%flag(iff)<=lev_files(iff)) THEN
2027	zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+
2028	$d_q_dyn(1:klon,1:klev)
2029	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2030	$o_tnhus%name,itau_w,zx_tmp_fi3d)
2031	ENDIF
2032
2033	IF (o_tnhusscpbl%flag(iff)<=lev_files(iff)) THEN
2034	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+
2035	$ d_q_eva(1:klon,1:klev)/pdtphys
2036	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2037	$o_tnhusscpbl%name,itau_w,zx_tmp_fi3d)
2038	ENDIF
2039
2040	IF (o_evu%flag(iff)<=lev_files(iff)) THEN
2041	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2042	$o_evu%name,itau_w,coefm)
2043	ENDIF
2044
2045	IF (o_h2o%flag(iff)<=lev_files(iff)) THEN
2046	zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+
2047	$ ql_seri(1:klon,1:klev)
2048	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2049	$o_h2o%name,itau_w,zx_tmp_fi3d)
2050	ENDIF
2051
2052	IF (o_mcd%flag(iff)<=lev_files(iff)) THEN
2053	zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+
2054	$ dnwd0(1:klon,1:klev))
2055	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2056	$o_mcd%name,itau_w,zx_tmp_fi3d)
2057	ENDIF
2058
2059	IF (o_dmc%flag(iff)<=lev_files(iff)) THEN
2060	zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) +
2061	$ dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
2062	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2063	$o_dmc%name,itau_w,zx_tmp_fi3d)
2064	ENDIF
2065
2066	IF (o_ref_liq%flag(iff)<=lev_files(iff)) THEN
2067	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2068	$o_ref_liq%name,itau_w,ref_liq)
2069	ENDIF
2070
2071	IF (o_ref_ice%flag(iff)<=lev_files(iff)) THEN
2072	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2073	$o_ref_ice%name,itau_w,ref_ice)
2074	ENDIF
2075
2076	if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR.
2077	$ RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR.
2078	$ RCFC12_per.NE.RCFC12_act) THEN
2079
2080	IF (o_rsut4co2%flag(iff)<=lev_files(iff)) THEN
2081	zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 )
2082	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2083	$o_rsut4co2%name,itau_w,zx_tmp_fi2d)
2084	ENDIF
2085
2086	IF (o_rlut4co2%flag(iff)<=lev_files(iff)) THEN
2087	zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 )
2088	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2089	$o_rlut4co2%name,itau_w,zx_tmp_fi2d)
2090	ENDIF
2091
2092	IF (o_rsutcs4co2%flag(iff)<=lev_files(iff)) THEN
2093	zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 )
2094	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2095	$o_rsutcs4co2%name,itau_w,zx_tmp_fi2d)
2096	ENDIF
2097
2098	IF (o_rlutcs4co2%flag(iff)<=lev_files(iff)) THEN
2099	zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 )
2100	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2101	$o_rlutcs4co2%name,itau_w,zx_tmp_fi2d)
2102	ENDIF
2103
2104	IF (o_rsu4co2%flag(iff)<=lev_files(iff)) THEN
2105	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2106	$o_rsu4co2%name,itau_w,swupp)
2107	ENDIF
2108
2109	IF (o_rlu4co2%flag(iff)<=lev_files(iff)) THEN
2110	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2111	$o_rlu4co2%name,itau_w,lwupp)
2112	ENDIF
2113
2114	IF (o_rsucs4co2%flag(iff)<=lev_files(iff)) THEN
2115	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2116	$o_rsucs4co2%name,itau_w,swup0p)
2117	ENDIF
2118
2119	IF (o_rlucs4co2%flag(iff)<=lev_files(iff)) THEN
2120	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2121	$o_rlucs4co2%name,itau_w,lwup0p)
2122	ENDIF
2123
2124	IF (o_rsd4co2%flag(iff)<=lev_files(iff)) THEN
2125	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2126	$o_rsd4co2%name,itau_w,swdnp)
2127	ENDIF
2128
2129	IF (o_rld4co2%flag(iff)<=lev_files(iff)) THEN
2130	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2131	$o_rld4co2%name,itau_w,lwdnp)
2132	ENDIF
2133
2134	IF (o_rsdcs4co2%flag(iff)<=lev_files(iff)) THEN
2135	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2136	$o_rsdcs4co2%name,itau_w,swdn0p)
2137	ENDIF
2138
2139	IF (o_rldcs4co2%flag(iff)<=lev_files(iff)) THEN
2140	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2141	$o_rldcs4co2%name,itau_w,lwdn0p)
2142	ENDIF
2143
2144	endif
2145
2146	if (nqtot.GE.3) THEN
2147	DO iq=3,nqtot
2148	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
2149	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2150	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
2151	ENDIF
2152	ENDDO
2153	endif
2154
2155	if (ok_sync) then
2156	c$OMP MASTER
2157	call histsync(nid_files(iff))
2158	c$OMP END MASTER
2159	endif
2160
2161	ENDIF ! clef_files
2162
2163	ENDDO ! iff=1,nfiles

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

Download in other formats: