Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

write_histmth.h @ 1332

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

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

Context Navigation

source: LMDZ4/tags/LMDZ4_V3_1/libf/phylmd/write_histmth.h @ 1332

Download in other formats: