Context Navigation

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

Last change on this file since 1645 was 1638, checked in by idelkadi, 12 years ago
Introduction du declenchement stochastique de la convection
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 76.6 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_files(iff).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>=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)) THEN
680	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
681	$o_prw%name,itau_w,prw)
682	ENDIF
683
684	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
685	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
686	$o_s_pblh%name,itau_w,s_pblh)
687	ENDIF
688
689	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
690	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
691	$o_s_pblt%name,itau_w,s_pblt)
692	ENDIF
693
694	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
695	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
696	$o_s_lcl%name,itau_w,s_lcl)
697	ENDIF
698
699	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
700	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
701	$o_s_therm%name,itau_w,s_therm)
702	ENDIF
703
704	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
705	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
706	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
707	! $o_s_capCL%name,itau_w,s_capCL)
708	! ENDIF
709
710	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
711	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
712	! $o_s_oliqCL%name,itau_w,s_oliqCL)
713	! ENDIF
714
715	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
716	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
717	! $o_s_cteiCL%name,itau_w,s_cteiCL)
718	! ENDIF
719
720	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
721	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
722	! $o_s_trmb1%name,itau_w,s_trmb1)
723	! ENDIF
724
725	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
726	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
727	! $o_s_trmb2%name,itau_w,s_trmb2)
728	! ENDIF
729
730	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
731	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
732	! $o_s_trmb3%name,itau_w,s_trmb3)
733	! ENDIF
734
735	! Champs interpolles sur des niveaux de pression
736
737	ll=0
738	DO k=1, nlevSTD
739	bb2=clevSTD(k)
740	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
741	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
742	$ bb2.EQ."100".OR.
743	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
744
745	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
746	ll=ll+1
747	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
748	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
749	$o_uSTDlevs(ll)%name,
750	& itau_w,uwriteSTD(:,k,iff))
751	ENDIF
752
753	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
754	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
755	$o_vSTDlevs(ll)%name,
756	& itau_w,vwriteSTD(:,k,iff))
757	ENDIF
758
759	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
760	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
761	$o_wSTDlevs(ll)%name,
762	& itau_w,wwriteSTD(:,k,iff))
763	ENDIF
764
765	IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
766	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
767	$o_zSTDlevs(ll)%name,
768	& itau_w,phiwriteSTD(:,k,iff))
769	ENDIF
770
771	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
772	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
773	$o_qSTDlevs(ll)%name,
774	& itau_w, qwriteSTD(:,k,iff))
775	ENDIF
776
777	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
778	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
779	$o_tSTDlevs(ll)%name,
780	& itau_w, twriteSTD(:,k,iff))
781	ENDIF
782
783	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
784	ENDDO
785
786	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
787	DO i=1, klon
788	IF (pctsrf(i,is_oce).GT.epsfra.OR.
789	$ pctsrf(i,is_sic).GT.epsfra) THEN
790	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
791	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
792	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
793	ELSE
794	zx_tmp_fi2d(i) = 273.15
795	ENDIF
796	ENDDO
797	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
798	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
799	ENDIF
800
801	! Couplage convection-couche limite
802	IF (iflag_con.GE.3) THEN
803	IF (iflag_coupl>=1) THEN
804	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
805	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
806	$o_ale_bl%name,itau_w,ale_bl)
807	ENDIF
808	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
809	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
810	$o_alp_bl%name,itau_w,alp_bl)
811	ENDIF
812	ENDIF !iflag_coupl>=1
813	ENDIF !(iflag_con.GE.3)
814
815	! Wakes
816	IF (iflag_con.EQ.3) THEN
817	IF (iflag_wake>=1) THEN
818	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
819	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
820	$o_ale_wk%name,itau_w,ale_wake)
821	ENDIF
822	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
823	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
824	$o_alp_wk%name,itau_w,alp_wake)
825	ENDIF
826
827	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
828	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
829	$o_ale%name,itau_w,ale)
830	ENDIF
831	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
832	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
833	$o_alp%name,itau_w,alp)
834	ENDIF
835	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
836	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
837	$o_cin%name,itau_w,cin)
838	ENDIF
839	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
840	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
841	$o_WAPE%name,itau_w,wake_pe)
842	ENDIF
843	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
844	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
845	$o_wake_h%name,itau_w,wake_h)
846	ENDIF
847
848	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
849	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
850	$o_wake_s%name,itau_w,wake_s)
851	ENDIF
852
853	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
854	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
855	$o_wake_deltat%name,
856	$ itau_w,wake_deltat)
857	ENDIF
858
859	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
860	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
861	$o_wake_deltaq%name,
862	$ itau_w,wake_deltaq)
863	ENDIF
864
865	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
866	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
867	s o_wake_omg%name,itau_w,wake_omg)
868	ENDIF
869
870	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
871	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
872	& /pdtphys
873	CALL histwrite_phy(nid_files(iff),
874	$clef_stations(iff),
875	& o_dtwak%name,itau_w,zx_tmp_fi3d)
876	ENDIF
877
878	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
879	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
880	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
881	& o_dqwak%name,itau_w,zx_tmp_fi3d)
882	ENDIF
883	ENDIF ! iflag_wake>=1
884
885	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
886	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
887	$o_Vprecip%name,itau_w,Vprecip)
888	ENDIF
889
890	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
891	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
892	$o_ftd%name,itau_w,ftd)
893	ENDIF
894
895	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
896	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
897	$o_fqd%name,itau_w,fqd)
898	ENDIF
899	ENDIF !(iflag_con.EQ.3)
900
901	!!! nrlmd le 10/04/2012
902
903	IF (o_zlcl_th%flag(iff)<=lev_files(iff)) THEN
904	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
905	s o_zlcl_th%name,itau_w,zlcl_th)
906	ENDIF
907
908	IF (o_fraca0%flag(iff)<=lev_files(iff)) THEN
909	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
910	s o_fraca0%name,itau_w,fraca0)
911	ENDIF
912
913	IF (o_w0%flag(iff)<=lev_files(iff)) THEN
914	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
915	s o_w0%name,itau_w,w0)
916	ENDIF
917
918	IF (o_w_conv%flag(iff)<=lev_files(iff)) THEN
919	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
920	s o_w_conv%name,itau_w,w_conv)
921	ENDIF
922
923	IF (o_therm_tke_max0%flag(iff)<=lev_files(iff)) THEN
924	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
925	s o_therm_tke_max0%name,itau_w,therm_tke_max0)
926	ENDIF
927
928	IF (o_env_tke_max0%flag(iff)<=lev_files(iff)) THEN
929	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
930	s o_env_tke_max0%name,itau_w,env_tke_max0)
931	ENDIF
932
933	IF (o_n2%flag(iff)<=lev_files(iff)) THEN
934	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
935	s o_n2%name,itau_w,n2)
936	ENDIF
937
938	IF (o_s2%flag(iff)<=lev_files(iff)) THEN
939	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
940	s o_s2%name,itau_w,s2)
941	ENDIF
942
943	IF (o_proba_notrig%flag(iff)<=lev_files(iff)) THEN
944	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
945	s o_proba_notrig%name,itau_w,proba_notrig)
946	ENDIF
947
948	IF (o_random_notrig%flag(iff)<=lev_files(iff)) THEN
949	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
950	s o_random_notrig%name,itau_w,random_notrig)
951	ENDIF
952
953	IF (o_ale_bl_stat%flag(iff)<=lev_files(iff)) THEN
954	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
955	s o_ale_bl_stat%name,itau_w,ale_bl_stat)
956	ENDIF
957
958	IF (o_ale_bl_trig%flag(iff)<=lev_files(iff)) THEN
959	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
960	s o_ale_bl_trig%name,itau_w,ale_bl_trig)
961	ENDIF
962
963	IF (o_therm_tke_max%flag(iff)<=lev_files(iff)) THEN
964	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
965	s o_therm_tke_max%name,itau_w,therm_tke_max)
966	ENDIF
967
968	IF (o_env_tke_max%flag(iff)<=lev_files(iff)) THEN
969	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
970	s o_env_tke_max%name,itau_w,env_tke_max)
971	ENDIF
972
973	IF (o_alp_bl_det%flag(iff)<=lev_files(iff)) THEN
974	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
975	s o_alp_bl_det%name,itau_w,alp_bl_det)
976	ENDIF
977
978	IF (o_alp_bl_fluct_m%flag(iff)<=lev_files(iff)) THEN
979	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
980	s o_alp_bl_fluct_m%name,itau_w,alp_bl_fluct_m)
981	ENDIF
982
983	IF (o_alp_bl_fluct_tke%flag(iff)<=lev_files(iff)) THEN
984	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
985	s o_alp_bl_fluct_tke%name,itau_w,alp_bl_fluct_tke)
986	ENDIF
987
988	IF (o_alp_bl_conv%flag(iff)<=lev_files(iff)) THEN
989	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
990	s o_alp_bl_conv%name,itau_w,alp_bl_conv)
991	ENDIF
992
993	IF (o_alp_bl_stat%flag(iff)<=lev_files(iff)) THEN
994	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
995	s o_alp_bl_stat%name,itau_w,alp_bl_stat)
996	ENDIF
997
998	!!! fin nrlmd le 10/04/2012
999
1000	IF (type_ocean=='slab ') THEN
1001	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
1002	$ CALL histwrite_phy(
1003	$ nid_files(iff),clef_stations(iff),
1004	$o_slab_bils%name,itau_w,slab_wfbils)
1005
1006	ENDIF !type_ocean == force/slab
1007
1008	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
1009	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1010	s o_weakinv%name,itau_w,weak_inversion)
1011	ENDIF
1012
1013	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
1014	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1015	$o_dthmin%name,itau_w,dthmin)
1016	ENDIF
1017
1018	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
1019	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1020	$o_cldtau%name,itau_w,cldtau)
1021	ENDIF
1022
1023	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
1024	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1025	$o_cldemi%name,itau_w,cldemi)
1026	ENDIF
1027
1028	IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
1029	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1030	s o_pr_con_l%name,itau_w,pmflxr(:,1:klev))
1031	ENDIF
1032
1033	IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
1034	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1035	s o_pr_con_i%name,itau_w,pmflxs(:,1:klev))
1036	ENDIF
1037
1038	IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
1039	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1040	s o_pr_lsc_l%name,itau_w,prfl(:,1:klev))
1041	ENDIF
1042
1043	IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
1044	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1045	s o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
1046	ENDIF
1047
1048	IF (o_re%flag(iff)<=lev_files(iff)) THEN
1049	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1050	$o_re%name,itau_w,re)
1051	ENDIF
1052
1053	IF (o_fl%flag(iff)<=lev_files(iff)) THEN
1054	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1055	$o_fl%name,itau_w,fl)
1056	ENDIF
1057
1058
1059
1060	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
1061	DO i=1, klon
1062	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
1063	ENDDO
1064	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1065	$o_rh2m%name,itau_w,zx_tmp_fi2d)
1066	ENDIF
1067
1068	IF (.NOT.clef_stations(iff)) THEN
1069	IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN
1070	DO i=1, klon
1071	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
1072	ENDDO
1073	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1074	$o_rh2m_min%name,
1075	s itau_w,zx_tmp_fi2d)
1076	ENDIF
1077
1078	IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN
1079	DO i=1, klon
1080	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
1081	ENDDO
1082	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1083	$o_rh2m_max%name,
1084	s itau_w,zx_tmp_fi2d)
1085	ENDIF
1086	ENDIF
1087
1088
1089	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
1090	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1091	$o_qsat2m%name,itau_w,qsat2m)
1092	ENDIF
1093
1094	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
1095	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1096	$o_tpot%name,itau_w,tpot)
1097	ENDIF
1098
1099	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
1100	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1101	$o_tpote%name,itau_w,tpote)
1102	ENDIF
1103
1104	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
1105	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
1106	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1107	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
1108	ENDIF
1109
1110	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
1111	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
1112	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1113	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
1114	ENDIF
1115
1116	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
1117	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1118	s o_LWdownOR%name,itau_w,sollwdown)
1119	ENDIF
1120
1121	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
1122	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1123	$o_snowl%name,itau_w,snow_lsc)
1124	ENDIF
1125
1126	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
1127	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1128	s o_solldown%name,itau_w,sollwdown)
1129	ENDIF
1130
1131	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
1132	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1133	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
1134	ENDIF
1135
1136	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
1137	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1138	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
1139	ENDIF
1140
1141	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
1142	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1143	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
1144	ENDIF
1145
1146	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
1147	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1148	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
1149	ENDIF
1150
1151	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
1152	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1153	$o_rugs%name,itau_w,zxrugs)
1154	ENDIF
1155
1156	! OD550 per species
1157	IF (new_aod .and. (.not. aerosol_couple)) THEN
1158	IF (ok_ade.OR.ok_aie) THEN
1159
1160	IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN
1161	CALL histwrite_phy(nid_files(iff),
1162	$clef_stations(iff),
1163	$o_od550aer%name,itau_w,
1164	$ od550aer)
1165	ENDIF
1166	IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN
1167	CALL histwrite_phy(nid_files(iff),
1168	$clef_stations(iff),
1169	$o_od865aer%name,itau_w,
1170	$ od865aer)
1171	ENDIF
1172	IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN
1173	CALL histwrite_phy(nid_files(iff),
1174	$clef_stations(iff),
1175	$o_absvisaer%name,itau_w,
1176	$ absvisaer)
1177	ENDIF
1178	IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN
1179	CALL histwrite_phy(nid_files(iff),
1180	$clef_stations(iff),
1181	$o_od550lt1aer%name,itau_w,
1182	$ od550lt1aer)
1183	ENDIF
1184
1185	IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN
1186	CALL histwrite_phy(nid_files(iff),
1187	$clef_stations(iff),
1188	$o_sconcso4%name,itau_w,
1189	$ sconcso4)
1190	ENDIF
1191	IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN
1192	CALL histwrite_phy(nid_files(iff),
1193	$clef_stations(iff),
1194	$o_sconcoa%name,itau_w,
1195	$ sconcoa)
1196	ENDIF
1197	IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN
1198	CALL histwrite_phy(nid_files(iff),
1199	$clef_stations(iff),
1200	$o_sconcbc%name,itau_w,
1201	$ sconcbc)
1202	ENDIF
1203	IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN
1204	CALL histwrite_phy(nid_files(iff),
1205	$clef_stations(iff),
1206	$o_sconcss%name,itau_w,
1207	$ sconcss)
1208	ENDIF
1209	IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN
1210	CALL histwrite_phy(nid_files(iff),
1211	$clef_stations(iff),
1212	$o_sconcdust%name,itau_w,
1213	$ sconcdust)
1214	ENDIF
1215
1216	IF (o_concso4%flag(iff)<=lev_files(iff)) THEN
1217	CALL histwrite_phy(nid_files(iff),
1218	$clef_stations(iff),
1219	$o_concso4%name,itau_w,
1220	$ concso4)
1221	ENDIF
1222	IF (o_concoa%flag(iff)<=lev_files(iff)) THEN
1223	CALL histwrite_phy(nid_files(iff),
1224	$clef_stations(iff),
1225	$o_concoa%name,itau_w,
1226	$ concoa)
1227	ENDIF
1228	IF (o_concbc%flag(iff)<=lev_files(iff)) THEN
1229	CALL histwrite_phy(nid_files(iff),
1230	$clef_stations(iff),
1231	$o_concbc%name,itau_w,
1232	$ concbc)
1233	ENDIF
1234	IF (o_concss%flag(iff)<=lev_files(iff)) THEN
1235	CALL histwrite_phy(nid_files(iff),
1236	$clef_stations(iff),
1237	$o_concss%name,itau_w,
1238	$ concss)
1239	ENDIF
1240	IF (o_concdust%flag(iff)<=lev_files(iff)) THEN
1241	CALL histwrite_phy(nid_files(iff),
1242	$clef_stations(iff),
1243	$o_concdust%name,itau_w,
1244	$ concdust)
1245	ENDIF
1246
1247	IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN
1248	CALL histwrite_phy(nid_files(iff),
1249	$clef_stations(iff),
1250	$o_loadso4%name,itau_w,
1251	$ loadso4)
1252	ENDIF
1253	IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN
1254	CALL histwrite_phy(nid_files(iff),
1255	$clef_stations(iff),
1256	$o_loadoa%name,itau_w,
1257	$ loadoa)
1258	ENDIF
1259	IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN
1260	CALL histwrite_phy(nid_files(iff),
1261	$clef_stations(iff),
1262	$o_loadbc%name,itau_w,
1263	$ loadbc)
1264	ENDIF
1265	IF (o_loadss%flag(iff)<=lev_files(iff)) THEN
1266	CALL histwrite_phy(nid_files(iff),
1267	$clef_stations(iff),
1268	$o_loadss%name,itau_w,
1269	$ loadss)
1270	ENDIF
1271	IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN
1272	CALL histwrite_phy(nid_files(iff),
1273	$clef_stations(iff),
1274	$o_loaddust%name,itau_w,
1275	$ loaddust)
1276	ENDIF
1277
1278	DO naero = 1, naero_spc
1279	IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
1280	CALL histwrite_phy(nid_files(iff),
1281	$clef_stations(iff),
1282	$ o_tausumaero(naero)%name,itau_w,
1283	$ tausum_aero(:,2,naero) )
1284	ENDIF
1285	END DO
1286	endif
1287	ENDIF
1288
1289	IF (ok_ade) THEN
1290	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
1291	CALL histwrite_phy(nid_files(iff),
1292	$clef_stations(iff),
1293	$o_topswad%name,itau_w,
1294	$ topswad_aero)
1295	ENDIF
1296	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
1297	CALL histwrite_phy(nid_files(iff),
1298	$clef_stations(iff),
1299	$o_solswad%name,itau_w,
1300	$ solswad_aero)
1301	ENDIF
1302
1303	!====MS forcing diagnostics
1304	if (new_aod) then
1305	IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
1306	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1307	$o_swtoaas_nat%name,itau_w,
1308	$ topsw_aero(:,1))
1309	ENDIF
1310
1311	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
1312	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1313	$o_swsrfas_nat%name,itau_w,
1314	$ solsw_aero(:,1))
1315	ENDIF
1316
1317	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
1318	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1319	$o_swtoacs_nat%name,itau_w,
1320	$ topsw0_aero(:,1))
1321	ENDIF
1322
1323	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
1324	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1325	$o_swsrfcs_nat%name,itau_w,
1326	$ solsw0_aero(:,1))
1327	ENDIF
1328
1329	!ant
1330	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
1331	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1332	$o_swtoaas_ant%name,itau_w,
1333	$ topsw_aero(:,2))
1334	ENDIF
1335
1336	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
1337	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1338	$o_swsrfas_ant%name,itau_w,
1339	$ solsw_aero(:,2))
1340	ENDIF
1341
1342	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
1343	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1344	$o_swtoacs_ant%name,itau_w,
1345	$ topsw0_aero(:,2))
1346	ENDIF
1347
1348	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
1349	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1350	$o_swsrfcs_ant%name,itau_w,
1351	$ solsw0_aero(:,2))
1352	ENDIF
1353
1354	!cf
1355
1356	if (.not. aerosol_couple) then
1357	IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
1358	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1359	$o_swtoacf_nat%name,itau_w,
1360	$ topswcf_aero(:,1))
1361	ENDIF
1362
1363	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
1364	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1365	$o_swsrfcf_nat%name,itau_w,
1366	$ solswcf_aero(:,1))
1367	ENDIF
1368
1369	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
1370	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1371	$o_swtoacf_ant%name,itau_w,
1372	$ topswcf_aero(:,2))
1373	ENDIF
1374
1375	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
1376	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1377	$o_swsrfcf_ant%name,itau_w,
1378	$ solswcf_aero(:,2))
1379	ENDIF
1380
1381	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
1382	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1383	$o_swtoacf_zero%name,itau_w,
1384	$ topswcf_aero(:,3))
1385	ENDIF
1386
1387	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
1388	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1389	$o_swsrfcf_zero%name,itau_w,
1390	$ solswcf_aero(:,3))
1391	ENDIF
1392	endif
1393
1394	endif ! new_aod
1395	!====MS forcing diagnostics
1396
1397	ENDIF
1398
1399	IF (ok_aie) THEN
1400	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
1401	CALL histwrite_phy(nid_files(iff),
1402	$clef_stations(iff),
1403	$o_topswai%name,itau_w,
1404	$ topswai_aero)
1405	ENDIF
1406	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
1407	CALL histwrite_phy(nid_files(iff),
1408	$clef_stations(iff),
1409	$o_solswai%name,itau_w,
1410	$ solswai_aero)
1411	ENDIF
1412	IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN
1413	CALL histwrite_phy(nid_files(iff),
1414	$clef_stations(iff),
1415	$o_scdnc%name,itau_w,
1416	$ scdnc)
1417	ENDIF
1418	IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN
1419	CALL histwrite_phy(nid_files(iff),
1420	$clef_stations(iff),
1421	$o_cldncl%name,itau_w,
1422	$ cldncl)
1423	ENDIF
1424	IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN
1425	CALL histwrite_phy(nid_files(iff),
1426	$clef_stations(iff),
1427	$o_reffclws%name,itau_w,
1428	$ reffclws)
1429	ENDIF
1430	IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN
1431	CALL histwrite_phy(nid_files(iff),
1432	$clef_stations(iff),
1433	$o_reffclwc%name,itau_w,
1434	$ reffclwc)
1435	ENDIF
1436	IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN
1437	CALL histwrite_phy(nid_files(iff),
1438	$clef_stations(iff),
1439	$o_cldnvi%name,itau_w,
1440	$ cldnvi)
1441	ENDIF
1442	IF (o_lcc%flag(iff)<=lev_files(iff)) THEN
1443	CALL histwrite_phy(nid_files(iff),
1444	$clef_stations(iff),
1445	$o_lcc%name,itau_w,
1446	$ lcc)
1447	ENDIF
1448	IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN
1449	CALL histwrite_phy(nid_files(iff),
1450	$clef_stations(iff),
1451	$o_lcc3d%name,itau_w,
1452	$ lcc3d)
1453	ENDIF
1454	IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN
1455	CALL histwrite_phy(nid_files(iff),
1456	$clef_stations(iff),
1457	$o_lcc3dcon%name,itau_w,
1458	$ lcc3dcon)
1459	ENDIF
1460	IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN
1461	CALL histwrite_phy(nid_files(iff),
1462	$clef_stations(iff),
1463	$o_lcc3dstra%name,itau_w,
1464	$ lcc3dstra)
1465	ENDIF
1466	IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN
1467	CALL histwrite_phy(nid_files(iff),
1468	$clef_stations(iff),
1469	$o_reffclwtop%name,itau_w,
1470	$ reffclwtop)
1471	ENDIF
1472	ENDIF
1473
1474	! Champs 3D:
1475	IF (ok_ade .OR. ok_aie) then
1476	IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN
1477	CALL histwrite_phy(nid_files(iff),
1478	$clef_stations(iff),
1479	$o_ec550aer%name,itau_w,
1480	& ec550aer)
1481	ENDIF
1482	ENDIF
1483
1484	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
1485	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1486	$o_lwcon%name,itau_w,flwc)
1487	ENDIF
1488
1489	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
1490	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1491	$o_iwcon%name,itau_w,fiwc)
1492	ENDIF
1493
1494	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
1495	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1496	$o_temp%name,itau_w,t_seri)
1497	ENDIF
1498
1499	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
1500	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1501	$o_theta%name,itau_w,theta)
1502	ENDIF
1503
1504	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
1505	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1506	$o_ovapinit%name,itau_w,
1507	$ qx(:,:,ivap))
1508	ENDIF
1509
1510	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
1511	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1512	$ o_ovap%name,itau_w,q_seri)
1513	ENDIF
1514
1515	IF (o_oliq%flag(iff)<=lev_files(iff)) THEN
1516	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1517	$ o_oliq%name,itau_w,ql_seri)
1518	ENDIF
1519
1520	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
1521	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1522	$o_geop%name,itau_w,zphi)
1523	ENDIF
1524
1525	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
1526	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1527	$o_vitu%name,itau_w,u_seri)
1528	ENDIF
1529
1530	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
1531	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1532	$o_vitv%name,itau_w,v_seri)
1533	ENDIF
1534
1535	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
1536	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1537	$o_vitw%name,itau_w,omega)
1538	ENDIF
1539
1540	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
1541	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1542	$o_pres%name,itau_w,pplay)
1543	ENDIF
1544
1545	IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
1546	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1547	$o_paprs%name,
1548	s itau_w,paprs(:,1:klev))
1549	ENDIF
1550
1551	IF (o_zfull%flag(iff)<=lev_files(iff)) THEN
1552	DO i=1, klon
1553	zx_tmp_fi3d1(i,1)= pphis(i)/RG
1554	!020611 zx_tmp_fi3d(i,1)= pphis(i)/RG
1555	ENDDO
1556	DO k=1, klev
1557	!020611 DO k=1, klev-1
1558	DO i=1, klon
1559	!020611 zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (t_seri(i,k) RD
1560	zx_tmp_fi3d1(i,k+1)= zx_tmp_fi3d1(i,k) - (t_seri(i,k) RD
1561	$ (paprs(i,k+1) - paprs(i,k))) / ( pplay(i,k) * RG )
1562	ENDDO
1563	ENDDO
1564	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1565	$o_zfull%name,itau_w,zx_tmp_fi3d1(:,2:klevp1))
1566	!020611 $o_zfull%name,itau_w,zx_tmp_fi3d)
1567	ENDIF
1568
1569	IF (o_zhalf%flag(iff)<=lev_files(iff)) THEN
1570	DO i=1, klon
1571	zx_tmp_fi3d(i,1)= pphis(i)/RG - (
1572	$ (t_seri(i,1)+zxtsol(i))/2. RD
1573	$ (pplay(i,1) - paprs(i,1)))/( (paprs(i,1)+pplay(i,1))/2. * RG)
1574	ENDDO
1575	DO k=1, klev-1
1576	DO i=1, klon
1577	zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (
1578	$ (t_seri(i,k)+t_seri(i,k+1))/2. RD
1579	$ (pplay(i,k+1) - pplay(i,k))) / ( paprs(i,k) * RG )
1580	ENDDO
1581	ENDDO
1582	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1583	$o_zhalf%name,itau_w,zx_tmp_fi3d)
1584	ENDIF
1585
1586	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
1587	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1588	$o_rneb%name,itau_w,cldfra)
1589	ENDIF
1590
1591	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
1592	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1593	$o_rnebcon%name,itau_w,rnebcon)
1594	ENDIF
1595
1596	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
1597	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1598	$o_rhum%name,itau_w,zx_rh)
1599	ENDIF
1600
1601	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
1602	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1603	$o_ozone%name, itau_w,
1604	$ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1605	ENDIF
1606
1607	IF (o_ozone_light%flag(iff)<=lev_files(iff) .and.
1608	$ read_climoz == 2) THEN
1609	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1610	$o_ozone_light%name, itau_w,
1611	$ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1612	ENDIF
1613
1614	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
1615	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1616	$o_dtphy%name,itau_w,d_t)
1617	ENDIF
1618
1619	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
1620	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1621	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
1622	ENDIF
1623
1624	DO nsrf=1, nbsrf
1625	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1626	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1627	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1628	s o_albe_srf(nsrf)%name,itau_w,
1629	$ zx_tmp_fi2d)
1630	ENDIF
1631
1632	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1633	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1634	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1635	s o_rugs_srf(nsrf)%name,itau_w,
1636	$ zx_tmp_fi2d)
1637	ENDIF
1638
1639	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1640	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1641	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1642	s o_ages_srf(nsrf)%name,itau_w
1643	$ ,zx_tmp_fi2d)
1644	ENDIF
1645	ENDDO !nsrf=1, nbsrf
1646
1647	IF (o_alb1%flag(iff)<=lev_files(iff)) THEN
1648	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1649	$o_alb1%name,itau_w,albsol1)
1650	ENDIF
1651
1652	IF (o_alb2%flag(iff)<=lev_files(iff)) THEN
1653	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1654	$o_alb2%name,itau_w,albsol2)
1655	ENDIF
1656
1657	!FH Sorties pour la couche limite
1658	if (iflag_pbl>1) then
1659	zx_tmp_fi3d=0.
1660	do nsrf=1,nbsrf
1661	do k=1,klev
1662	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
1663	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1664	enddo
1665	enddo
1666	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
1667	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1668	$o_tke%name,itau_w,zx_tmp_fi3d)
1669	ENDIF
1670
1671	IF (.NOT.clef_stations(iff)) THEN
1672	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
1673	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1674	s o_tke_max%name,itau_w,zx_tmp_fi3d)
1675	ENDIF
1676	ENDIF
1677	endif
1678
1679	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
1680	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1681	$o_kz%name,itau_w,coefh)
1682	ENDIF
1683
1684	IF (.NOT.clef_stations(iff)) THEN
1685	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
1686	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1687	$o_kz_max%name,itau_w,coefh)
1688	ENDIF
1689	ENDIF
1690
1691	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
1692	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1693	$o_clwcon%name,itau_w,clwcon0)
1694	ENDIF
1695
1696	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
1697	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1698	$o_dtdyn%name,itau_w,d_t_dyn)
1699	ENDIF
1700
1701	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
1702	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1703	$o_dqdyn%name,itau_w,d_q_dyn)
1704	ENDIF
1705
1706	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
1707	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1708	$o_dudyn%name,itau_w,d_u_dyn)
1709	ENDIF
1710
1711	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
1712	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1713	$o_dvdyn%name,itau_w,d_v_dyn)
1714	ENDIF
1715
1716	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
1717	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1718	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1719	$o_dtcon%name,itau_w,zx_tmp_fi3d)
1720	ENDIF
1721
1722	if(iflag_thermals.eq.1)then
1723	IF (o_tntc%flag(iff)<=lev_files(iff)) THEN
1724	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys +
1725	$ d_t_ajsb(1:klon,1:klev)/pdtphys
1726	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1727	$o_tntc%name,itau_w,zx_tmp_fi3d)
1728	ENDIF
1729	else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then
1730	IF (o_tntc%flag(iff)<=lev_files(iff)) THEN
1731	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys +
1732	$ d_t_ajs(1:klon,1:klev)/pdtphys +
1733	$ d_t_wake(1:klon,1:klev)/pdtphys
1734	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1735	$o_tntc%name,itau_w,zx_tmp_fi3d)
1736	ENDIF
1737	endif
1738
1739	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
1740	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1741	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1742	$o_ducon%name,itau_w,zx_tmp_fi3d)
1743	ENDIF
1744
1745	IF (o_dvcon%flag(iff)<=lev_files(iff)) THEN
1746	zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1747	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1748	$o_dvcon%name,itau_w,zx_tmp_fi3d)
1749	ENDIF
1750
1751	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
1752	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1753	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1754	$o_dqcon%name,itau_w,zx_tmp_fi3d)
1755	ENDIF
1756
1757	if(iflag_thermals.eq.1)then
1758	IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN
1759	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1760	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1761	$o_tnhusc%name,itau_w,zx_tmp_fi3d)
1762	ENDIF
1763	else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then
1764	IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN
1765	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys +
1766	$ d_q_ajs(1:klon,1:klev)/pdtphys +
1767	$ d_q_wake(1:klon,1:klev)/pdtphys
1768	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1769	$o_tnhusc%name,itau_w,zx_tmp_fi3d)
1770	ENDIF
1771	endif
1772
1773	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
1774	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1775	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1776	$o_dtlsc%name,itau_w,zx_tmp_fi3d)
1777	ENDIF
1778
1779	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
1780	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1781	$ d_t_eva(1:klon,1:klev))/pdtphys
1782	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1783	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
1784	ENDIF
1785
1786	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
1787	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1788	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1789	$o_dqlsc%name,itau_w,zx_tmp_fi3d)
1790	ENDIF
1791
1792	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1793	! Sorties specifiques a la separation thermiques/non thermiques
1794	if (iflag_thermals>=1) then
1795
1796	IF (o_dtlscth%flag(iff)<=lev_files(iff)) THEN
1797	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1798	CALL
1799	s histwrite_phy(nid_files(iff),clef_stations(iff),
1800	$o_dtlscth%name,itau_w,zx_tmp_fi3d)
1801	ENDIF
1802
1803	IF (o_dtlscst%flag(iff)<=lev_files(iff)) THEN
1804	zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1805	CALL
1806	s histwrite_phy(nid_files(iff),clef_stations(iff),
1807	$o_dtlscst%name,itau_w,zx_tmp_fi3d)
1808	ENDIF
1809
1810	IF (o_dqlscth%flag(iff)<=lev_files(iff)) THEN
1811	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1812	CALL
1813	s histwrite_phy(nid_files(iff),clef_stations(iff),
1814	$o_dqlscth%name,itau_w,zx_tmp_fi3d)
1815	ENDIF
1816
1817	IF (o_dqlscst%flag(iff)<=lev_files(iff)) THEN
1818	zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1819	CALL
1820	s histwrite_phy(nid_files(iff),clef_stations(iff),
1821	$o_dqlscst%name,itau_w,zx_tmp_fi3d)
1822	ENDIF
1823
1824	IF (o_plulth%flag(iff)<=lev_files(iff)) THEN
1825	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1826	$o_plulth%name,itau_w,plul_th)
1827	ENDIF
1828
1829	IF (o_plulst%flag(iff)<=lev_files(iff)) THEN
1830	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1831	$o_plulst%name,itau_w,plul_st)
1832	ENDIF
1833
1834	do k=1,klev
1835	do i=1,klon
1836	if (ptconvth(i,k)) then
1837	zx_tmp_fi3d(i,k)=1.
1838	else
1839	zx_tmp_fi3d(i,k)=0.
1840	endif
1841	enddo
1842	enddo
1843	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1844	CALL
1845	s histwrite_phy(nid_files(iff),clef_stations(iff),
1846	$o_ptconvth%name,itau_w,zx_tmp_fi3d)
1847	ENDIF
1848
1849	do i=1,klon
1850	zx_tmp_fi2d(1:klon)=lmax_th(:)
1851	enddo
1852	IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN
1853	CALL
1854	s histwrite_phy(nid_files(iff),clef_stations(iff),
1855	$o_lmaxth%name,itau_w,zx_tmp_fi2d)
1856	ENDIF
1857
1858	endif ! iflag_thermals>=1
1859
1860	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1861
1862	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
1863	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1864	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1865	$o_dtvdf%name,itau_w,zx_tmp_fi3d)
1866	ENDIF
1867
1868	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
1869	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1870	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1871	$o_dqvdf%name,itau_w,zx_tmp_fi3d)
1872	ENDIF
1873
1874	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
1875	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1876	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1877	$o_dteva%name,itau_w,zx_tmp_fi3d)
1878	ENDIF
1879
1880	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
1881	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1882	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1883	$o_dqeva%name,itau_w,zx_tmp_fi3d)
1884	ENDIF
1885
1886	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
1887	zpt_conv = 0.
1888	where (ptconv) zpt_conv = 1.
1889	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1890	$o_ptconv%name,itau_w,zpt_conv)
1891	ENDIF
1892
1893	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
1894	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1895	$o_ratqs%name,itau_w,ratqs)
1896	ENDIF
1897
1898	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
1899	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys -
1900	$ d_t_ajsb(1:klon,1:klev)/pdtphys
1901	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1902	$o_dtthe%name,itau_w,zx_tmp_fi3d)
1903	ENDIF
1904
1905	IF (iflag_thermals>=1) THEN
1906	IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN
1907	! Pour l instant 0 a y reflichir pour les thermiques
1908	zx_tmp_fi2d=0.
1909	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1910	$o_ftime_th%name,
1911	s itau_w,zx_tmp_fi2d)
1912	ENDIF
1913
1914	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
1915	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1916	$o_f_th%name,itau_w,fm_therm)
1917	ENDIF
1918
1919	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
1920	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1921	$o_e_th%name,itau_w,entr_therm)
1922	ENDIF
1923
1924	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
1925	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1926	$o_w_th%name,itau_w,zw2)
1927	ENDIF
1928
1929	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
1930	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1931	$o_q_th%name,itau_w,zqasc)
1932	ENDIF
1933
1934
1935	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
1936	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1937	$o_a_th%name,itau_w,fraca)
1938	ENDIF
1939
1940	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
1941	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1942	$o_d_th%name,itau_w,detr_therm)
1943	ENDIF
1944
1945	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1946	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1947	$o_f0_th%name,itau_w,f0)
1948	ENDIF
1949
1950	IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN
1951	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1952	$o_zmax_th%name,itau_w,zmax_th)
1953	ENDIF
1954
1955	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
1956	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys -
1957	$ d_q_ajsb(1:klon,1:klev)/pdtphys
1958	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1959	$o_dqthe%name,itau_w,zx_tmp_fi3d)
1960	ENDIF
1961
1962	ENDIF !iflag_thermals
1963
1964	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
1965	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1966	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1967	$o_dtajs%name,itau_w,zx_tmp_fi3d)
1968	ENDIF
1969
1970	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
1971	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1972	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1973	$o_dqajs%name,itau_w,zx_tmp_fi3d)
1974	ENDIF
1975
1976	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
1977	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1978	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1979	$o_dtswr%name,itau_w,zx_tmp_fi3d)
1980	ENDIF
1981
1982	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
1983	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1984	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1985	$o_dtsw0%name,itau_w,zx_tmp_fi3d)
1986	ENDIF
1987
1988	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
1989	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1990	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1991	$o_dtlwr%name,itau_w,zx_tmp_fi3d)
1992	ENDIF
1993
1994	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
1995	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1996	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1997	$o_dtlw0%name,itau_w,zx_tmp_fi3d)
1998	ENDIF
1999
2000	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
2001	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
2002	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2003	$o_dtec%name,itau_w,zx_tmp_fi3d)
2004	ENDIF
2005
2006	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
2007	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
2008	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2009	$o_duvdf%name,itau_w,zx_tmp_fi3d)
2010	ENDIF
2011
2012	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
2013	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
2014	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2015	$o_dvvdf%name,itau_w,zx_tmp_fi3d)
2016	ENDIF
2017
2018	IF (ok_orodr) THEN
2019	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
2020	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
2021	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2022	$o_duoro%name,itau_w,zx_tmp_fi3d)
2023	ENDIF
2024
2025	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
2026	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
2027	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2028	$o_dvoro%name,itau_w,zx_tmp_fi3d)
2029	ENDIF
2030
2031	IF (o_dtoro%flag(iff)<=lev_files(iff)) THEN
2032	zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
2033	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2034	$o_dtoro%name,itau_w,zx_tmp_fi3d)
2035	ENDIF
2036	ENDIF
2037
2038	IF (ok_orolf) THEN
2039	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
2040	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
2041	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2042	$o_dulif%name,itau_w,zx_tmp_fi3d)
2043	ENDIF
2044
2045	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
2046	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
2047	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2048	$o_dvlif%name,itau_w,zx_tmp_fi3d)
2049	ENDIF
2050
2051	IF (o_dtlif%flag(iff)<=lev_files(iff)) THEN
2052	zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
2053	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2054	$o_dtlif%name,itau_w,zx_tmp_fi3d)
2055	ENDIF
2056	ENDIF
2057
2058	IF (ok_hines) THEN
2059	IF (o_duhin%flag(iff)<=lev_files(iff)) THEN
2060	zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys
2061	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2062	$o_duhin%name,itau_w,zx_tmp_fi3d)
2063	ENDIF
2064
2065
2066	IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN
2067	zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys
2068	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2069	$o_dvhin%name,itau_w,zx_tmp_fi3d)
2070	ENDIF
2071
2072	IF (o_dthin%flag(iff)<=lev_files(iff)) THEN
2073	zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
2074	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2075	$o_dthin%name,itau_w,zx_tmp_fi3d)
2076	ENDIF
2077	ENDIF
2078
2079	IF (o_rsu%flag(iff)<=lev_files(iff)) THEN
2080	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2081	$o_rsu%name,itau_w,swup)
2082	ENDIF
2083	IF (o_rsd%flag(iff)<=lev_files(iff)) THEN
2084	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2085	$o_rsd%name,itau_w,swdn)
2086	ENDIF
2087	IF (o_rlu%flag(iff)<=lev_files(iff)) THEN
2088	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2089	$o_rlu%name,itau_w,lwup)
2090	ENDIF
2091	IF (o_rld%flag(iff)<=lev_files(iff)) THEN
2092	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2093	$o_rld%name,itau_w,lwdn)
2094	ENDIF
2095
2096	IF (o_rsucs%flag(iff)<=lev_files(iff)) THEN
2097	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2098	$o_rsucs%name,itau_w,swup0)
2099	ENDIF
2100	IF (o_rsdcs%flag(iff)<=lev_files(iff)) THEN
2101	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2102	$o_rsdcs%name,itau_w,swdn0)
2103	ENDIF
2104	IF (o_rlucs%flag(iff)<=lev_files(iff)) THEN
2105	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2106	$o_rlucs%name,itau_w,lwup0)
2107	ENDIF
2108	IF (o_rldcs%flag(iff)<=lev_files(iff)) THEN
2109	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2110	$o_rldcs%name,itau_w,lwdn0)
2111	ENDIF
2112
2113	IF (o_tnt%flag(iff)<=lev_files(iff)) THEN
2114	zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+
2115	$d_t_dyn(1:klon,1:klev)
2116	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2117	$o_tnt%name,itau_w,zx_tmp_fi3d)
2118	ENDIF
2119
2120	IF (o_tntr%flag(iff)<=lev_files(iff)) THEN
2121	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY -
2122	$cool(1:klon,1:klev)/RDAY
2123	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2124	$o_tntr%name,itau_w,zx_tmp_fi3d)
2125	ENDIF
2126
2127	IF (o_tntscpbl%flag(iff)<=lev_files(iff)) THEN
2128	zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+
2129	$ d_t_eva(1:klon,1:klev)+
2130	$ d_t_vdf(1:klon,1:klev))/pdtphys
2131	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2132	$o_tntscpbl%name,itau_w,zx_tmp_fi3d)
2133	ENDIF
2134
2135	IF (o_tnhus%flag(iff)<=lev_files(iff)) THEN
2136	zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+
2137	$d_q_dyn(1:klon,1:klev)
2138	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2139	$o_tnhus%name,itau_w,zx_tmp_fi3d)
2140	ENDIF
2141
2142	IF (o_tnhusscpbl%flag(iff)<=lev_files(iff)) THEN
2143	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+
2144	$ d_q_eva(1:klon,1:klev)/pdtphys
2145	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2146	$o_tnhusscpbl%name,itau_w,zx_tmp_fi3d)
2147	ENDIF
2148
2149	IF (o_evu%flag(iff)<=lev_files(iff)) THEN
2150	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2151	$o_evu%name,itau_w,coefm)
2152	ENDIF
2153
2154	IF (o_h2o%flag(iff)<=lev_files(iff)) THEN
2155	zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+
2156	$ ql_seri(1:klon,1:klev)
2157	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2158	$o_h2o%name,itau_w,zx_tmp_fi3d)
2159	ENDIF
2160
2161	IF (o_mcd%flag(iff)<=lev_files(iff)) THEN
2162	zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+
2163	$ dnwd0(1:klon,1:klev))
2164	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2165	$o_mcd%name,itau_w,zx_tmp_fi3d)
2166	ENDIF
2167
2168	IF (o_dmc%flag(iff)<=lev_files(iff)) THEN
2169	zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) +
2170	$ dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
2171	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2172	$o_dmc%name,itau_w,zx_tmp_fi3d)
2173	ENDIF
2174
2175	IF (o_ref_liq%flag(iff)<=lev_files(iff)) THEN
2176	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2177	$o_ref_liq%name,itau_w,ref_liq)
2178	ENDIF
2179
2180	IF (o_ref_ice%flag(iff)<=lev_files(iff)) THEN
2181	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2182	$o_ref_ice%name,itau_w,ref_ice)
2183	ENDIF
2184
2185	if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR.
2186	$ RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR.
2187	$ RCFC12_per.NE.RCFC12_act) THEN
2188
2189	IF (o_rsut4co2%flag(iff)<=lev_files(iff)) THEN
2190	zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 )
2191	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2192	$o_rsut4co2%name,itau_w,zx_tmp_fi2d)
2193	ENDIF
2194
2195	IF (o_rlut4co2%flag(iff)<=lev_files(iff)) THEN
2196	zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 )
2197	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2198	$o_rlut4co2%name,itau_w,zx_tmp_fi2d)
2199	ENDIF
2200
2201	IF (o_rsutcs4co2%flag(iff)<=lev_files(iff)) THEN
2202	zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 )
2203	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2204	$o_rsutcs4co2%name,itau_w,zx_tmp_fi2d)
2205	ENDIF
2206
2207	IF (o_rlutcs4co2%flag(iff)<=lev_files(iff)) THEN
2208	zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 )
2209	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2210	$o_rlutcs4co2%name,itau_w,zx_tmp_fi2d)
2211	ENDIF
2212
2213	IF (o_rsu4co2%flag(iff)<=lev_files(iff)) THEN
2214	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2215	$o_rsu4co2%name,itau_w,swupp)
2216	ENDIF
2217
2218	IF (o_rlu4co2%flag(iff)<=lev_files(iff)) THEN
2219	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2220	$o_rlu4co2%name,itau_w,lwupp)
2221	ENDIF
2222
2223	IF (o_rsucs4co2%flag(iff)<=lev_files(iff)) THEN
2224	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2225	$o_rsucs4co2%name,itau_w,swup0p)
2226	ENDIF
2227
2228	IF (o_rlucs4co2%flag(iff)<=lev_files(iff)) THEN
2229	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2230	$o_rlucs4co2%name,itau_w,lwup0p)
2231	ENDIF
2232
2233	IF (o_rsd4co2%flag(iff)<=lev_files(iff)) THEN
2234	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2235	$o_rsd4co2%name,itau_w,swdnp)
2236	ENDIF
2237
2238	IF (o_rld4co2%flag(iff)<=lev_files(iff)) THEN
2239	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2240	$o_rld4co2%name,itau_w,lwdnp)
2241	ENDIF
2242
2243	IF (o_rsdcs4co2%flag(iff)<=lev_files(iff)) THEN
2244	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2245	$o_rsdcs4co2%name,itau_w,swdn0p)
2246	ENDIF
2247
2248	IF (o_rldcs4co2%flag(iff)<=lev_files(iff)) THEN
2249	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2250	$o_rldcs4co2%name,itau_w,lwdn0p)
2251	ENDIF
2252
2253	endif
2254
2255	if (nqtot.GE.3) THEN
2256	DO iq=3,nqtot
2257	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
2258	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
2259	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
2260	ENDIF
2261	ENDDO
2262	endif
2263
2264	if (ok_sync) then
2265	c$OMP MASTER
2266	call histsync(nid_files(iff))
2267	c$OMP END MASTER
2268	endif
2269
2270	ENDIF ! clef_files
2271
2272	ENDDO ! iff=1,nfiles

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

Download in other formats: