Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

write_histmth.h

Last change on this file was 735, checked in by Laurent Fairhead, 18 years ago
Modifications version parallele YM/LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 66.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
[704]	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)
[704]	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.
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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
[704]	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
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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 )
[704]	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
[704]	145	c
[524]	146	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
[704]	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 )
[704]	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 )
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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 )
[704]	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 )
[704]	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 )
[704]	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 )
[704]	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)
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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)
[704]	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.
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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
[704]	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
[704]	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)
[704]	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)
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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.
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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)
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[704]	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
[644]	821	IF (nqmax.GE.3) THEN
[653]	822	DO iq=3,nqmax
[704]	823	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
	824	CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
[644]	825	ENDDO
	826	ENDIF
[524]	827	c
[644]	828	ENDIF !lev_histmth.GE.4
[524]	829	c
[644]	830	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
	831	c
	832	IF(lev_histmth.GE.1) THEN
	833	c
	834	c Champs 2D:
	835	c
[704]	836	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d)
	837	CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG)
[644]	838	C
	839	cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
[704]	840	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
	841	CALL histwrite_phy(nid_mth,"aire",itau_w,airephy)
[644]	842	c
	843	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100.
[704]	844	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	845	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w,
	846	$ zx_tmp_fi2d)
[644]	847	c
	848	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)
[704]	849	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	850	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w,
	851	$ zx_tmp_fi2d)
[644]	852	c
	853	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100.
[704]	854	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	855	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w,
	856	$ zx_tmp_fi2d)
[644]	857	c
	858	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)
[704]	859	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d)
	860	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w,
	861	$ zx_tmp_fi2d)
[644]	862	c
	863	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100.
[704]	864	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	865	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w,
	866	$ zx_tmp_fi2d)
[644]	867	c
	868	zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)
[704]	869	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	870	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w,
	871	$ zx_tmp_fi2d)
[644]	872	c
[704]	873	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
	874	CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat)
[644]	875	c
	876	zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon)
[704]	877	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	878	CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d)
[644]	879	c
	880	c IF(itap.EQ.ecrit_mth) THEN
	881	c DO i=1, klon
	882	c IF(nday_clr(i).GT.0.) THEN
	883	c tops01(i)=tops01(i)/nday_clr(i)
	884	c ELSE
	885	c tops01(i)=1.e+20
	886	c ENDIF
	887	c ENDDO
	888	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d)
[704]	889	c CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d,
[644]	890	c . iim*jjmp1,ndex2d)
	891	c ENDIF !itap.EQ.ecrit_mth
	892	c
	893	c220404 BEG
	894	c
	895	c temperature tendency due to total diabatic heating
	896	DO l=1, klev
	897	DO i=1, klon
	898	d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime
	899	ENDDO !i
	900	ENDDO !l
	901	c
[704]	902	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d)
	903	CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t)
[644]	904	c
	905	c temperature tendency due to SW radiation
	906	DO l=1, klev
	907	DO i=1, klon
	908	zx_tmp_fi3d(i,l)=heat(i,l)/RDAY
	909	ENDDO !i
	910	ENDDO !l
	911	c
[704]	912	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	913	CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d)
[644]	914	c
	915	c temperature tendency due to LW radiation
	916	DO l=1, klev
	917	DO i=1, klon
	918	zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY
	919	ENDDO !i
	920	ENDDO !l
	921	c
[704]	922	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	923	CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d)
[644]	924	c
	925	c temperature tendency due to moist convective processes
	926	DO l=1, klev
	927	DO i=1, klon
	928	zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys
	929	ENDDO !i
	930	ENDDO !l
	931	c
[704]	932	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	933	CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d)
[644]	934	c
	935	c temperature tendency due to dry convective processes
	936	DO l=1, klev
	937	DO i=1, klon
	938	zx_tmp_fi3d(i,l)=d_t_ajs(i,l)/pdtphys
	939	ENDDO !i
	940	ENDDO !l
	941	c
[704]	942	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	943	CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d)
[644]	944	c
	945	c temperature tendency due to large scale precipitation
	946	DO l=1, klev
	947	DO i=1, klon
	948	zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys
	949	ENDDO !i
	950	ENDDO !l
	951	c
[704]	952	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	953	CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d)
[644]	954	c
	955	c total moisture tendency due to diabatic processes
	956	DO l=1, klev
	957	DO i=1, klon
	958	d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime
	959	zx_tmp_fi3d(i,l)=d_qx(i,l,ivap)
	960	ENDDO !i
	961	ENDDO !l
	962	c
[704]	963	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d)
	964	CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d)
[644]	965	c
	966	cBEG 210404 BEG
	967	c
	968	IF(ok_orodr.AND.ok_orolf) THEN
	969	DO k = 1, klev
	970	DO i = 1, klon
	971	d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys
	972	d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys
	973	ENDDO
	974	ENDDO
	975	ENDIF !(ok_orodr.AND.ok_orolf) THEN
	976	c
[704]	977	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d)
	978	CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli)
[644]	979	c
[704]	980	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d)
	981	CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli)
[644]	982	c
	983	zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
[704]	984	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	985	CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d)
[644]	986	c
	987	zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
[704]	988	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	989	CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d)
[644]	990	c
[704]	991	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d)
	992	CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u)
[644]	993	c
[704]	994	cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d)
	995	CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v)
[644]	996
	997	cEND 210404 END
	998	c
	999	c cldtrue
	1000	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d)
[704]	1001	c CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d,
[644]	1002	c $ iimjjmp1klev,ndex3d)
	1003	c
	1004	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d)
[704]	1005	c CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d,
[644]	1006	c $ iimjjmp1klev,ndex3d)
	1007	c
	1008	c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d)
[704]	1009	c CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d,
[644]	1010	c $ iimjjmp1klev,ndex3d)
	1011	c
	1012	DO i=1, klon
	1013	cIM 151004 zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1)
	1014	zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i)
	1015	ENDDO
[704]	1016	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1017	CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d)
[644]	1018	c
	1019	DO i=1, klon
	1020	cIM 151004 zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1)
	1021	zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i)
	1022	ENDDO
[704]	1023	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1024	CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d)
[644]	1025	c
	1026	c220404 END
	1027	cAMIP2 END
	1028	c
[704]	1029	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
	1030	CALL histwrite_phy(nid_mth,"slp",itau_w,slp)
[644]	1031
[704]	1032	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
	1033	CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol)
[644]	1034	c
[704]	1035	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d)
	1036	CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m)
[644]	1037	c
	1038	c ENSEMBLES BEG
[704]	1039	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	1040	CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m)
[644]	1041	c
[704]	1042	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	1043	CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m)
[644]	1044	c
	1045	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d)
	1046	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d)
	1047	c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d)
[704]	1048	c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d,
[644]	1049	c . iim*jjmp1,ndex2d)
	1050	c
	1051	DO i=1, klon
	1052	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	1053	ENDDO
[704]	1054	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1055	CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d)
[644]	1056	c
	1057	DO i = 1, klon
	1058	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	1059	ENDDO
[704]	1060	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1061	CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d)
[644]	1062	c
	1063	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
[704]	1064	c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d,
[644]	1065	c . iimjjmp1klev,ndex3d)
	1066	c
	1067	c ENSEMBLES END
	1068	c
[704]	1069	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d)
	1070	CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m)
[644]	1071	c
[704]	1072	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d)
	1073	CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m)
[644]	1074	C
[704]	1075	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d)
	1076	CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m)
[644]	1077	C
	1078	DO i = 1, klon
	1079	zx_tmp_fi2d(i) = paprs(i,1)
	1080	ENDDO
[704]	1081	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1082	CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d)
[644]	1083	c
[704]	1084	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d)
	1085	CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf)
[644]	1086	c
	1087	if (.not. ok_veget) then
[704]	1088	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d)
	1089	CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol)
[644]	1090	endif
	1091	c
	1092	DO i = 1, klon
	1093	zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
	1094	ENDDO
[704]	1095	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1096	CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d)
[644]	1097	c
[704]	1098	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d)
	1099	CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain)
[644]	1100	c
	1101	DO i = 1, klon
	1102	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	1103	ENDDO
[704]	1104	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1105	CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d)
[644]	1106	c
	1107	DO i = 1, klon
	1108	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	1109	ENDDO
[704]	1110	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1111	CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d)
[644]	1112	c
[704]	1113	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
	1114	CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall)
[644]	1115	c
	1116	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
[704]	1117	c CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
[644]	1118	c . ndex2d)
	1119	c
[704]	1120	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
	1121	CALL histwrite_phy(nid_mth,"evap",itau_w,evap)
[644]	1122	c
[704]	1123	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
	1124	CALL histwrite_phy(nid_mth,"tops",itau_w,topsw)
[644]	1125	c
[704]	1126	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d)
	1127	CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0)
[644]	1128	c
[704]	1129	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
	1130	CALL histwrite_phy(nid_mth,"topl",itau_w,toplw)
[644]	1131	c
[704]	1132	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d)
	1133	CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0)
[644]	1134	c
	1135	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 )
[704]	1136	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1137	CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d)
[644]	1138	c
	1139	c
	1140	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 )
[704]	1141	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1142	CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d)
[644]	1143	c
	1144	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 )
[704]	1145	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1146	CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d)
[644]	1147	c
	1148	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 )
[704]	1149	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1150	CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d)
[644]	1151	c
[704]	1152	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d)
	1153	CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200)
[644]	1154	c
[704]	1155	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d)
	1156	CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr)
[644]	1157	c
[704]	1158	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d)
	1159	CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200)
[644]	1160	c
[704]	1161	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d)
	1162	CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr)
[644]	1163	c
[704]	1164	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d)
	1165	CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200)
[644]	1166	c
[704]	1167	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d)
	1168	CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr)
[644]	1169	c
	1170	zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon)
[704]	1171	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1172	CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d)
[644]	1173	c
	1174	zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon)
[704]	1175	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1176	CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d)
[644]	1177	c
	1178	c200 END
[704]	1179	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
	1180	CALL histwrite_phy(nid_mth,"sols",itau_w,solsw)
[644]	1181	c
[704]	1182	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d)
	1183	CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0)
[644]	1184	c
[704]	1185	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
	1186	CALL histwrite_phy(nid_mth,"soll",itau_w,sollw)
[644]	1187	c
[704]	1188	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d)
	1189	CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol)
[684]	1190	c
[704]	1191	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d)
	1192	CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0)
[644]	1193	c
	1194	c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
	1195	zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 )
[704]	1196	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1197	CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d)
[644]	1198	c
	1199	c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
	1200	zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 )
[704]	1201	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1202	CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d)
[644]	1203	c
	1204	c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
	1205	zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 )
[704]	1206	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1207	CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d)
[644]	1208	c
	1209	c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
	1210	zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 )
[704]	1211	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1212	CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d)
[644]	1213	c
	1214	cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon)
	1215	cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon)
	1216	zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon)
[704]	1217	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	1218	CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d)
[644]	1219	c
	1220	DO i=1, klon
	1221	sollwdownclr(i)=-1.*lwdn0(i,1)
	1222	ENDDO !i=1, klon
	1223	c
	1224	cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon)
	1225	cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon)
	1226	zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
[704]	1227	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	1228	CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d)
[644]	1229	c
[704]	1230	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
	1231	CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown)
[644]	1232	c
	1233	c DO i=1, klon
	1234	cIM 08004 sollwdownclr(i)=lwdn0(i,1)
	1235	c sollwdownclr(i)=-1.*lwdn0(i,1)
	1236	c ENDDO !i=1, klon
	1237	c
[704]	1238	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d)
	1239	CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr)
[644]	1240	c
[704]	1241	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
	1242	CALL histwrite_phy(nid_mth,"bils",itau_w,bils)
[644]	1243	c
	1244	zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
	1245	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
[704]	1246	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1247	CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d)
[644]	1248	c
[704]	1249	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
	1250	CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
[644]	1251	c
	1252	c
	1253	c DO i = 1, klon
	1254	c zx_tmp_fi2d(i) = fluxu(i,1)
	1255	c ENDDO
	1256	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
[704]	1257	c CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[644]	1258	c
	1259	c DO i = 1, klon
	1260	c zx_tmp_fi2d(i) = fluxv(i,1)
	1261	c ENDDO
	1262	c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
[704]	1263	c CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[644]	1264	c
	1265	c
[704]	1266	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d)
	1267	CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte)
[644]	1268	c
[704]	1269	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d)
	1270	CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving)
[644]	1271	c
	1272	DO nsrf = 1, nbsrf
	1273	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
[704]	1274	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1275	CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w,
	1276	$ zx_tmp_fi2d)
[644]	1277	C
	1278	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
[704]	1279	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1280	CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w,
	1281	$ zx_tmp_fi2d)
[644]	1282	ENDDO
	1283	c
	1284	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100.
[704]	1285	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1286	CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w,
	1287	$ zx_tmp_fi2d)
[644]	1288	c
	1289	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)
[704]	1290	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	1291	CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w,
	1292	$ zx_tmp_fi2d)
[644]	1293	c
	1294	DO nsrf = 1, nbsrf
	1295	CYYY
	1296	C
	1297	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
[704]	1298	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1299	CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w,
	1300	$ zx_tmp_fi2d)
[644]	1301	C
	1302	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
[704]	1303	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1304	CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w,
	1305	$ zx_tmp_fi2d)
[644]	1306	C
	1307	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
[704]	1308	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1309	CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w,
	1310	$ zx_tmp_fi2d)
[644]	1311	C
	1312	zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
[704]	1313	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1314	CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w,
	1315	$ zx_tmp_fi2d)
[644]	1316	c
	1317	zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
[704]	1318	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1319	CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w,
	1320	$ zx_tmp_fi2d)
[644]	1321	C
	1322	zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
[704]	1323	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1324	CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w,
	1325	$ zx_tmp_fi2d)
[644]	1326	C
[684]	1327	zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
[704]	1328	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1329	CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w,
	1330	$ zx_tmp_fi2d)
[684]	1331	C
[644]	1332	IF(1.EQ.0) THEN
	1333	C
	1334	zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf)
[704]	1335	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1336	CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w,
	1337	$ zx_tmp_fi2d)
[644]	1338	C
	1339	zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon))
	1340	$ *pctsrf(1:klon,nsrf)
[704]	1341	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1342	CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w,
	1343	$ zx_tmp_fi2d)
[644]	1344	C
	1345	ENDIF !(1.EQ.0) THEN
	1346	C
	1347	END DO
	1348	c
[704]	1349	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
	1350	CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm)
[644]	1351	c
[704]	1352	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
	1353	CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh)
[644]	1354	c
[704]	1355	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
	1356	CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl)
[644]	1357	c
[704]	1358	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
	1359	CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm)
[644]	1360	c
[704]	1361	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
	1362	CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh)
[644]	1363	c
[704]	1364	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d)
	1365	CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.)
[644]	1366	c
[704]	1367	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
	1368	CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq)
[644]	1369	c
	1370	zx_tmp_fi2d(1:klon) = flwp(1:klon)
[704]	1371	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1372	CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d)
[644]	1373	c
	1374	zx_tmp_fi2d(1:klon) = fiwp(1:klon)
[704]	1375	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1376	CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d)
[644]	1377	c
[704]	1378	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
	1379	CALL histwrite_phy(nid_mth,"ue",itau_w,ue)
[644]	1380	c
[704]	1381	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
	1382	CALL histwrite_phy(nid_mth,"ve",itau_w,ve)
[644]	1383	c
[704]	1384	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d)
	1385	CALL histwrite_phy(nid_mth,"uq",itau_w,uq)
[644]	1386	c
[704]	1387	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d)
	1388	CALL histwrite_phy(nid_mth,"vq",itau_w,vq)
[644]	1389	cKE43
	1390	IF (iflag_con.GE.3) THEN ! sb
	1391	c
[704]	1392	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d)
	1393	CALL histwrite_phy(nid_mth,"cape",itau_w,cape)
[644]	1394	c
[704]	1395	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d)
	1396	CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase)
[644]	1397	c
[704]	1398	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d)
	1399	CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct)
[644]	1400	c
[704]	1401	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d)
	1402	CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf)
[644]	1403	c
	1404	c34EK
	1405	c
[704]	1406	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
	1407	CALL histwrite_phy(nid_mth,"prw",itau_w,prw)
[644]	1408	c
	1409	ENDIF !iflag_con.GE.3
	1410	c
	1411	cIM cf. AM 081204 BEG
	1412	cHBTM2
	1413	c
	1414	c
[704]	1415	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
	1416	CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh)
[644]	1417	c
[704]	1418	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
	1419	CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt)
[644]	1420	c
[704]	1421	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
	1422	CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl)
[644]	1423	c
[704]	1424	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
	1425	CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL)
[644]	1426	c
[704]	1427	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
	1428	CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL)
[644]	1429	c
[704]	1430	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
	1431	CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL)
[644]	1432	c
[704]	1433	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
	1434	CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm)
[644]	1435	c
[704]	1436	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
	1437	CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1)
[644]	1438	c
[704]	1439	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
	1440	CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2)
[644]	1441	c
[704]	1442	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
	1443	CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3)
[644]	1444	c
	1445	cIM cf. AM 081204 END
	1446	c
	1447	c Champs interpolles sur des niveaux de pression
	1448	c DO k=1, nlevENS
	1449	DO k=1, nlevSTD
	1450	c
	1451	c bb=clevSTD(k)
	1452	c
	1453	c IF(k.GE.2) THEN
	1454	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	1455	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	1456	c aa=clevSTD(k)
	1457	c bb=aa(1:lnblnk1(aa))
	1458	c ENDIF
	1459	c
	1460	IF(bb2.EQ."850".OR.bb2.EQ."700"
	1461	$ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN
	1462	c
[704]	1463	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d)
	1464	CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2))
[644]	1465	c
[704]	1466	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d)
	1467	CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2))
[644]	1468	c
[704]	1469	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d)
	1470	CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2))
[644]	1471	c
[704]	1472	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d)
	1473	CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2))
[644]	1474	c
	1475	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700"
	1476	ENDDO !nlevSTD
	1477	c
	1478	DO i=1, klon
	1479	IF (pctsrf(i,is_oce).GT.epsfra) THEN
	1480	zx_tmp_fi2d(i) = fluxo(i)
	1481	ELSE
[684]	1482	zx_tmp_fi2d(i) = 0.
[644]	1483	ENDIF
	1484	ENDDO
	1485	c
[704]	1486	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1487	CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d)
[644]	1488	c
	1489	DO i=1, klon
	1490	IF (pctsrf(i,is_sic).GT.epsfra) THEN
	1491	zx_tmp_fi2d(i) = fluxg(i)
	1492	ELSE
[684]	1493	zx_tmp_fi2d(i) = 0.
[644]	1494	ENDIF
	1495	ENDDO
	1496	c
[704]	1497	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1498	CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d)
[644]	1499	c
[684]	1500	DO i=1, klon
	1501	IF (pctsrf(i,is_oce).GT.epsfra.OR.
	1502	. pctsrf(i,is_sic).GT.epsfra) THEN
	1503	zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+
	1504	. ftsol(i, is_sic) * pctsrf(i,is_sic))/
	1505	. (pctsrf(i,is_oce)+pctsrf(i,is_sic))
	1506	ELSE
	1507	zx_tmp_fi2d(i) = 273.15
	1508	ENDIF
	1509	ENDDO
	1510	c
[704]	1511	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1512	CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d)
[684]	1513	c
[644]	1514	IF (OCEAN.EQ.'force ') THEN
	1515	c
	1516	DO i=1, klon
	1517	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	1518	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	1519	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	1520	. fluxg(i)*pctsrf(i,is_sic)
	1521	ELSE
	1522	zx_tmp_fi2d(i) = 1.E+20
	1523	ENDIF
	1524	ENDDO
[704]	1525	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1526	CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,
	1527	. zx_tmp_fi2d)
[644]	1528	c
	1529	ELSE IF (OCEAN.EQ.'slab ') THEN
	1530	c
	1531	c
	1532	DO i=1, klon
	1533	IF((pctsrf(i,is_oce).GT.epsfra).OR.
	1534	. (pctsrf(i,is_sic).GT.epsfra)) THEN
	1535	zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+
	1536	. fluxg(i)*pctsrf(i,is_sic)
	1537	ELSE
	1538	zx_tmp_fi2d(i) = 1.E+20
	1539	ENDIF
	1540	ENDDO
[704]	1541	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1542	CALL histwrite_phy(nid_mth,"slab_bils",itau_w,
	1543	. zx_tmp_fi2d)
[644]	1544	c
	1545	DO i=1, klon
	1546	IF(pctsrf(i,is_oce).GT.epsfra.OR.
	1547	$ pctsrf(i,is_sic).GT.epsfra) THEN
	1548	zx_tmp_fi2d(i)=tslab(i)
	1549	ELSE
	1550	zx_tmp_fi2d(i) = 1.E+20
	1551	ENDIF
	1552	ENDDO !i=1, klon
[704]	1553	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	1554	CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d)
[644]	1555	c
[704]	1556	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d)
	1557	CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice)
[644]	1558	c
[704]	1559	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d)
	1560	CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.)
[644]	1561	c
	1562	ENDIF !OCEAN.EQ.force/slab
	1563	c
	1564	ENDIF !lev_histmth.GE.1
	1565	c
	1566	IF(lev_histmth.GE.2) THEN
	1567	c
	1568	c Champs 3D:
	1569	C
	1570	c
[704]	1571	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d)
	1572	CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc)
[644]	1573	c
[704]	1574	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d)
	1575	CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc)
[644]	1576	c
[704]	1577	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
	1578	CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri)
[524]	1579	c
[644]	1580	c#ifdef histmthNMC
	1581	c#include "write_histmthNMC.h"
	1582	c#endif
	1583	c
[704]	1584	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
	1585	CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap))
[524]	1586	c
[644]	1587	c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/
	1588	c . (1-qx(1:klon,1:klev,ivap))
	1589	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
[704]	1590	c CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d,
[644]	1591	c . iimjjmp1klev,ndex3d)
	1592	c
[704]	1593	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
	1594	CALL histwrite_phy(nid_mth,"geop",itau_w,zphi)
[524]	1595	c
[704]	1596	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
	1597	CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri)
[524]	1598	c
[704]	1599	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
	1600	CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri)
[644]	1601	c
[704]	1602	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
	1603	CALL histwrite_phy(nid_mth,"vitw",itau_w,omega)
[644]	1604	c
[704]	1605	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
	1606	CALL histwrite_phy(nid_mth,"pres",itau_w,pplay)
[644]	1607	c
[704]	1608	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d)
	1609	CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.)
[644]	1610	c
[704]	1611	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d)
	1612	CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon)
[644]	1613	c
[704]	1614	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d)
	1615	CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh)
[644]	1616	c
	1617	DO k=1, klev
	1618	DO i=1, klon
	1619	cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968
	1620	c $ /(paprs(i,k)-paprs(i,k+1))
	1621	c $ *(paprs(i,1)/101325.0)
	1622	c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)1.0e+6,1.0e-6)RG/46.6968
	1623	zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968
	1624	$ /(paprs(i,k)-paprs(i,k+1))
	1625	$ *(paprs(i,1)/101325.0)
	1626	ENDDO !i
	1627	ENDDO !k
	1628	c
	1629	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d)
	1630	c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d)
[704]	1631	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1632	CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d)
[644]	1633	c
[704]	1634	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
	1635	CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t)
[644]	1636	c
[704]	1637	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d)
	1638	CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap))
[644]	1639	c
	1640	ENDIF !lev_histmth.GE.2
	1641	c
	1642	IF(lev_histmth.GE.3) THEN
	1643	c
	1644	DO nsrf=1, nbsrf
	1645	c
	1646	zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
[704]	1647	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1648	CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w,
	1649	$ zx_tmp_fi2d)
[644]	1650	C
	1651	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
[704]	1652	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	1653	CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w,
	1654	$ zx_tmp_fi2d)
[644]	1655	c
	1656	zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
[704]	1657	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d)
	1658	CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w
	1659	$ ,agesno)
[644]	1660	c
	1661	ENDDO !nsrf=1, nbsrf
	1662	c
[704]	1663	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
	1664	CALL histwrite_phy(nid_mth,"albs",itau_w,albsol)
[644]	1665	c
[704]	1666	cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
	1667	CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw)
[644]	1668	c
	1669	ENDIF !lev_histmth.GE.3
	1670	c
	1671	c
	1672	IF(lev_histmth.GE.4) THEN
	1673	c
[704]	1674	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
	1675	CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0)
[644]	1676	c
[653]	1677	IF (iflag_con.GE.3) THEN
	1678	c
[704]	1679	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d)
	1680	CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma)
[644]	1681	c
[704]	1682	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d)
	1683	CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd)
[644]	1684	c
[704]	1685	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d)
	1686	CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd)
[644]	1687	c
[704]	1688	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d)
	1689	CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0)
[644]	1690	c
[653]	1691	ENDIF !iflag_con.GE.3
[644]	1692	c
[704]	1693	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
	1694	CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn)
[644]	1695	c
[704]	1696	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d)
	1697	CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn)
[644]	1698	c
	1699	cIM: 101003 : K/30min ==> K/s
	1700	c
	1701	zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
[704]	1702	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1703	CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d)
[644]	1704	c
	1705	cIM: 101003 : K/30min ==> K/s
	1706	zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+
	1707	. d_t_eva(1:klon,1:klev))/pdtphys
[704]	1708	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d)
	1709	CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d)
[644]	1710	c
	1711	zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
[704]	1712	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1713	CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d)
[644]	1714	c
	1715	cIM: 101003 : K/30min ==> K/s
	1716	zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
[704]	1717	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1718	CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d)
[644]	1719	c
	1720	zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
[704]	1721	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1722	CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d)
[644]	1723	c
	1724	cIM: 101003 : K/30min ==> K/s
	1725	zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
[704]	1726	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1727	CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d)
[644]	1728	c
	1729	zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
[704]	1730	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1731	CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d)
[644]	1732	c
	1733	zpt_conv = 0.
	1734	where (ptconv) zpt_conv = 1.
[704]	1735	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
	1736	CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv)
[644]	1737	c
[704]	1738	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d)
	1739	CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs)
[644]	1740	c
	1741	zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys
[704]	1742	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1743	CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d)
[644]	1744	c
	1745	cIM: 101003 : K/day ==> K/s
	1746	zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
[704]	1747	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1748	CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d)
[644]	1749	c
	1750	cIM: 101003 : K/day ==> K/s
	1751	zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
[704]	1752	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1753	CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d)
[644]	1754	c
	1755	cIM: 101003 : K/day ==> K/s
	1756	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
[704]	1757	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1758	CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d)
[644]	1759	c
	1760	cIM: 101003 : K/day ==> K/s
	1761	zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
[704]	1762	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1763	CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d)
[644]	1764	c
	1765	cIM: 101003 : deja en K/s
	1766	zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
[704]	1767	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1768	CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d)
[644]	1769	c
	1770	zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
[704]	1771	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1772	CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d)
[644]	1773	c
	1774	zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
[704]	1775	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1776	CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d)
[644]	1777	c
	1778	IF (ok_orodr) THEN
	1779	zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
[704]	1780	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1781	CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d)
[644]	1782	c
	1783	zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
[704]	1784	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1785	CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d)
[644]	1786	c
	1787	ENDIF
	1788	C
	1789	IF (ok_orolf) THEN
	1790	zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
[704]	1791	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1792	CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d)
[644]	1793	c
	1794	zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
[704]	1795	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
	1796	CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d)
[644]	1797	ENDIF
	1798	C
[524]	1799	IF (nqmax.GE.3) THEN
[653]	1800	DO iq=3,nqmax
[704]	1801	cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d)
	1802	CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq))
[524]	1803	ENDDO
	1804	ENDIF
	1805	c
	1806	ENDIF !lev_histmth.GE.4
	1807	c
[644]	1808	ENDIF !(type_run.EQ."AMIP")
	1809	c
[524]	1810	if (ok_sync) then
[735]	1811	c$OMP MASTER
[524]	1812	call histsync(nid_mth)
[735]	1813	c$OMP END MASTER
[524]	1814	endif
[644]	1815	c
	1816	ENDIF !ok_mensuel

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

Context Navigation

source: LMDZ4/branches/V3_test/libf/phylmd/write_histmth.h

Download in other formats: