Context Navigation

source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_write.h @ 1419

Last change on this file since 1419 was 1419, checked in by musat, 14 years ago
Modify levels for some outputs to diminsh the CMIP5 storage (JLD) Add Hines GWD tendencies duhin, dvhin (F.Lott) IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 52.6 KB

Rev	Line
[907]	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 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1146]	10	IF (o_phis%flag(iff)<=lev_files(iff)) THEN
[907]	11	CALL histwrite_phy(nid_files(iff),
[1146]	12	$ o_phis%name,itau_w,pphis)
[907]	13	ENDIF
	14
[1146]	15	IF (o_aire%flag(iff)<=lev_files(iff)) THEN
	16	CALL histwrite_phy(nid_files(iff),o_aire%name,itau_w,airephy)
[907]	17	ENDIF
	18
[1146]	19	IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN
[907]	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),
[1146]	24	$ o_contfracATM%name,itau_w,zx_tmp_fi2d)
[907]	25	ENDIF
	26
[1146]	27	IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN
	28	CALL histwrite_phy(nid_files(iff),o_contfracOR%name,itau_w,
[996]	29	$ pctsrf(:,is_ter))
[907]	30	ENDIF
	31
[1146]	32	IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN
[907]	33	CALL histwrite_phy(nid_files(iff),
[1146]	34	$ o_aireTER%name,itau_w,paire_ter)
[907]	35	ENDIF
	36
	37	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	38
[1146]	39	IF (o_flat%flag(iff)<=lev_files(iff)) THEN
	40	CALL histwrite_phy(nid_files(iff),o_flat%name,itau_w,zxfluxlat)
[907]	41	ENDIF
	42
[1146]	43	IF (o_slp%flag(iff)<=lev_files(iff)) THEN
	44	CALL histwrite_phy(nid_files(iff),o_slp%name,itau_w,slp)
[907]	45	ENDIF
	46
[1146]	47	IF (o_tsol%flag(iff)<=lev_files(iff)) THEN
	48	CALL histwrite_phy(nid_files(iff),o_tsol%name,itau_w,zxtsol)
[907]	49	ENDIF
	50
[1146]	51	IF (o_t2m%flag(iff)<=lev_files(iff)) THEN
	52	CALL histwrite_phy(nid_files(iff),o_t2m%name,itau_w,zt2m)
[907]	53	ENDIF
	54
[1146]	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)
[907]	57	ENDIF
	58
[1146]	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)
[907]	61	ENDIF
	62
[1146]	63	IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN
[907]	64	DO i=1, klon
	65	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	66	ENDDO
[1146]	67	CALL histwrite_phy(nid_files(iff),
	68	s o_wind10m%name,itau_w,zx_tmp_fi2d)
[907]	69	ENDIF
	70
[1146]	71	IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN
[907]	72	DO i=1, klon
	73	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	74	ENDDO
[1146]	75	CALL histwrite_phy(nid_files(iff),o_wind10max%name,
[907]	76	$ itau_w,zx_tmp_fi2d)
	77	ENDIF
	78
[1146]	79	IF (o_sicf%flag(iff)<=lev_files(iff)) THEN
[907]	80	DO i = 1, klon
	81	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	82	ENDDO
[1146]	83	CALL histwrite_phy(nid_files(iff),
	84	$ o_sicf%name,itau_w,zx_tmp_fi2d)
[907]	85	ENDIF
	86
[1146]	87	IF (o_q2m%flag(iff)<=lev_files(iff)) THEN
	88	CALL histwrite_phy(nid_files(iff),o_q2m%name,itau_w,zq2m)
[907]	89	ENDIF
	90
[1146]	91	IF (o_u10m%flag(iff)<=lev_files(iff)) THEN
	92	CALL histwrite_phy(nid_files(iff),o_u10m%name,itau_w,zu10m)
[907]	93	ENDIF
	94
[1146]	95	IF (o_v10m%flag(iff)<=lev_files(iff)) THEN
	96	CALL histwrite_phy(nid_files(iff),o_v10m%name,itau_w,zv10m)
[907]	97	ENDIF
	98
[1146]	99	IF (o_psol%flag(iff)<=lev_files(iff)) THEN
[907]	100	DO i = 1, klon
	101	zx_tmp_fi2d(i) = paprs(i,1)
	102	ENDDO
[1146]	103	CALL histwrite_phy(nid_files(iff),
	104	s o_psol%name,itau_w,zx_tmp_fi2d)
[907]	105	ENDIF
	106
[1146]	107	IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN
	108	CALL histwrite_phy(nid_files(iff),o_qsurf%name,itau_w,zxqsurf)
[907]	109	ENDIF
	110
	111	if (.not. ok_veget) then
[1146]	112	IF (o_qsol%flag(iff)<=lev_files(iff)) THEN
	113	CALL histwrite_phy(nid_files(iff),o_qsol%name,itau_w,qsol)
[907]	114	ENDIF
	115	endif
	116
[1146]	117	IF (o_precip%flag(iff)<=lev_files(iff)) THEN
[907]	118	DO i = 1, klon
	119	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
	120	ENDDO
[1146]	121	CALL histwrite_phy(nid_files(iff),o_precip%name,
	122	s itau_w,zx_tmp_fi2d)
[907]	123	ENDIF
	124
[1146]	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)
[907]	128	ENDIF
	129
[1146]	130	IF (o_plul%flag(iff)<=lev_files(iff)) THEN
[907]	131	DO i = 1, klon
	132	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	133	ENDDO
[1146]	134	CALL histwrite_phy(nid_files(iff),o_plul%name,itau_w,zx_tmp_fi2d)
[907]	135	ENDIF
	136
[1146]	137	IF (o_pluc%flag(iff)<=lev_files(iff)) THEN
[907]	138	DO i = 1, klon
	139	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	140	ENDDO
[1146]	141	CALL histwrite_phy(nid_files(iff),o_pluc%name,itau_w,zx_tmp_fi2d)
[907]	142	ENDIF
	143
[1146]	144	IF (o_snow%flag(iff)<=lev_files(iff)) THEN
	145	CALL histwrite_phy(nid_files(iff),o_snow%name,itau_w,snow_fall)
[907]	146	ENDIF
	147
[1334]	148	IF (o_msnow%flag(iff)<=lev_files(iff)) THEN
	149	CALL histwrite_phy(nid_files(iff),o_msnow%name,itau_w,snow_o)
	150	ENDIF
	151
	152	IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN
	153	CALL histwrite_phy(nid_files(iff),o_fsnow%name,itau_w,zfra_o)
	154	ENDIF
	155
[1146]	156	IF (o_evap%flag(iff)<=lev_files(iff)) THEN
	157	CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap)
[907]	158	ENDIF
	159
[1146]	160	IF (o_tops%flag(iff)<=lev_files(iff)) THEN
	161	CALL histwrite_phy(nid_files(iff),o_tops%name,itau_w,topsw)
[907]	162	ENDIF
	163
[1146]	164	IF (o_tops0%flag(iff)<=lev_files(iff)) THEN
	165	CALL histwrite_phy(nid_files(iff),o_tops0%name,itau_w,topsw0)
[907]	166	ENDIF
	167
[1146]	168	IF (o_topl%flag(iff)<=lev_files(iff)) THEN
	169	CALL histwrite_phy(nid_files(iff),o_topl%name,itau_w,toplw)
[907]	170	ENDIF
	171
[1146]	172	IF (o_topl0%flag(iff)<=lev_files(iff)) THEN
	173	CALL histwrite_phy(nid_files(iff),o_topl0%name,itau_w,toplw0)
[907]	174	ENDIF
	175
[1146]	176	IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN
[907]	177	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
[1146]	178	CALL histwrite_phy(nid_files(iff),o_SWupTOA%name,
	179	s itau_w,zx_tmp_fi2d)
[907]	180	ENDIF
	181
[1146]	182	IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN
[907]	183	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
	184	CALL histwrite_phy(nid_files(iff),
[1146]	185	$ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d)
[907]	186	ENDIF
	187
[1146]	188	IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN
[907]	189	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
[1146]	190	CALL histwrite_phy(nid_files(iff),
	191	s o_SWdnTOA%name,itau_w,zx_tmp_fi2d)
[907]	192	ENDIF
	193
[1146]	194	IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN
[907]	195	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
	196	CALL histwrite_phy(nid_files(iff),
[1146]	197	$ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d)
[907]	198	ENDIF
	199
[1334]	200	IF (o_nettop%flag(iff)<=lev_files(iff)) THEN
	201	zx_tmp_fi2d(:) = topsw(:)-toplw(:)
	202	CALL histwrite_phy(nid_files(iff),
	203	$ o_nettop%name,itau_w,zx_tmp_fi2d)
	204	ENDIF
	205
[1146]	206	IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN
	207	CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200)
[907]	208	ENDIF
	209
[1146]	210	IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN
	211	CALL histwrite_phy(nid_files(iff),
	212	s o_SWup200clr%name,itau_w,SWup200clr)
[907]	213	ENDIF
	214
[1146]	215	IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN
	216	CALL histwrite_phy(nid_files(iff),o_SWdn200%name,itau_w,SWdn200)
[907]	217	ENDIF
	218
[1146]	219	IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN
	220	CALL histwrite_phy(nid_files(iff),
	221	s o_SWdn200clr%name,itau_w,SWdn200clr)
[907]	222	ENDIF
	223
[1146]	224	IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN
	225	CALL histwrite_phy(nid_files(iff),o_LWup200%name,itau_w,LWup200)
[907]	226	ENDIF
	227
[1146]	228	IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN
	229	CALL histwrite_phy(nid_files(iff),
	230	s o_LWup200clr%name,itau_w,LWup200clr)
[907]	231	ENDIF
	232
[1146]	233	IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN
	234	CALL histwrite_phy(nid_files(iff),
[1279]	235	s o_LWdn200%name,itau_w,LWdn200)
[907]	236	ENDIF
	237
[1146]	238	IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN
	239	CALL histwrite_phy(nid_files(iff),
[1279]	240	s o_LWdn200clr%name,itau_w,LWdn200clr)
[907]	241	ENDIF
	242
[1146]	243	IF (o_sols%flag(iff)<=lev_files(iff)) THEN
	244	CALL histwrite_phy(nid_files(iff),o_sols%name,itau_w,solsw)
[907]	245	ENDIF
	246
[1146]	247	IF (o_sols0%flag(iff)<=lev_files(iff)) THEN
	248	CALL histwrite_phy(nid_files(iff),o_sols0%name,itau_w,solsw0)
[907]	249	ENDIF
	250
[1146]	251	IF (o_soll%flag(iff)<=lev_files(iff)) THEN
	252	CALL histwrite_phy(nid_files(iff),o_soll%name,itau_w,sollw)
[907]	253	ENDIF
	254
[1146]	255	IF (o_radsol%flag(iff)<=lev_files(iff)) THEN
	256	CALL histwrite_phy(nid_files(iff),o_radsol%name,itau_w,radsol)
[907]	257	ENDIF
	258
[1146]	259	IF (o_soll0%flag(iff)<=lev_files(iff)) THEN
	260	CALL histwrite_phy(nid_files(iff),o_soll0%name,itau_w,sollw0)
[907]	261	ENDIF
	262
[1146]	263	IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN
[907]	264	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
[1146]	265	CALL histwrite_phy(nid_files(iff),
	266	s o_SWupSFC%name,itau_w,zx_tmp_fi2d)
[907]	267	ENDIF
	268
[1146]	269	IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	270	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
	271	CALL histwrite_phy(nid_files(iff),
[1146]	272	$ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	273	ENDIF
	274
[1146]	275	IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN
[907]	276	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
	277	CALL histwrite_phy(nid_files(iff),
[1146]	278	$ o_SWdnSFC%name,itau_w,zx_tmp_fi2d)
[907]	279	ENDIF
	280
[1146]	281	IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	282	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
	283	CALL histwrite_phy(nid_files(iff),
[1146]	284	$ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	285	ENDIF
	286
[1146]	287	IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN
[907]	288	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
	289	CALL histwrite_phy(nid_files(iff),
[1146]	290	$ o_LWupSFC%name,itau_w,zx_tmp_fi2d)
[907]	291	ENDIF
	292
[1146]	293	IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN
[907]	294	CALL histwrite_phy(nid_files(iff),
[1146]	295	$ o_LWdnSFC%name,itau_w,sollwdown)
[907]	296	ENDIF
	297
[958]	298	sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
[1146]	299	IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	300	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
	301	CALL histwrite_phy(nid_files(iff),
[1146]	302	$ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d)
[907]	303	ENDIF
	304
[1146]	305	IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN
[907]	306	CALL histwrite_phy(nid_files(iff),
[1146]	307	$ o_LWdnSFCclr%name,itau_w,sollwdownclr)
[907]	308	ENDIF
	309
[1146]	310	IF (o_bils%flag(iff)<=lev_files(iff)) THEN
	311	CALL histwrite_phy(nid_files(iff),o_bils%name,itau_w,bils)
[907]	312	ENDIF
	313
[1146]	314	IF (o_sens%flag(iff)<=lev_files(iff)) THEN
[907]	315	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
[1146]	316	CALL histwrite_phy(nid_files(iff),o_sens%name,itau_w,zx_tmp_fi2d)
[907]	317	ENDIF
	318
[1146]	319	IF (o_fder%flag(iff)<=lev_files(iff)) THEN
	320	CALL histwrite_phy(nid_files(iff),o_fder%name,itau_w,fder)
[907]	321	ENDIF
	322
[1146]	323	IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN
	324	CALL histwrite_phy(nid_files(iff),o_ffonte%name,itau_w,zxffonte)
[907]	325	ENDIF
	326
[1146]	327	IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN
[907]	328	CALL histwrite_phy(nid_files(iff),
[1146]	329	$ o_fqcalving%name,itau_w,zxfqcalving)
[907]	330	ENDIF
	331
[1146]	332	IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN
[907]	333	CALL histwrite_phy(nid_files(iff),
[1146]	334	$ o_fqfonte%name,itau_w,zxfqfonte)
[907]	335	ENDIF
	336
[1334]	337	IF (o_taux%flag(iff)<=lev_files(iff)) THEN
	338	zx_tmp_fi2d=0.
	339	do nsrf=1,nbsrf
	340	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
	341	enddo
	342	CALL histwrite_phy(nid_files(iff),
	343	$ o_taux%name,itau_w,zx_tmp_fi2d)
	344	ENDIF
	345
	346	IF (o_tauy%flag(iff)<=lev_files(iff)) THEN
	347	zx_tmp_fi2d=0.
	348	do nsrf=1,nbsrf
	349	zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
	350	enddo
	351	CALL histwrite_phy(nid_files(iff),
	352	$ o_tauy%name,itau_w,zx_tmp_fi2d)
	353	ENDIF
	354
[907]	355	DO nsrf = 1, nbsrf
[1146]	356	! IF(nsrf.GE.2) THEN
	357	IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	358	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
	359	CALL histwrite_phy(nid_files(iff),
[1146]	360	$ o_pourc_srf(nsrf)%name,itau_w,
[907]	361	$ zx_tmp_fi2d)
	362	ENDIF
	363
[1146]	364	IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	365	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
	366	CALL histwrite_phy(nid_files(iff),
[1146]	367	$ o_fract_srf(nsrf)%name,itau_w,
	368	$ zx_tmp_fi2d)
[907]	369	ENDIF
[1146]	370	! ENDIF !nsrf.GT.2
[907]	371
[1146]	372	IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	373	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
	374	CALL histwrite_phy(nid_files(iff),
[1146]	375	$ o_taux_srf(nsrf)%name,itau_w,
	376	$ zx_tmp_fi2d)
[907]	377	ENDIF
	378
[1146]	379	IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	380	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
	381	CALL histwrite_phy(nid_files(iff),
[1146]	382	$ o_tauy_srf(nsrf)%name,itau_w,
	383	$ zx_tmp_fi2d)
[907]	384	ENDIF
	385
[1146]	386	IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	387	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
	388	CALL histwrite_phy(nid_files(iff),
[1146]	389	$ o_tsol_srf(nsrf)%name,itau_w,
[907]	390	$ zx_tmp_fi2d)
	391	ENDIF
	392
[1146]	393	IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	394	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
[1146]	395	CALL histwrite_phy(nid_files(iff),o_u10m_srf(nsrf)%name,
[907]	396	$ itau_w,zx_tmp_fi2d)
	397	ENDIF
	398
[1146]	399	IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	400	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
[1146]	401	CALL histwrite_phy(nid_files(iff),o_v10m_srf(nsrf)%name,
[907]	402	$ itau_w,zx_tmp_fi2d)
	403	ENDIF
	404
[1146]	405	IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	406	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
[1146]	407	CALL histwrite_phy(nid_files(iff),o_t2m_srf(nsrf)%name,
[907]	408	$ itau_w,zx_tmp_fi2d)
	409	ENDIF
	410
[1334]	411	IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
	412	zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
	413	CALL histwrite_phy(nid_files(iff),o_evap_srf(nsrf)%name,
	414	$ itau_w,zx_tmp_fi2d)
	415	ENDIF
	416
[1146]	417	IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	418	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
	419	CALL histwrite_phy(nid_files(iff),
[1146]	420	$ o_sens_srf(nsrf)%name,itau_w,
[907]	421	$ zx_tmp_fi2d)
	422	ENDIF
	423
[1146]	424	IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	425	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
	426	CALL histwrite_phy(nid_files(iff),
[1146]	427	$ o_lat_srf(nsrf)%name,itau_w,
[907]	428	$ zx_tmp_fi2d)
	429	ENDIF
	430
[1146]	431	IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	432	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
	433	CALL histwrite_phy(nid_files(iff),
[1146]	434	$ o_flw_srf(nsrf)%name,itau_w,
[907]	435	$ zx_tmp_fi2d)
	436	ENDIF
	437
[1146]	438	IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	439	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
	440	CALL histwrite_phy(nid_files(iff),
[1146]	441	$ o_fsw_srf(nsrf)%name,itau_w,
[907]	442	$ zx_tmp_fi2d)
	443	ENDIF
	444
[1146]	445	IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	446	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
	447	CALL histwrite_phy(nid_files(iff),
[1146]	448	$ o_wbils_srf(nsrf)%name,itau_w,
[907]	449	$ zx_tmp_fi2d)
	450	ENDIF
	451
[1146]	452	IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	453	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
	454	CALL histwrite_phy(nid_files(iff),
[1146]	455	$ o_wbilo_srf(nsrf)%name,itau_w,
[907]	456	$ zx_tmp_fi2d)
	457	ENDIF
	458
	459	if (iflag_pbl>1 .and. lev_histday.gt.10 ) then
[1146]	460	IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	461	CALL histwrite_phy(nid_files(iff),
[1146]	462	$ o_tke_srf(nsrf)%name,itau_w,
	463	$ pbl_tke(:,1:klev,nsrf))
[907]	464	ENDIF
	465
[1146]	466	IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	467	CALL histwrite_phy(nid_files(iff),
[1146]	468	$ o_tke_max_srf(nsrf)%name,itau_w,
[907]	469	$ pbl_tke(:,1:klev,nsrf))
	470	ENDIF
	471	endif
	472	ENDDO
	473
[1146]	474	IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN
	475	CALL histwrite_phy(nid_files(iff),o_cdrm%name,itau_w,cdragm)
[907]	476	ENDIF
	477
[1146]	478	IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN
	479	CALL histwrite_phy(nid_files(iff),o_cdrh%name,itau_w,cdragh)
[907]	480	ENDIF
	481
[1146]	482	IF (o_cldl%flag(iff)<=lev_files(iff)) THEN
	483	CALL histwrite_phy(nid_files(iff),o_cldl%name,itau_w,cldl)
[907]	484	ENDIF
	485
[1146]	486	IF (o_cldm%flag(iff)<=lev_files(iff)) THEN
	487	CALL histwrite_phy(nid_files(iff),o_cldm%name,itau_w,cldm)
[907]	488	ENDIF
	489
[1146]	490	IF (o_cldh%flag(iff)<=lev_files(iff)) THEN
	491	CALL histwrite_phy(nid_files(iff),o_cldh%name,itau_w,cldh)
[907]	492	ENDIF
	493
[1146]	494	IF (o_cldt%flag(iff)<=lev_files(iff)) THEN
	495	CALL histwrite_phy(nid_files(iff),o_cldt%name,
[1285]	496	& itau_w,cldt)
[907]	497	ENDIF
	498
[1146]	499	IF (o_cldq%flag(iff)<=lev_files(iff)) THEN
	500	CALL histwrite_phy(nid_files(iff),o_cldq%name,itau_w,cldq)
[907]	501	ENDIF
	502
[1146]	503	IF (o_lwp%flag(iff)<=lev_files(iff)) THEN
[907]	504	zx_tmp_fi2d(1:klon) = flwp(1:klon)
[1146]	505	CALL histwrite_phy(nid_files(iff),
	506	s o_lwp%name,itau_w,zx_tmp_fi2d)
[907]	507	ENDIF
	508
[1146]	509	IF (o_iwp%flag(iff)<=lev_files(iff)) THEN
[907]	510	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
[1146]	511	CALL histwrite_phy(nid_files(iff),
	512	s o_iwp%name,itau_w,zx_tmp_fi2d)
[907]	513	ENDIF
	514
[1146]	515	IF (o_ue%flag(iff)<=lev_files(iff)) THEN
	516	CALL histwrite_phy(nid_files(iff),o_ue%name,itau_w,ue)
[907]	517	ENDIF
	518
[1146]	519	IF (o_ve%flag(iff)<=lev_files(iff)) THEN
	520	CALL histwrite_phy(nid_files(iff),o_ve%name,itau_w,ve)
[907]	521	ENDIF
	522
[1146]	523	IF (o_uq%flag(iff)<=lev_files(iff)) THEN
	524	CALL histwrite_phy(nid_files(iff),o_uq%name,itau_w,uq)
[907]	525	ENDIF
	526
[1146]	527	IF (o_vq%flag(iff)<=lev_files(iff)) THEN
	528	CALL histwrite_phy(nid_files(iff),o_vq%name,itau_w,vq)
[907]	529	ENDIF
	530
	531	IF(iflag_con.GE.3) THEN ! sb
[1146]	532	IF (o_cape%flag(iff)<=lev_files(iff)) THEN
	533	CALL histwrite_phy(nid_files(iff),o_cape%name,itau_w,cape)
[907]	534	ENDIF
	535
[1146]	536	IF (o_pbase%flag(iff)<=lev_files(iff)) THEN
[1334]	537	CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w,ema_pcb)
[907]	538	ENDIF
	539
[1146]	540	IF (o_ptop%flag(iff)<=lev_files(iff)) THEN
	541	CALL histwrite_phy(nid_files(iff),o_ptop%name,itau_w,ema_pct)
[907]	542	ENDIF
	543
[1146]	544	IF (o_fbase%flag(iff)<=lev_files(iff)) THEN
	545	CALL histwrite_phy(nid_files(iff),o_fbase%name,itau_w,ema_cbmf)
[907]	546	ENDIF
	547
[1146]	548	IF (o_prw%flag(iff)<=lev_files(iff)) THEN
	549	CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw)
[907]	550	ENDIF
[1100]	551
[1146]	552	IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
	553	CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)
[1100]	554	ENDIF
	555
[1146]	556	IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
	557	CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)
	558	ENDIF
	559
	560	IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
	561	CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)
	562	ENDIF
	563
	564	IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
	565	CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)
	566	ENDIF
	567
	568	IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
	569	CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)
	570	ENDIF
	571
[1334]	572	IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN
	573	zx_tmp_fi2d=float(itau_con)/float(itap)
	574	CALL histwrite_phy(nid_files(iff),o_ftime_con%name,
	575	s itau_w,zx_tmp_fi2d)
	576	ENDIF
	577
	578	IF (o_mc%flag(iff)<=lev_files(iff)) THEN
	579	if(iflag_thermals.gt.1)then
	580	zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm
	581	else
	582	zx_tmp_fi3d=dnwd+dnwd0+upwd
	583	endif
	584	CALL histwrite_phy(nid_files(iff),o_mc%name,itau_w,zx_tmp_fi3d)
	585	ENDIF
	586
[907]	587	ENDIF !iflag_con .GE. 3
	588
[1146]	589	IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN
	590	CALL histwrite_phy(nid_files(iff),o_s_pblh%name,itau_w,s_pblh)
[907]	591	ENDIF
	592
[1146]	593	IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN
	594	CALL histwrite_phy(nid_files(iff),o_s_pblt%name,itau_w,s_pblt)
[907]	595	ENDIF
	596
[1146]	597	IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN
	598	CALL histwrite_phy(nid_files(iff),o_s_lcl%name,itau_w,s_lcl)
[907]	599	ENDIF
	600
[1398]	601	IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN
	602	CALL histwrite_phy(nid_files(iff),o_s_therm%name,itau_w,s_therm)
[907]	603	ENDIF
	604
[1398]	605	!IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
	606	! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
	607	! CALL histwrite_phy(nid_files(iff),o_s_capCL%name,itau_w,s_capCL)
	608	! ENDIF
[907]	609
[1398]	610	! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
	611	! CALL histwrite_phy(nid_files(iff),o_s_oliqCL%name,itau_w,s_oliqCL)
	612	! ENDIF
[907]	613
[1398]	614	! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
	615	! CALL histwrite_phy(nid_files(iff),o_s_cteiCL%name,itau_w,s_cteiCL)
	616	! ENDIF
[907]	617
[1398]	618	! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
	619	! CALL histwrite_phy(nid_files(iff),o_s_trmb1%name,itau_w,s_trmb1)
	620	! ENDIF
[907]	621
[1398]	622	! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
	623	! CALL histwrite_phy(nid_files(iff),o_s_trmb2%name,itau_w,s_trmb2)
	624	! ENDIF
[907]	625
[1398]	626	! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
	627	! CALL histwrite_phy(nid_files(iff),o_s_trmb3%name,itau_w,s_trmb3)
	628	! ENDIF
[907]	629
	630	! Champs interpolles sur des niveaux de pression
	631
[1146]	632	ll=0
[907]	633	DO k=1, nlevSTD
[1279]	634	! IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	635	! IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	636	bb2=clevSTD(k)
[907]	637	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
[1279]	638	$ bb2.EQ."500".OR.bb2.EQ."200".OR.
	639	$ bb2.EQ."50".OR.bb2.EQ."10") THEN
[907]	640
	641	! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1146]	642	ll=ll+1
	643	IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	644	CALL histwrite_phy(nid_files(iff),o_uSTDlevs(ll)%name,
	645	& itau_w,uwriteSTD(:,k,iff))
[907]	646	ENDIF
	647
[1146]	648	IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	649	CALL histwrite_phy(nid_files(iff),o_vSTDlevs(ll)%name,
	650	& itau_w,vwriteSTD(:,k,iff))
[907]	651	ENDIF
	652
[1146]	653	IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	654	CALL histwrite_phy(nid_files(iff),o_wSTDlevs(ll)%name,
	655	& itau_w,wwriteSTD(:,k,iff))
[907]	656	ENDIF
	657
[1398]	658	IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	659	CALL histwrite_phy(nid_files(iff),o_zSTDlevs(ll)%name,
[1055]	660	& itau_w,phiwriteSTD(:,k,iff))
[907]	661	ENDIF
	662
[1146]	663	IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	664	CALL histwrite_phy(nid_files(iff),o_qSTDlevs(ll)%name,
	665	& itau_w, qwriteSTD(:,k,iff))
[907]	666	ENDIF
	667
[1146]	668	IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN
	669	CALL histwrite_phy(nid_files(iff),o_tSTDlevs(ll)%name,
	670	& itau_w, twriteSTD(:,k,iff))
[907]	671	ENDIF
	672
	673	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	674	ENDDO
	675
[1146]	676	IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN
[907]	677	DO i=1, klon
	678	IF (pctsrf(i,is_oce).GT.epsfra.OR.
	679	$ pctsrf(i,is_sic).GT.epsfra) THEN
	680	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
	681	$ ftsol(i, is_sic) * pctsrf(i,is_sic))/
	682	$ (pctsrf(i,is_oce)+pctsrf(i,is_sic))
	683	ELSE
	684	zx_tmp_fi2d(i) = 273.15
	685	ENDIF
	686	ENDDO
[1146]	687	CALL histwrite_phy(nid_files(iff),
	688	s o_t_oce_sic%name,itau_w,zx_tmp_fi2d)
[907]	689	ENDIF
	690
[1146]	691	! Couplage convection-couche limite
	692	IF (iflag_con.GE.3) THEN
[951]	693	IF (iflag_coupl.EQ.1) THEN
[1146]	694	IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
	695	CALL histwrite_phy(nid_files(iff),o_ale_bl%name,itau_w,ale_bl)
[951]	696	ENDIF
[1146]	697	IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN
	698	CALL histwrite_phy(nid_files(iff),o_alp_bl%name,itau_w,alp_bl)
[951]	699	ENDIF
	700	ENDIF !iflag_coupl.EQ.1
[1146]	701	ENDIF !(iflag_con.GE.3)
[951]	702
[1146]	703	! Wakes
	704	IF (iflag_con.EQ.3) THEN
[951]	705	IF (iflag_wake.EQ.1) THEN
[1146]	706	IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
	707	CALL histwrite_phy(nid_files(iff),o_ale_wk%name,itau_w,ale_wake)
[951]	708	ENDIF
[1146]	709	IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN
	710	CALL histwrite_phy(nid_files(iff),o_alp_wk%name,itau_w,alp_wake)
[951]	711	ENDIF
[973]	712
[1146]	713	IF (o_ale%flag(iff)<=lev_files(iff)) THEN
	714	CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)
[973]	715	ENDIF
[1146]	716	IF (o_alp%flag(iff)<=lev_files(iff)) THEN
	717	CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)
[1063]	718	ENDIF
[1146]	719	IF (o_cin%flag(iff)<=lev_files(iff)) THEN
	720	CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)
	721	ENDIF
	722	IF (o_wape%flag(iff)<=lev_files(iff)) THEN
	723	CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe)
	724	ENDIF
	725	IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN
	726	CALL histwrite_phy(nid_files(iff),o_wake_h%name,itau_w,wake_h)
	727	ENDIF
[1063]	728
[1146]	729	IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN
	730	CALL histwrite_phy(nid_files(iff),o_wake_s%name,itau_w,wake_s)
[1063]	731	ENDIF
	732
[1146]	733	IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN
	734	CALL histwrite_phy(nid_files(iff),o_wake_deltat%name,
[1063]	735	$ itau_w,wake_deltat)
	736	ENDIF
	737
[1146]	738	IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN
	739	CALL histwrite_phy(nid_files(iff),o_wake_deltaq%name,
[1063]	740	$ itau_w,wake_deltaq)
	741	ENDIF
	742
[1146]	743	IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN
	744	CALL histwrite_phy(nid_files(iff),
	745	s o_wake_omg%name,itau_w,wake_omg)
[1063]	746	ENDIF
[1146]	747
	748	IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
	749	zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)
[1100]	750	& /pdtphys
[1146]	751	CALL histwrite_phy(nid_files(iff),
	752	& o_dtwak%name,itau_w,zx_tmp_fi3d)
	753	ENDIF
	754
	755	IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
	756	zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
	757	CALL histwrite_phy(nid_files(iff),
	758	& o_dqwak%name,itau_w,zx_tmp_fi3d)
[1100]	759	ENDIF
[1146]	760	ENDIF ! iflag_wake.EQ.1
[1063]	761
[1146]	762	IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN
	763	CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip)
[1063]	764	ENDIF
	765
[1146]	766	IF (o_ftd%flag(iff)<=lev_files(iff)) THEN
	767	CALL histwrite_phy(nid_files(iff),o_ftd%name,itau_w,ftd)
[1063]	768	ENDIF
	769
[1146]	770	IF (o_fqd%flag(iff)<=lev_files(iff)) THEN
	771	CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd)
[1063]	772	ENDIF
[1100]	773	ENDIF !(iflag_con.EQ.3)
[996]	774
[1146]	775	IF (type_ocean=='slab ') THEN
	776	IF ( o_slab_bils%flag(iff)<=lev_files(iff))
[996]	777	$ CALL histwrite_phy(
[1146]	778	$ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils)
[907]	779
[996]	780	ENDIF !type_ocean == force/slab
[907]	781
[1146]	782	IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN
	783	CALL histwrite_phy(nid_files(iff),
	784	s o_weakinv%name,itau_w,weak_inversion)
[907]	785	ENDIF
	786
[1146]	787	IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN
	788	CALL histwrite_phy(nid_files(iff),o_dthmin%name,itau_w,dthmin)
[907]	789	ENDIF
	790
[1146]	791	IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN
	792	CALL histwrite_phy(nid_files(iff),o_cldtau%name,itau_w,cldtau)
[907]	793	ENDIF
	794
[1146]	795	IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN
	796	CALL histwrite_phy(nid_files(iff),o_cldemi%name,itau_w,cldemi)
[907]	797	ENDIF
	798
[1279]	799	IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN
	800	CALL histwrite_phy(nid_files(iff),
	801	s o_pr_con_l%name,itau_w,pmflxr(:,1:klev))
	802	ENDIF
[929]	803
[1279]	804	IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN
	805	CALL histwrite_phy(nid_files(iff),
	806	s o_pr_con_i%name,itau_w,pmflxs(:,1:klev))
	807	ENDIF
[929]	808
[1279]	809	IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN
	810	CALL histwrite_phy(nid_files(iff),
	811	s o_pr_lsc_l%name,itau_w,prfl(:,1:klev))
	812	ENDIF
[929]	813
[1279]	814	IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN
	815	CALL histwrite_phy(nid_files(iff),
	816	s o_pr_lsc_i%name,itau_w,psfl(:,1:klev))
	817	ENDIF
[929]	818
[1284]	819	IF (o_re%flag(iff)<=lev_files(iff)) THEN
	820	CALL histwrite_phy(nid_files(iff),o_re%name,itau_w,re)
	821	ENDIF
[1279]	822
[1284]	823	IF (o_fl%flag(iff)<=lev_files(iff)) THEN
	824	CALL histwrite_phy(nid_files(iff),o_fl%name,itau_w,fl)
	825	ENDIF
	826
	827
	828
[1146]	829	IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN
[907]	830	DO i=1, klon
	831	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	832	ENDDO
[1146]	833	CALL histwrite_phy(nid_files(iff),o_rh2m%name,itau_w,zx_tmp_fi2d)
[907]	834	ENDIF
	835
[1334]	836	IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN
	837	DO i=1, klon
	838	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	839	ENDDO
	840	CALL histwrite_phy(nid_files(iff),o_rh2m_min%name,
	841	s itau_w,zx_tmp_fi2d)
	842	ENDIF
	843
	844	IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN
	845	DO i=1, klon
	846	zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
	847	ENDDO
	848	CALL histwrite_phy(nid_files(iff),o_rh2m_max%name,
	849	s itau_w,zx_tmp_fi2d)
	850	ENDIF
	851
	852
[1146]	853	IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN
	854	CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m)
[907]	855	ENDIF
	856
[1146]	857	IF (o_tpot%flag(iff)<=lev_files(iff)) THEN
	858	CALL histwrite_phy(nid_files(iff),o_tpot%name,itau_w,tpot)
[907]	859	ENDIF
	860
[1146]	861	IF (o_tpote%flag(iff)<=lev_files(iff)) THEN
	862	CALL histwrite_phy(nid_files(iff),o_tpote%name,itau_w,tpote)
[907]	863	ENDIF
	864
[1146]	865	IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN
[907]	866	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
[1146]	867	CALL histwrite_phy(nid_files(iff),
	868	s o_SWnetOR%name,itau_w, zx_tmp_fi2d)
[907]	869	ENDIF
	870
[1146]	871	IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN
[907]	872	zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
[1146]	873	CALL histwrite_phy(nid_files(iff),
	874	s o_SWdownOR%name,itau_w, zx_tmp_fi2d)
[907]	875	ENDIF
	876
[1146]	877	IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN
	878	CALL histwrite_phy(nid_files(iff),
	879	s o_LWdownOR%name,itau_w,sollwdown)
[907]	880	ENDIF
	881
[1146]	882	IF (o_snowl%flag(iff)<=lev_files(iff)) THEN
	883	CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc)
[907]	884	ENDIF
	885
[1146]	886	IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
	887	CALL histwrite_phy(nid_files(iff),
	888	s o_solldown%name,itau_w,sollwdown)
[907]	889	ENDIF
	890
[1146]	891	IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN
	892	CALL histwrite_phy(nid_files(iff),
	893	s o_dtsvdfo%name,itau_w,d_ts(:,is_oce))
[907]	894	ENDIF
	895
[1146]	896	IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN
	897	CALL histwrite_phy(nid_files(iff),
	898	s o_dtsvdft%name,itau_w,d_ts(:,is_ter))
[907]	899	ENDIF
	900
[1146]	901	IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN
[907]	902	CALL histwrite_phy(nid_files(iff),
[1146]	903	$ o_dtsvdfg%name,itau_w, d_ts(:,is_lic))
[907]	904	ENDIF
	905
[1146]	906	IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN
	907	CALL histwrite_phy(nid_files(iff),
	908	s o_dtsvdfi%name,itau_w,d_ts(:,is_sic))
[907]	909	ENDIF
	910
[1146]	911	IF (o_rugs%flag(iff)<=lev_files(iff)) THEN
	912	CALL histwrite_phy(nid_files(iff),o_rugs%name,itau_w,zxrugs)
[907]	913	ENDIF
	914
[1279]	915	! OD550 per species
	916	IF (new_aod .and. (.not. aerosol_couple)) THEN
[1345]	917	IF (ok_ade.OR.ok_aie) THEN
[1337]	918
	919	IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN
	920	CALL histwrite_phy(nid_files(iff),o_od550aer%name,itau_w,
	921	$ od550aer)
	922	ENDIF
	923	IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN
	924	CALL histwrite_phy(nid_files(iff),o_od865aer%name,itau_w,
	925	$ od865aer)
	926	ENDIF
	927	IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN
	928	CALL histwrite_phy(nid_files(iff),o_absvisaer%name,itau_w,
	929	$ absvisaer)
	930	ENDIF
	931	IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN
	932	CALL histwrite_phy(nid_files(iff),o_od550lt1aer%name,itau_w,
	933	$ od550lt1aer)
	934	ENDIF
	935
[1345]	936	IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN
	937	CALL histwrite_phy(nid_files(iff),o_sconcso4%name,itau_w,
	938	$ sconcso4)
[1279]	939	ENDIF
[1345]	940	IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN
	941	CALL histwrite_phy(nid_files(iff),o_sconcoa%name,itau_w,
	942	$ sconcoa)
	943	ENDIF
	944	IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN
	945	CALL histwrite_phy(nid_files(iff),o_sconcbc%name,itau_w,
	946	$ sconcbc)
	947	ENDIF
	948	IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN
	949	CALL histwrite_phy(nid_files(iff),o_sconcss%name,itau_w,
	950	$ sconcss)
	951	ENDIF
	952	IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN
	953	CALL histwrite_phy(nid_files(iff),o_sconcdust%name,itau_w,
	954	$ sconcdust)
	955	ENDIF
	956
	957	IF (o_concso4%flag(iff)<=lev_files(iff)) THEN
	958	CALL histwrite_phy(nid_files(iff),o_concso4%name,itau_w,
	959	$ concso4)
	960	ENDIF
	961	IF (o_concoa%flag(iff)<=lev_files(iff)) THEN
	962	CALL histwrite_phy(nid_files(iff),o_concoa%name,itau_w,
	963	$ concoa)
	964	ENDIF
	965	IF (o_concbc%flag(iff)<=lev_files(iff)) THEN
	966	CALL histwrite_phy(nid_files(iff),o_concbc%name,itau_w,
	967	$ concbc)
	968	ENDIF
	969	IF (o_concss%flag(iff)<=lev_files(iff)) THEN
	970	CALL histwrite_phy(nid_files(iff),o_concss%name,itau_w,
	971	$ concss)
	972	ENDIF
	973	IF (o_concdust%flag(iff)<=lev_files(iff)) THEN
	974	CALL histwrite_phy(nid_files(iff),o_concdust%name,itau_w,
	975	$ concdust)
	976	ENDIF
	977
	978	IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN
	979	CALL histwrite_phy(nid_files(iff),o_loadso4%name,itau_w,
	980	$ loadso4)
	981	ENDIF
	982	IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN
	983	CALL histwrite_phy(nid_files(iff),o_loadoa%name,itau_w,
	984	$ loadoa)
	985	ENDIF
	986	IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN
	987	CALL histwrite_phy(nid_files(iff),o_loadbc%name,itau_w,
	988	$ loadbc)
	989	ENDIF
	990	IF (o_loadss%flag(iff)<=lev_files(iff)) THEN
	991	CALL histwrite_phy(nid_files(iff),o_loadss%name,itau_w,
	992	$ loadss)
	993	ENDIF
	994	IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN
	995	CALL histwrite_phy(nid_files(iff),o_loaddust%name,itau_w,
	996	$ loaddust)
	997	ENDIF
	998
	999	DO naero = 1, naero_spc
	1000	IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
	1001	CALL histwrite_phy(nid_files(iff),
	1002	$ o_tausumaero(naero)%name,itau_w,
	1003	$ tausum_aero(:,2,naero) )
	1004	ENDIF
	1005	END DO
	1006	endif
[1279]	1007	ENDIF
[1345]	1008
[1146]	1009	IF (ok_ade) THEN
[1279]	1010	IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
	1011	CALL histwrite_phy(nid_files(iff),o_topswad%name,itau_w,
	1012	$ topswad_aero)
	1013	ENDIF
	1014	IF (o_solswad%flag(iff)<=lev_files(iff)) THEN
	1015	CALL histwrite_phy(nid_files(iff),o_solswad%name,itau_w,
	1016	$ solswad_aero)
	1017	ENDIF
	1018
	1019	!====MS forcing diagnostics
	1020	if (new_aod) then
	1021	IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
	1022	CALL histwrite_phy(nid_files(iff),o_swtoaas_nat%name,itau_w,
	1023	$ topsw_aero(:,1))
[1146]	1024	ENDIF
[1279]	1025
	1026	IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
	1027	CALL histwrite_phy(nid_files(iff),o_swsrfas_nat%name,itau_w,
	1028	$ solsw_aero(:,1))
[1146]	1029	ENDIF
	1030
[1279]	1031	IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
	1032	CALL histwrite_phy(nid_files(iff),o_swtoacs_nat%name,itau_w,
	1033	$ topsw0_aero(:,1))
[1146]	1034	ENDIF
[1279]	1035
	1036	IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
	1037	CALL histwrite_phy(nid_files(iff),o_swsrfcs_nat%name,itau_w,
	1038	$ solsw0_aero(:,1))
[1146]	1039	ENDIF
[1279]	1040
	1041	!ant
	1042	IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
	1043	CALL histwrite_phy(nid_files(iff),o_swtoaas_ant%name,itau_w,
	1044	$ topsw_aero(:,2))
	1045	ENDIF
	1046
	1047	IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
	1048	CALL histwrite_phy(nid_files(iff),o_swsrfas_ant%name,itau_w,
	1049	$ solsw_aero(:,2))
	1050	ENDIF
	1051
	1052	IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
	1053	CALL histwrite_phy(nid_files(iff),o_swtoacs_ant%name,itau_w,
	1054	$ topsw0_aero(:,2))
	1055	ENDIF
	1056
	1057	IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
	1058	CALL histwrite_phy(nid_files(iff),o_swsrfcs_ant%name,itau_w,
	1059	$ solsw0_aero(:,2))
	1060	ENDIF
	1061
	1062	!cf
	1063
	1064	if (.not. aerosol_couple) then
	1065	IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
	1066	CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
	1067	$ topswcf_aero(:,1))
	1068	ENDIF
	1069
	1070	IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
	1071	CALL histwrite_phy(nid_files(iff),o_swsrfcf_nat%name,itau_w,
	1072	$ solswcf_aero(:,1))
	1073	ENDIF
	1074
	1075	IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
	1076	CALL histwrite_phy(nid_files(iff),o_swtoacf_ant%name,itau_w,
	1077	$ topswcf_aero(:,2))
	1078	ENDIF
	1079
	1080	IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
	1081	CALL histwrite_phy(nid_files(iff),o_swsrfcf_ant%name,itau_w,
	1082	$ solswcf_aero(:,2))
	1083	ENDIF
	1084
	1085	IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
	1086	CALL histwrite_phy(nid_files(iff),o_swtoacf_zero%name,itau_w,
	1087	$ topswcf_aero(:,3))
	1088	ENDIF
	1089
	1090	IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
	1091	CALL histwrite_phy(nid_files(iff),o_swsrfcf_zero%name,itau_w,
	1092	$ solswcf_aero(:,3))
	1093	ENDIF
	1094	endif
	1095
	1096	endif ! new_aod
	1097	!====MS forcing diagnostics
	1098
[1146]	1099	ENDIF
	1100
[1279]	1101	IF (ok_aie) THEN
	1102	IF (o_topswai%flag(iff)<=lev_files(iff)) THEN
	1103	CALL histwrite_phy(nid_files(iff),o_topswai%name,itau_w,
	1104	$ topswai_aero)
	1105	ENDIF
	1106	IF (o_solswai%flag(iff)<=lev_files(iff)) THEN
	1107	CALL histwrite_phy(nid_files(iff),o_solswai%name,itau_w,
	1108	$ solswai_aero)
	1109	ENDIF
[1337]	1110	IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN
	1111	CALL histwrite_phy(nid_files(iff),o_scdnc%name,itau_w,
	1112	$ scdnc)
	1113	ENDIF
	1114	IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN
	1115	CALL histwrite_phy(nid_files(iff),o_cldncl%name,itau_w,
	1116	$ cldncl)
	1117	ENDIF
	1118	IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN
	1119	CALL histwrite_phy(nid_files(iff),o_reffclws%name,itau_w,
	1120	$ reffclws)
	1121	ENDIF
	1122	IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN
	1123	CALL histwrite_phy(nid_files(iff),o_reffclwc%name,itau_w,
	1124	$ reffclwc)
	1125	ENDIF
	1126	IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN
	1127	CALL histwrite_phy(nid_files(iff),o_cldnvi%name,itau_w,
	1128	$ cldnvi)
	1129	ENDIF
	1130	IF (o_lcc%flag(iff)<=lev_files(iff)) THEN
	1131	CALL histwrite_phy(nid_files(iff),o_lcc%name,itau_w,
	1132	$ lcc)
	1133	ENDIF
	1134	IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN
	1135	CALL histwrite_phy(nid_files(iff),o_lcc3d%name,itau_w,
	1136	$ lcc3d)
	1137	ENDIF
	1138	IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN
	1139	CALL histwrite_phy(nid_files(iff),o_lcc3dcon%name,itau_w,
	1140	$ lcc3dcon)
	1141	ENDIF
	1142	IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN
	1143	CALL histwrite_phy(nid_files(iff),o_lcc3dstra%name,itau_w,
	1144	$ lcc3dstra)
	1145	ENDIF
	1146	IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN
	1147	CALL histwrite_phy(nid_files(iff),o_reffclwtop%name,itau_w,
	1148	$ reffclwtop)
	1149	ENDIF
[1279]	1150	ENDIF
	1151
[907]	1152	! Champs 3D:
[1364]	1153	IF (ok_ade .OR. ok_aie) then
	1154	IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN
[1367]	1155	CALL histwrite_phy(nid_files(iff),o_ec550aer%name,itau_w,
	1156	& ec550aer)
[1364]	1157	ENDIF
[1337]	1158	ENDIF
	1159
[1146]	1160	IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN
	1161	CALL histwrite_phy(nid_files(iff),o_lwcon%name,itau_w,flwc)
[907]	1162	ENDIF
	1163
[1146]	1164	IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN
	1165	CALL histwrite_phy(nid_files(iff),o_iwcon%name,itau_w,fiwc)
[907]	1166	ENDIF
	1167
[1146]	1168	IF (o_temp%flag(iff)<=lev_files(iff)) THEN
	1169	CALL histwrite_phy(nid_files(iff),o_temp%name,itau_w,t_seri)
[907]	1170	ENDIF
	1171
[1146]	1172	IF (o_theta%flag(iff)<=lev_files(iff)) THEN
	1173	CALL histwrite_phy(nid_files(iff),o_theta%name,itau_w,theta)
[945]	1174	ENDIF
	1175
[1398]	1176	IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN
	1177	CALL histwrite_phy(nid_files(iff),o_ovapinit%name,itau_w,
	1178	$ qx(:,:,ivap))
[907]	1179	ENDIF
	1180
[1398]	1181	IF (o_ovap%flag(iff)<=lev_files(iff)) THEN
[1146]	1182	CALL histwrite_phy(nid_files(iff),
[1398]	1183	$ o_ovap%name,itau_w,q_seri)
[907]	1184	ENDIF
	1185
[1146]	1186	IF (o_geop%flag(iff)<=lev_files(iff)) THEN
	1187	CALL histwrite_phy(nid_files(iff),o_geop%name,itau_w,zphi)
[907]	1188	ENDIF
	1189
[1146]	1190	IF (o_vitu%flag(iff)<=lev_files(iff)) THEN
	1191	CALL histwrite_phy(nid_files(iff),o_vitu%name,itau_w,u_seri)
[907]	1192	ENDIF
	1193
[1146]	1194	IF (o_vitv%flag(iff)<=lev_files(iff)) THEN
	1195	CALL histwrite_phy(nid_files(iff),o_vitv%name,itau_w,v_seri)
[907]	1196	ENDIF
	1197
[1146]	1198	IF (o_vitw%flag(iff)<=lev_files(iff)) THEN
	1199	CALL histwrite_phy(nid_files(iff),o_vitw%name,itau_w,omega)
	1200	ENDIF
	1201
	1202	IF (o_pres%flag(iff)<=lev_files(iff)) THEN
	1203	CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay)
[907]	1204	ENDIF
	1205
[1334]	1206	IF (o_paprs%flag(iff)<=lev_files(iff)) THEN
	1207	CALL histwrite_phy(nid_files(iff),o_paprs%name,
	1208	s itau_w,paprs(:,1:klev))
	1209	ENDIF
	1210
[1146]	1211	IF (o_rneb%flag(iff)<=lev_files(iff)) THEN
	1212	CALL histwrite_phy(nid_files(iff),o_rneb%name,itau_w,cldfra)
[907]	1213	ENDIF
	1214
[1146]	1215	IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN
	1216	CALL histwrite_phy(nid_files(iff),o_rnebcon%name,itau_w,rnebcon)
[907]	1217	ENDIF
	1218
[1146]	1219	IF (o_rhum%flag(iff)<=lev_files(iff)) THEN
	1220	CALL histwrite_phy(nid_files(iff),o_rhum%name,itau_w,zx_rh)
[907]	1221	ENDIF
	1222
[1146]	1223	IF (o_ozone%flag(iff)<=lev_files(iff)) THEN
[1279]	1224	CALL histwrite_phy(nid_files(iff), o_ozone%name, itau_w,
	1225	$ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
[907]	1226	ENDIF
	1227
[1279]	1228	IF (o_ozone_light%flag(iff)<=lev_files(iff) .and.
	1229	$ read_climoz == 2) THEN
	1230	CALL histwrite_phy(nid_files(iff), o_ozone_light%name, itau_w,
	1231	$ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)
	1232	ENDIF
	1233
[1146]	1234	IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
	1235	CALL histwrite_phy(nid_files(iff),o_dtphy%name,itau_w,d_t)
[907]	1236	ENDIF
	1237
[1146]	1238	IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN
	1239	CALL histwrite_phy(nid_files(iff),
	1240	s o_dqphy%name,itau_w, d_qx(:,:,ivap))
[907]	1241	ENDIF
	1242
	1243	DO nsrf=1, nbsrf
[1146]	1244	IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1245	zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
[1146]	1246	CALL histwrite_phy(nid_files(iff),
	1247	s o_albe_srf(nsrf)%name,itau_w,
	1248	$ zx_tmp_fi2d)
[907]	1249	ENDIF
	1250
[1146]	1251	IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1252	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
[1146]	1253	CALL histwrite_phy(nid_files(iff),
	1254	s o_rugs_srf(nsrf)%name,itau_w,
[907]	1255	$ zx_tmp_fi2d)
	1256	ENDIF
	1257
[1146]	1258	IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN
[907]	1259	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
[1146]	1260	CALL histwrite_phy(nid_files(iff),
	1261	s o_ages_srf(nsrf)%name,itau_w
[907]	1262	$ ,zx_tmp_fi2d)
	1263	ENDIF
	1264	ENDDO !nsrf=1, nbsrf
	1265
[1398]	1266	IF (o_alb1%flag(iff)<=lev_files(iff)) THEN
	1267	CALL histwrite_phy(nid_files(iff),o_alb1%name,itau_w,albsol1)
[907]	1268	ENDIF
	1269
[1398]	1270	IF (o_alb2%flag(iff)<=lev_files(iff)) THEN
	1271	CALL histwrite_phy(nid_files(iff),o_alb2%name,itau_w,albsol2)
[907]	1272	ENDIF
	1273
	1274	!FH Sorties pour la couche limite
	1275	if (iflag_pbl>1) then
	1276	zx_tmp_fi3d=0.
	1277	do nsrf=1,nbsrf
	1278	do k=1,klev
	1279	zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k)
	1280	$ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
	1281	enddo
	1282	enddo
[1146]	1283	IF (o_tke%flag(iff)<=lev_files(iff)) THEN
	1284	CALL histwrite_phy(nid_files(iff),o_tke%name,itau_w,zx_tmp_fi3d)
[907]	1285	ENDIF
	1286
[1146]	1287	IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN
	1288	CALL histwrite_phy(nid_files(iff),
	1289	s o_tke_max%name,itau_w,zx_tmp_fi3d)
[907]	1290	ENDIF
	1291	endif
	1292
[1146]	1293	IF (o_kz%flag(iff)<=lev_files(iff)) THEN
	1294	CALL histwrite_phy(nid_files(iff),o_kz%name,itau_w,coefh)
[907]	1295	ENDIF
	1296
[1146]	1297	IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN
	1298	CALL histwrite_phy(nid_files(iff),o_kz_max%name,itau_w,coefh)
[907]	1299	ENDIF
	1300
[1146]	1301	IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN
	1302	CALL histwrite_phy(nid_files(iff),o_clwcon%name,itau_w,clwcon0)
[907]	1303	ENDIF
	1304
[1146]	1305	IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
	1306	CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn)
[907]	1307	ENDIF
	1308
[1146]	1309	IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN
	1310	CALL histwrite_phy(nid_files(iff),o_dqdyn%name,itau_w,d_q_dyn)
[907]	1311	ENDIF
	1312
[1146]	1313	IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN
	1314	CALL histwrite_phy(nid_files(iff),o_dudyn%name,itau_w,d_u_dyn)
[1054]	1315	ENDIF
	1316
[1146]	1317	IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN
	1318	CALL histwrite_phy(nid_files(iff),o_dvdyn%name,itau_w,d_v_dyn)
[1054]	1319	ENDIF
	1320
[1146]	1321	IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN
[907]	1322	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
[1146]	1323	CALL histwrite_phy(nid_files(iff),o_dtcon%name,itau_w,zx_tmp_fi3d)
[907]	1324	ENDIF
	1325
[1146]	1326	IF (o_ducon%flag(iff)<=lev_files(iff)) THEN
[907]	1327	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
[1146]	1328	CALL histwrite_phy(nid_files(iff),o_ducon%name,itau_w,zx_tmp_fi3d)
[907]	1329	ENDIF
	1330
[1146]	1331	IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN
[907]	1332	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
[1146]	1333	CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d)
[907]	1334	ENDIF
	1335
[1146]	1336	IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
[907]	1337	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
[1146]	1338	CALL histwrite_phy(nid_files(iff),o_dtlsc%name,itau_w,zx_tmp_fi3d)
[907]	1339	ENDIF
	1340
[1146]	1341	IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN
[907]	1342	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
	1343	$ d_t_eva(1:klon,1:klev))/pdtphys
[1146]	1344	CALL histwrite_phy(nid_files(iff),
	1345	s o_dtlschr%name,itau_w,zx_tmp_fi3d)
[907]	1346	ENDIF
	1347
[1146]	1348	IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN
[907]	1349	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
[1146]	1350	CALL histwrite_phy(nid_files(iff),o_dqlsc%name,itau_w,zx_tmp_fi3d)
[907]	1351	ENDIF
	1352
[1146]	1353	IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1354	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
[1146]	1355	CALL histwrite_phy(nid_files(iff),o_dtvdf%name,itau_w,zx_tmp_fi3d)
[907]	1356	ENDIF
	1357
[1146]	1358	IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1359	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
[1146]	1360	CALL histwrite_phy(nid_files(iff),o_dqvdf%name,itau_w,zx_tmp_fi3d)
[907]	1361	ENDIF
	1362
[1146]	1363	IF (o_dteva%flag(iff)<=lev_files(iff)) THEN
[907]	1364	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
[1146]	1365	CALL histwrite_phy(nid_files(iff),o_dteva%name,itau_w,zx_tmp_fi3d)
[907]	1366	ENDIF
	1367
[1146]	1368	IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN
[907]	1369	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
[1146]	1370	CALL histwrite_phy(nid_files(iff),o_dqeva%name,itau_w,zx_tmp_fi3d)
[907]	1371	ENDIF
	1372
[1146]	1373	IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN
[907]	1374	zpt_conv = 0.
	1375	where (ptconv) zpt_conv = 1.
[1146]	1376	CALL histwrite_phy(nid_files(iff),o_ptconv%name,itau_w,zpt_conv)
[907]	1377	ENDIF
	1378
[1146]	1379	IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN
	1380	CALL histwrite_phy(nid_files(iff),o_ratqs%name,itau_w,ratqs)
[907]	1381	ENDIF
	1382
[1146]	1383	IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN
[907]	1384	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
[1146]	1385	CALL histwrite_phy(nid_files(iff),o_dtthe%name,itau_w,zx_tmp_fi3d)
[907]	1386	ENDIF
	1387
[1036]	1388	IF (iflag_thermals.gt.1) THEN
[1334]	1389	IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN
	1390	! Pour l instant 0 a y reflichir pour les thermiques
	1391	zx_tmp_fi2d=0.
	1392	CALL histwrite_phy(nid_files(iff),o_ftime_th%name,
	1393	s itau_w,zx_tmp_fi2d)
	1394	ENDIF
	1395
[1146]	1396	IF (o_f_th%flag(iff)<=lev_files(iff)) THEN
	1397	CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm)
[951]	1398	ENDIF
[945]	1399
[1146]	1400	IF (o_e_th%flag(iff)<=lev_files(iff)) THEN
	1401	CALL histwrite_phy(nid_files(iff),o_e_th%name,itau_w,entr_therm)
[951]	1402	ENDIF
[945]	1403
[1146]	1404	IF (o_w_th%flag(iff)<=lev_files(iff)) THEN
	1405	CALL histwrite_phy(nid_files(iff),o_w_th%name,itau_w,zw2)
[1028]	1406	ENDIF
	1407
[1146]	1408	IF (o_q_th%flag(iff)<=lev_files(iff)) THEN
	1409	CALL histwrite_phy(nid_files(iff),o_q_th%name,itau_w,zqasc)
[1028]	1410	ENDIF
	1411
[1146]	1412	IF (o_lambda_th%flag(iff)<=lev_files(iff)) THEN
	1413	CALL histwrite_phy(nid_files(iff),
	1414	s o_lambda_th%name,itau_w,lambda_th)
[1028]	1415	ENDIF
	1416
[1146]	1417	IF (o_a_th%flag(iff)<=lev_files(iff)) THEN
	1418	CALL histwrite_phy(nid_files(iff),o_a_th%name,itau_w,fraca)
[1028]	1419	ENDIF
	1420
[1146]	1421	IF (o_d_th%flag(iff)<=lev_files(iff)) THEN
	1422	CALL histwrite_phy(nid_files(iff),o_d_th%name,itau_w,detr_therm)
[1028]	1423	ENDIF
	1424
[1398]	1425	!IM ENDIF !iflag_thermals
[945]	1426
[1146]	1427	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
	1428	CALL histwrite_phy(nid_files(iff),o_f0_th%name,itau_w,f0)
[945]	1429	ENDIF
	1430
[1146]	1431	IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
	1432	CALL histwrite_phy(nid_files(iff),o_zmax_th%name,itau_w,zmax0)
[945]	1433	ENDIF
	1434
[1146]	1435	IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN
[907]	1436	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
[1146]	1437	CALL histwrite_phy(nid_files(iff),o_dqthe%name,itau_w,zx_tmp_fi3d)
[907]	1438	ENDIF
	1439
[1398]	1440	ENDIF !iflag_thermals
	1441
[1146]	1442	IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN
[907]	1443	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
[1146]	1444	CALL histwrite_phy(nid_files(iff),o_dtajs%name,itau_w,zx_tmp_fi3d)
[907]	1445	ENDIF
	1446
[1146]	1447	IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN
[907]	1448	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
[1146]	1449	CALL histwrite_phy(nid_files(iff),o_dqajs%name,itau_w,zx_tmp_fi3d)
[907]	1450	ENDIF
	1451
[1146]	1452	IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN
[907]	1453	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
[1146]	1454	CALL histwrite_phy(nid_files(iff),o_dtswr%name,itau_w,zx_tmp_fi3d)
[907]	1455	ENDIF
	1456
[1146]	1457	IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN
[907]	1458	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
[1146]	1459	CALL histwrite_phy(nid_files(iff),o_dtsw0%name,itau_w,zx_tmp_fi3d)
[907]	1460	ENDIF
	1461
[1146]	1462	IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN
[907]	1463	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
[1146]	1464	CALL histwrite_phy(nid_files(iff),o_dtlwr%name,itau_w,zx_tmp_fi3d)
[907]	1465	ENDIF
	1466
[1146]	1467	IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN
[907]	1468	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
[1146]	1469	CALL histwrite_phy(nid_files(iff),o_dtlw0%name,itau_w,zx_tmp_fi3d)
[907]	1470	ENDIF
	1471
[1146]	1472	IF (o_dtec%flag(iff)<=lev_files(iff)) THEN
[907]	1473	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
[1146]	1474	CALL histwrite_phy(nid_files(iff),o_dtec%name,itau_w,zx_tmp_fi3d)
[907]	1475	ENDIF
	1476
[1146]	1477	IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1478	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[1146]	1479	CALL histwrite_phy(nid_files(iff),o_duvdf%name,itau_w,zx_tmp_fi3d)
[907]	1480	ENDIF
	1481
[1146]	1482	IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN
[907]	1483	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[1146]	1484	CALL histwrite_phy(nid_files(iff),o_dvvdf%name,itau_w,zx_tmp_fi3d)
[907]	1485	ENDIF
	1486
	1487	IF (ok_orodr) THEN
[1146]	1488	IF (o_duoro%flag(iff)<=lev_files(iff)) THEN
[907]	1489	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[1146]	1490	CALL histwrite_phy(nid_files(iff),o_duoro%name,itau_w,zx_tmp_fi3d)
[907]	1491	ENDIF
	1492
[1146]	1493	IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN
[907]	1494	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[1146]	1495	CALL histwrite_phy(nid_files(iff),o_dvoro%name,itau_w,zx_tmp_fi3d)
[907]	1496	ENDIF
	1497	ENDIF
	1498
	1499	IF (ok_orolf) THEN
[1146]	1500	IF (o_dulif%flag(iff)<=lev_files(iff)) THEN
[907]	1501	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[1146]	1502	CALL histwrite_phy(nid_files(iff),o_dulif%name,itau_w,zx_tmp_fi3d)
[907]	1503	ENDIF
	1504
[1146]	1505	IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN
[907]	1506	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[1146]	1507	CALL histwrite_phy(nid_files(iff),o_dvlif%name,itau_w,zx_tmp_fi3d)
[907]	1508	ENDIF
	1509	ENDIF
	1510
[1419]	1511	IF (ok_hines) THEN
	1512	IF (o_duhin%flag(iff)<=lev_files(iff)) THEN
	1513	zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys
	1514	CALL histwrite_phy(nid_files(iff),o_duhin%name,itau_w,zx_tmp_fi3d)
	1515	ENDIF
[907]	1516
[1419]	1517	IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN
	1518	zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys
	1519	CALL histwrite_phy(nid_files(iff),o_dvhin%name,itau_w,zx_tmp_fi3d)
	1520	ENDIF
	1521	ENDIF
	1522
[1405]	1523	IF (nqtot.GE.3 .AND. o_trac(1)%flag(iff)<=lev_files(iff)) THEN
	1524	CALL histwrite_phy(nid_files(iff),
	1525	s o_trac(1)%name,itau_w,qx(:,:,3))
	1526	ENDIF
	1527
	1528	IF (nqtot.GE.4 .AND. o_trac(2)%flag(iff)<=lev_files(iff)) THEN
	1529	CALL histwrite_phy(nid_files(iff),
	1530	s o_trac(2)%name,itau_w,qx(:,:,4))
	1531	ENDIF
	1532
	1533
	1534
[907]	1535	if (ok_sync) then
	1536	c$OMP MASTER
	1537	call histsync(nid_files(iff))
	1538	c$OMP END MASTER
	1539	endif
	1540
	1541	ENDIF ! clef_files
	1542
	1543	ENDDO ! iff=1,nfiles

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

Download in other formats: