Context Navigation

source: LMDZ5/branches/testing/libf/phylmd/ini_histrac.h @ 2283

Last change on this file since 2283 was 1910, checked in by Laurent Fairhead, 11 years ago
Merged trunk changes r1860:1909 into testing branch
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 11.0 KB

Rev	Line
[524]	1	!
[1279]	2	! $Id $
[524]	3	!
[1664]	4	IF (ecrit_tra>0.) THEN
[1279]	5	!$OMP MASTER
[1750]	6	!!! CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
	7	! correction pour l heure initiale !jyg
	8	! !jyg
	9	CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian) !jyg
	10
[1279]	11	CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
[1750]	12	CALL histvert(nid_tra, "presnivs", "Vertical levels", "Pa",klev, presnivs, nvert,"down")
[524]	13
[1279]	14	zsto = pdtphys
	15	zout = ecrit_tra
	16	CALL histdef(nid_tra, "phis", "Surface geop. height", "-", &
	17	iim,jj_nb,nhori, 1,1,1, -99, 32,"once", zsto,zout)
	18	CALL histdef(nid_tra, "aire", "Grid area", "-", &
	19	iim,jj_nb,nhori, 1,1,1, -99, 32,"once", zsto,zout)
[1403]	20	CALL histdef(nid_tra, "zmasse", "column density of air in cell", &
	21	"kg m-2", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
	22	zsto,zout)
[1750]	23	! RomP >>>
	24	CALL histdef(nid_tra, "sourceBE", "source 7Be", &
	25	"at/kgA/s", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
	26	zsto,zout)
	27	! RomP <<<
[524]	28
[1279]	29	!TRACEURS
	30	!----------------
	31	DO it = 1,nbtr
	32	iiq = niadv(it+2)
[524]	33
[1279]	34	! CONCENTRATIONS
	35	CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga", &
	36	iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
[549]	37
[1279]	38	! TD LESSIVAGE
	39	IF (lessivage .AND. aerosol(it)) THEN
	40	CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
[1750]	41	"at/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
[1279]	42	"ave(X)", zsto,zout)
[1750]	43	CALL histdef(nid_tra, "d_tr_ls_"//tname(iiq), &
	44	"tendance lessivage large scale"// ttext(iiq), "?",&
	45	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	46	"ave(X)", zsto,zout)
	47	CALL histdef(nid_tra, "d_tr_insc_"//tname(iiq), &
	48	"tendance lessivage large scale"// ttext(iiq), "?",&
	49	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	50	"ave(X)", zsto,zout)
	51	CALL histdef(nid_tra, "d_tr_bcscav_"//tname(iiq), &
	52	"tendance lessivage large scale"// ttext(iiq), "?",&
	53	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	54	"ave(X)", zsto,zout)
	55	CALL histdef(nid_tra, "d_tr_evls_"//tname(iiq), &
	56	"tendance lessivage large scale"// ttext(iiq), "?",&
	57	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	58	"ave(X)", zsto,zout)
	59	! Tracer concentration in LS precipitation at surface
	60	CALL histdef(nid_tra, "qpr_ls_"//tname(iiq), &
	61	"concentration in LS precip"// ttext(iiq), "at/kgw", &
	62	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	63	"ave(X)", zsto,zout)
	64	END IF
[549]	65
[1279]	66	! TD THERMIQUES
	67	IF (iflag_thermals.gt.0) THEN
	68	CALL histdef(nid_tra, "d_tr_th_"//tname(iiq), &
	69	"tendance thermique"// ttext(iiq), "?", &
	70	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	71	"ave(X)", zsto,zout)
	72	ENDIF
[549]	73
[1279]	74	! TD CONVECTION
	75	IF (iflag_con.GE.2) THEN
	76	CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq), &
	77	"tendance convection"// ttext(iiq), "?", &
	78	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	79	"ave(X)", zsto,zout)
	80	ENDIF
[524]	81
[1750]	82	! RomP >>>
	83	IF (iflag_con.EQ.30) THEN
	84	CALL histdef(nid_tra, "d_tr_cvMA_"//tname(iiq), &
	85	"tendance convection"// ttext(iiq), "?",&
	86	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	87	"ave(X)", zsto,zout)
	88	CALL histdef(nid_tra, "d_tr_trsp_"//tname(iiq), &
	89	"tendance transport "// ttext(iiq), "at/kga", &
	90	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	91	"ave(X)", zsto,zout)
	92	CALL histdef(nid_tra, "d_tr_sscav_"//tname(iiq), &
	93	"tendance lessivage flux satures "// ttext(iiq), "at/kga", &
	94	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	95	"ave(X)", zsto,zout)
	96	CALL histdef(nid_tra, "d_tr_sat_"//tname(iiq), &
	97	"tendance flux satures "// ttext(iiq), "at/kga", &
	98	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	99	"ave(X)", zsto,zout)
	100	CALL histdef(nid_tra, "d_tr_uscav_"//tname(iiq), &
	101	"tendance flux insatures "// ttext(iiq), "at/kga", &
	102	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	103	"ave(X)", zsto,zout)
	104	CALL histdef(nid_tra, "tr_pr_"//tname(iiq), &
	105	"concentration dans precip"// ttext(iiq), "at/kga", &
	106	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	107	"ave(X)", zsto,zout)
	108	CALL histdef(nid_tra, "tr_aa_"//tname(iiq), &
	109	"concentration precip issu AA"// ttext(iiq), "at/kga", &
	110	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	111	"ave(X)", zsto,zout)
	112	CALL histdef(nid_tra, "tr_mel_"//tname(iiq), &
	113	"concentration precip issu melange"// ttext(iiq), "at/kga", &
	114	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	115	"ave(X)", zsto,zout)
	116	CALL histdef(nid_tra, "tr_di_"//tname(iiq), &
	117	"concentration dans descente insaturee"// ttext(iiq), "at/kga", &
	118	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	119	"ave(X)", zsto,zout)
	120	CALL histdef(nid_tra, "tr_trspdi_"//tname(iiq), &
	121	"conc descente insaturee MA"// ttext(iiq), "at/kga", &
	122	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	123	"ave(X)", zsto,zout)
	124	CALL histdef(nid_tra, "zmfd1a_"//tname(iiq), &
	125	"zmfd1a"// ttext(iiq), "_", &
	126	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	127	"ave(X)", zsto,zout)
	128	CALL histdef(nid_tra, "zmfphi2_"//tname(iiq), &
	129	"zmfphi2"// ttext(iiq), "_", &
	130	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	131	"ave(X)", zsto,zout)
	132	CALL histdef(nid_tra, "zmfdam_"//tname(iiq), &
	133	"zmfdam"// ttext(iiq), "_", &
	134	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	135	"ave(X)", zsto,zout)
	136	ENDIF
	137	! RomP <<<
	138	CALL histdef(nid_tra, "dtrdyn_"//tname(iiq), &
	139	"td dyn tra"// ttext(iiq), "at/kga", &
	140	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	141	"ave(X)", zsto,zout)
	142	! TD decroissance radioactive
	143	CALL histdef(nid_tra, "d_tr_dec_"//tname(iiq), &
	144	"tendance decroi radio "// ttext(iiq), "", &
	145	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	146	"ave(X)", zsto,zout)
	147
[1279]	148	! TD COUCHE-LIMITE
[1864]	149	IF (iflag_vdf_trac>=0) THEN
[1279]	150	CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq), &
	151	"tendance couche limite"// ttext(iiq), "?", &
	152	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	153	"ave(X)", zsto,zout)
[1750]	154	! Dry deposit (1st layer and surface)
	155	CALL histdef(nid_tra, "d_tr_dry_"//tname(iiq), &
	156	"tendancy dry deposit"// ttext(iiq), "at/kga/step", &
	157	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	158	"ave(X)", zsto,zout)
	159	CALL histdef(nid_tra, "flux_tr_dry_"//tname(iiq), &
	160	"dry deposit at surf (downward)"// ttext(iiq), "at/m2/step", &
	161	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	162	"ave(X)", zsto,zout)
	163	ENDIF
[1279]	164	ENDDO
[1750]	165
	166	CALL histdef(nid_tra, "Mint", "Mint","", &
	167	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	168	"inst(X)", zout,zout)
	169	CALL histdef(nid_tra, "frac_impa", "frac_impa","", &
	170	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	171	"inst(X)", zout,zout)
	172	CALL histdef(nid_tra, "frac_nucl", "frac_nucl","", &
	173	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	174	"inst(X)", zout,zout)
[1279]	175	!---------------
	176	!
	177	! VENT (niveau 1)
	178	CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", &
	179	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	180	"inst(X)", zout,zout)
	181	CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", &
	182	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	183	"inst(X)", zout,zout)
[524]	184
[1279]	185	! TEMPERATURE DU SOL
	186	CALL histdef(nid_tra, "ftsol1", "temper sol", "-", &
	187	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	188	"inst(X)", zout,zout)
	189	CALL histdef(nid_tra, "ftsol2", "temper sol", "-", &
	190	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	191	"inst(X)", zout,zout)
	192	CALL histdef(nid_tra, "ftsol3", "temper sol", "-", &
	193	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	194	"inst", zout,zout)
	195	CALL histdef(nid_tra, "ftsol4", "temper sol", "-", &
	196	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	197	"inst(X)", zout,zout)
	198
	199	! NATURE DU SOL
	200	CALL histdef(nid_tra, "psrf1", "nature sol", "-", &
	201	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	202	"inst(X)", zout,zout)
	203	CALL histdef(nid_tra, "psrf2", "nature sol", "-", &
	204	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	205	"inst(X)", zout,zout)
	206	CALL histdef(nid_tra, "psrf3", "nature sol", "-", &
	207	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	208	"inst(X)", zout,zout)
	209	CALL histdef(nid_tra, "psrf4", "nature sol", "-", &
	210	iim,jj_nb,nhori, 1,1,1, -99, 32, &
	211	"inst(X)", zout,zout)
	212	! DIVERS
[1403]	213	CALL histdef(nid_tra, "pplay", "pressure","-", &
[1279]	214	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	215	"inst(X)", zout,zout)
[1403]	216	CALL histdef(nid_tra, "T", "temperature","K", &
[1279]	217	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	218	"inst(X)", zout,zout)
	219	CALL histdef(nid_tra, "mfu", "flux u mont","-", &
	220	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	221	"ave(X)", zsto,zout)
	222	CALL histdef(nid_tra, "mfd", "flux u decen","-", &
	223	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	224	"ave(X)", zsto,zout)
	225	CALL histdef(nid_tra, "en_u", "flux u mont","-", &
	226	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	227	"ave(X)", zsto,zout)
	228	CALL histdef(nid_tra, "en_d", "flux u mont","-", &
	229	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	230	"ave(X)", zsto,zout)
	231	CALL histdef(nid_tra, "de_d", "flux u mont","-", &
	232	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	233	"ave(X)", zsto,zout)
	234	CALL histdef(nid_tra, "de_u", "flux u decen","-", &
	235	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	236	"ave(X)", zsto,zout)
	237	CALL histdef(nid_tra, "coefh", "turbulent coef","-", &
	238	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
	239	"ave(X)", zsto,zout)
	240
	241	CALL histend(nid_tra)
	242	!$OMP END MASTER
[1664]	243	END IF ! ecrit_tra>0.
[1279]	244

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

Download in other formats: