Context Navigation

phys_output_write.h @ 1114

Last change on this file since 1114 was 1114, checked in by jghattas, 15 years ago

Creation du module infotrac:

contient les variables de advtrac.h
contient la subroutine iniadvtrac renommer en infotrac_init
le nombre des traceurs est lu dans tracer.def en dynamique (ou par default ou recu par INCA)
ce module est utilise dans la dynamique et la physique
contient aussi la variable nbtr qui avant etait stockee dans dimphy

Le fichier advtrac.h n'existe plus.
La compilation ne prend plus en compte le nombre de traceur.

/JG

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

File size: 41.1 KB

Line
1	itau_w = itau_phy + itap
2
3	DO iff=1,nfiles
4
5	IF (clef_files(iff)) THEN
6	ndex2d = 0
7	ndex3d = 0
8
9	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
10	IF (o_phis%flag(iff)<=lev_files(iff)) THEN
11	CALL histwrite_phy(nid_files(iff),
12	$ o_phis%name,itau_w,pphis)
13	ENDIF
14
15	IF (o_aire%flag(iff)<=lev_files(iff)) THEN
16	CALL histwrite_phy(nid_files(iff),o_aire%name,itau_w,airephy)
17	ENDIF
18
19	IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN
20	DO i=1, klon
21	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
22	ENDDO
23	CALL histwrite_phy(nid_files(iff),
24	$ o_contfracATM%name,itau_w,zx_tmp_fi2d)
25	ENDIF
26
27	IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN
28	CALL histwrite_phy(nid_files(iff),o_contfracOR%name,itau_w,
29	$ pctsrf(:,is_ter))
30	ENDIF
31
32	IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN
33	CALL histwrite_phy(nid_files(iff),
34	$ o_aireTER%name,itau_w,paire_ter)
35	ENDIF
36
37	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
38
39	IF (o_flat%flag(iff)<=lev_files(iff)) THEN
40	CALL histwrite_phy(nid_files(iff),o_flat%name,itau_w,zxfluxlat)
41	ENDIF
42
43	IF (o_slp%flag(iff)<=lev_files(iff)) THEN
44	CALL histwrite_phy(nid_files(iff),o_slp%name,itau_w,slp)
45	ENDIF
46
47	IF (o_tsol%flag(iff)<=lev_files(iff)) THEN
48	CALL histwrite_phy(nid_files(iff),o_tsol%name,itau_w,zxtsol)
49	ENDIF
50
51	IF (o_t2m%flag(iff)<=lev_files(iff)) THEN
52	CALL histwrite_phy(nid_files(iff),o_t2m%name,itau_w,zt2m)
53	ENDIF
54
55	IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN
56	CALL histwrite_phy(nid_files(iff),o_t2m_min%name,itau_w,zt2m)
57	ENDIF
58
59	IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN
60	CALL histwrite_phy(nid_files(iff),o_t2m_max%name,itau_w,zt2m)
61	ENDIF
62
63	IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN
64	DO i=1, klon
65	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
66	ENDDO
67	CALL histwrite_phy(nid_files(iff),
68	s o_wind10m%name,itau_w,zx_tmp_fi2d)
69	ENDIF
70
71	IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN
72	DO i=1, klon
73	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
74	ENDDO
75	CALL histwrite_phy(nid_files(iff),o_wind10max%name,
76	$ itau_w,zx_tmp_fi2d)
77	ENDIF
78
79	IF (o_sicf%flag(iff)<=lev_files(iff)) THEN
80	DO i = 1, klon
81	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
82	ENDDO
83	CALL histwrite_phy(nid_files(iff),
84	$ o_sicf%name,itau_w,zx_tmp_fi2d)
85	ENDIF
86
87	IF (o_q2m%flag(iff)<=lev_files(iff)) THEN
88	CALL histwrite_phy(nid_files(iff),o_q2m%name,itau_w,zq2m)
89	ENDIF
90
91	IF (o_u10m%flag(iff)<=lev_files(iff)) THEN
92	CALL histwrite_phy(nid_files(iff),o_u10m%name,itau_w,zu10m)
93	ENDIF
94
95	IF (o_v10m%flag(iff)<=lev_files(iff)) THEN
96	CALL histwrite_phy(nid_files(iff),o_v10m%name,itau_w,zv10m)
97	ENDIF
98
99	IF (o_psol%flag(iff)<=lev_files(iff)) THEN
100	DO i = 1, klon
101	zx_tmp_fi2d(i) = paprs(i,1)
102	ENDDO
103	CALL histwrite_phy(nid_files(iff),
104	s o_psol%name,itau_w,zx_tmp_fi2d)
105	ENDIF
106
107	IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN
108	CALL histwrite_phy(nid_files(iff),o_qsurf%name,itau_w,zxqsurf)
109	ENDIF
110
111	if (.not. ok_veget) then
112	IF (o_qsol%flag(iff)<=lev_files(iff)) THEN
113	CALL histwrite_phy(nid_files(iff),o_qsol%name,itau_w,qsol)
114	ENDIF
115	endif
116
117	IF (o_precip%flag(iff)<=lev_files(iff)) THEN
118	DO i = 1, klon
119	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
120	ENDDO
121	CALL histwrite_phy(nid_files(iff),o_precip%name,
122	s itau_w,zx_tmp_fi2d)
123	ENDIF
124
125	IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN
126	CALL histwrite_phy(nid_files(iff),o_ndayrain%name,
127	s itau_w,nday_rain)
128	ENDIF
129
130	IF (o_plul%flag(iff)<=lev_files(iff)) THEN
131	DO i = 1, klon
132	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
133	ENDDO
134	CALL histwrite_phy(nid_files(iff),o_plul%name,itau_w,zx_tmp_fi2d)
135	ENDIF
136
137	IF (o_pluc%flag(iff)<=lev_files(iff)) THEN
138	DO i = 1, klon
139	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
140	ENDDO
141	CALL histwrite_phy(nid_files(iff),o_pluc%name,itau_w,zx_tmp_fi2d)
142	ENDIF
143
144	IF (o_snow%flag(iff)<=lev_files(iff)) THEN
145	CALL histwrite_phy(nid_files(iff),o_snow%name,itau_w,snow_fall)
146	ENDIF
147
148	IF (o_evap%flag(iff)<=lev_files(iff)) THEN
149	CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap)
150	ENDIF
151
152	IF (o_tops%flag(iff)<=lev_files(iff)) THEN
153	CALL histwrite_phy(nid_files(iff),o_tops%name,itau_w,topsw)
154	ENDIF
155
156	IF (o_tops0%flag(iff)<=lev_files(iff)) THEN
157	CALL histwrite_phy(nid_files(iff),o_tops0%name,itau_w,topsw0)
158	ENDIF
159
160	IF (o_topl%flag(iff)<=lev_files(iff)) THEN
161	CALL histwrite_phy(nid_files(iff),o_topl%name,itau_w,toplw)
162	ENDIF
163
164	IF (o_topl0%flag(iff)<=lev_files(iff)) THEN
165	CALL histwrite_phy(nid_files(iff),o_topl0%name,itau_w,toplw0)
166	ENDIF
167
168	IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN
169	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
170	CALL histwrite_phy(nid_files(iff),o_SWupTOA%name,
171	s itau_w,zx_tmp_fi2d)
172	ENDIF
173
174	IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN
175	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
176	CALL histwrite_phy(nid_files(iff),
177	$ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d)
178	ENDIF
179
180	IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN
181	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
182	CALL histwrite_phy(nid_files(iff),
183	s o_SWdnTOA%name,itau_w,zx_tmp_fi2d)
184	ENDIF
185
186	IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN
187	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
188	CALL histwrite_phy(nid_files(iff),
189	$ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d)
190	ENDIF
191
192	IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
193	CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200)
194	ENDIF
195
196	IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN
197	CALL histwrite_phy(nid_files(iff),
198	s o_SWup200clr%name,itau_w,SWup200clr)
199	ENDIF
200
201	IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN
202	CALL histwrite_phy(nid_files(iff),o_SWdn200%name,itau_w,SWdn200)
203	ENDIF
204
205	IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN
206	CALL histwrite_phy(nid_files(iff),
207	s o_SWdn200clr%name,itau_w,SWdn200clr)
208	ENDIF
209
210	IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN
211	CALL histwrite_phy(nid_files(iff),o_LWup200%name,itau_w,LWup200)
212	ENDIF
213
214	IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN
215	CALL histwrite_phy(nid_files(iff),
216	s o_LWup200clr%name,itau_w,LWup200clr)
217	ENDIF
218
219	IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN
220	CALL histwrite_phy(nid_files(iff),
221	s o_LWdn200%name,itau_w,zx_tmp_fi2d)
222	ENDIF
223
224	IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN
225	CALL histwrite_phy(nid_files(iff),
226	s o_LWdn200clr%name,itau_w,zx_tmp_fi2d)
227	ENDIF
228
229	IF (o_sols%flag(iff)<=lev_files(iff)) THEN
230	CALL histwrite_phy(nid_files(iff),o_sols%name,itau_w,solsw)
231	ENDIF
232
233	IF (o_sols0%flag(iff)<=lev_files(iff)) THEN
234	CALL histwrite_phy(nid_files(iff),o_sols0%name,itau_w,solsw0)
235	ENDIF
236
237	IF (o_soll%flag(iff)<=lev_files(iff)) THEN
238	CALL histwrite_phy(nid_files(iff),o_soll%name,itau_w,sollw)
239	ENDIF
240
241	IF (o_radsol%flag(iff)<=lev_files(iff)) THEN
242	CALL histwrite_phy(nid_files(iff),o_radsol%name,itau_w,radsol)
243	ENDIF
244
245	IF (o_soll0%flag(iff)<=lev_files(iff)) THEN
246	CALL histwrite_phy(nid_files(iff),o_soll0%name,itau_w,sollw0)
247	ENDIF
248
249	IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN
250	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
251	CALL histwrite_phy(nid_files(iff),
252	s o_SWupSFC%name,itau_w,zx_tmp_fi2d)
253	ENDIF
254
255	IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN
256	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
257	CALL histwrite_phy(nid_files(iff),
258	$ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d)
259	ENDIF
260
261	IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN
262	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
263	CALL histwrite_phy(nid_files(iff),
264	$ o_SWdnSFC%name,itau_w,zx_tmp_fi2d)
265	ENDIF
266
267	IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
268	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
269	CALL histwrite_phy(nid_files(iff),
270	$ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d)
271	ENDIF
272
273	IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN
274	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
275	CALL histwrite_phy(nid_files(iff),
276	$ o_LWupSFC%name,itau_w,zx_tmp_fi2d)
277	ENDIF
278
279	IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN
280	CALL histwrite_phy(nid_files(iff),
281	$ o_LWdnSFC%name,itau_w,sollwdown)
282	ENDIF
283
284	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
285	IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN
286	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
287	CALL histwrite_phy(nid_files(iff),
288	$ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d)
289	ENDIF
290
291	IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
292	CALL histwrite_phy(nid_files(iff),
293	$ o_LWdnSFCclr%name,itau_w,sollwdownclr)
294	ENDIF
295
296	IF (o_bils%flag(iff)<=lev_files(iff)) THEN
297	CALL histwrite_phy(nid_files(iff),o_bils%name,itau_w,bils)
298	ENDIF
299
300	IF (o_sens%flag(iff)<=lev_files(iff)) THEN
301	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
302	CALL histwrite_phy(nid_files(iff),o_sens%name,itau_w,zx_tmp_fi2d)
303	ENDIF
304
305	IF (o_fder%flag(iff)<=lev_files(iff)) THEN
306	CALL histwrite_phy(nid_files(iff),o_fder%name,itau_w,fder)
307	ENDIF
308
309	IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN
310	CALL histwrite_phy(nid_files(iff),o_ffonte%name,itau_w,zxffonte)
311	ENDIF
312
313	IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN
314	CALL histwrite_phy(nid_files(iff),
315	$ o_fqcalving%name,itau_w,zxfqcalving)
316	ENDIF
317
318	IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN
319	CALL histwrite_phy(nid_files(iff),
320	$ o_fqfonte%name,itau_w,zxfqfonte)
321	ENDIF
322
323	DO nsrf = 1, nbsrf
324	! IF(nsrf.GE.2) THEN
325	IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
326	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
327	CALL histwrite_phy(nid_files(iff),
328	$ o_pourc_srf(nsrf)%name,itau_w,
329	$ zx_tmp_fi2d)
330	ENDIF
331
332	IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
333	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
334	CALL histwrite_phy(nid_files(iff),
335	$ o_fract_srf(nsrf)%name,itau_w,
336	$ zx_tmp_fi2d)
337	ENDIF
338	! ENDIF !nsrf.GT.2
339
340	IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
341	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
342	CALL histwrite_phy(nid_files(iff),
343	$ o_taux_srf(nsrf)%name,itau_w,
344	$ zx_tmp_fi2d)
345	ENDIF
346
347	IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
348	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
349	CALL histwrite_phy(nid_files(iff),
350	$ o_tauy_srf(nsrf)%name,itau_w,
351	$ zx_tmp_fi2d)
352	ENDIF
353
354	IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
355	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
356	CALL histwrite_phy(nid_files(iff),
357	$ o_tsol_srf(nsrf)%name,itau_w,
358	$ zx_tmp_fi2d)
359	ENDIF
360
361	IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
362	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
363	CALL histwrite_phy(nid_files(iff),o_u10m_srf(nsrf)%name,
364	$ itau_w,zx_tmp_fi2d)
365	ENDIF
366
367	IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
368	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
369	CALL histwrite_phy(nid_files(iff),o_v10m_srf(nsrf)%name,
370	$ itau_w,zx_tmp_fi2d)
371	ENDIF
372
373	IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
374	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
375	CALL histwrite_phy(nid_files(iff),o_t2m_srf(nsrf)%name,
376	$ itau_w,zx_tmp_fi2d)
377	ENDIF
378
379	IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
380	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
381	CALL histwrite_phy(nid_files(iff),
382	$ o_sens_srf(nsrf)%name,itau_w,
383	$ zx_tmp_fi2d)
384	ENDIF
385
386	IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
387	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
388	CALL histwrite_phy(nid_files(iff),
389	$ o_lat_srf(nsrf)%name,itau_w,
390	$ zx_tmp_fi2d)
391	ENDIF
392
393	IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
394	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
395	CALL histwrite_phy(nid_files(iff),
396	$ o_flw_srf(nsrf)%name,itau_w,
397	$ zx_tmp_fi2d)
398	ENDIF
399
400	IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
401	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
402	CALL histwrite_phy(nid_files(iff),
403	$ o_fsw_srf(nsrf)%name,itau_w,
404	$ zx_tmp_fi2d)
405	ENDIF
406
407	IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
408	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
409	CALL histwrite_phy(nid_files(iff),
410	$ o_wbils_srf(nsrf)%name,itau_w,
411	$ zx_tmp_fi2d)
412	ENDIF
413
414	IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
415	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
416	CALL histwrite_phy(nid_files(iff),
417	$ o_wbilo_srf(nsrf)%name,itau_w,
418	$ zx_tmp_fi2d)
419	ENDIF
420
421	if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
422	IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
423	CALL histwrite_phy(nid_files(iff),
424	$ o_tke_srf(nsrf)%name,itau_w,
425	$ pbl_tke(:,1:klev,nsrf))
426	ENDIF
427
428	IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
429	CALL histwrite_phy(nid_files(iff),
430	$ o_tke_max_srf(nsrf)%name,itau_w,
431	$ pbl_tke(:,1:klev,nsrf))
432	ENDIF
433	endif
434	ENDDO
435
436	IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN
437	CALL histwrite_phy(nid_files(iff),o_cdrm%name,itau_w,cdragm)
438	ENDIF
439
440	IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN
441	CALL histwrite_phy(nid_files(iff),o_cdrh%name,itau_w,cdragh)
442	ENDIF
443
444	IF (o_cldl%flag(iff)<=lev_files(iff)) THEN
445	CALL histwrite_phy(nid_files(iff),o_cldl%name,itau_w,cldl)
446	ENDIF
447
448	IF (o_cldm%flag(iff)<=lev_files(iff)) THEN
449	CALL histwrite_phy(nid_files(iff),o_cldm%name,itau_w,cldm)
450	ENDIF
451
452	IF (o_cldh%flag(iff)<=lev_files(iff)) THEN
453	CALL histwrite_phy(nid_files(iff),o_cldh%name,itau_w,cldh)
454	ENDIF
455
456	IF (o_cldt%flag(iff)<=lev_files(iff)) THEN
457	CALL histwrite_phy(nid_files(iff),o_cldt%name,
458	& itau_w,cldt*100)
459	ENDIF
460
461	IF (o_cldq%flag(iff)<=lev_files(iff)) THEN
462	CALL histwrite_phy(nid_files(iff),o_cldq%name,itau_w,cldq)
463	ENDIF
464
465	IF (o_lwp%flag(iff)<=lev_files(iff)) THEN
466	zx_tmp_fi2d(1:klon) = flwp(1:klon)
467	CALL histwrite_phy(nid_files(iff),
468	s o_lwp%name,itau_w,zx_tmp_fi2d)
469	ENDIF
470
471	IF (o_iwp%flag(iff)<=lev_files(iff)) THEN
472	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
473	CALL histwrite_phy(nid_files(iff),
474	s o_iwp%name,itau_w,zx_tmp_fi2d)
475	ENDIF
476
477	IF (o_ue%flag(iff)<=lev_files(iff)) THEN
478	CALL histwrite_phy(nid_files(iff),o_ue%name,itau_w,ue)
479	ENDIF
480
481	IF (o_ve%flag(iff)<=lev_files(iff)) THEN
482	CALL histwrite_phy(nid_files(iff),o_ve%name,itau_w,ve)
483	ENDIF
484
485	IF (o_uq%flag(iff)<=lev_files(iff)) THEN
486	CALL histwrite_phy(nid_files(iff),o_uq%name,itau_w,uq)
487	ENDIF
488
489	IF (o_vq%flag(iff)<=lev_files(iff)) THEN
490	CALL histwrite_phy(nid_files(iff),o_vq%name,itau_w,vq)
491	ENDIF
492
493	IF(iflag_con.GE.3) THEN ! sb
494	IF (o_cape%flag(iff)<=lev_files(iff)) THEN
495	CALL histwrite_phy(nid_files(iff),o_cape%name,itau_w,cape)
496	ENDIF
497
498	IF (o_pbase%flag(iff)<=lev_files(iff)) THEN
499	CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w,pbase)
500	ENDIF
501
502	IF (o_ptop%flag(iff)<=lev_files(iff)) THEN
503	CALL histwrite_phy(nid_files(iff),o_ptop%name,itau_w,ema_pct)
504	ENDIF
505
506	IF (o_fbase%flag(iff)<=lev_files(iff)) THEN
507	CALL histwrite_phy(nid_files(iff),o_fbase%name,itau_w,ema_cbmf)
508	ENDIF
509
510	IF (o_prw%flag(iff)<=lev_files(iff)) THEN
511	CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw)
512	ENDIF
513	ENDIF !iflag_con .GE. 3
514
515	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
516	CALL histwrite_phy(nid_files(iff),o_s_pblh%name,itau_w,s_pblh)
517	ENDIF
518
519	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
520	CALL histwrite_phy(nid_files(iff),o_s_pblt%name,itau_w,s_pblt)
521	ENDIF
522
523	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
524	CALL histwrite_phy(nid_files(iff),o_s_lcl%name,itau_w,s_lcl)
525	ENDIF
526
527	IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
528	CALL histwrite_phy(nid_files(iff),o_s_capCL%name,itau_w,s_capCL)
529	ENDIF
530
531	IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
532	CALL histwrite_phy(nid_files(iff),o_s_oliqCL%name,itau_w,s_oliqCL)
533	ENDIF
534
535	IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
536	CALL histwrite_phy(nid_files(iff),o_s_cteiCL%name,itau_w,s_cteiCL)
537	ENDIF
538
539	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
540	CALL histwrite_phy(nid_files(iff),o_s_therm%name,itau_w,s_therm)
541	ENDIF
542
543	IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
544	CALL histwrite_phy(nid_files(iff),o_s_trmb1%name,itau_w,s_trmb1)
545	ENDIF
546
547	IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
548	CALL histwrite_phy(nid_files(iff),o_s_trmb2%name,itau_w,s_trmb2)
549	ENDIF
550
551	IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
552	CALL histwrite_phy(nid_files(iff),o_s_trmb3%name,itau_w,s_trmb3)
553	ENDIF
554
555	! Champs interpolles sur des niveaux de pression
556	! if=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
557	! if=2 on ecrit w et ph 500 seulement au niv 1
558	! et u v sur 850 700 500 200
559	! if=3 on ecrit ph a 500 seulement au niv 1
560	! on ecrit u v t q a 850 700 500 200 au niv 3
561
562	ll=0
563	DO k=1, nlevSTD
564	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
565	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
566	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
567	$ bb2.EQ."500".OR.bb2.EQ."200") THEN
568
569	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
570	ll=ll+1
571	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
572	CALL histwrite_phy(nid_files(iff),o_uSTDlevs(ll)%name,
573	& itau_w,uwriteSTD(:,k,iff))
574	ENDIF
575
576	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
577	CALL histwrite_phy(nid_files(iff),o_vSTDlevs(ll)%name,
578	& itau_w,vwriteSTD(:,k,iff))
579	ENDIF
580
581	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
582	CALL histwrite_phy(nid_files(iff),o_wSTDlevs(ll)%name,
583	& itau_w,wwriteSTD(:,k,iff))
584	ENDIF
585
586	IF (o_phiSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
587	CALL histwrite_phy(nid_files(iff),o_phiSTDlevs(ll)%name,
588	& itau_w,phiwriteSTD(:,k,iff))
589	ENDIF
590
591	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
592	CALL histwrite_phy(nid_files(iff),o_qSTDlevs(ll)%name,
593	& itau_w, qwriteSTD(:,k,iff))
594	ENDIF
595
596	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
597	CALL histwrite_phy(nid_files(iff),o_tSTDlevs(ll)%name,
598	& itau_w, twriteSTD(:,k,iff))
599	ENDIF
600
601	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
602	ENDDO
603
604	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
605	DO i=1, klon
606	IF (pctsrf(i,is_oce).GT.epsfra.OR.
607	$ pctsrf(i,is_sic).GT.epsfra) THEN
608	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
609	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
610	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
611	ELSE
612	zx_tmp_fi2d(i) = 273.15
613	ENDIF
614	ENDDO
615	CALL histwrite_phy(nid_files(iff),
616	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
617	ENDIF
618
619	IF (type_ocean=='force ') THEN
620
621	IF (iflag_con.GE.3) THEN
622	IF (iflag_coupl.EQ.1) THEN
623	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
624	CALL histwrite_phy(nid_files(iff),o_ale_bl%name,itau_w,ale_bl)
625	ENDIF
626	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
627	CALL histwrite_phy(nid_files(iff),o_alp_bl%name,itau_w,alp_bl)
628	ENDIF
629	ENDIF !iflag_coupl.EQ.1
630	ENDIF !(iflag_con.GE.3)
631	IF (iflag_con.EQ.3) THEN
632	IF (iflag_wake.EQ.1) THEN
633	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
634	CALL histwrite_phy(nid_files(iff),o_ale_wk%name,itau_w,ale_wake)
635	ENDIF
636	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
637	CALL histwrite_phy(nid_files(iff),o_alp_wk%name,itau_w,alp_wake)
638	ENDIF
639
640	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
641	CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)
642	ENDIF
643	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
644	CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)
645	ENDIF
646	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
647	CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)
648	ENDIF
649	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
650	CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe)
651	ENDIF
652	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
653	CALL histwrite_phy(nid_files(iff),o_wake_h%name,itau_w,wake_h)
654	ENDIF
655
656	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
657	CALL histwrite_phy(nid_files(iff),o_wake_s%name,itau_w,wake_s)
658	ENDIF
659
660	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
661	CALL histwrite_phy(nid_files(iff),o_wake_deltat%name,
662	$ itau_w,wake_deltat)
663	ENDIF
664
665	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
666	CALL histwrite_phy(nid_files(iff),o_wake_deltaq%name,
667	$ itau_w,wake_deltaq)
668	ENDIF
669
670	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
671	CALL histwrite_phy(nid_files(iff),
672	s o_wake_omg%name,itau_w,wake_omg)
673	ENDIF
674	ENDIF
675
676	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
677	CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip)
678	ENDIF
679
680	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
681	CALL histwrite_phy(nid_files(iff),o_ftd%name,itau_w,ftd)
682	ENDIF
683
684	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
685	CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd)
686	ENDIF
687
688	ENDIF !(iflag_con.EQ.3)
689
690	ELSE IF (type_ocean=='slab ') THEN
691
692	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
693	$ CALL histwrite_phy(
694	$ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils)
695
696	ENDIF !type_ocean == force/slab
697
698	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
699	CALL histwrite_phy(nid_files(iff),
700	s o_weakinv%name,itau_w,weak_inversion)
701	ENDIF
702
703	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
704	CALL histwrite_phy(nid_files(iff),o_dthmin%name,itau_w,dthmin)
705	ENDIF
706
707	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
708	CALL histwrite_phy(nid_files(iff),o_cldtau%name,itau_w,cldtau)
709	ENDIF
710
711	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
712	CALL histwrite_phy(nid_files(iff),o_cldemi%name,itau_w,cldemi)
713	ENDIF
714
715	! IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
716	! CALL histwrite_phy(nid_files(iff),o_pmflxr%name,itau_w,pmflxr)
717	! ENDIF
718
719	! IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
720	! CALL histwrite_phy(nid_files(iff),o_pmflxs%name,itau_w,pmflxs)
721	! ENDIF
722
723	! IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
724	! CALL histwrite_phy(nid_files(iff),o_prfl%name,itau_w,prfl)
725	! ENDIF
726
727	! IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
728	! CALL histwrite_phy(nid_files(iff),o_psfl%name,itau_w,psfl)
729	! ENDIF
730
731	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
732	DO i=1, klon
733	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
734	ENDDO
735	CALL histwrite_phy(nid_files(iff),o_rh2m%name,itau_w,zx_tmp_fi2d)
736	ENDIF
737
738	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
739	CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m)
740	ENDIF
741
742	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
743	CALL histwrite_phy(nid_files(iff),o_tpot%name,itau_w,tpot)
744	ENDIF
745
746	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
747	CALL histwrite_phy(nid_files(iff),o_tpote%name,itau_w,tpote)
748	ENDIF
749
750	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
751	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
752	CALL histwrite_phy(nid_files(iff),
753	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
754	ENDIF
755
756	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
757	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
758	CALL histwrite_phy(nid_files(iff),
759	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
760	ENDIF
761
762	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
763	CALL histwrite_phy(nid_files(iff),
764	s o_LWdownOR%name,itau_w,sollwdown)
765	ENDIF
766
767	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
768	CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc)
769	ENDIF
770
771	IF (iflag_con.GE.3) THEN
772	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
773	CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)
774	ENDIF
775	ENDIF !(iflag_con.GE.3)
776
777	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
778	CALL histwrite_phy(nid_files(iff),
779	s o_solldown%name,itau_w,sollwdown)
780	ENDIF
781
782	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
783	CALL histwrite_phy(nid_files(iff),
784	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
785	ENDIF
786
787	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
788	CALL histwrite_phy(nid_files(iff),
789	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
790	ENDIF
791
792	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
793	CALL histwrite_phy(nid_files(iff),
794	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
795	ENDIF
796
797	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
798	CALL histwrite_phy(nid_files(iff),
799	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
800	ENDIF
801
802	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
803	CALL histwrite_phy(nid_files(iff),o_rugs%name,itau_w,zxrugs)
804	ENDIF
805
806	IF (ok_ade) THEN
807	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
808	CALL histwrite_phy(nid_files(iff),o_topswad%name,itau_w,topswad)
809	ENDIF
810	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
811	CALL histwrite_phy(nid_files(iff),o_solswad%name,itau_w,solswad)
812	ENDIF
813	ENDIF
814
815	IF (ok_aie) THEN
816	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
817	CALL histwrite_phy(nid_files(iff),o_topswai%name,itau_w,topswai)
818	ENDIF
819	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
820	CALL histwrite_phy(nid_files(iff),o_solswai%name,itau_w,solswai)
821	ENDIF
822	ENDIF
823
824	! Champs 3D:
825	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
826	CALL histwrite_phy(nid_files(iff),o_lwcon%name,itau_w,flwc)
827	ENDIF
828
829	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
830	CALL histwrite_phy(nid_files(iff),o_iwcon%name,itau_w,fiwc)
831	ENDIF
832
833	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
834	CALL histwrite_phy(nid_files(iff),o_temp%name,itau_w,t_seri)
835	ENDIF
836
837	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
838	CALL histwrite_phy(nid_files(iff),o_theta%name,itau_w,theta)
839	ENDIF
840
841	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
842	CALL histwrite_phy(nid_files(iff),o_ovap%name,itau_w,qx(:,:,ivap))
843	ENDIF
844
845	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
846	CALL histwrite_phy(nid_files(iff),o_geop%name,itau_w,zphi)
847	ENDIF
848
849	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
850	CALL histwrite_phy(nid_files(iff),o_vitu%name,itau_w,u_seri)
851	ENDIF
852
853	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
854	CALL histwrite_phy(nid_files(iff),o_vitv%name,itau_w,v_seri)
855	ENDIF
856
857	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
858	CALL histwrite_phy(nid_files(iff),o_vitw%name,itau_w,omega)
859	ENDIF
860
861	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
862	CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay)
863	ENDIF
864
865	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
866	CALL histwrite_phy(nid_files(iff),o_rneb%name,itau_w,cldfra)
867	ENDIF
868
869	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
870	CALL histwrite_phy(nid_files(iff),o_rnebcon%name,itau_w,rnebcon)
871	ENDIF
872
873	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
874	CALL histwrite_phy(nid_files(iff),o_rhum%name,itau_w,zx_rh)
875	ENDIF
876
877	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
878	DO k=1, klev
879	DO i=1, klon
880	zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
881	$ /(paprs(i,k)-paprs(i,k+1))
882	$ *(paprs(i,1)/101325.0)
883	ENDDO !i
884	ENDDO !k
885	CALL histwrite_phy(nid_files(iff),o_ozone%name,itau_w,zx_tmp_fi3d)
886	ENDIF
887
888	IF(iflag_con.GE.3) THEN
889	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
890	CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)
891	ENDIF
892	ENDIF !(iflag_con.GE.3)
893
894	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
895	CALL histwrite_phy(nid_files(iff),o_dtphy%name,itau_w,d_t)
896	ENDIF
897
898	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
899	CALL histwrite_phy(nid_files(iff),
900	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
901	ENDIF
902
903	DO nsrf=1, nbsrf
904	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
905	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
906	CALL histwrite_phy(nid_files(iff),
907	s o_albe_srf(nsrf)%name,itau_w,
908	$ zx_tmp_fi2d)
909	ENDIF
910
911	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
912	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
913	CALL histwrite_phy(nid_files(iff),
914	s o_rugs_srf(nsrf)%name,itau_w,
915	$ zx_tmp_fi2d)
916	ENDIF
917
918	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
919	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
920	CALL histwrite_phy(nid_files(iff),
921	s o_ages_srf(nsrf)%name,itau_w
922	$ ,zx_tmp_fi2d)
923	ENDIF
924	ENDDO !nsrf=1, nbsrf
925
926	IF (o_albs%flag(iff)<=lev_files(iff)) THEN
927	CALL histwrite_phy(nid_files(iff),o_albs%name,itau_w,albsol1)
928	ENDIF
929
930	IF (o_albslw%flag(iff)<=lev_files(iff)) THEN
931	CALL histwrite_phy(nid_files(iff),o_albslw%name,itau_w,albsol2)
932	ENDIF
933
934	!FH Sorties pour la couche limite
935	if (iflag_pbl>1) then
936	zx_tmp_fi3d=0.
937	do nsrf=1,nbsrf
938	do k=1,klev
939	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
940	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
941	enddo
942	enddo
943	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
944	CALL histwrite_phy(nid_files(iff),o_tke%name,itau_w,zx_tmp_fi3d)
945	ENDIF
946
947	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
948	CALL histwrite_phy(nid_files(iff),
949	s o_tke_max%name,itau_w,zx_tmp_fi3d)
950	ENDIF
951	endif
952
953	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
954	CALL histwrite_phy(nid_files(iff),o_kz%name,itau_w,coefh)
955	ENDIF
956
957	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
958	CALL histwrite_phy(nid_files(iff),o_kz_max%name,itau_w,coefh)
959	ENDIF
960
961	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
962	CALL histwrite_phy(nid_files(iff),o_clwcon%name,itau_w,clwcon0)
963	ENDIF
964
965	IF(iflag_con.GE.3) THEN
966	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
967	CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)
968	ENDIF
969
970	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
971	CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)
972	ENDIF
973
974	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
975	CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)
976	ENDIF
977	ENDIF !(iflag_con.GE.3)
978
979	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
980	CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn)
981	ENDIF
982
983	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
984	CALL histwrite_phy(nid_files(iff),o_dqdyn%name,itau_w,d_q_dyn)
985	ENDIF
986
987	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
988	CALL histwrite_phy(nid_files(iff),o_dudyn%name,itau_w,d_u_dyn)
989	ENDIF
990
991	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
992	CALL histwrite_phy(nid_files(iff),o_dvdyn%name,itau_w,d_v_dyn)
993	ENDIF
994
995	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
996	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
997	CALL histwrite_phy(nid_files(iff),o_dtcon%name,itau_w,zx_tmp_fi3d)
998	ENDIF
999
1000	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
1001	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1002	CALL histwrite_phy(nid_files(iff),o_ducon%name,itau_w,zx_tmp_fi3d)
1003	ENDIF
1004
1005	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
1006	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1007	CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d)
1008	ENDIF
1009
1010	IF(iflag_con.GE.3) THEN
1011	IF (iflag_wake == 1) THEN
1012	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
1013	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
1014	& /pdtphys
1015	CALL histwrite_phy(nid_files(iff),
1016	& o_dtwak%name,itau_w,zx_tmp_fi3d)
1017	ENDIF
1018
1019	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
1020	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
1021	CALL histwrite_phy(nid_files(iff),
1022	& o_dqwak%name,itau_w,zx_tmp_fi3d)
1023	ENDIF
1024	ENDIF
1025	ENDIF !(iflag_con.GE.3)
1026
1027	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
1028	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1029	CALL histwrite_phy(nid_files(iff),o_dtlsc%name,itau_w,zx_tmp_fi3d)
1030	ENDIF
1031
1032	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
1033	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1034	$ d_t_eva(1:klon,1:klev))/pdtphys
1035	CALL histwrite_phy(nid_files(iff),
1036	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
1037	ENDIF
1038
1039	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
1040	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1041	CALL histwrite_phy(nid_files(iff),o_dqlsc%name,itau_w,zx_tmp_fi3d)
1042	ENDIF
1043
1044	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
1045	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1046	CALL histwrite_phy(nid_files(iff),o_dtvdf%name,itau_w,zx_tmp_fi3d)
1047	ENDIF
1048
1049	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
1050	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1051	CALL histwrite_phy(nid_files(iff),o_dqvdf%name,itau_w,zx_tmp_fi3d)
1052	ENDIF
1053
1054	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
1055	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1056	CALL histwrite_phy(nid_files(iff),o_dteva%name,itau_w,zx_tmp_fi3d)
1057	ENDIF
1058
1059	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
1060	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1061	CALL histwrite_phy(nid_files(iff),o_dqeva%name,itau_w,zx_tmp_fi3d)
1062	ENDIF
1063
1064	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
1065	zpt_conv = 0.
1066	where (ptconv) zpt_conv = 1.
1067	CALL histwrite_phy(nid_files(iff),o_ptconv%name,itau_w,zpt_conv)
1068	ENDIF
1069
1070	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
1071	CALL histwrite_phy(nid_files(iff),o_ratqs%name,itau_w,ratqs)
1072	ENDIF
1073
1074	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
1075	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
1076	CALL histwrite_phy(nid_files(iff),o_dtthe%name,itau_w,zx_tmp_fi3d)
1077	ENDIF
1078
1079	IF (iflag_thermals.gt.1) THEN
1080	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
1081	CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm)
1082	ENDIF
1083
1084	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
1085	CALL histwrite_phy(nid_files(iff),o_e_th%name,itau_w,entr_therm)
1086	ENDIF
1087
1088	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
1089	CALL histwrite_phy(nid_files(iff),o_w_th%name,itau_w,zw2)
1090	ENDIF
1091
1092	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
1093	CALL histwrite_phy(nid_files(iff),o_q_th%name,itau_w,zqasc)
1094	ENDIF
1095
1096	IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN
1097	CALL histwrite_phy(nid_files(iff),
1098	s o_lambda_th%name,itau_w,lambda_th)
1099	ENDIF
1100
1101	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
1102	CALL histwrite_phy(nid_files(iff),o_a_th%name,itau_w,fraca)
1103	ENDIF
1104
1105	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
1106	CALL histwrite_phy(nid_files(iff),o_d_th%name,itau_w,detr_therm)
1107	ENDIF
1108
1109	ENDIF !iflag_thermals
1110
1111	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1112	CALL histwrite_phy(nid_files(iff),o_f0_th%name,itau_w,f0)
1113	ENDIF
1114
1115	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1116	CALL histwrite_phy(nid_files(iff),o_zmax_th%name,itau_w,zmax0)
1117	ENDIF
1118
1119	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
1120	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
1121	CALL histwrite_phy(nid_files(iff),o_dqthe%name,itau_w,zx_tmp_fi3d)
1122	ENDIF
1123
1124	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
1125	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1126	CALL histwrite_phy(nid_files(iff),o_dtajs%name,itau_w,zx_tmp_fi3d)
1127	ENDIF
1128
1129	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
1130	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1131	CALL histwrite_phy(nid_files(iff),o_dqajs%name,itau_w,zx_tmp_fi3d)
1132	ENDIF
1133
1134	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
1135	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1136	CALL histwrite_phy(nid_files(iff),o_dtswr%name,itau_w,zx_tmp_fi3d)
1137	ENDIF
1138
1139	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
1140	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1141	CALL histwrite_phy(nid_files(iff),o_dtsw0%name,itau_w,zx_tmp_fi3d)
1142	ENDIF
1143
1144	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
1145	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1146	CALL histwrite_phy(nid_files(iff),o_dtlwr%name,itau_w,zx_tmp_fi3d)
1147	ENDIF
1148
1149	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
1150	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1151	CALL histwrite_phy(nid_files(iff),o_dtlw0%name,itau_w,zx_tmp_fi3d)
1152	ENDIF
1153
1154	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
1155	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1156	CALL histwrite_phy(nid_files(iff),o_dtec%name,itau_w,zx_tmp_fi3d)
1157	ENDIF
1158
1159	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
1160	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1161	CALL histwrite_phy(nid_files(iff),o_duvdf%name,itau_w,zx_tmp_fi3d)
1162	ENDIF
1163
1164	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
1165	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1166	CALL histwrite_phy(nid_files(iff),o_dvvdf%name,itau_w,zx_tmp_fi3d)
1167	ENDIF
1168
1169	IF (ok_orodr) THEN
1170	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
1171	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1172	CALL histwrite_phy(nid_files(iff),o_duoro%name,itau_w,zx_tmp_fi3d)
1173	ENDIF
1174
1175	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
1176	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1177	CALL histwrite_phy(nid_files(iff),o_dvoro%name,itau_w,zx_tmp_fi3d)
1178	ENDIF
1179	ENDIF
1180
1181	IF (ok_orolf) THEN
1182	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
1183	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1184	CALL histwrite_phy(nid_files(iff),o_dulif%name,itau_w,zx_tmp_fi3d)
1185	ENDIF
1186
1187	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
1188	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1189	CALL histwrite_phy(nid_files(iff),o_dvlif%name,itau_w,zx_tmp_fi3d)
1190	ENDIF
1191	ENDIF
1192
1193	! IF (o_trac%flag(iff)<=lev_files(iff)) THEN
1194	if (nqtot.GE.3) THEN
1195	! DO iq=3,nqtot
1196	DO iq=3,4
1197	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
1198	CALL histwrite_phy(nid_files(iff),
1199	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
1200	ENDIF
1201	ENDDO
1202	endif
1203
1204	if (ok_sync) then
1205	c$OMP MASTER
1206	call histsync(nid_files(iff))
1207	c$OMP END MASTER
1208	endif
1209
1210	ENDIF ! clef_files
1211
1212	ENDDO ! iff=1,nfiles

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

Context Navigation

source: LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h @ 1114

Download in other formats: