Context Navigation

phys_output_write.h @ 1255

Last change on this file since 1255 was 1246, checked in by Laurent Fairhead, 15 years ago

En deconnectant les aérosols (ok_ade=ok_aie=n) on a les mêmes

résultats avant et après les modifs.

preindustrial readin fields are used to compute natural aerosol fields

to allow for clean double calls to radiation

full forcing diagnostics (NAT, ANT, ZERO, Cloud forcing, CS,AS) are

activated with lev_histmth 4, If lev_histmth is not 4, the call to the
radiation is minimized, for efficiency, but ade and aie are computed and
applied (however for species wise forcing one would need to do
difference runs) (still quite a bit new forcing info, requires probably

some more explanation)

there is a hardcoded key in sw_aeroAR4.F90 which lets you choose to use the zero aerosol, or natural aerosol perturbation acting on the meteorology, but still would put out the full forcing diagnostics.
aod fields from offline aerosol fields are also output in histmth for

all aerosol tracers read in and available for evaluation

aeropt contains the ss humidity correction from nicolas&yves

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

File size: 43.6 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,LWdn200)
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,LWdn200clr)
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
514	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
515	CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)
516	ENDIF
517
518	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
519	CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)
520	ENDIF
521
522	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
523	CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)
524	ENDIF
525
526	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
527	CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)
528	ENDIF
529
530	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
531	CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)
532	ENDIF
533
534	ENDIF !iflag_con .GE. 3
535
536	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
537	CALL histwrite_phy(nid_files(iff),o_s_pblh%name,itau_w,s_pblh)
538	ENDIF
539
540	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
541	CALL histwrite_phy(nid_files(iff),o_s_pblt%name,itau_w,s_pblt)
542	ENDIF
543
544	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
545	CALL histwrite_phy(nid_files(iff),o_s_lcl%name,itau_w,s_lcl)
546	ENDIF
547
548	IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
549	CALL histwrite_phy(nid_files(iff),o_s_capCL%name,itau_w,s_capCL)
550	ENDIF
551
552	IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
553	CALL histwrite_phy(nid_files(iff),o_s_oliqCL%name,itau_w,s_oliqCL)
554	ENDIF
555
556	IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
557	CALL histwrite_phy(nid_files(iff),o_s_cteiCL%name,itau_w,s_cteiCL)
558	ENDIF
559
560	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
561	CALL histwrite_phy(nid_files(iff),o_s_therm%name,itau_w,s_therm)
562	ENDIF
563
564	IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
565	CALL histwrite_phy(nid_files(iff),o_s_trmb1%name,itau_w,s_trmb1)
566	ENDIF
567
568	IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
569	CALL histwrite_phy(nid_files(iff),o_s_trmb2%name,itau_w,s_trmb2)
570	ENDIF
571
572	IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
573	CALL histwrite_phy(nid_files(iff),o_s_trmb3%name,itau_w,s_trmb3)
574	ENDIF
575
576	! Champs interpolles sur des niveaux de pression
577
578	ll=0
579	DO k=1, nlevSTD
580	! IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
581	! IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
582	bb2=clevSTD(k)
583	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
584	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
585	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
586
587	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
588	ll=ll+1
589	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
590	CALL histwrite_phy(nid_files(iff),o_uSTDlevs(ll)%name,
591	& itau_w,uwriteSTD(:,k,iff))
592	ENDIF
593
594	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
595	CALL histwrite_phy(nid_files(iff),o_vSTDlevs(ll)%name,
596	& itau_w,vwriteSTD(:,k,iff))
597	ENDIF
598
599	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
600	CALL histwrite_phy(nid_files(iff),o_wSTDlevs(ll)%name,
601	& itau_w,wwriteSTD(:,k,iff))
602	ENDIF
603
604	IF (o_phiSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
605	CALL histwrite_phy(nid_files(iff),o_phiSTDlevs(ll)%name,
606	& itau_w,phiwriteSTD(:,k,iff))
607	ENDIF
608
609	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
610	CALL histwrite_phy(nid_files(iff),o_qSTDlevs(ll)%name,
611	& itau_w, qwriteSTD(:,k,iff))
612	ENDIF
613
614	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
615	CALL histwrite_phy(nid_files(iff),o_tSTDlevs(ll)%name,
616	& itau_w, twriteSTD(:,k,iff))
617	ENDIF
618
619	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
620	ENDDO
621
622	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
623	DO i=1, klon
624	IF (pctsrf(i,is_oce).GT.epsfra.OR.
625	$ pctsrf(i,is_sic).GT.epsfra) THEN
626	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
627	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
628	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
629	ELSE
630	zx_tmp_fi2d(i) = 273.15
631	ENDIF
632	ENDDO
633	CALL histwrite_phy(nid_files(iff),
634	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
635	ENDIF
636
637	! Couplage convection-couche limite
638	IF (iflag_con.GE.3) THEN
639	IF (iflag_coupl.EQ.1) THEN
640	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
641	CALL histwrite_phy(nid_files(iff),o_ale_bl%name,itau_w,ale_bl)
642	ENDIF
643	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
644	CALL histwrite_phy(nid_files(iff),o_alp_bl%name,itau_w,alp_bl)
645	ENDIF
646	ENDIF !iflag_coupl.EQ.1
647	ENDIF !(iflag_con.GE.3)
648
649	! Wakes
650	IF (iflag_con.EQ.3) THEN
651	IF (iflag_wake.EQ.1) THEN
652	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
653	CALL histwrite_phy(nid_files(iff),o_ale_wk%name,itau_w,ale_wake)
654	ENDIF
655	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
656	CALL histwrite_phy(nid_files(iff),o_alp_wk%name,itau_w,alp_wake)
657	ENDIF
658
659	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
660	CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)
661	ENDIF
662	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
663	CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)
664	ENDIF
665	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
666	CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)
667	ENDIF
668	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
669	CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe)
670	ENDIF
671	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
672	CALL histwrite_phy(nid_files(iff),o_wake_h%name,itau_w,wake_h)
673	ENDIF
674
675	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
676	CALL histwrite_phy(nid_files(iff),o_wake_s%name,itau_w,wake_s)
677	ENDIF
678
679	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
680	CALL histwrite_phy(nid_files(iff),o_wake_deltat%name,
681	$ itau_w,wake_deltat)
682	ENDIF
683
684	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
685	CALL histwrite_phy(nid_files(iff),o_wake_deltaq%name,
686	$ itau_w,wake_deltaq)
687	ENDIF
688
689	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
690	CALL histwrite_phy(nid_files(iff),
691	s o_wake_omg%name,itau_w,wake_omg)
692	ENDIF
693
694	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
695	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
696	& /pdtphys
697	CALL histwrite_phy(nid_files(iff),
698	& o_dtwak%name,itau_w,zx_tmp_fi3d)
699	ENDIF
700
701	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
702	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
703	CALL histwrite_phy(nid_files(iff),
704	& o_dqwak%name,itau_w,zx_tmp_fi3d)
705	ENDIF
706	ENDIF ! iflag_wake.EQ.1
707
708	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
709	CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip)
710	ENDIF
711
712	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
713	CALL histwrite_phy(nid_files(iff),o_ftd%name,itau_w,ftd)
714	ENDIF
715
716	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
717	CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd)
718	ENDIF
719	ENDIF !(iflag_con.EQ.3)
720
721	IF (type_ocean=='slab ') THEN
722	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
723	$ CALL histwrite_phy(
724	$ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils)
725
726	ENDIF !type_ocean == force/slab
727
728	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
729	CALL histwrite_phy(nid_files(iff),
730	s o_weakinv%name,itau_w,weak_inversion)
731	ENDIF
732
733	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
734	CALL histwrite_phy(nid_files(iff),o_dthmin%name,itau_w,dthmin)
735	ENDIF
736
737	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
738	CALL histwrite_phy(nid_files(iff),o_cldtau%name,itau_w,cldtau)
739	ENDIF
740
741	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
742	CALL histwrite_phy(nid_files(iff),o_cldemi%name,itau_w,cldemi)
743	ENDIF
744
745	! IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
746	! CALL histwrite_phy(nid_files(iff),o_pmflxr%name,itau_w,pmflxr)
747	! ENDIF
748
749	! IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
750	! CALL histwrite_phy(nid_files(iff),o_pmflxs%name,itau_w,pmflxs)
751	! ENDIF
752
753	! IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
754	! CALL histwrite_phy(nid_files(iff),o_prfl%name,itau_w,prfl)
755	! ENDIF
756
757	! IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
758	! CALL histwrite_phy(nid_files(iff),o_psfl%name,itau_w,psfl)
759	! ENDIF
760
761	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
762	DO i=1, klon
763	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
764	ENDDO
765	CALL histwrite_phy(nid_files(iff),o_rh2m%name,itau_w,zx_tmp_fi2d)
766	ENDIF
767
768	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
769	CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m)
770	ENDIF
771
772	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
773	CALL histwrite_phy(nid_files(iff),o_tpot%name,itau_w,tpot)
774	ENDIF
775
776	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
777	CALL histwrite_phy(nid_files(iff),o_tpote%name,itau_w,tpote)
778	ENDIF
779
780	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
781	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
782	CALL histwrite_phy(nid_files(iff),
783	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
784	ENDIF
785
786	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
787	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
788	CALL histwrite_phy(nid_files(iff),
789	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
790	ENDIF
791
792	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
793	CALL histwrite_phy(nid_files(iff),
794	s o_LWdownOR%name,itau_w,sollwdown)
795	ENDIF
796
797	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
798	CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc)
799	ENDIF
800
801	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
802	CALL histwrite_phy(nid_files(iff),
803	s o_solldown%name,itau_w,sollwdown)
804	ENDIF
805
806	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
807	CALL histwrite_phy(nid_files(iff),
808	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
809	ENDIF
810
811	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
812	CALL histwrite_phy(nid_files(iff),
813	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
814	ENDIF
815
816	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
817	CALL histwrite_phy(nid_files(iff),
818	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
819	ENDIF
820
821	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
822	CALL histwrite_phy(nid_files(iff),
823	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
824	ENDIF
825
826	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
827	CALL histwrite_phy(nid_files(iff),o_rugs%name,itau_w,zxrugs)
828	ENDIF
829
830	! OD550 per species
831	IF (new_aod) THEN
832	DO naero = 1, naero_spc
833	IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
834	CALL histwrite_phy(nid_files(iff),
835	$ o_tausumaero(naero)%name,itau_w,
836	$ tausum_aero(:,2,naero) )
837	ENDIF
838	END DO
839	ENDIF
840
841	IF (ok_ade) THEN
842	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
843	CALL histwrite_phy(nid_files(iff),o_topswad%name,itau_w,
844	$ topswad_aero)
845	ENDIF
846	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
847	CALL histwrite_phy(nid_files(iff),o_solswad%name,itau_w,
848	$ solswad_aero)
849	ENDIF
850
851	!====MS forcing diagnostics
852	if (new_aod) then
853	IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
854	CALL histwrite_phy(nid_files(iff),o_swtoaas_nat%name,itau_w,
855	$ topsw_aero(:,1))
856	ENDIF
857
858	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
859	CALL histwrite_phy(nid_files(iff),o_swsrfas_nat%name,itau_w,
860	$ solsw_aero(:,1))
861	ENDIF
862
863	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
864	CALL histwrite_phy(nid_files(iff),o_swtoacs_nat%name,itau_w,
865	$ topsw0_aero(:,1))
866	ENDIF
867
868	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
869	CALL histwrite_phy(nid_files(iff),o_swsrfcs_nat%name,itau_w,
870	$ solsw0_aero(:,1))
871	ENDIF
872
873	!ant
874	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
875	CALL histwrite_phy(nid_files(iff),o_swtoaas_ant%name,itau_w,
876	$ topsw_aero(:,2))
877	ENDIF
878
879	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
880	CALL histwrite_phy(nid_files(iff),o_swsrfas_ant%name,itau_w,
881	$ solsw_aero(:,2))
882	ENDIF
883
884	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
885	CALL histwrite_phy(nid_files(iff),o_swtoacs_ant%name,itau_w,
886	$ topsw0_aero(:,2))
887	ENDIF
888
889	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
890	CALL histwrite_phy(nid_files(iff),o_swsrfcs_ant%name,itau_w,
891	$ solsw0_aero(:,2))
892	ENDIF
893
894	!cf
895	IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
896	CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
897	$ topswcf_aero(:,1))
898	ENDIF
899
900	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
901	CALL histwrite_phy(nid_files(iff),o_swsrfcf_nat%name,itau_w,
902	$ solswcf_aero(:,1))
903	ENDIF
904
905	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
906	CALL histwrite_phy(nid_files(iff),o_swtoacf_ant%name,itau_w,
907	$ topswcf_aero(:,2))
908	ENDIF
909
910	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
911	CALL histwrite_phy(nid_files(iff),o_swsrfcf_ant%name,itau_w,
912	$ solswcf_aero(:,2))
913	ENDIF
914
915	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
916	CALL histwrite_phy(nid_files(iff),o_swtoacf_zero%name,itau_w,
917	$ topswcf_aero(:,3))
918	ENDIF
919
920	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
921	CALL histwrite_phy(nid_files(iff),o_swsrfcf_zero%name,itau_w,
922	$ solswcf_aero(:,3))
923	ENDIF
924
925	endif ! new_aod
926	!====MS forcing diagnostics
927
928	ENDIF
929
930	IF (ok_aie) THEN
931	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
932	CALL histwrite_phy(nid_files(iff),o_topswai%name,itau_w,
933	$ topswai_aero)
934	ENDIF
935	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
936	CALL histwrite_phy(nid_files(iff),o_solswai%name,itau_w,
937	$ solswai_aero)
938	ENDIF
939	ENDIF
940
941	! Champs 3D:
942	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
943	CALL histwrite_phy(nid_files(iff),o_lwcon%name,itau_w,flwc)
944	ENDIF
945
946	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
947	CALL histwrite_phy(nid_files(iff),o_iwcon%name,itau_w,fiwc)
948	ENDIF
949
950	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
951	CALL histwrite_phy(nid_files(iff),o_temp%name,itau_w,t_seri)
952	ENDIF
953
954	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
955	CALL histwrite_phy(nid_files(iff),o_theta%name,itau_w,theta)
956	ENDIF
957
958	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
959	CALL histwrite_phy(nid_files(iff),o_ovap%name,itau_w,qx(:,:,ivap))
960	ENDIF
961
962	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
963	CALL histwrite_phy(nid_files(iff),
964	$ o_ovapinit%name,itau_w,q_seri)
965	ENDIF
966
967	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
968	CALL histwrite_phy(nid_files(iff),o_geop%name,itau_w,zphi)
969	ENDIF
970
971	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
972	CALL histwrite_phy(nid_files(iff),o_vitu%name,itau_w,u_seri)
973	ENDIF
974
975	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
976	CALL histwrite_phy(nid_files(iff),o_vitv%name,itau_w,v_seri)
977	ENDIF
978
979	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
980	CALL histwrite_phy(nid_files(iff),o_vitw%name,itau_w,omega)
981	ENDIF
982
983	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
984	CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay)
985	ENDIF
986
987	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
988	CALL histwrite_phy(nid_files(iff),o_rneb%name,itau_w,cldfra)
989	ENDIF
990
991	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
992	CALL histwrite_phy(nid_files(iff),o_rnebcon%name,itau_w,rnebcon)
993	ENDIF
994
995	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
996	CALL histwrite_phy(nid_files(iff),o_rhum%name,itau_w,zx_rh)
997	ENDIF
998
999	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
1000	CALL histwrite_phy(nid_files(iff), o_ozone%name, itau_w,
1001	$ wo * dobson_u * 1e3 / zmasse / rmo3 * rmd)
1002	ENDIF
1003
1004	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
1005	CALL histwrite_phy(nid_files(iff),o_dtphy%name,itau_w,d_t)
1006	ENDIF
1007
1008	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
1009	CALL histwrite_phy(nid_files(iff),
1010	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
1011	ENDIF
1012
1013	DO nsrf=1, nbsrf
1014	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1015	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1016	CALL histwrite_phy(nid_files(iff),
1017	s o_albe_srf(nsrf)%name,itau_w,
1018	$ zx_tmp_fi2d)
1019	ENDIF
1020
1021	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1022	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
1023	CALL histwrite_phy(nid_files(iff),
1024	s o_rugs_srf(nsrf)%name,itau_w,
1025	$ zx_tmp_fi2d)
1026	ENDIF
1027
1028	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
1029	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1030	CALL histwrite_phy(nid_files(iff),
1031	s o_ages_srf(nsrf)%name,itau_w
1032	$ ,zx_tmp_fi2d)
1033	ENDIF
1034	ENDDO !nsrf=1, nbsrf
1035
1036	IF (o_albs%flag(iff)<=lev_files(iff)) THEN
1037	CALL histwrite_phy(nid_files(iff),o_albs%name,itau_w,albsol1)
1038	ENDIF
1039
1040	IF (o_albslw%flag(iff)<=lev_files(iff)) THEN
1041	CALL histwrite_phy(nid_files(iff),o_albslw%name,itau_w,albsol2)
1042	ENDIF
1043
1044	!FH Sorties pour la couche limite
1045	if (iflag_pbl>1) then
1046	zx_tmp_fi3d=0.
1047	do nsrf=1,nbsrf
1048	do k=1,klev
1049	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
1050	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1051	enddo
1052	enddo
1053	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
1054	CALL histwrite_phy(nid_files(iff),o_tke%name,itau_w,zx_tmp_fi3d)
1055	ENDIF
1056
1057	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
1058	CALL histwrite_phy(nid_files(iff),
1059	s o_tke_max%name,itau_w,zx_tmp_fi3d)
1060	ENDIF
1061	endif
1062
1063	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
1064	CALL histwrite_phy(nid_files(iff),o_kz%name,itau_w,coefh)
1065	ENDIF
1066
1067	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
1068	CALL histwrite_phy(nid_files(iff),o_kz_max%name,itau_w,coefh)
1069	ENDIF
1070
1071	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
1072	CALL histwrite_phy(nid_files(iff),o_clwcon%name,itau_w,clwcon0)
1073	ENDIF
1074
1075	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
1076	CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn)
1077	ENDIF
1078
1079	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
1080	CALL histwrite_phy(nid_files(iff),o_dqdyn%name,itau_w,d_q_dyn)
1081	ENDIF
1082
1083	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
1084	CALL histwrite_phy(nid_files(iff),o_dudyn%name,itau_w,d_u_dyn)
1085	ENDIF
1086
1087	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
1088	CALL histwrite_phy(nid_files(iff),o_dvdyn%name,itau_w,d_v_dyn)
1089	ENDIF
1090
1091	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
1092	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1093	CALL histwrite_phy(nid_files(iff),o_dtcon%name,itau_w,zx_tmp_fi3d)
1094	ENDIF
1095
1096	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
1097	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1098	CALL histwrite_phy(nid_files(iff),o_ducon%name,itau_w,zx_tmp_fi3d)
1099	ENDIF
1100
1101	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
1102	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1103	CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d)
1104	ENDIF
1105
1106	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
1107	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1108	CALL histwrite_phy(nid_files(iff),o_dtlsc%name,itau_w,zx_tmp_fi3d)
1109	ENDIF
1110
1111	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
1112	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
1113	$ d_t_eva(1:klon,1:klev))/pdtphys
1114	CALL histwrite_phy(nid_files(iff),
1115	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
1116	ENDIF
1117
1118	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
1119	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1120	CALL histwrite_phy(nid_files(iff),o_dqlsc%name,itau_w,zx_tmp_fi3d)
1121	ENDIF
1122
1123	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
1124	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1125	CALL histwrite_phy(nid_files(iff),o_dtvdf%name,itau_w,zx_tmp_fi3d)
1126	ENDIF
1127
1128	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
1129	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1130	CALL histwrite_phy(nid_files(iff),o_dqvdf%name,itau_w,zx_tmp_fi3d)
1131	ENDIF
1132
1133	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
1134	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1135	CALL histwrite_phy(nid_files(iff),o_dteva%name,itau_w,zx_tmp_fi3d)
1136	ENDIF
1137
1138	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
1139	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1140	CALL histwrite_phy(nid_files(iff),o_dqeva%name,itau_w,zx_tmp_fi3d)
1141	ENDIF
1142
1143	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
1144	zpt_conv = 0.
1145	where (ptconv) zpt_conv = 1.
1146	CALL histwrite_phy(nid_files(iff),o_ptconv%name,itau_w,zpt_conv)
1147	ENDIF
1148
1149	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
1150	CALL histwrite_phy(nid_files(iff),o_ratqs%name,itau_w,ratqs)
1151	ENDIF
1152
1153	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
1154	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
1155	CALL histwrite_phy(nid_files(iff),o_dtthe%name,itau_w,zx_tmp_fi3d)
1156	ENDIF
1157
1158	IF (iflag_thermals.gt.1) THEN
1159	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
1160	CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm)
1161	ENDIF
1162
1163	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
1164	CALL histwrite_phy(nid_files(iff),o_e_th%name,itau_w,entr_therm)
1165	ENDIF
1166
1167	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
1168	CALL histwrite_phy(nid_files(iff),o_w_th%name,itau_w,zw2)
1169	ENDIF
1170
1171	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
1172	CALL histwrite_phy(nid_files(iff),o_q_th%name,itau_w,zqasc)
1173	ENDIF
1174
1175	IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN
1176	CALL histwrite_phy(nid_files(iff),
1177	s o_lambda_th%name,itau_w,lambda_th)
1178	ENDIF
1179
1180	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
1181	CALL histwrite_phy(nid_files(iff),o_a_th%name,itau_w,fraca)
1182	ENDIF
1183
1184	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
1185	CALL histwrite_phy(nid_files(iff),o_d_th%name,itau_w,detr_therm)
1186	ENDIF
1187
1188	ENDIF !iflag_thermals
1189
1190	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1191	CALL histwrite_phy(nid_files(iff),o_f0_th%name,itau_w,f0)
1192	ENDIF
1193
1194	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
1195	CALL histwrite_phy(nid_files(iff),o_zmax_th%name,itau_w,zmax0)
1196	ENDIF
1197
1198	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
1199	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
1200	CALL histwrite_phy(nid_files(iff),o_dqthe%name,itau_w,zx_tmp_fi3d)
1201	ENDIF
1202
1203	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
1204	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1205	CALL histwrite_phy(nid_files(iff),o_dtajs%name,itau_w,zx_tmp_fi3d)
1206	ENDIF
1207
1208	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
1209	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1210	CALL histwrite_phy(nid_files(iff),o_dqajs%name,itau_w,zx_tmp_fi3d)
1211	ENDIF
1212
1213	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
1214	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
1215	CALL histwrite_phy(nid_files(iff),o_dtswr%name,itau_w,zx_tmp_fi3d)
1216	ENDIF
1217
1218	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
1219	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
1220	CALL histwrite_phy(nid_files(iff),o_dtsw0%name,itau_w,zx_tmp_fi3d)
1221	ENDIF
1222
1223	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
1224	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
1225	CALL histwrite_phy(nid_files(iff),o_dtlwr%name,itau_w,zx_tmp_fi3d)
1226	ENDIF
1227
1228	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
1229	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
1230	CALL histwrite_phy(nid_files(iff),o_dtlw0%name,itau_w,zx_tmp_fi3d)
1231	ENDIF
1232
1233	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
1234	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
1235	CALL histwrite_phy(nid_files(iff),o_dtec%name,itau_w,zx_tmp_fi3d)
1236	ENDIF
1237
1238	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
1239	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1240	CALL histwrite_phy(nid_files(iff),o_duvdf%name,itau_w,zx_tmp_fi3d)
1241	ENDIF
1242
1243	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
1244	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1245	CALL histwrite_phy(nid_files(iff),o_dvvdf%name,itau_w,zx_tmp_fi3d)
1246	ENDIF
1247
1248	IF (ok_orodr) THEN
1249	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
1250	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1251	CALL histwrite_phy(nid_files(iff),o_duoro%name,itau_w,zx_tmp_fi3d)
1252	ENDIF
1253
1254	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
1255	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1256	CALL histwrite_phy(nid_files(iff),o_dvoro%name,itau_w,zx_tmp_fi3d)
1257	ENDIF
1258	ENDIF
1259
1260	IF (ok_orolf) THEN
1261	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
1262	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1263	CALL histwrite_phy(nid_files(iff),o_dulif%name,itau_w,zx_tmp_fi3d)
1264	ENDIF
1265
1266	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
1267	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1268	CALL histwrite_phy(nid_files(iff),o_dvlif%name,itau_w,zx_tmp_fi3d)
1269	ENDIF
1270	ENDIF
1271
1272	! IF (o_trac%flag(iff)<=lev_files(iff)) THEN
1273	if (nqtot.GE.3) THEN
1274	! DO iq=3,nqtot
1275	DO iq=3,4
1276	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
1277	CALL histwrite_phy(nid_files(iff),
1278	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
1279	ENDIF
1280	ENDDO
1281	endif
1282
1283	if (ok_sync) then
1284	c$OMP MASTER
1285	call histsync(nid_files(iff))
1286	c$OMP END MASTER
1287	endif
1288
1289	ENDIF ! clef_files
1290
1291	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 @ 1255

Download in other formats: