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

Rev	Line
[1577]	1	itau_w = itau_phy + itap + start_time * day_step / iphysiq
[907]	2
	3	DO iff=1,nfiles
	4
	5	IF (clef_files(iff)) THEN
	6	ndex2d = 0
	7	ndex3d = 0
	8
	9	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1146]	10	IF (o_phis%flag(iff)<=lev_files(iff)) THEN
[1539]	11	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	12	$ o_phis%name,itau_w,pphis)
[907]	13	ENDIF
	14
[1539]	15	IF (.NOT.clef_stations(iff)) THEN
[1146]	16	IF (o_aire%flag(iff)<=lev_files(iff)) THEN
[1539]	17	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	18	$o_aire%name,itau_w,airephy)
[907]	19	ENDIF
	20
[1146]	21	IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN
[907]	22	DO i=1, klon
	23	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
	24	ENDDO
[1539]	25	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	26	$ o_contfracATM%name,itau_w,zx_tmp_fi2d)
[907]	27	ENDIF
[1539]	28	ENDIF
[907]	29
[1146]	30	IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN
[1539]	31	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	32	$ o_contfracOR%name,itau_w,pctsrf(:,is_ter))
[907]	33	ENDIF
	34
[1146]	35	IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN
[1539]	36	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	37	$ o_aireTER%name,itau_w,paire_ter)
[907]	38	ENDIF
	39
	40	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	41
[1146]	42	IF (o_flat%flag(iff)<=lev_files(iff)) THEN
[1539]	43	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	44	$ o_flat%name,itau_w,zxfluxlat)
[907]	45	ENDIF
	46
[1146]	47	IF (o_slp%flag(iff)<=lev_files(iff)) THEN
[1539]	48	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	49	$ o_slp%name,itau_w,slp)
[907]	50	ENDIF
	51
[1146]	52	IF (o_tsol%flag(iff)<=lev_files(iff)) THEN
[1539]	53	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	54	$ o_tsol%name,itau_w,zxtsol)
[907]	55	ENDIF
	56
[1146]	57	IF (o_t2m%flag(iff)<=lev_files(iff)) THEN
[1539]	58	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	59	$ o_t2m%name,itau_w,zt2m)
[907]	60	ENDIF
	61
[1539]	62	IF (.NOT.clef_stations(iff)) THEN
[1146]	63	IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN
[1539]	64	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	65	$o_t2m_min%name,itau_w,zt2m)
[907]	66	ENDIF
	67
[1146]	68	IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN
[1539]	69	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	70	$o_t2m_max%name,itau_w,zt2m)
[907]	71	ENDIF
[1539]	72	ENDIF
[907]	73
[1146]	74	IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN
[907]	75	DO i=1, klon
	76	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	77	ENDDO
[1539]	78	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	79	s o_wind10m%name,itau_w,zx_tmp_fi2d)
[907]	80	ENDIF
	81
[1539]	82	IF (.NOT.clef_stations(iff)) THEN
[1146]	83	IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN
[907]	84	DO i=1, klon
	85	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	86	ENDDO
[1539]	87	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	88	$o_wind10max%name,itau_w,zx_tmp_fi2d)
[907]	89	ENDIF
[1539]	90	ENDIF
[907]	91
[1146]	92	IF (o_sicf%flag(iff)<=lev_files(iff)) THEN
[907]	93	DO i = 1, klon
	94	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	95	ENDDO
[1539]	96	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	97	$ o_sicf%name,itau_w,zx_tmp_fi2d)
[907]	98	ENDIF
	99
[1146]	100	IF (o_q2m%flag(iff)<=lev_files(iff)) THEN
[1539]	101	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	102	$o_q2m%name,itau_w,zq2m)
[907]	103	ENDIF
	104
[1146]	105	IF (o_u10m%flag(iff)<=lev_files(iff)) THEN
[1539]	106	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	107	$o_u10m%name,itau_w,zu10m)
[907]	108	ENDIF
	109
[1146]	110	IF (o_v10m%flag(iff)<=lev_files(iff)) THEN
[1539]	111	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	112	$o_v10m%name,itau_w,zv10m)
[907]	113	ENDIF
	114
[1146]	115	IF (o_psol%flag(iff)<=lev_files(iff)) THEN
[907]	116	DO i = 1, klon
	117	zx_tmp_fi2d(i) = paprs(i,1)
	118	ENDDO
[1539]	119	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	120	s o_psol%name,itau_w,zx_tmp_fi2d)
[907]	121	ENDIF
	122
[1492]	123	IF (o_mass%flag(iff)<=lev_files(iff)) THEN
[1539]	124	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	125	s o_mass%name,itau_w,zmasse)
[1492]	126	ENDIF
	127
	128
[1146]	129	IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN
[1539]	130	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	131	$o_qsurf%name,itau_w,zxqsurf)
[907]	132	ENDIF
	133
	134	if (.not. ok_veget) then
[1146]	135	IF (o_qsol%flag(iff)<=lev_files(iff)) THEN
[1539]	136	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	137	$o_qsol%name,itau_w,qsol)
[907]	138	ENDIF
	139	endif
	140
[1146]	141	IF (o_precip%flag(iff)<=lev_files(iff)) THEN
[907]	142	DO i = 1, klon
	143	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
	144	ENDDO
[1539]	145	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	146	$o_precip%name,itau_w,zx_tmp_fi2d)
[907]	147	ENDIF
	148
[1146]	149	IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN
[1539]	150	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	151	$o_ndayrain%name,itau_w,nday_rain)
[907]	152	ENDIF
	153
[1146]	154	IF (o_plul%flag(iff)<=lev_files(iff)) THEN
[907]	155	DO i = 1, klon
	156	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	157	ENDDO
[1539]	158	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	159	$o_plul%name,itau_w,zx_tmp_fi2d)
[907]	160	ENDIF
	161
[1146]	162	IF (o_pluc%flag(iff)<=lev_files(iff)) THEN
[907]	163	DO i = 1, klon
	164	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	165	ENDDO
[1539]	166	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	167	$o_pluc%name,itau_w,zx_tmp_fi2d)
[907]	168	ENDIF
	169
[1146]	170	IF (o_snow%flag(iff)<=lev_files(iff)) THEN
[1539]	171	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	172	$o_snow%name,itau_w,snow_fall)
[907]	173	ENDIF
	174
[1334]	175	IF (o_msnow%flag(iff)<=lev_files(iff)) THEN
[1539]	176	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	177	$o_msnow%name,itau_w,snow_o)
[1334]	178	ENDIF
	179
	180	IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN
[1539]	181	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	182	$o_fsnow%name,itau_w,zfra_o)
[1334]	183	ENDIF
	184
[1146]	185	IF (o_evap%flag(iff)<=lev_files(iff)) THEN
[1539]	186	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	187	$o_evap%name,itau_w,evap)
[907]	188	ENDIF
	189
[1146]	190	IF (o_tops%flag(iff)<=lev_files(iff)) THEN
[1539]	191	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	192	$o_tops%name,itau_w,topsw)
[907]	193	ENDIF
	194
[1146]	195	IF (o_tops0%flag(iff)<=lev_files(iff)) THEN
[1539]	196	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	197	$o_tops0%name,itau_w,topsw0)
[907]	198	ENDIF
	199
[1146]	200	IF (o_topl%flag(iff)<=lev_files(iff)) THEN
[1539]	201	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	202	$o_topl%name,itau_w,toplw)
[907]	203	ENDIF
	204
[1146]	205	IF (o_topl0%flag(iff)<=lev_files(iff)) THEN
[1539]	206	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	207	$o_topl0%name,itau_w,toplw0)
[907]	208	ENDIF
	209
[1146]	210	IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN
[907]	211	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
[1539]	212	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	213	$o_SWupTOA%name,itau_w,zx_tmp_fi2d)
[907]	214	ENDIF
	215
[1146]	216	IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN
[907]	217	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
[1539]	218	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	219	$ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d)
[907]	220	ENDIF
	221
[1146]	222	IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN
[907]	223	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
[1539]	224	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	225	s o_SWdnTOA%name,itau_w,zx_tmp_fi2d)
[907]	226	ENDIF
	227
[1146]	228	IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN
[907]	229	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
[1539]	230	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	231	$ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d)
[907]	232	ENDIF
	233
[1334]	234	IF (o_nettop%flag(iff)<=lev_files(iff)) THEN
	235	zx_tmp_fi2d(:) = topsw(:)-toplw(:)
[1539]	236	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1334]	237	$ o_nettop%name,itau_w,zx_tmp_fi2d)
	238	ENDIF
	239
[1146]	240	IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
[1539]	241	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	242	$o_SWup200%name,itau_w,SWup200)
[907]	243	ENDIF
	244
[1146]	245	IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN
[1539]	246	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	247	s o_SWup200clr%name,itau_w,SWup200clr)
[907]	248	ENDIF
	249
[1146]	250	IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN
[1539]	251	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	252	$o_SWdn200%name,itau_w,SWdn200)
[907]	253	ENDIF
	254
[1146]	255	IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN
[1539]	256	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	257	s o_SWdn200clr%name,itau_w,SWdn200clr)
[907]	258	ENDIF
	259
[1146]	260	IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN
[1539]	261	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	262	$o_LWup200%name,itau_w,LWup200)
[907]	263	ENDIF
	264
[1146]	265	IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN
[1539]	266	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	267	s o_LWup200clr%name,itau_w,LWup200clr)
[907]	268	ENDIF
	269
[1146]	270	IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN
[1539]	271	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	272	s o_LWdn200%name,itau_w,LWdn200)
[907]	273	ENDIF
	274
[1146]	275	IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN
[1539]	276	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	277	s o_LWdn200clr%name,itau_w,LWdn200clr)
[907]	278	ENDIF
	279
[1146]	280	IF (o_sols%flag(iff)<=lev_files(iff)) THEN
[1539]	281	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	282	$o_sols%name,itau_w,solsw)
[907]	283	ENDIF
	284
[1146]	285	IF (o_sols0%flag(iff)<=lev_files(iff)) THEN
[1539]	286	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	287	$o_sols0%name,itau_w,solsw0)
[907]	288	ENDIF
	289
[1146]	290	IF (o_soll%flag(iff)<=lev_files(iff)) THEN
[1539]	291	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	292	$o_soll%name,itau_w,sollw)
[907]	293	ENDIF
	294
[1146]	295	IF (o_radsol%flag(iff)<=lev_files(iff)) THEN
[1539]	296	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	297	$o_radsol%name,itau_w,radsol)
[907]	298	ENDIF
	299
[1146]	300	IF (o_soll0%flag(iff)<=lev_files(iff)) THEN
[1539]	301	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	302	$o_soll0%name,itau_w,sollw0)
[907]	303	ENDIF
	304
[1146]	305	IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN
[907]	306	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
[1539]	307	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	308	s o_SWupSFC%name,itau_w,zx_tmp_fi2d)
[907]	309	ENDIF
	310
[1146]	311	IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	312	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
[1539]	313	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	314	$ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	315	ENDIF
	316
[1146]	317	IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN
[907]	318	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
[1539]	319	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	320	$ o_SWdnSFC%name,itau_w,zx_tmp_fi2d)
[907]	321	ENDIF
	322
[1146]	323	IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	324	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
[1539]	325	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	326	$ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	327	ENDIF
	328
[1146]	329	IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN
[907]	330	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
[1539]	331	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	332	$ o_LWupSFC%name,itau_w,zx_tmp_fi2d)
[907]	333	ENDIF
	334
[1146]	335	IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN
[1539]	336	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	337	$ o_LWdnSFC%name,itau_w,sollwdown)
[907]	338	ENDIF
	339
[958]	340	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
[1146]	341	IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	342	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
[1539]	343	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	344	$ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	345	ENDIF
	346
[1146]	347	IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
[1539]	348	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	349	$ o_LWdnSFCclr%name,itau_w,sollwdownclr)
[907]	350	ENDIF
	351
[1146]	352	IF (o_bils%flag(iff)<=lev_files(iff)) THEN
[1539]	353	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	354	$o_bils%name,itau_w,bils)
[907]	355	ENDIF
	356
[1146]	357	IF (o_sens%flag(iff)<=lev_files(iff)) THEN
[907]	358	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
[1539]	359	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	360	$o_sens%name,itau_w,zx_tmp_fi2d)
[907]	361	ENDIF
	362
[1146]	363	IF (o_fder%flag(iff)<=lev_files(iff)) THEN
[1539]	364	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	365	$o_fder%name,itau_w,fder)
[907]	366	ENDIF
	367
[1146]	368	IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN
[1539]	369	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	370	$o_ffonte%name,itau_w,zxffonte)
[907]	371	ENDIF
	372
[1146]	373	IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN
[1539]	374	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	375	$ o_fqcalving%name,itau_w,zxfqcalving)
[907]	376	ENDIF
	377
[1146]	378	IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN
[1539]	379	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	380	$ o_fqfonte%name,itau_w,zxfqfonte)
[907]	381	ENDIF
	382
[1334]	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
[1539]	388	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1334]	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
[1539]	397	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1334]	398	$ o_tauy%name,itau_w,zx_tmp_fi2d)
	399	ENDIF
	400
[907]	401	DO nsrf = 1, nbsrf
[1146]	402	! IF(nsrf.GE.2) THEN
	403	IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	404	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
	405	CALL histwrite_phy(nid_files(iff),
[1539]	406	$ clef_stations(iff),
[1146]	407	$ o_pourc_srf(nsrf)%name,itau_w,
[907]	408	$ zx_tmp_fi2d)
	409	ENDIF
	410
[1146]	411	IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	412	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
[1539]	413	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	414	$ o_fract_srf(nsrf)%name,itau_w,
	415	$ zx_tmp_fi2d)
[907]	416	ENDIF
[1146]	417	! ENDIF !nsrf.GT.2
[907]	418
[1146]	419	IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	420	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
[1539]	421	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	422	$ o_taux_srf(nsrf)%name,itau_w,
	423	$ zx_tmp_fi2d)
[907]	424	ENDIF
	425
[1146]	426	IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	427	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
[1539]	428	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	429	$ o_tauy_srf(nsrf)%name,itau_w,
	430	$ zx_tmp_fi2d)
[907]	431	ENDIF
	432
[1146]	433	IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	434	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
[1539]	435	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	436	$ o_tsol_srf(nsrf)%name,itau_w,
[907]	437	$ zx_tmp_fi2d)
	438	ENDIF
	439
[1146]	440	IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	441	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
[1539]	442	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	443	$o_u10m_srf(nsrf)%name,
[907]	444	$ itau_w,zx_tmp_fi2d)
	445	ENDIF
	446
[1146]	447	IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	448	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
[1539]	449	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	450	$o_v10m_srf(nsrf)%name,
[907]	451	$ itau_w,zx_tmp_fi2d)
	452	ENDIF
	453
[1146]	454	IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	455	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
[1539]	456	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	457	$o_t2m_srf(nsrf)%name,
[907]	458	$ itau_w,zx_tmp_fi2d)
	459	ENDIF
	460
[1334]	461	IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
	462	zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
[1539]	463	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	464	$o_evap_srf(nsrf)%name,
[1334]	465	$ itau_w,zx_tmp_fi2d)
	466	ENDIF
	467
[1146]	468	IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	469	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
[1539]	470	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	471	$ o_sens_srf(nsrf)%name,itau_w,
[907]	472	$ zx_tmp_fi2d)
	473	ENDIF
	474
[1146]	475	IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	476	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
[1539]	477	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	478	$ o_lat_srf(nsrf)%name,itau_w,
[907]	479	$ zx_tmp_fi2d)
	480	ENDIF
	481
[1146]	482	IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	483	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
[1539]	484	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	485	$ o_flw_srf(nsrf)%name,itau_w,
[907]	486	$ zx_tmp_fi2d)
	487	ENDIF
	488
[1146]	489	IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	490	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
[1539]	491	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	492	$ o_fsw_srf(nsrf)%name,itau_w,
[907]	493	$ zx_tmp_fi2d)
	494	ENDIF
	495
[1146]	496	IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	497	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
[1539]	498	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	499	$ o_wbils_srf(nsrf)%name,itau_w,
[907]	500	$ zx_tmp_fi2d)
	501	ENDIF
	502
[1146]	503	IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	504	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
[1539]	505	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	506	$ o_wbilo_srf(nsrf)%name,itau_w,
[907]	507	$ zx_tmp_fi2d)
	508	ENDIF
	509
[1638]	510	if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then
[1146]	511	IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[1539]	512	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	513	$ o_tke_srf(nsrf)%name,itau_w,
	514	$ pbl_tke(:,1:klev,nsrf))
[907]	515	ENDIF
	516
[1539]	517	IF (.NOT.clef_stations(iff)) THEN
[1146]	518	IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[1539]	519	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	520	$ o_tke_max_srf(nsrf)%name,itau_w,
[907]	521	$ pbl_tke(:,1:klev,nsrf))
	522	ENDIF
[1539]	523	ENDIF
[907]	524	endif
	525	ENDDO
	526
[1146]	527	IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN
[1539]	528	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	529	$o_cdrm%name,itau_w,cdragm)
[907]	530	ENDIF
	531
[1146]	532	IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN
[1539]	533	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	534	$o_cdrh%name,itau_w,cdragh)
[907]	535	ENDIF
	536
[1146]	537	IF (o_cldl%flag(iff)<=lev_files(iff)) THEN
[1539]	538	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	539	$o_cldl%name,itau_w,cldl)
[907]	540	ENDIF
	541
[1146]	542	IF (o_cldm%flag(iff)<=lev_files(iff)) THEN
[1539]	543	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	544	$o_cldm%name,itau_w,cldm)
[907]	545	ENDIF
	546
[1146]	547	IF (o_cldh%flag(iff)<=lev_files(iff)) THEN
[1539]	548	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	549	$o_cldh%name,itau_w,cldh)
[907]	550	ENDIF
	551
[1146]	552	IF (o_cldt%flag(iff)<=lev_files(iff)) THEN
[1539]	553	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	554	$o_cldt%name,
[1285]	555	& itau_w,cldt)
[907]	556	ENDIF
	557
[1146]	558	IF (o_cldq%flag(iff)<=lev_files(iff)) THEN
[1539]	559	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	560	$o_cldq%name,itau_w,cldq)
[907]	561	ENDIF
	562
[1146]	563	IF (o_lwp%flag(iff)<=lev_files(iff)) THEN
[907]	564	zx_tmp_fi2d(1:klon) = flwp(1:klon)
[1539]	565	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	566	s o_lwp%name,itau_w,zx_tmp_fi2d)
[907]	567	ENDIF
	568
[1146]	569	IF (o_iwp%flag(iff)<=lev_files(iff)) THEN
[907]	570	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
[1539]	571	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	572	s o_iwp%name,itau_w,zx_tmp_fi2d)
[907]	573	ENDIF
	574
[1146]	575	IF (o_ue%flag(iff)<=lev_files(iff)) THEN
[1539]	576	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	577	$o_ue%name,itau_w,ue)
[907]	578	ENDIF
	579
[1146]	580	IF (o_ve%flag(iff)<=lev_files(iff)) THEN
[1539]	581	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	582	$o_ve%name,itau_w,ve)
[907]	583	ENDIF
	584
[1146]	585	IF (o_uq%flag(iff)<=lev_files(iff)) THEN
[1539]	586	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	587	$o_uq%name,itau_w,uq)
[907]	588	ENDIF
	589
[1146]	590	IF (o_vq%flag(iff)<=lev_files(iff)) THEN
[1539]	591	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	592	$o_vq%name,itau_w,vq)
[907]	593	ENDIF
	594
	595	IF(iflag_con.GE.3) THEN ! sb
[1146]	596	IF (o_cape%flag(iff)<=lev_files(iff)) THEN
[1539]	597	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	598	$o_cape%name,itau_w,cape)
[907]	599	ENDIF
	600
[1146]	601	IF (o_pbase%flag(iff)<=lev_files(iff)) THEN
[1539]	602	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	603	$o_pbase%name,itau_w,ema_pcb)
[907]	604	ENDIF
	605
[1146]	606	IF (o_ptop%flag(iff)<=lev_files(iff)) THEN
[1539]	607	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	608	$o_ptop%name,itau_w,ema_pct)
[907]	609	ENDIF
	610
[1146]	611	IF (o_fbase%flag(iff)<=lev_files(iff)) THEN
[1539]	612	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	613	$o_fbase%name,itau_w,ema_cbmf)
[907]	614	ENDIF
	615
[1562]	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
[1518]	621
[1562]	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
[1518]	626
[1562]	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
[1518]	632
[1539]	633	IF (.NOT.clef_stations(iff)) THEN
[1146]	634	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
[1539]	635	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	636	$o_cape_max%name,itau_w,cape)
[1100]	637	ENDIF
[1539]	638	ENDIF
[1100]	639
[1146]	640	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
[1539]	641	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	642	$o_upwd%name,itau_w,upwd)
[1146]	643	ENDIF
	644
	645	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
[1539]	646	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	647	$o_Ma%name,itau_w,Ma)
[1146]	648	ENDIF
	649
	650	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
[1539]	651	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	652	$o_dnwd%name,itau_w,dnwd)
[1146]	653	ENDIF
	654
	655	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
[1539]	656	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	657	$o_dnwd0%name,itau_w,dnwd0)
[1146]	658	ENDIF
	659
[1334]	660	IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN
	661	zx_tmp_fi2d=float(itau_con)/float(itap)
[1539]	662	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	663	$o_ftime_con%name,
[1334]	664	s itau_w,zx_tmp_fi2d)
	665	ENDIF
	666
	667	IF (o_mc%flag(iff)<=lev_files(iff)) THEN
[1638]	668	if(iflag_thermals>=1)then
[1507]	669	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
[1334]	670	else
	671	zx_tmp_fi3d=dnwd+dnwd0+upwd
	672	endif
[1539]	673	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	674	$o_mc%name,itau_w,zx_tmp_fi3d)
[1334]	675	ENDIF
	676
[907]	677	ENDIF !iflag_con .GE. 3
	678
[1638]	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
[1633]	683
[1146]	684	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
[1539]	685	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	686	$o_s_pblh%name,itau_w,s_pblh)
[907]	687	ENDIF
	688
[1146]	689	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
[1539]	690	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	691	$o_s_pblt%name,itau_w,s_pblt)
[907]	692	ENDIF
	693
[1146]	694	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
[1539]	695	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	696	$o_s_lcl%name,itau_w,s_lcl)
[907]	697	ENDIF
	698
[1398]	699	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
[1539]	700	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	701	$o_s_therm%name,itau_w,s_therm)
[907]	702	ENDIF
	703
[1398]	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
[1539]	706	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	707	! $o_s_capCL%name,itau_w,s_capCL)
[1398]	708	! ENDIF
[907]	709
[1398]	710	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
[1539]	711	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	712	! $o_s_oliqCL%name,itau_w,s_oliqCL)
[1398]	713	! ENDIF
[907]	714
[1398]	715	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
[1539]	716	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	717	! $o_s_cteiCL%name,itau_w,s_cteiCL)
[1398]	718	! ENDIF
[907]	719
[1398]	720	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
[1539]	721	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	722	! $o_s_trmb1%name,itau_w,s_trmb1)
[1398]	723	! ENDIF
[907]	724
[1398]	725	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
[1539]	726	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	727	! $o_s_trmb2%name,itau_w,s_trmb2)
[1398]	728	! ENDIF
[907]	729
[1398]	730	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
[1539]	731	! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	732	! $o_s_trmb3%name,itau_w,s_trmb3)
[1398]	733	! ENDIF
[907]	734
	735	! Champs interpolles sur des niveaux de pression
	736
[1146]	737	ll=0
[907]	738	DO k=1, nlevSTD
[1279]	739	bb2=clevSTD(k)
[907]	740	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
[1279]	741	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
[1539]	742	$ bb2.EQ."100".OR.
[1279]	743	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
[907]	744
	745	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1146]	746	ll=ll+1
	747	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	748	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	749	$o_uSTDlevs(ll)%name,
[1146]	750	& itau_w,uwriteSTD(:,k,iff))
[907]	751	ENDIF
	752
[1146]	753	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	754	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	755	$o_vSTDlevs(ll)%name,
[1146]	756	& itau_w,vwriteSTD(:,k,iff))
[907]	757	ENDIF
	758
[1146]	759	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	760	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	761	$o_wSTDlevs(ll)%name,
[1146]	762	& itau_w,wwriteSTD(:,k,iff))
[907]	763	ENDIF
	764
[1398]	765	IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	766	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	767	$o_zSTDlevs(ll)%name,
[1055]	768	& itau_w,phiwriteSTD(:,k,iff))
[907]	769	ENDIF
	770
[1146]	771	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	772	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	773	$o_qSTDlevs(ll)%name,
[1146]	774	& itau_w, qwriteSTD(:,k,iff))
[907]	775	ENDIF
	776
[1146]	777	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
[1539]	778	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	779	$o_tSTDlevs(ll)%name,
[1146]	780	& itau_w, twriteSTD(:,k,iff))
[907]	781	ENDIF
	782
	783	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	784	ENDDO
	785
[1146]	786	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
[907]	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
[1539]	797	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	798	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
[907]	799	ENDIF
	800
[1146]	801	! Couplage convection-couche limite
	802	IF (iflag_con.GE.3) THEN
[1638]	803	IF (iflag_coupl>=1) THEN
[1146]	804	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
[1539]	805	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	806	$o_ale_bl%name,itau_w,ale_bl)
[951]	807	ENDIF
[1146]	808	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
[1539]	809	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	810	$o_alp_bl%name,itau_w,alp_bl)
[951]	811	ENDIF
[1638]	812	ENDIF !iflag_coupl>=1
[1146]	813	ENDIF !(iflag_con.GE.3)
[951]	814
[1146]	815	! Wakes
	816	IF (iflag_con.EQ.3) THEN
[1507]	817	IF (iflag_wake>=1) THEN
[1146]	818	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
[1539]	819	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	820	$o_ale_wk%name,itau_w,ale_wake)
[951]	821	ENDIF
[1146]	822	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
[1539]	823	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	824	$o_alp_wk%name,itau_w,alp_wake)
[951]	825	ENDIF
[973]	826
[1638]	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
[1146]	839	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
[1539]	840	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	841	$o_WAPE%name,itau_w,wake_pe)
[1146]	842	ENDIF
	843	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
[1539]	844	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	845	$o_wake_h%name,itau_w,wake_h)
[1146]	846	ENDIF
[1063]	847
[1146]	848	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
[1539]	849	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	850	$o_wake_s%name,itau_w,wake_s)
[1063]	851	ENDIF
	852
[1146]	853	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
[1539]	854	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	855	$o_wake_deltat%name,
[1063]	856	$ itau_w,wake_deltat)
	857	ENDIF
	858
[1146]	859	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
[1539]	860	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	861	$o_wake_deltaq%name,
[1063]	862	$ itau_w,wake_deltaq)
	863	ENDIF
	864
[1146]	865	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
[1539]	866	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	867	s o_wake_omg%name,itau_w,wake_omg)
[1063]	868	ENDIF
[1146]	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)
[1100]	872	& /pdtphys
[1146]	873	CALL histwrite_phy(nid_files(iff),
[1539]	874	$clef_stations(iff),
[1146]	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
[1539]	880	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	881	& o_dqwak%name,itau_w,zx_tmp_fi3d)
[1100]	882	ENDIF
[1507]	883	ENDIF ! iflag_wake>=1
[1063]	884
[1146]	885	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
[1539]	886	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	887	$o_Vprecip%name,itau_w,Vprecip)
[1063]	888	ENDIF
	889
[1146]	890	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
[1539]	891	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	892	$o_ftd%name,itau_w,ftd)
[1063]	893	ENDIF
	894
[1146]	895	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
[1539]	896	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	897	$o_fqd%name,itau_w,fqd)
[1063]	898	ENDIF
[1100]	899	ENDIF !(iflag_con.EQ.3)
[996]	900
[1638]	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
[1146]	1000	IF (type_ocean=='slab ') THEN
	1001	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
[996]	1002	$ CALL histwrite_phy(
[1539]	1003	$ nid_files(iff),clef_stations(iff),
	1004	$o_slab_bils%name,itau_w,slab_wfbils)
[907]	1005
[996]	1006	ENDIF !type_ocean == force/slab
[907]	1007
[1146]	1008	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
[1539]	1009	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1010	s o_weakinv%name,itau_w,weak_inversion)
[907]	1011	ENDIF
	1012
[1146]	1013	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
[1539]	1014	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1015	$o_dthmin%name,itau_w,dthmin)
[907]	1016	ENDIF
	1017
[1146]	1018	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
[1539]	1019	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1020	$o_cldtau%name,itau_w,cldtau)
[907]	1021	ENDIF
	1022
[1146]	1023	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
[1539]	1024	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1025	$o_cldemi%name,itau_w,cldemi)
[907]	1026	ENDIF
	1027
[1279]	1028	IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
[1539]	1029	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	1030	s o_pr_con_l%name,itau_w,pmflxr(:,1:klev))
	1031	ENDIF
[929]	1032
[1279]	1033	IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
[1539]	1034	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	1035	s o_pr_con_i%name,itau_w,pmflxs(:,1:klev))
	1036	ENDIF
[929]	1037
[1279]	1038	IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
[1539]	1039	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	1040	s o_pr_lsc_l%name,itau_w,prfl(:,1:klev))
	1041	ENDIF
[929]	1042
[1279]	1043	IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
[1539]	1044	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1279]	1045	s o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
	1046	ENDIF
[929]	1047
[1284]	1048	IF (o_re%flag(iff)<=lev_files(iff)) THEN
[1539]	1049	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1050	$o_re%name,itau_w,re)
[1284]	1051	ENDIF
[1279]	1052
[1284]	1053	IF (o_fl%flag(iff)<=lev_files(iff)) THEN
[1539]	1054	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1055	$o_fl%name,itau_w,fl)
[1284]	1056	ENDIF
	1057
	1058
	1059
[1146]	1060	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
[907]	1061	DO i=1, klon
	1062	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	1063	ENDDO
[1539]	1064	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1065	$o_rh2m%name,itau_w,zx_tmp_fi2d)
[907]	1066	ENDIF
	1067
[1539]	1068	IF (.NOT.clef_stations(iff)) THEN
[1334]	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
[1539]	1073	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1074	$o_rh2m_min%name,
[1334]	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
[1539]	1082	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1083	$o_rh2m_max%name,
[1334]	1084	s itau_w,zx_tmp_fi2d)
	1085	ENDIF
[1539]	1086	ENDIF
[1334]	1087
	1088
[1146]	1089	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
[1539]	1090	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1091	$o_qsat2m%name,itau_w,qsat2m)
[907]	1092	ENDIF
	1093
[1146]	1094	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
[1539]	1095	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1096	$o_tpot%name,itau_w,tpot)
[907]	1097	ENDIF
	1098
[1146]	1099	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
[1539]	1100	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1101	$o_tpote%name,itau_w,tpote)
[907]	1102	ENDIF
	1103
[1146]	1104	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
[907]	1105	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
[1539]	1106	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1107	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
[907]	1108	ENDIF
	1109
[1146]	1110	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
[907]	1111	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
[1539]	1112	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1113	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
[907]	1114	ENDIF
	1115
[1146]	1116	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
[1539]	1117	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1118	s o_LWdownOR%name,itau_w,sollwdown)
[907]	1119	ENDIF
	1120
[1146]	1121	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
[1539]	1122	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1123	$o_snowl%name,itau_w,snow_lsc)
[907]	1124	ENDIF
	1125
[1146]	1126	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
[1539]	1127	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1128	s o_solldown%name,itau_w,sollwdown)
[907]	1129	ENDIF
	1130
[1146]	1131	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
[1539]	1132	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1133	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
[907]	1134	ENDIF
	1135
[1146]	1136	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
[1539]	1137	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1138	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
[907]	1139	ENDIF
	1140
[1146]	1141	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
[1539]	1142	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1143	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
[907]	1144	ENDIF
	1145
[1146]	1146	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
[1539]	1147	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1148	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
[907]	1149	ENDIF
	1150
[1146]	1151	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
[1539]	1152	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1153	$o_rugs%name,itau_w,zxrugs)
[907]	1154	ENDIF
	1155
[1279]	1156	! OD550 per species
	1157	IF (new_aod .and. (.not. aerosol_couple)) THEN
[1345]	1158	IF (ok_ade.OR.ok_aie) THEN
[1337]	1159
	1160	IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN
[1539]	1161	CALL histwrite_phy(nid_files(iff),
	1162	$clef_stations(iff),
	1163	$o_od550aer%name,itau_w,
[1337]	1164	$ od550aer)
	1165	ENDIF
	1166	IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN
[1539]	1167	CALL histwrite_phy(nid_files(iff),
	1168	$clef_stations(iff),
	1169	$o_od865aer%name,itau_w,
[1337]	1170	$ od865aer)
	1171	ENDIF
	1172	IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN
[1539]	1173	CALL histwrite_phy(nid_files(iff),
	1174	$clef_stations(iff),
	1175	$o_absvisaer%name,itau_w,
[1337]	1176	$ absvisaer)
	1177	ENDIF
	1178	IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN
[1539]	1179	CALL histwrite_phy(nid_files(iff),
	1180	$clef_stations(iff),
	1181	$o_od550lt1aer%name,itau_w,
[1337]	1182	$ od550lt1aer)
	1183	ENDIF
	1184
[1345]	1185	IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN
[1539]	1186	CALL histwrite_phy(nid_files(iff),
	1187	$clef_stations(iff),
	1188	$o_sconcso4%name,itau_w,
[1345]	1189	$ sconcso4)
[1279]	1190	ENDIF
[1345]	1191	IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN
[1539]	1192	CALL histwrite_phy(nid_files(iff),
	1193	$clef_stations(iff),
	1194	$o_sconcoa%name,itau_w,
[1345]	1195	$ sconcoa)
	1196	ENDIF
	1197	IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN
[1539]	1198	CALL histwrite_phy(nid_files(iff),
	1199	$clef_stations(iff),
	1200	$o_sconcbc%name,itau_w,
[1345]	1201	$ sconcbc)
	1202	ENDIF
	1203	IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN
[1539]	1204	CALL histwrite_phy(nid_files(iff),
	1205	$clef_stations(iff),
	1206	$o_sconcss%name,itau_w,
[1345]	1207	$ sconcss)
	1208	ENDIF
	1209	IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN
[1539]	1210	CALL histwrite_phy(nid_files(iff),
	1211	$clef_stations(iff),
	1212	$o_sconcdust%name,itau_w,
[1345]	1213	$ sconcdust)
	1214	ENDIF
	1215
	1216	IF (o_concso4%flag(iff)<=lev_files(iff)) THEN
[1539]	1217	CALL histwrite_phy(nid_files(iff),
	1218	$clef_stations(iff),
	1219	$o_concso4%name,itau_w,
[1345]	1220	$ concso4)
	1221	ENDIF
	1222	IF (o_concoa%flag(iff)<=lev_files(iff)) THEN
[1539]	1223	CALL histwrite_phy(nid_files(iff),
	1224	$clef_stations(iff),
	1225	$o_concoa%name,itau_w,
[1345]	1226	$ concoa)
	1227	ENDIF
	1228	IF (o_concbc%flag(iff)<=lev_files(iff)) THEN
[1539]	1229	CALL histwrite_phy(nid_files(iff),
	1230	$clef_stations(iff),
	1231	$o_concbc%name,itau_w,
[1345]	1232	$ concbc)
	1233	ENDIF
	1234	IF (o_concss%flag(iff)<=lev_files(iff)) THEN
[1539]	1235	CALL histwrite_phy(nid_files(iff),
	1236	$clef_stations(iff),
	1237	$o_concss%name,itau_w,
[1345]	1238	$ concss)
	1239	ENDIF
	1240	IF (o_concdust%flag(iff)<=lev_files(iff)) THEN
[1539]	1241	CALL histwrite_phy(nid_files(iff),
	1242	$clef_stations(iff),
	1243	$o_concdust%name,itau_w,
[1345]	1244	$ concdust)
	1245	ENDIF
	1246
	1247	IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN
[1539]	1248	CALL histwrite_phy(nid_files(iff),
	1249	$clef_stations(iff),
	1250	$o_loadso4%name,itau_w,
[1345]	1251	$ loadso4)
	1252	ENDIF
	1253	IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN
[1539]	1254	CALL histwrite_phy(nid_files(iff),
	1255	$clef_stations(iff),
	1256	$o_loadoa%name,itau_w,
[1345]	1257	$ loadoa)
	1258	ENDIF
	1259	IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN
[1539]	1260	CALL histwrite_phy(nid_files(iff),
	1261	$clef_stations(iff),
	1262	$o_loadbc%name,itau_w,
[1345]	1263	$ loadbc)
	1264	ENDIF
	1265	IF (o_loadss%flag(iff)<=lev_files(iff)) THEN
[1539]	1266	CALL histwrite_phy(nid_files(iff),
	1267	$clef_stations(iff),
	1268	$o_loadss%name,itau_w,
[1345]	1269	$ loadss)
	1270	ENDIF
	1271	IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN
[1539]	1272	CALL histwrite_phy(nid_files(iff),
	1273	$clef_stations(iff),
	1274	$o_loaddust%name,itau_w,
[1345]	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),
[1539]	1281	$clef_stations(iff),
[1345]	1282	$ o_tausumaero(naero)%name,itau_w,
	1283	$ tausum_aero(:,2,naero) )
	1284	ENDIF
	1285	END DO
	1286	endif
[1279]	1287	ENDIF
[1345]	1288
[1146]	1289	IF (ok_ade) THEN
[1279]	1290	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
[1539]	1291	CALL histwrite_phy(nid_files(iff),
	1292	$clef_stations(iff),
	1293	$o_topswad%name,itau_w,
[1279]	1294	$ topswad_aero)
	1295	ENDIF
	1296	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
[1539]	1297	CALL histwrite_phy(nid_files(iff),
	1298	$clef_stations(iff),
	1299	$o_solswad%name,itau_w,
[1279]	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
[1539]	1306	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1307	$o_swtoaas_nat%name,itau_w,
[1279]	1308	$ topsw_aero(:,1))
[1146]	1309	ENDIF
[1279]	1310
	1311	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
[1539]	1312	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1313	$o_swsrfas_nat%name,itau_w,
[1279]	1314	$ solsw_aero(:,1))
[1146]	1315	ENDIF
	1316
[1279]	1317	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
[1539]	1318	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1319	$o_swtoacs_nat%name,itau_w,
[1279]	1320	$ topsw0_aero(:,1))
[1146]	1321	ENDIF
[1279]	1322
	1323	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
[1539]	1324	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1325	$o_swsrfcs_nat%name,itau_w,
[1279]	1326	$ solsw0_aero(:,1))
[1146]	1327	ENDIF
[1279]	1328
	1329	!ant
	1330	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1331	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1332	$o_swtoaas_ant%name,itau_w,
[1279]	1333	$ topsw_aero(:,2))
	1334	ENDIF
	1335
	1336	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1337	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1338	$o_swsrfas_ant%name,itau_w,
[1279]	1339	$ solsw_aero(:,2))
	1340	ENDIF
	1341
	1342	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1343	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1344	$o_swtoacs_ant%name,itau_w,
[1279]	1345	$ topsw0_aero(:,2))
	1346	ENDIF
	1347
	1348	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1349	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1350	$o_swsrfcs_ant%name,itau_w,
[1279]	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
[1539]	1358	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1359	$o_swtoacf_nat%name,itau_w,
[1279]	1360	$ topswcf_aero(:,1))
	1361	ENDIF
	1362
	1363	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
[1539]	1364	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1365	$o_swsrfcf_nat%name,itau_w,
[1279]	1366	$ solswcf_aero(:,1))
	1367	ENDIF
	1368
	1369	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1370	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1371	$o_swtoacf_ant%name,itau_w,
[1279]	1372	$ topswcf_aero(:,2))
	1373	ENDIF
	1374
	1375	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
[1539]	1376	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1377	$o_swsrfcf_ant%name,itau_w,
[1279]	1378	$ solswcf_aero(:,2))
	1379	ENDIF
	1380
	1381	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
[1539]	1382	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1383	$o_swtoacf_zero%name,itau_w,
[1279]	1384	$ topswcf_aero(:,3))
	1385	ENDIF
	1386
	1387	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
[1539]	1388	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1389	$o_swsrfcf_zero%name,itau_w,
[1279]	1390	$ solswcf_aero(:,3))
	1391	ENDIF
	1392	endif
	1393
	1394	endif ! new_aod
	1395	!====MS forcing diagnostics
	1396
[1146]	1397	ENDIF
	1398
[1279]	1399	IF (ok_aie) THEN
	1400	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
[1539]	1401	CALL histwrite_phy(nid_files(iff),
	1402	$clef_stations(iff),
	1403	$o_topswai%name,itau_w,
[1279]	1404	$ topswai_aero)
	1405	ENDIF
	1406	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
[1539]	1407	CALL histwrite_phy(nid_files(iff),
	1408	$clef_stations(iff),
	1409	$o_solswai%name,itau_w,
[1279]	1410	$ solswai_aero)
	1411	ENDIF
[1337]	1412	IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN
[1539]	1413	CALL histwrite_phy(nid_files(iff),
	1414	$clef_stations(iff),
	1415	$o_scdnc%name,itau_w,
[1337]	1416	$ scdnc)
	1417	ENDIF
	1418	IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN
[1539]	1419	CALL histwrite_phy(nid_files(iff),
	1420	$clef_stations(iff),
	1421	$o_cldncl%name,itau_w,
[1337]	1422	$ cldncl)
	1423	ENDIF
[1539]	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,
[1337]	1428	$ reffclws)
[1539]	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,
[1337]	1434	$ reffclwc)
[1539]	1435	ENDIF
[1337]	1436	IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN
[1539]	1437	CALL histwrite_phy(nid_files(iff),
	1438	$clef_stations(iff),
	1439	$o_cldnvi%name,itau_w,
[1337]	1440	$ cldnvi)
	1441	ENDIF
	1442	IF (o_lcc%flag(iff)<=lev_files(iff)) THEN
[1539]	1443	CALL histwrite_phy(nid_files(iff),
	1444	$clef_stations(iff),
	1445	$o_lcc%name,itau_w,
[1337]	1446	$ lcc)
	1447	ENDIF
	1448	IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN
[1539]	1449	CALL histwrite_phy(nid_files(iff),
	1450	$clef_stations(iff),
	1451	$o_lcc3d%name,itau_w,
[1337]	1452	$ lcc3d)
	1453	ENDIF
	1454	IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN
[1539]	1455	CALL histwrite_phy(nid_files(iff),
	1456	$clef_stations(iff),
	1457	$o_lcc3dcon%name,itau_w,
[1337]	1458	$ lcc3dcon)
	1459	ENDIF
	1460	IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN
[1539]	1461	CALL histwrite_phy(nid_files(iff),
	1462	$clef_stations(iff),
	1463	$o_lcc3dstra%name,itau_w,
[1337]	1464	$ lcc3dstra)
	1465	ENDIF
	1466	IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN
[1539]	1467	CALL histwrite_phy(nid_files(iff),
	1468	$clef_stations(iff),
	1469	$o_reffclwtop%name,itau_w,
[1337]	1470	$ reffclwtop)
	1471	ENDIF
[1279]	1472	ENDIF
	1473
[907]	1474	! Champs 3D:
[1364]	1475	IF (ok_ade .OR. ok_aie) then
	1476	IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN
[1539]	1477	CALL histwrite_phy(nid_files(iff),
	1478	$clef_stations(iff),
	1479	$o_ec550aer%name,itau_w,
[1367]	1480	& ec550aer)
[1364]	1481	ENDIF
[1337]	1482	ENDIF
	1483
[1146]	1484	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
[1539]	1485	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1486	$o_lwcon%name,itau_w,flwc)
[907]	1487	ENDIF
	1488
[1146]	1489	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
[1539]	1490	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1491	$o_iwcon%name,itau_w,fiwc)
[907]	1492	ENDIF
	1493
[1146]	1494	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
[1539]	1495	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1496	$o_temp%name,itau_w,t_seri)
[907]	1497	ENDIF
	1498
[1146]	1499	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
[1539]	1500	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1501	$o_theta%name,itau_w,theta)
[945]	1502	ENDIF
	1503
[1398]	1504	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
[1539]	1505	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1506	$o_ovapinit%name,itau_w,
[1398]	1507	$ qx(:,:,ivap))
[907]	1508	ENDIF
	1509
[1398]	1510	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
[1539]	1511	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1398]	1512	$ o_ovap%name,itau_w,q_seri)
[907]	1513	ENDIF
[1638]	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
[907]	1519
[1146]	1520	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
[1539]	1521	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1522	$o_geop%name,itau_w,zphi)
[907]	1523	ENDIF
	1524
[1146]	1525	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
[1539]	1526	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1527	$o_vitu%name,itau_w,u_seri)
[907]	1528	ENDIF
	1529
[1146]	1530	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
[1539]	1531	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1532	$o_vitv%name,itau_w,v_seri)
[907]	1533	ENDIF
	1534
[1146]	1535	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
[1539]	1536	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1537	$o_vitw%name,itau_w,omega)
[1146]	1538	ENDIF
	1539
	1540	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
[1539]	1541	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1542	$o_pres%name,itau_w,pplay)
[907]	1543	ENDIF
	1544
[1334]	1545	IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
[1539]	1546	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1547	$o_paprs%name,
[1334]	1548	s itau_w,paprs(:,1:klev))
	1549	ENDIF
	1550
[1539]	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
[1146]	1586	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
[1539]	1587	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1588	$o_rneb%name,itau_w,cldfra)
[907]	1589	ENDIF
	1590
[1146]	1591	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
[1539]	1592	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1593	$o_rnebcon%name,itau_w,rnebcon)
[907]	1594	ENDIF
	1595
[1146]	1596	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
[1539]	1597	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1598	$o_rhum%name,itau_w,zx_rh)
[907]	1599	ENDIF
	1600
[1146]	1601	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
[1539]	1602	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1603	$o_ozone%name, itau_w,
[1279]	1604	$ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
[907]	1605	ENDIF
	1606
[1279]	1607	IF (o_ozone_light%flag(iff)<=lev_files(iff) .and.
	1608	$ read_climoz == 2) THEN
[1539]	1609	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1610	$o_ozone_light%name, itau_w,
[1279]	1611	$ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
	1612	ENDIF
	1613
[1146]	1614	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
[1539]	1615	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1616	$o_dtphy%name,itau_w,d_t)
[907]	1617	ENDIF
	1618
[1146]	1619	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
[1539]	1620	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1621	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
[907]	1622	ENDIF
	1623
	1624	DO nsrf=1, nbsrf
[1146]	1625	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1626	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
[1539]	1627	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1628	s o_albe_srf(nsrf)%name,itau_w,
	1629	$ zx_tmp_fi2d)
[907]	1630	ENDIF
	1631
[1146]	1632	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1633	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
[1539]	1634	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1635	s o_rugs_srf(nsrf)%name,itau_w,
[907]	1636	$ zx_tmp_fi2d)
	1637	ENDIF
	1638
[1146]	1639	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1640	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
[1539]	1641	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1642	s o_ages_srf(nsrf)%name,itau_w
[907]	1643	$ ,zx_tmp_fi2d)
	1644	ENDIF
	1645	ENDDO !nsrf=1, nbsrf
	1646
[1398]	1647	IF (o_alb1%flag(iff)<=lev_files(iff)) THEN
[1539]	1648	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1649	$o_alb1%name,itau_w,albsol1)
[907]	1650	ENDIF
	1651
[1398]	1652	IF (o_alb2%flag(iff)<=lev_files(iff)) THEN
[1539]	1653	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1654	$o_alb2%name,itau_w,albsol2)
[907]	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
[1146]	1666	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
[1539]	1667	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1668	$o_tke%name,itau_w,zx_tmp_fi3d)
[907]	1669	ENDIF
	1670
[1539]	1671	IF (.NOT.clef_stations(iff)) THEN
[1146]	1672	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
[1539]	1673	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1674	s o_tke_max%name,itau_w,zx_tmp_fi3d)
[907]	1675	ENDIF
[1539]	1676	ENDIF
[907]	1677	endif
	1678
[1146]	1679	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
[1539]	1680	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1681	$o_kz%name,itau_w,coefh)
[907]	1682	ENDIF
	1683
[1539]	1684	IF (.NOT.clef_stations(iff)) THEN
[1146]	1685	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
[1539]	1686	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1687	$o_kz_max%name,itau_w,coefh)
[907]	1688	ENDIF
[1539]	1689	ENDIF
[907]	1690
[1146]	1691	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
[1539]	1692	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1693	$o_clwcon%name,itau_w,clwcon0)
[907]	1694	ENDIF
	1695
[1146]	1696	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
[1539]	1697	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1698	$o_dtdyn%name,itau_w,d_t_dyn)
[907]	1699	ENDIF
	1700
[1146]	1701	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
[1539]	1702	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1703	$o_dqdyn%name,itau_w,d_q_dyn)
[907]	1704	ENDIF
	1705
[1146]	1706	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
[1539]	1707	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1708	$o_dudyn%name,itau_w,d_u_dyn)
[1054]	1709	ENDIF
	1710
[1146]	1711	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
[1539]	1712	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1713	$o_dvdyn%name,itau_w,d_v_dyn)
[1054]	1714	ENDIF
	1715
[1146]	1716	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
[907]	1717	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
[1539]	1718	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1719	$o_dtcon%name,itau_w,zx_tmp_fi3d)
	1720	ENDIF
[907]	1721
[1539]	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
[1146]	1739	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
[907]	1740	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
[1539]	1741	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1742	$o_ducon%name,itau_w,zx_tmp_fi3d)
[907]	1743	ENDIF
	1744
[1638]	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
[1146]	1751	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
[907]	1752	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
[1539]	1753	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1754	$o_dqcon%name,itau_w,zx_tmp_fi3d)
[907]	1755	ENDIF
	1756
[1539]	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
[1146]	1773	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
[907]	1774	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
[1539]	1775	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1776	$o_dtlsc%name,itau_w,zx_tmp_fi3d)
[907]	1777	ENDIF
	1778
[1146]	1779	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
[907]	1780	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
	1781	$ d_t_eva(1:klon,1:klev))/pdtphys
[1539]	1782	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	1783	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
[907]	1784	ENDIF
	1785
[1146]	1786	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
[907]	1787	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
[1539]	1788	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1789	$o_dqlsc%name,itau_w,zx_tmp_fi3d)
[907]	1790	ENDIF
	1791
[1507]	1792	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1793	! Sorties specifiques a la separation thermiques/non thermiques
[1638]	1794	if (iflag_thermals>=1) then
[1507]	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
[1539]	1799	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1800	$o_dtlscth%name,itau_w,zx_tmp_fi3d)
[1507]	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
[1539]	1806	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1807	$o_dtlscst%name,itau_w,zx_tmp_fi3d)
[1507]	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
[1539]	1813	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1814	$o_dqlscth%name,itau_w,zx_tmp_fi3d)
[1507]	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
[1539]	1820	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1821	$o_dqlscst%name,itau_w,zx_tmp_fi3d)
[1507]	1822	ENDIF
	1823
	1824	IF (o_plulth%flag(iff)<=lev_files(iff)) THEN
[1539]	1825	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1826	$o_plulth%name,itau_w,plul_th)
[1507]	1827	ENDIF
	1828
	1829	IF (o_plulst%flag(iff)<=lev_files(iff)) THEN
[1539]	1830	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1831	$o_plulst%name,itau_w,plul_st)
[1507]	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
[1539]	1845	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1846	$o_ptconvth%name,itau_w,zx_tmp_fi3d)
[1507]	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
[1539]	1854	s histwrite_phy(nid_files(iff),clef_stations(iff),
	1855	$o_lmaxth%name,itau_w,zx_tmp_fi2d)
[1507]	1856	ENDIF
	1857
[1638]	1858	endif ! iflag_thermals>=1
[1507]	1859
	1860	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1861
[1146]	1862	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1863	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
[1539]	1864	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1865	$o_dtvdf%name,itau_w,zx_tmp_fi3d)
[907]	1866	ENDIF
	1867
[1146]	1868	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1869	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
[1539]	1870	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1871	$o_dqvdf%name,itau_w,zx_tmp_fi3d)
[907]	1872	ENDIF
	1873
[1146]	1874	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
[907]	1875	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
[1539]	1876	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1877	$o_dteva%name,itau_w,zx_tmp_fi3d)
[907]	1878	ENDIF
	1879
[1146]	1880	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
[907]	1881	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
[1539]	1882	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1883	$o_dqeva%name,itau_w,zx_tmp_fi3d)
[907]	1884	ENDIF
	1885
[1146]	1886	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
[907]	1887	zpt_conv = 0.
	1888	where (ptconv) zpt_conv = 1.
[1539]	1889	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1890	$o_ptconv%name,itau_w,zpt_conv)
[907]	1891	ENDIF
	1892
[1146]	1893	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
[1539]	1894	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1895	$o_ratqs%name,itau_w,ratqs)
[907]	1896	ENDIF
	1897
[1146]	1898	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
[1539]	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)
[907]	1903	ENDIF
	1904
[1638]	1905	IF (iflag_thermals>=1) THEN
[1334]	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.
[1539]	1909	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1910	$o_ftime_th%name,
[1334]	1911	s itau_w,zx_tmp_fi2d)
	1912	ENDIF
	1913
[1146]	1914	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1915	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1916	$o_f_th%name,itau_w,fm_therm)
[951]	1917	ENDIF
[945]	1918
[1146]	1919	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1920	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1921	$o_e_th%name,itau_w,entr_therm)
[951]	1922	ENDIF
[945]	1923
[1146]	1924	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1925	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1926	$o_w_th%name,itau_w,zw2)
[1028]	1927	ENDIF
	1928
[1146]	1929	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1930	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1931	$o_q_th%name,itau_w,zqasc)
[1028]	1932	ENDIF
	1933
	1934
[1146]	1935	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1936	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1937	$o_a_th%name,itau_w,fraca)
[1028]	1938	ENDIF
	1939
[1146]	1940	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1941	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1942	$o_d_th%name,itau_w,detr_therm)
[1028]	1943	ENDIF
	1944
[1146]	1945	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1946	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1947	$o_f0_th%name,itau_w,f0)
[945]	1948	ENDIF
	1949
[1507]	1950	IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN
[1539]	1951	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1952	$o_zmax_th%name,itau_w,zmax_th)
[945]	1953	ENDIF
	1954
[1146]	1955	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
[1539]	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)
[907]	1960	ENDIF
	1961
[1398]	1962	ENDIF !iflag_thermals
	1963
[1146]	1964	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
[907]	1965	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
[1539]	1966	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1967	$o_dtajs%name,itau_w,zx_tmp_fi3d)
[907]	1968	ENDIF
	1969
[1146]	1970	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
[907]	1971	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
[1539]	1972	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1973	$o_dqajs%name,itau_w,zx_tmp_fi3d)
[907]	1974	ENDIF
	1975
[1146]	1976	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
[907]	1977	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
[1539]	1978	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1979	$o_dtswr%name,itau_w,zx_tmp_fi3d)
[907]	1980	ENDIF
	1981
[1146]	1982	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
[907]	1983	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
[1539]	1984	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1985	$o_dtsw0%name,itau_w,zx_tmp_fi3d)
[907]	1986	ENDIF
	1987
[1146]	1988	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
[907]	1989	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
[1539]	1990	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1991	$o_dtlwr%name,itau_w,zx_tmp_fi3d)
[907]	1992	ENDIF
	1993
[1146]	1994	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
[907]	1995	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
[1539]	1996	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	1997	$o_dtlw0%name,itau_w,zx_tmp_fi3d)
[907]	1998	ENDIF
	1999
[1146]	2000	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
[907]	2001	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
[1539]	2002	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2003	$o_dtec%name,itau_w,zx_tmp_fi3d)
[907]	2004	ENDIF
	2005
[1146]	2006	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
[907]	2007	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[1539]	2008	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2009	$o_duvdf%name,itau_w,zx_tmp_fi3d)
[907]	2010	ENDIF
	2011
[1146]	2012	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
[907]	2013	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[1539]	2014	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2015	$o_dvvdf%name,itau_w,zx_tmp_fi3d)
[907]	2016	ENDIF
	2017
	2018	IF (ok_orodr) THEN
[1146]	2019	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
[907]	2020	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[1539]	2021	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2022	$o_duoro%name,itau_w,zx_tmp_fi3d)
[907]	2023	ENDIF
	2024
[1146]	2025	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
[907]	2026	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[1539]	2027	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2028	$o_dvoro%name,itau_w,zx_tmp_fi3d)
[907]	2029	ENDIF
[1539]	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
[907]	2036	ENDIF
	2037
	2038	IF (ok_orolf) THEN
[1146]	2039	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
[907]	2040	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[1539]	2041	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2042	$o_dulif%name,itau_w,zx_tmp_fi3d)
[907]	2043	ENDIF
	2044
[1146]	2045	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
[907]	2046	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[1539]	2047	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
	2048	$o_dvlif%name,itau_w,zx_tmp_fi3d)
[907]	2049	ENDIF
[1539]	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
[907]	2056	ENDIF
	2057
[1539]	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
[1403]	2255	if (nqtot.GE.3) THEN
	2256	DO iq=3,nqtot
[1146]	2257	IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN
[1539]	2258	CALL histwrite_phy(nid_files(iff),clef_stations(iff),
[1146]	2259	s o_trac(iq-2)%name,itau_w,qx(:,:,iq))
	2260	ENDIF
[1403]	2261	ENDDO
	2262	endif
[907]	2263
	2264	if (ok_sync) then
	2265	c$OMP MASTER
[1539]	2266	call histsync(nid_files(iff))
[907]	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: