Context Navigation

write_histmth.h @ 1076

Last change on this file since 1076 was 925, checked in by Laurent Fairhead, 16 years ago
Erreur lors de la reintegration des sorties aerosols SD LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 69.8 KB

Rev	Line
[644]	1	c
	2	c $Header$
	3	c
[524]	4	IF (ok_mensuel) THEN
	5	c
	6	ndex2d = 0
	7	ndex3d = 0
	8
	9	itau_w = itau_phy + itap
	10
	11	c
[644]	12	IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
	13	c
[524]	14	IF(lev_histmth.GE.1) THEN
	15	c
	16	c Champs 2D:
	17	c
[766]	18	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
	19	CALL histwrite_phy(nid_mth,"phis",itau_w,pphis)
[524]	20	C
[644]	21	cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
[766]	22	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
	23	CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
[524]	24	c
[644]	25	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
[766]	26	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	27	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
	28	$ zx_tmp_fi2d)
[524]	29	c
[644]	30	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
[766]	31	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	32	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
	33	$ zx_tmp_fi2d)
[644]	34	c
[766]	35	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
	36	CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
[524]	37
[766]	38	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
	39	CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
[524]	40	c
[766]	41	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
	42	CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
[524]	43	c
	44	c ENSEMBLES BEG
[766]	45	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	46	CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
[524]	47	c
[766]	48	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	49	CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
[524]	50	c
	51	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
	52	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
	53	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
[766]	54	c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
[524]	55	c . iim*jjmp1,ndex2d)
	56	c
	57	DO i=1, klon
	58	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	59	ENDDO
[766]	60	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	61	CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
[524]	62	c
	63	DO i = 1, klon
	64	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	65	ENDDO
[766]	66	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	67	CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
[524]	68	c
	69	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
[766]	70	c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
[524]	71	c . iimjjmp1klev,ndex3d)
	72	c
	73	c ENSEMBLES END
	74	c
[766]	75	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
	76	CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
[524]	77	c
[766]	78	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
	79	CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
[524]	80	C
[766]	81	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
	82	CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
[524]	83	C
	84	DO i = 1, klon
	85	zx_tmp_fi2d(i) = paprs(i,1)
	86	ENDDO
[766]	87	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	88	CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
[524]	89	c
[766]	90	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
	91	CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
[524]	92	c
	93	if (.not. ok_veget) then
[766]	94	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
	95	CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
[524]	96	endif
	97	c
	98	DO i = 1, klon
	99	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
	100	ENDDO
[766]	101	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	102	CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
[524]	103	c
[766]	104	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
	105	CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
[524]	106	c
	107	DO i = 1, klon
	108	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	109	ENDDO
[766]	110	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	111	CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
[524]	112	c
	113	DO i = 1, klon
	114	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	115	ENDDO
[766]	116	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	117	CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
[524]	118	c
[766]	119	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
	120	CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall)
[524]	121	c
[644]	122	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
	123	c CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
	124	c . ndex2d)
[524]	125	c
[766]	126	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
	127	CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
[524]	128	c
[766]	129	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
	130	CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
[524]	131	c
[766]	132	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
	133	CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
[524]	134	c
[766]	135	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
	136	CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
[524]	137	c
[766]	138	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
	139	CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
[524]	140	c
	141	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
[766]	142	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	143	CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
[524]	144	c
[766]	145	c
[524]	146	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
[766]	147	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	148	CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
[524]	149	c
	150	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
[766]	151	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	152	CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
[524]	153	c
	154	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
[766]	155	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	156	CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
[524]	157	c
[766]	158	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
	159	CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
[524]	160	c
[766]	161	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
	162	CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
[524]	163	c
[766]	164	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
	165	CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
[524]	166	c
[766]	167	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
	168	CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
[524]	169	c
[766]	170	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
	171	CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
[524]	172	c
[766]	173	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
	174	CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
[524]	175	c
	176	zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
[766]	177	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	178	CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
[524]	179	c
	180	zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
[766]	181	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	182	CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
[524]	183	c
	184	c200 END
[766]	185	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
	186	CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
[524]	187	c
[766]	188	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
	189	CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
[524]	190	c
[766]	191	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
	192	CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
[524]	193	c
[766]	194	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
	195	CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
[684]	196	c
[766]	197	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
	198	CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
[524]	199	c
	200	c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
	201	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
[766]	202	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	203	CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
[524]	204	c
	205	c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
	206	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
[766]	207	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	208	CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
[524]	209	c
	210	c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
	211	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
[766]	212	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	213	CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
[524]	214	c
	215	c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
	216	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
[766]	217	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	218	CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
[524]	219	c
[644]	220	cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
	221	cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
	222	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
[766]	223	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	224	CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
[524]	225	c
[644]	226	DO i=1, klon
	227	sollwdownclr(i)=-1.*lwdn0(i,1)
	228	ENDDO !i=1, klon
	229
	230	cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
	231	cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
	232	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
[766]	233	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	234	CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
[524]	235	c
[766]	236	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
	237	CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
[524]	238	c
[644]	239	c DO i=1, klon
	240	cIM 080904 sollwdownclr(i)=lwdn0(i,1)
	241	c sollwdownclr(i)=-1.*lwdn0(i,1)
	242	c ENDDO !i=1, klon
[524]	243	c
[766]	244	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
	245	CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
[524]	246	c
[896]	247	c
	248	c effets des aerosols
	249	c
	250	IF (ok_ade.OR.ok_aie) THEN
	251	zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
[925]	252	CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
[896]	253	c
	254	zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
[925]	255	CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
[896]	256	c
	257	zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
[925]	258	CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
[896]	259	c
	260	zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
[925]	261	CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
[896]	262	endif
	263
	264	c
[766]	265	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
	266	CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
[524]	267	c
	268	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
	269	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
[766]	270	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	271	CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
[524]	272	c
[766]	273	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
	274	CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
[524]	275	c
	276	c
	277	c DO i = 1, klon
	278	c zx_tmp_fi2d(i) = fluxu(i,1)
	279	c ENDDO
	280	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	281	c CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	282	c
	283	c DO i = 1, klon
	284	c zx_tmp_fi2d(i) = fluxv(i,1)
	285	c ENDDO
	286	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	287	c CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	288	c
	289	c
[766]	290	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
	291	CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
[524]	292	c
[766]	293	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
	294	CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
	295
	296	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqfonte,zx_tmp_2d)
	297	CALL histwrite_phy(nid_mth,"fqfonte",itau_w,zxfqfonte)
	298
[524]	299	DO nsrf = 1, nbsrf
	300	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
[766]	301	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	302	CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
	303	$ zx_tmp_fi2d)
[524]	304	C
	305	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
[766]	306	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	307	CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
	308	$ zx_tmp_fi2d)
[524]	309	ENDDO
	310	c
	311	DO nsrf = 1, nbsrf
	312	CYYY
[892]	313	IF(nsrf.GE.2) THEN
[644]	314	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
[766]	315	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	316	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w,
	317	& zx_tmp_fi2d)
[644]	318	c
	319	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
[766]	320	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	321	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w,
	322	$ zx_tmp_fi2d)
[897]	323	ENDIF !nsrf.GE.2
[524]	324	C
	325	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
[766]	326	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	327	CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
	328	$ zx_tmp_fi2d)
[524]	329	C
	330	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
[766]	331	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	332	CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
	333	$ zx_tmp_fi2d)
[524]	334	C
	335	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
[766]	336	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	337	CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
	338	$ zx_tmp_fi2d)
[524]	339	C
	340	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
[766]	341	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	342	CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
	343	$ zx_tmp_fi2d)
[524]	344	c
	345	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
[766]	346	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	347	CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
	348	$ zx_tmp_fi2d)
[524]	349	C
	350	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
[766]	351	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	352	CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
	353	$ zx_tmp_fi2d)
[524]	354	C
[684]	355	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
[766]	356	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	357	CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
	358	$ zx_tmp_fi2d)
[684]	359	C
[644]	360	IF(1.EQ.0) THEN
	361	C
	362	zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
[766]	363	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	364	CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
	365	$ zx_tmp_fi2d)
[644]	366	C
	367	zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
	368	$ *pctsrf(1:klon,nsrf)
[766]	369	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	370	CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
	371	$ zx_tmp_fi2d)
[644]	372	C
	373	ENDIF !(1.EQ.0) THEN
	374	C
[524]	375	END DO
	376	c
[766]	377	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
	378	CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
[524]	379	c
[766]	380	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
	381	CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
[524]	382	c
[766]	383	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
	384	CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
[524]	385	c
[766]	386	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
	387	CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
[524]	388	c
[766]	389	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
	390	CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
[524]	391	c
[766]	392	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
	393	CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt)
[524]	394	c
[766]	395	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
	396	CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
[524]	397	c
	398	zx_tmp_fi2d(1:klon) = flwp(1:klon)
[766]	399	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	400	CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
[524]	401	c
	402	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
[766]	403	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	404	CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
[524]	405	c
[766]	406	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
	407	CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
[524]	408	c
[766]	409	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
	410	CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
[524]	411	c
[766]	412	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
	413	CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
[524]	414	c
[766]	415	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
	416	CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
[524]	417	cKE43
	418	IF (iflag_con.GE.3) THEN ! sb
	419	c
[766]	420	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
	421	CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
[524]	422	c
[766]	423	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
	424	CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
[524]	425	c
[766]	426	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
	427	CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
[524]	428	c
[766]	429	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
	430	CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
[524]	431	c
	432	c34EK
	433	c
[766]	434	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
	435	CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
[524]	436	c
[644]	437	ENDIF !iflag_con.GE.3
[524]	438	c
	439	c
[644]	440	cIM cf. AM 081204 BEG
	441	cHBTM2
[524]	442	c
[644]	443	c
[766]	444	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
	445	CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
[644]	446	c
[766]	447	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
	448	CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
[644]	449	c
[766]	450	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
	451	CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
[644]	452	c
[766]	453	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
	454	CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
[644]	455	c
[766]	456	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
	457	CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
[644]	458	c
[766]	459	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
	460	CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
[644]	461	c
[766]	462	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
	463	CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
[644]	464	c
[766]	465	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
	466	CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
[644]	467	c
[766]	468	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
	469	CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
[644]	470	c
[766]	471	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
	472	CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
[644]	473	c
	474	cIM cf. AM 081204 END
	475	c
[524]	476	c Champs interpolles sur des niveaux de pression
[644]	477	DO k=1, nlevSTD
[524]	478	c
[644]	479	c bb=clevSTD(k)
	480	c
	481	c IF(k.GE.2) THEN
	482	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	483	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	484	c aa=clevSTD(k)
	485	c bb=aa(1:lnblnk1(aa))
	486	c ENDIF
	487	c
	488	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
	489	$ bb2.EQ."500".OR.bb2.EQ."200") THEN
	490	c
[766]	491	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
	492	CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
[524]	493	c
[766]	494	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
	495	CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
[524]	496	c
[766]	497	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
	498	CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
[524]	499	c
[766]	500	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
	501	CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
[524]	502	c
[644]	503	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
	504	ENDDO !nlevSTD
	505	c
	506	DO i=1, klon
	507	IF (pctsrf(i,is_oce).GT.epsfra) THEN
	508	zx_tmp_fi2d(i) = fluxo(i)
	509	ELSE
[684]	510	zx_tmp_fi2d(i) = 0.
[644]	511	ENDIF
	512	ENDDO
	513	c
[766]	514	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	515	CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
[644]	516	c
	517	DO i=1, klon
	518	IF (pctsrf(i,is_sic).GT.epsfra) THEN
	519	zx_tmp_fi2d(i) = fluxg(i)
	520	ELSE
[684]	521	zx_tmp_fi2d(i) = 0.
[644]	522	ENDIF
	523	ENDDO
	524	c
[766]	525	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	526	CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
[644]	527	c
[684]	528	DO i=1, klon
	529	IF (pctsrf(i,is_oce).GT.epsfra.OR.
	530	. pctsrf(i,is_sic).GT.epsfra) THEN
	531	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
	532	. ftsol(i, is_sic) * pctsrf(i,is_sic))/
	533	. (pctsrf(i,is_oce)+pctsrf(i,is_sic))
	534	ELSE
	535	zx_tmp_fi2d(i) = 273.15
	536	ENDIF
	537	ENDDO
	538	c
[766]	539	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	540	CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
[684]	541	c
[644]	542	IF (OCEAN.EQ.'force ') THEN
	543	c
	544	DO i=1, klon
	545	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	546	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	547	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	548	. fluxg(i)*pctsrf(i,is_sic)
	549	ELSE
	550	zx_tmp_fi2d(i) = 1.E+20
	551	ENDIF
	552	ENDDO
[766]	553	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	554	CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,zx_tmp_fi2d)
[644]	555	c
	556	ELSE IF (OCEAN.EQ.'slab ') THEN
	557	c
	558	DO i=1, klon
	559	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	560	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	561	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	562	. fluxg(i)*pctsrf(i,is_sic)
	563	ELSE
	564	zx_tmp_fi2d(i) = 1.E+20
	565	ENDIF
	566	ENDDO
[766]	567	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	568	CALL histwrite_phy(nid_mth,"slab_bils",itau_w,zx_tmp_fi2d)
[644]	569	c
	570	DO i=1, klon
	571	IF(pctsrf(i,is_oce).GT.epsfra.OR.
	572	$ pctsrf(i,is_sic).GT.epsfra) THEN
	573	zx_tmp_fi2d(i)=tslab(i)
	574	ELSE
	575	zx_tmp_fi2d(i) = 1.E+20
	576	ENDIF
	577	ENDDO !i=1, klon
[766]	578	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	579	CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
[644]	580	c
[766]	581	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
	582	CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
[644]	583	c
[766]	584	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
	585	CALL histwrite_phy(nid_mth,"siceh",itau_w, seaice/1000.)
[644]	586	c
	587	ENDIF !OCEAN.EQ.force/slab
	588	c
[524]	589	ENDIF !lev_histmth.GE.1
	590	c
	591	IF(lev_histmth.GE.2) THEN
	592	c
	593	c Champs 3D:
	594	C
	595	c
[766]	596	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
	597	CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
[524]	598	c
[766]	599	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
	600	CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
[524]	601	c
[766]	602	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
	603	CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
[524]	604	c
	605	c#ifdef histmthNMC
	606	c#include "write_histmthNMC.h"
	607	c#endif
	608	c
[766]	609	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
	610	CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap))
[524]	611	c
	612	c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
	613	c . (1-qx(1:klon,1:klev,ivap))
	614	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
	615	c CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d,
	616	c . iimjjmp1klev,ndex3d)
	617	c
[766]	618	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
	619	CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
[524]	620	c
[766]	621	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
	622	CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
[524]	623	c
[766]	624	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
	625	CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
[524]	626	c
[766]	627	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
	628	CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
[524]	629	c
[766]	630	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
	631	CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
[524]	632	c
[766]	633	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d)
	634	CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra)
[524]	635	c
[766]	636	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
	637	CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
[524]	638	c
[766]	639	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
	640	CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
	641
[644]	642	DO k=1, klev
	643	DO i=1, klon
	644	cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
	645	c $ /(paprs(i,k)-paprs(i,k+1))
	646	c $ *(paprs(i,1)/101325.0)
	647	c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)1.0e+6,1.0e-6)RG/46.6968
	648	zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
	649	$ /(paprs(i,k)-paprs(i,k+1))
	650	$ *(paprs(i,1)/101325.0)
	651	ENDDO !i
	652	ENDDO !k
	653	c
	654	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
	655	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
[766]	656	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	657	CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
[524]	658	c
[766]	659	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
	660	CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
[524]	661	c
[766]	662	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),
	663	cym . zx_tmp_3d)
	664	CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap))
	665	cIM 22.03.04 END
[524]	666	c
	667	ENDIF !lev_histmth.GE.2
	668	c
	669	IF(lev_histmth.GE.3) THEN
	670	c
	671	DO nsrf=1, nbsrf
	672	c
	673	zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
[766]	674	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	675	CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
	676	$ zx_tmp_fi2d)
[524]	677	C
	678	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
[766]	679	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	680	CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
	681	$ zx_tmp_fi2d)
[524]	682	c
	683	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
[766]	684	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
	685	CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
	686	$ ,zx_tmp_fi2d)
[524]	687	c
	688	ENDDO !nsrf=1, nbsrf
	689	c
[766]	690	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
	691	CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
[524]	692	c
[766]	693	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
	694	CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
[524]	695	c
	696	ENDIF !lev_histmth.GE.3
	697	c
	698	c
	699	IF(lev_histmth.GE.4) THEN
	700	c
[766]	701	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
	702	CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
	703
[653]	704	IF (iflag_con.GE.3) THEN
	705	c
[766]	706	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
	707	CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
[524]	708	c
[766]	709	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
	710	CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
[684]	711	c
[766]	712	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
	713	CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
[524]	714	c
[766]	715	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
	716	CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
[524]	717	c
[653]	718	ENDIF !iflag_con.GE.3
[524]	719	c
[766]	720	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
	721	CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
[524]	722	c
[766]	723	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
	724	CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
[524]	725	c
	726	cIM: 101003 : K/30min ==> K/s
	727	zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
[766]	728	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	729	CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
[524]	730	c
[644]	731	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
[766]	732	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	733	CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
[524]	734	c
	735	cIM: 101003 : K/30min ==> K/s
	736	zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
[766]	737	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	738	CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
	739	cIM: 071003
[524]	740	cIM: 101003 : K/30min ==> K/s
	741	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
	742	. d_t_eva(1:klon,1:klev))/pdtphys
[766]	743	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
	744	CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
[524]	745	c
[644]	746	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
[766]	747	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	748	CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
[524]	749	c
	750	cIM: 101003 : K/30min ==> K/s
	751	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
[766]	752	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	753	CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
[524]	754	c
[644]	755	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
[766]	756	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	757	CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
[524]	758	c
	759	cIM: 101003 : K/30min ==> K/s
	760	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
[766]	761	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	762	CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
[524]	763	c
[644]	764	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
[766]	765	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	766	CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
[524]	767	c
	768	zpt_conv = 0.
	769	where (ptconv) zpt_conv = 1.
[766]	770	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
	771	CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
[524]	772	c
[766]	773	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
	774	CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
[524]	775	c
	776	cIM: 101003 : K/30min ==> K/s
	777	zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys
[766]	778	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	779	CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
[524]	780	c
[644]	781	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
[766]	782	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	783	CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
[524]	784	c
	785	cIM: 101003 : K/day ==> K/s
	786	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
[766]	787	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	788	CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
[524]	789	c
	790	cIM: 101003 : K/day ==> K/s
	791	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
[766]	792	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	793	CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
[524]	794	c
	795	cIM: 101003 : K/day ==> K/s
	796	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
[766]	797	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	798	CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
[524]	799	c
	800	cIM: 101003 : K/day ==> K/s
	801	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
[766]	802	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	803	CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
[524]	804	c
	805	cIM: 101003 : deja en K/s
	806	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
[766]	807	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	808	CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
[524]	809	c
[644]	810	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[766]	811	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	812	CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
[524]	813	c
[644]	814	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[766]	815	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	816	CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
[524]	817	c
	818	IF (ok_orodr) THEN
[644]	819	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[766]	820	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	821	CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
[524]	822	c
[644]	823	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[766]	824	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	825	CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
[524]	826	c
	827	ENDIF
	828	C
	829	IF (ok_orolf) THEN
[644]	830	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[766]	831	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	832	CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
[524]	833	c
[644]	834	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[766]	835	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	836	CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
[524]	837	ENDIF
	838	C
[862]	839
	840	#if defined(INCA_AER) && defined(CPP_COUPLE)
	841	CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
	842	CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
	843	CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
	844	CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
	845	CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
	846	CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
	847	CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
	848	CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
	849	CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
	850	CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
	851	CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
	852	CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
	853	CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
	854	CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
	855	#endif
	856	#ifndef INCA
[644]	857	IF (nqmax.GE.3) THEN
[653]	858	DO iq=3,nqmax
[766]	859	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
	860	CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
[644]	861	ENDDO
	862	ENDIF
[862]	863	#endif
[644]	864	ENDIF !lev_histmth.GE.4
[524]	865	c
[644]	866	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
	867	c
	868	IF(lev_histmth.GE.1) THEN
	869	c
	870	c Champs 2D:
	871	c
[766]	872	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
	873	CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG)
[644]	874	C
	875	cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
[766]	876	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
	877	CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
[644]	878	c
	879	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
[766]	880	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	881	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
	882	$ zx_tmp_fi2d)
[644]	883	c
	884	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
[766]	885	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	886	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
	887	$ zx_tmp_fi2d)
[644]	888	c
	889	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100.
[766]	890	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	891	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
	892	$ zx_tmp_fi2d)
[644]	893	c
	894	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)
[766]	895	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	896	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
	897	$ zx_tmp_fi2d)
[644]	898	c
	899	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100.
[766]	900	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	901	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
	902	$ zx_tmp_fi2d)
[644]	903	c
	904	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)
[766]	905	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	906	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
	907	$ zx_tmp_fi2d)
[644]	908	c
[766]	909	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
	910	CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat)
[644]	911	c
	912	zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon)
[766]	913	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	914	CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d)
[644]	915	c
	916	c IF(itap.EQ.ecrit_mth) THEN
	917	c DO i=1, klon
	918	c IF(nday_clr(i).GT.0.) THEN
	919	c tops01(i)=tops01(i)/nday_clr(i)
	920	c ELSE
	921	c tops01(i)=1.e+20
	922	c ENDIF
	923	c ENDDO
	924	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d)
[766]	925	c CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d,
[644]	926	c . iim*jjmp1,ndex2d)
	927	c ENDIF !itap.EQ.ecrit_mth
	928	c
	929	c220404 BEG
	930	c
	931	c temperature tendency due to total diabatic heating
	932	DO l=1, klev
	933	DO i=1, klon
	934	d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime
	935	ENDDO !i
	936	ENDDO !l
	937	c
[766]	938	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
	939	CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t)
[644]	940	c
	941	c temperature tendency due to SW radiation
	942	DO l=1, klev
	943	DO i=1, klon
	944	zx_tmp_fi3d(i,l)=heat(i,l)/RDAY
	945	ENDDO !i
	946	ENDDO !l
	947	c
[766]	948	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	949	CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d)
[644]	950	c
	951	c temperature tendency due to LW radiation
	952	DO l=1, klev
	953	DO i=1, klon
	954	zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY
	955	ENDDO !i
	956	ENDDO !l
	957	c
[766]	958	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	959	CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
[644]	960	c
	961	c temperature tendency due to moist convective processes
	962	DO l=1, klev
	963	DO i=1, klon
	964	zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
	965	ENDDO !i
	966	ENDDO !l
	967	c
[766]	968	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	969	CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
[644]	970	c
	971	c temperature tendency due to dry convective processes
	972	DO l=1, klev
	973	DO i=1, klon
	974	zx_tmp_fi3d(i,l)=d_t_ajs(i,l)/pdtphys
	975	ENDDO !i
	976	ENDDO !l
	977	c
[766]	978	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	979	CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
[644]	980	c
	981	c temperature tendency due to large scale precipitation
	982	DO l=1, klev
	983	DO i=1, klon
	984	zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys
	985	ENDDO !i
	986	ENDDO !l
	987	c
[766]	988	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	989	CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
[644]	990	c
	991	c total moisture tendency due to diabatic processes
	992	DO l=1, klev
	993	DO i=1, klon
	994	d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime
	995	zx_tmp_fi3d(i,l)=d_qx(i,l,ivap)
	996	ENDDO !i
	997	ENDDO !l
	998	c
[766]	999	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	1000	CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d)
[644]	1001	c
	1002	cBEG 210404 BEG
	1003	c
	1004	IF(ok_orodr.AND.ok_orolf) THEN
	1005	DO k = 1, klev
	1006	DO i = 1, klon
	1007	d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys
	1008	d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys
	1009	ENDDO
	1010	ENDDO
	1011	ENDIF !(ok_orodr.AND.ok_orolf) THEN
	1012	c
[766]	1013	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
	1014	CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli)
[644]	1015	c
[766]	1016	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
	1017	CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli)
[644]	1018	c
	1019	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
[766]	1020	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1021	CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d)
[644]	1022	c
	1023	zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
[766]	1024	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1025	CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d)
[644]	1026	c
[766]	1027	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
	1028	CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u)
[644]	1029	c
[766]	1030	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
	1031	CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v)
[644]	1032
	1033	cEND 210404 END
	1034	c
	1035	c cldtrue
	1036	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d)
[766]	1037	c CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
[644]	1038	c $ iimjjmp1klev,ndex3d)
	1039	c
	1040	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d)
[766]	1041	c CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
[644]	1042	c $ iimjjmp1klev,ndex3d)
	1043	c
	1044	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d)
[766]	1045	c CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
[644]	1046	c $ iimjjmp1klev,ndex3d)
	1047	c
	1048	DO i=1, klon
	1049	cIM 151004 zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1)
	1050	zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i)
	1051	ENDDO
[766]	1052	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1053	CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d)
[644]	1054	c
	1055	DO i=1, klon
	1056	cIM 151004 zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1)
	1057	zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i)
	1058	ENDDO
[766]	1059	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1060	CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d)
[644]	1061	c
	1062	c220404 END
	1063	cAMIP2 END
	1064	c
[766]	1065	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
	1066	CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
[644]	1067
[766]	1068	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
	1069	CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
[644]	1070	c
[766]	1071	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
	1072	CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
[644]	1073	c
	1074	c ENSEMBLES BEG
[766]	1075	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	1076	CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
[644]	1077	c
[766]	1078	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	1079	CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
[644]	1080	c
	1081	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
	1082	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
	1083	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
[766]	1084	c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
[644]	1085	c . iim*jjmp1,ndex2d)
	1086	c
	1087	DO i=1, klon
	1088	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	1089	ENDDO
[766]	1090	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1091	CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
[644]	1092	c
	1093	DO i = 1, klon
	1094	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	1095	ENDDO
[766]	1096	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1097	CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
[644]	1098	c
	1099	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
[766]	1100	c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
[644]	1101	c . iimjjmp1klev,ndex3d)
	1102	c
	1103	c ENSEMBLES END
	1104	c
[766]	1105	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
	1106	CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
[644]	1107	c
[766]	1108	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
	1109	CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
[644]	1110	C
[766]	1111	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
	1112	CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
[644]	1113	C
	1114	DO i = 1, klon
	1115	zx_tmp_fi2d(i) = paprs(i,1)
	1116	ENDDO
[766]	1117	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1118	CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
[644]	1119	c
[766]	1120	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
	1121	CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
[644]	1122	c
	1123	if (.not. ok_veget) then
[766]	1124	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
	1125	CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
[644]	1126	endif
	1127	c
	1128	DO i = 1, klon
	1129	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
	1130	ENDDO
[766]	1131	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1132	CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
[644]	1133	c
[766]	1134	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
	1135	CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
[644]	1136	c
	1137	DO i = 1, klon
	1138	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	1139	ENDDO
[766]	1140	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1141	CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
[644]	1142	c
	1143	DO i = 1, klon
	1144	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	1145	ENDDO
[766]	1146	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1147	CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
[644]	1148	c
[766]	1149	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
	1150	CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall)
[644]	1151	c
	1152	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
[766]	1153	c CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
[644]	1154	c . ndex2d)
	1155	c
[766]	1156	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
	1157	CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
[644]	1158	c
[766]	1159	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
	1160	CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
[644]	1161	c
[766]	1162	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
	1163	CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
[644]	1164	c
[766]	1165	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
	1166	CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
[644]	1167	c
[766]	1168	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
	1169	CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
[644]	1170	c
	1171	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
[766]	1172	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1173	CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
[644]	1174	c
	1175	c
	1176	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
[766]	1177	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1178	CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
[644]	1179	c
	1180	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
[766]	1181	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1182	CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
[644]	1183	c
	1184	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
[766]	1185	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1186	CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
[644]	1187	c
[766]	1188	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
	1189	CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
[644]	1190	c
[766]	1191	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
	1192	CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
[644]	1193	c
[766]	1194	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
	1195	CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
[644]	1196	c
[766]	1197	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
	1198	CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
[644]	1199	c
[766]	1200	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
	1201	CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
[644]	1202	c
[766]	1203	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
	1204	CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
[644]	1205	c
	1206	zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
[766]	1207	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1208	CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
[644]	1209	c
	1210	zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
[766]	1211	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1212	CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
[644]	1213	c
	1214	c200 END
[766]	1215	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
	1216	CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
[644]	1217	c
[766]	1218	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
	1219	CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
[644]	1220	c
[766]	1221	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
	1222	CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
[644]	1223	c
[766]	1224	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
	1225	CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
[684]	1226	c
[766]	1227	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
	1228	CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
[644]	1229	c
	1230	c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
	1231	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
[766]	1232	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1233	CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
[644]	1234	c
	1235	c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
	1236	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
[766]	1237	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1238	CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
[644]	1239	c
	1240	c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
	1241	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
[766]	1242	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1243	CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
[644]	1244	c
	1245	c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
	1246	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
[766]	1247	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1248	CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
[644]	1249	c
	1250	cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
	1251	cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
	1252	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
[766]	1253	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	1254	CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
[644]	1255	c
	1256	DO i=1, klon
	1257	sollwdownclr(i)=-1.*lwdn0(i,1)
	1258	ENDDO !i=1, klon
	1259	c
	1260	cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
	1261	cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
	1262	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
[766]	1263	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	1264	CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
[644]	1265	c
[766]	1266	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
	1267	CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
[644]	1268	c
	1269	c DO i=1, klon
	1270	cIM 08004 sollwdownclr(i)=lwdn0(i,1)
	1271	c sollwdownclr(i)=-1.*lwdn0(i,1)
	1272	c ENDDO !i=1, klon
	1273	c
[766]	1274	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
	1275	CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
[644]	1276	c
[766]	1277	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
	1278	CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
[644]	1279	c
	1280	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
	1281	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
[766]	1282	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1283	CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
[644]	1284	c
[766]	1285	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
	1286	CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
[644]	1287	c
[897]	1288	c Effets des aerosols
[644]	1289	c
[897]	1290	IF (ok_ade.OR.ok_aie) THEN
	1291	zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
	1292	CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
	1293	c
	1294	zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
	1295	CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
	1296	c
	1297	zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
	1298	CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
	1299	c
	1300	zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
	1301	CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
	1302	ENDIF
	1303	c
	1304	c
[644]	1305	c DO i = 1, klon
	1306	c zx_tmp_fi2d(i) = fluxu(i,1)
	1307	c ENDDO
	1308	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
[766]	1309	c CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[644]	1310	c
	1311	c DO i = 1, klon
	1312	c zx_tmp_fi2d(i) = fluxv(i,1)
	1313	c ENDDO
	1314	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
[766]	1315	c CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[644]	1316	c
	1317	c
[766]	1318	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
	1319	CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
[644]	1320	c
[766]	1321	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
	1322	CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
[644]	1323	c
	1324	DO nsrf = 1, nbsrf
	1325	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
[766]	1326	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1327	CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
	1328	$ zx_tmp_fi2d)
[644]	1329	C
	1330	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
[766]	1331	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1332	CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
	1333	$ zx_tmp_fi2d)
[644]	1334	ENDDO
	1335	c
	1336	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
[766]	1337	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1338	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
	1339	$ zx_tmp_fi2d)
[644]	1340	c
	1341	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
[766]	1342	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1343	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
	1344	$ zx_tmp_fi2d)
[644]	1345	c
	1346	DO nsrf = 1, nbsrf
	1347	CYYY
	1348	C
	1349	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
[766]	1350	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1351	CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
	1352	$ zx_tmp_fi2d)
[644]	1353	C
	1354	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
[766]	1355	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1356	CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
	1357	$ zx_tmp_fi2d)
[644]	1358	C
	1359	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
[766]	1360	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1361	CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
	1362	$ zx_tmp_fi2d)
[644]	1363	C
	1364	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
[766]	1365	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1366	CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
	1367	$ zx_tmp_fi2d)
[644]	1368	c
	1369	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
[766]	1370	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1371	CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
	1372	$ zx_tmp_fi2d)
[644]	1373	C
	1374	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
[766]	1375	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1376	CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
	1377	$ zx_tmp_fi2d)
[644]	1378	C
[684]	1379	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
[766]	1380	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1381	CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
	1382	$ zx_tmp_fi2d)
[684]	1383	C
[644]	1384	IF(1.EQ.0) THEN
	1385	C
	1386	zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
[766]	1387	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1388	CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
	1389	$ zx_tmp_fi2d)
[644]	1390	C
	1391	zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
	1392	$ *pctsrf(1:klon,nsrf)
[766]	1393	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1394	CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
	1395	$ zx_tmp_fi2d)
[644]	1396	C
	1397	ENDIF !(1.EQ.0) THEN
	1398	C
	1399	END DO
	1400	c
[766]	1401	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
	1402	CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
[644]	1403	c
[766]	1404	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
	1405	CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
[644]	1406	c
[766]	1407	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
	1408	CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
[644]	1409	c
[766]	1410	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
	1411	CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
[644]	1412	c
[766]	1413	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
	1414	CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
[644]	1415	c
[766]	1416	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
	1417	CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
[644]	1418	c
[766]	1419	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
	1420	CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
[644]	1421	c
	1422	zx_tmp_fi2d(1:klon) = flwp(1:klon)
[766]	1423	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1424	CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
[644]	1425	c
	1426	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
[766]	1427	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1428	CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
[644]	1429	c
[766]	1430	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
	1431	CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
[644]	1432	c
[766]	1433	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
	1434	CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
[644]	1435	c
[766]	1436	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
	1437	CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
[644]	1438	c
[766]	1439	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
	1440	CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
[644]	1441	cKE43
	1442	IF (iflag_con.GE.3) THEN ! sb
	1443	c
[766]	1444	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
	1445	CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
[644]	1446	c
[766]	1447	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
	1448	CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
[644]	1449	c
[766]	1450	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
	1451	CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
[644]	1452	c
[766]	1453	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
	1454	CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
[644]	1455	c
	1456	c34EK
	1457	c
[766]	1458	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
	1459	CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
[644]	1460	c
	1461	ENDIF !iflag_con.GE.3
	1462	c
	1463	cIM cf. AM 081204 BEG
	1464	cHBTM2
	1465	c
	1466	c
[766]	1467	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
	1468	CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
[644]	1469	c
[766]	1470	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
	1471	CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
[644]	1472	c
[766]	1473	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
	1474	CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
[644]	1475	c
[766]	1476	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
	1477	CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
[644]	1478	c
[766]	1479	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
	1480	CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
[644]	1481	c
[766]	1482	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
	1483	CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
[644]	1484	c
[766]	1485	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
	1486	CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
[644]	1487	c
[766]	1488	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
	1489	CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
[644]	1490	c
[766]	1491	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
	1492	CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
[644]	1493	c
[766]	1494	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
	1495	CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
[644]	1496	c
	1497	cIM cf. AM 081204 END
	1498	c
	1499	c Champs interpolles sur des niveaux de pression
	1500	c DO k=1, nlevENS
	1501	DO k=1, nlevSTD
	1502	c
	1503	c bb=clevSTD(k)
	1504	c
	1505	c IF(k.GE.2) THEN
	1506	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	1507	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	1508	c aa=clevSTD(k)
	1509	c bb=aa(1:lnblnk1(aa))
	1510	c ENDIF
	1511	c
	1512	IF(bb2.EQ."850".OR.bb2.EQ."700"
	1513	$ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	1514	c
[766]	1515	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
	1516	CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
[644]	1517	c
[766]	1518	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
	1519	CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
[644]	1520	c
[766]	1521	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
	1522	CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
[644]	1523	c
[766]	1524	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
	1525	CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
[644]	1526	c
	1527	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
	1528	ENDDO !nlevSTD
	1529	c
	1530	DO i=1, klon
	1531	IF (pctsrf(i,is_oce).GT.epsfra) THEN
	1532	zx_tmp_fi2d(i) = fluxo(i)
	1533	ELSE
[684]	1534	zx_tmp_fi2d(i) = 0.
[644]	1535	ENDIF
	1536	ENDDO
	1537	c
[766]	1538	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1539	CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
[644]	1540	c
	1541	DO i=1, klon
	1542	IF (pctsrf(i,is_sic).GT.epsfra) THEN
	1543	zx_tmp_fi2d(i) = fluxg(i)
	1544	ELSE
[684]	1545	zx_tmp_fi2d(i) = 0.
[644]	1546	ENDIF
	1547	ENDDO
	1548	c
[766]	1549	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1550	CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
[644]	1551	c
[684]	1552	DO i=1, klon
	1553	IF (pctsrf(i,is_oce).GT.epsfra.OR.
	1554	. pctsrf(i,is_sic).GT.epsfra) THEN
	1555	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
	1556	. ftsol(i, is_sic) * pctsrf(i,is_sic))/
	1557	. (pctsrf(i,is_oce)+pctsrf(i,is_sic))
	1558	ELSE
	1559	zx_tmp_fi2d(i) = 273.15
	1560	ENDIF
	1561	ENDDO
	1562	c
[766]	1563	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1564	CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
[684]	1565	c
[644]	1566	IF (OCEAN.EQ.'force ') THEN
	1567	c
	1568	DO i=1, klon
	1569	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	1570	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	1571	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	1572	. fluxg(i)*pctsrf(i,is_sic)
	1573	ELSE
	1574	zx_tmp_fi2d(i) = 1.E+20
	1575	ENDIF
	1576	ENDDO
[766]	1577	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1578	CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
	1579	. zx_tmp_fi2d)
[644]	1580	c
	1581	ELSE IF (OCEAN.EQ.'slab ') THEN
	1582	c
	1583	c
	1584	DO i=1, klon
	1585	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	1586	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	1587	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	1588	. fluxg(i)*pctsrf(i,is_sic)
	1589	ELSE
	1590	zx_tmp_fi2d(i) = 1.E+20
	1591	ENDIF
	1592	ENDDO
[766]	1593	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1594	CALL histwrite_phy(nid_mth,"slab_bils",itau_w,
	1595	. zx_tmp_fi2d)
[644]	1596	c
	1597	DO i=1, klon
	1598	IF(pctsrf(i,is_oce).GT.epsfra.OR.
	1599	$ pctsrf(i,is_sic).GT.epsfra) THEN
	1600	zx_tmp_fi2d(i)=tslab(i)
	1601	ELSE
	1602	zx_tmp_fi2d(i) = 1.E+20
	1603	ENDIF
	1604	ENDDO !i=1, klon
[766]	1605	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1606	CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
[644]	1607	c
[766]	1608	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
	1609	CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
[644]	1610	c
[766]	1611	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
	1612	CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
[644]	1613	c
	1614	ENDIF !OCEAN.EQ.force/slab
	1615	c
	1616	ENDIF !lev_histmth.GE.1
	1617	c
	1618	IF(lev_histmth.GE.2) THEN
	1619	c
	1620	c Champs 3D:
	1621	C
	1622	c
[766]	1623	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
	1624	CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
[644]	1625	c
[766]	1626	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
	1627	CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
[644]	1628	c
[766]	1629	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
	1630	CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
[524]	1631	c
[644]	1632	c#ifdef histmthNMC
	1633	c#include "write_histmthNMC.h"
	1634	c#endif
	1635	c
[766]	1636	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
	1637	CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
[524]	1638	c
[644]	1639	c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
	1640	c . (1-qx(1:klon,1:klev,ivap))
	1641	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
[766]	1642	c CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
[644]	1643	c . iimjjmp1klev,ndex3d)
	1644	c
[766]	1645	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
	1646	CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
[524]	1647	c
[766]	1648	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
	1649	CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
[524]	1650	c
[766]	1651	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
	1652	CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
[644]	1653	c
[766]	1654	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
	1655	CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
[644]	1656	c
[766]	1657	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
	1658	CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
[644]	1659	c
[766]	1660	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
	1661	CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
[644]	1662	c
[766]	1663	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
	1664	CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
[644]	1665	c
[766]	1666	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
	1667	CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
[644]	1668	c
	1669	DO k=1, klev
	1670	DO i=1, klon
	1671	cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
	1672	c $ /(paprs(i,k)-paprs(i,k+1))
	1673	c $ *(paprs(i,1)/101325.0)
	1674	c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)1.0e+6,1.0e-6)RG/46.6968
	1675	zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
	1676	$ /(paprs(i,k)-paprs(i,k+1))
	1677	$ *(paprs(i,1)/101325.0)
	1678	ENDDO !i
	1679	ENDDO !k
	1680	c
	1681	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
	1682	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
[766]	1683	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1684	CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
[644]	1685	c
[766]	1686	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
	1687	CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
[644]	1688	c
[766]	1689	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
	1690	CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
[644]	1691	c
	1692	ENDIF !lev_histmth.GE.2
	1693	c
	1694	IF(lev_histmth.GE.3) THEN
	1695	c
	1696	DO nsrf=1, nbsrf
	1697	c
	1698	zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
[766]	1699	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1700	CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
	1701	$ zx_tmp_fi2d)
[644]	1702	C
	1703	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
[766]	1704	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1705	CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
	1706	$ zx_tmp_fi2d)
[644]	1707	c
	1708	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
[766]	1709	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
	1710	CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
	1711	$ ,agesno)
[644]	1712	c
	1713	ENDDO !nsrf=1, nbsrf
	1714	c
[766]	1715	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
	1716	CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
[644]	1717	c
[766]	1718	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
	1719	CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
[644]	1720	c
	1721	ENDIF !lev_histmth.GE.3
	1722	c
	1723	c
	1724	IF(lev_histmth.GE.4) THEN
	1725	c
[766]	1726	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
	1727	CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
[644]	1728	c
[653]	1729	IF (iflag_con.GE.3) THEN
	1730	c
[766]	1731	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
	1732	CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
[644]	1733	c
[766]	1734	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
	1735	CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
[644]	1736	c
[766]	1737	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
	1738	CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
[644]	1739	c
[766]	1740	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
	1741	CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
[644]	1742	c
[653]	1743	ENDIF !iflag_con.GE.3
[644]	1744	c
[766]	1745	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
	1746	CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
[644]	1747	c
[766]	1748	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
	1749	CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
[644]	1750	c
	1751	cIM: 101003 : K/30min ==> K/s
	1752	c
	1753	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
[766]	1754	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1755	CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
[644]	1756	c
	1757	cIM: 101003 : K/30min ==> K/s
	1758	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
	1759	. d_t_eva(1:klon,1:klev))/pdtphys
[766]	1760	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
	1761	CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
[644]	1762	c
	1763	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
[766]	1764	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1765	CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
[644]	1766	c
	1767	cIM: 101003 : K/30min ==> K/s
	1768	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
[766]	1769	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1770	CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
[644]	1771	c
	1772	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
[766]	1773	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1774	CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
[644]	1775	c
	1776	cIM: 101003 : K/30min ==> K/s
	1777	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
[766]	1778	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1779	CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
[644]	1780	c
	1781	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
[766]	1782	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1783	CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
[644]	1784	c
	1785	zpt_conv = 0.
	1786	where (ptconv) zpt_conv = 1.
[766]	1787	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
	1788	CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
[644]	1789	c
[766]	1790	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
	1791	CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
[644]	1792	c
	1793	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
[766]	1794	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1795	CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
[644]	1796	c
	1797	cIM: 101003 : K/day ==> K/s
	1798	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
[766]	1799	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1800	CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
[644]	1801	c
	1802	cIM: 101003 : K/day ==> K/s
	1803	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
[766]	1804	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1805	CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
[644]	1806	c
	1807	cIM: 101003 : K/day ==> K/s
	1808	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
[766]	1809	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1810	CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
[644]	1811	c
	1812	cIM: 101003 : K/day ==> K/s
	1813	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
[766]	1814	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1815	CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
[644]	1816	c
	1817	cIM: 101003 : deja en K/s
	1818	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
[766]	1819	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1820	CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
[644]	1821	c
	1822	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[766]	1823	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1824	CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
[644]	1825	c
	1826	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[766]	1827	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1828	CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
[644]	1829	c
	1830	IF (ok_orodr) THEN
	1831	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[766]	1832	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1833	CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
[644]	1834	c
	1835	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[766]	1836	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1837	CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
[644]	1838	c
	1839	ENDIF
	1840	C
	1841	IF (ok_orolf) THEN
	1842	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[766]	1843	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1844	CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
[644]	1845	c
	1846	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[766]	1847	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1848	CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
[644]	1849	ENDIF
	1850	C
[862]	1851	#if defined(INCA_AER) && defined(CPP_COUPLE)
	1852	CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau)
	1853	CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1))
	1854	CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2))
	1855	CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi)
	1856	CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi)
	1857	CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1))
	1858	CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2))
	1859	CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3))
	1860	CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:))
	1861	CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:))
	1862	CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1))
	1863	CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2))
	1864	CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:))
	1865	CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3))
	1866	#endif
	1867	#ifndef INCA
[524]	1868	IF (nqmax.GE.3) THEN
[653]	1869	DO iq=3,nqmax
[766]	1870	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
	1871	CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
[524]	1872	ENDDO
	1873	ENDIF
	1874	c
[862]	1875	#endif
[524]	1876	ENDIF !lev_histmth.GE.4
	1877	c
[644]	1878	ENDIF !(type_run.EQ."AMIP")
	1879	c
[524]	1880	if (ok_sync) then
[766]	1881	c$OMP MASTER
[524]	1882	call histsync(nid_mth)
[766]	1883	c$OMP END MASTER
[524]	1884	endif
[644]	1885	c
	1886	ENDIF !ok_mensuel

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

Context Navigation

source: LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h @ 1076

Download in other formats: