Context Navigation

source: LMDZ4/trunk/libf/phylmd/write_histmth.h @ 895

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

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

Download in other formats: