Context Navigation

source: LMDZ4/branches/LMDZ4-dev/libf/phylmd/ini_histrac.h @ 1118

Last change on this file since 1118 was 1115, checked in by jghattas, 16 years ago

Si ecrit_tra=0, pas de creation ou ecriture de la fichier histrac.nc

JG

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 5.1 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
[1115]	4	IF (ecrit_tra>0. .AND. config_inca == 'none') THEN
[766]	5	c$OMP MASTER
[524]	6	CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
	7	c
[766]	8	CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,
[524]	9	. nhori, nid_tra)
	10	CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
	11	. klev, presnivs, nvert)
	12
	13
	14
	15	zsto = pdtphys
[956]	16	zout = ecrit_tra
[524]	17	c
	18	CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
[776]	19	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[524]	20	. "once", zsto,zout)
	21	c
	22	CALL histdef(nid_tra, "aire", "Grid area", "-",
[776]	23	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[524]	24	. "once", zsto,zout)
[1114]	25	DO it=1,nbtr
[524]	26	C champ 2D
	27	iq=it+2
	28	iiq=niadv(iq)
[1114]	29	CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",
[776]	30	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	31	. "ave(X)", zsto,zout)
	32	if (lessivage) THEN
[1114]	33	CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq),
[776]	34	. "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	35	. "ave(X)", zsto,zout)
	36	endif
[549]	37
	38	c---Ajout Olivia
[1114]	39	CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),
[549]	40	. "tendance thermique"// ttext(iiq), "?",
[776]	41	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]	42	. "ave(X)", zsto,zout)
[973]	43	c
	44	if(iflag_con.GE.2) then
[1114]	45	CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),
[549]	46	. "tendance convection"// ttext(iiq), "?",
[776]	47	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]	48	. "ave(X)", zsto,zout)
[973]	49	endif !(iflag_con.GE.2) then
[1114]	50	CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),
[549]	51	. "tendance couche limite"// ttext(iiq), "?",
[776]	52	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]	53	. "ave(X)", zsto,zout)
	54	c---fin Olivia
	55
[524]	56	ENDDO
[549]	57
[524]	58	CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
[776]	59	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[766]	60	. "inst(X)", zout,zout)
[524]	61
	62	CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
[776]	63	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	64	. "inst(X)", zout,zout)
[524]	65	CALL histdef(nid_tra, "psrf1", "nature sol", "-",
[776]	66	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	67	. "inst(X)", zout,zout)
[524]	68	CALL histdef(nid_tra, "psrf2", "nature sol", "-",
[776]	69	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	70	. "inst(X)", zout,zout)
[524]	71	CALL histdef(nid_tra, "psrf3", "nature sol", "-",
[776]	72	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	73	. "inst(X)", zout,zout)
[524]	74	CALL histdef(nid_tra, "psrf4", "nature sol", "-",
[776]	75	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	76	. "inst(X)", zout,zout)
[524]	77	CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
[776]	78	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	79	. "inst(X)", zout,zout)
[524]	80	CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
[776]	81	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	82	. "inst(X)", zout,zout)
[524]	83	CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
[776]	84	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[766]	85	. "inst", zout,zout)
[524]	86	CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
[776]	87	. iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]	88	. "inst(X)", zout,zout)
[524]	89	CALL histdef(nid_tra, "pplay", "flux u mont","-",
[776]	90	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[657]	91	. "inst(X)", zout,zout)
[524]	92	CALL histdef(nid_tra, "t", "flux u mont","-",
[776]	93	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[657]	94	. "inst(X)", zout,zout)
[524]	95	CALL histdef(nid_tra, "mfu", "flux u mont","-",
[776]	96	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	97	. "ave(X)", zsto,zout)
	98	CALL histdef(nid_tra, "mfd", "flux u decen","-",
[776]	99	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	100	. "ave(X)", zsto,zout)
	101	CALL histdef(nid_tra, "en_u", "flux u mont","-",
[776]	102	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	103	. "ave(X)", zsto,zout)
	104	CALL histdef(nid_tra, "en_d", "flux u mont","-",
[776]	105	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	106	. "ave(X)", zsto,zout)
	107	CALL histdef(nid_tra, "de_d", "flux u mont","-",
[776]	108	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	109	. "ave(X)", zsto,zout)
	110	CALL histdef(nid_tra, "de_u", "flux u decen","-",
[776]	111	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	112	. "ave(X)", zsto,zout)
	113	CALL histdef(nid_tra, "coefh", "turbulent coef","-",
[776]	114	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]	115	. "ave(X)", zsto,zout)
	116
	117	c
	118	CALL histend(nid_tra)
	119	ndex2d = 0
	120	ndex3d = 0
	121	ndex = 0
[766]	122	c$OMP END MASTER
[1115]	123	END IF ! ecrit_tra>0. .AND. config_inca == 'none'

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

Download in other formats: