Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histrac.h @ 646

Last change on this file since 646 was 646, checked in by lmdzadmin, 20 years ago
Nouvelles sorties pour INCA_AER. AC MAF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 18.8 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
	5	c
	6	CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon)
	7	DO i = 1, iim
	8	zx_lon(i,1) = xlon(i+1)
	9	zx_lon(i,jjm+1) = xlon(i+1)
	10	ENDDO
	11	CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat)
	12	CALL histbeg("histrac", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),
	13	. 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,
	14	. nhori, nid_tra)
	15	CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
	16	. klev, presnivs, nvert)
	17
	18	#ifdef INCA_AER
	19	CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,
	20	. 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,
	21	. nhori, nid_tra2)
	22
	23	CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,
	24	. 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,
	25	. nhori, nid_tra3)
	26
	27	call histvert(nid_tra2, "presnivs", "presnivs", "mb",
	28	. klev, presnivs, nvert)
	29	call histvert(nid_tra3, "presnivs", "presnivs", "mb",
	30	. klev, presnivs, nvert)
	31	#endif
	32
	33	#ifdef INCA
	34	! call histvert(nid_tra, "ap", "Hybrid A parameter", "-",
	35	! . klev+1, ap, nverta)
	36	! call histvert(nid_tra, "bp", "Hybrid B parameter", "-",
	37	! . klev+1, bp, nvertb)
	38	#endif
	39
	40	zsto = pdtphys
	41	zout = pdtphys * FLOAT(ecrit_tra)
	42	c
	43	CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
	44	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	45	. "once", zsto,zout)
	46	c
	47	CALL histdef(nid_tra, "aire", "Grid area", "-",
	48	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	49	. "once", zsto,zout)
	50	#ifdef INCA
	51	CALL histdef(nid_tra, "ps", "Surface pressure", "Pa",
	52	. iim,jjm+1,nhori, 1,1,1,-99, 32,
	53	. "ave(X)", zsto,zout)
	54
	55	CALL histdef(nid_tra, "ptrop", "Tropopause pressure", "Pa",
	56	. iim,jjm+1,nhori, 1,1,1,-99, 32,
	57	. "ave(X)", zsto,zout)
	58
	59	C 3d FIELDS
	60	CALL histdef(nid_tra, "temp", "Air temperature", "K",
	61	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	62	. "ave(X)", zsto,zout)
	63
	64	CALL histdef(nid_tra, "u", "zonal wind component", "m/s",
	65	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	66	. "ave(X)", zsto,zout)
	67
	68	CALL histdef(nid_tra, "v", "zonal wind component", "m/s",
	69	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	70	. "ave(X)", zsto,zout)
	71
	72	CALL histdef(nid_tra, "h2o", "Specific Humidity", "MMR",
	73	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	74	. "ave(X)", zsto,zout)
	75
	76	CALL histdef(nid_tra, "pmid", "Pressure", "Pa",
	77	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	78	. "ave(X)", zsto,zout)
	79
	80	CALL histdef(nid_tra, "pdel", "Delta Pressure", "Pa",
	81	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	82	. "ave(X)", zsto,zout)
	83
	84	! MS info
[635]	85	#ifdef INCA_AER
[524]	86	CALL histdef(nid_tra, "airm", "Air mass", "kg",
	87	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	88	. "ave(X)", zsto,zout)
[635]	89	#endif
[616]	90	#ifdef INCA
[524]	91	#ifdef INCAINFO
	92	DO it=1, phtcnt
	93	WRITE(str2,'(i2.2)') it
	94	CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
	95	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	96	. "ave(X)", zsto,zout)
	97	ENDDO
	98	DO it=1, hetcnt
	99	WRITE(str2,'(i2.2)') it
	100	CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
	101	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	102	. "ave(X)", zsto,zout)
	103	ENDDO
	104
	105	DO it=1, extcnt
	106	WRITE(str2,'(i2.2)') it
	107	CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
	108	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	109	. "ave(X)", zsto,zout)
	110	ENDDO
	111
	112	DO it=1, nfs
	113	WRITE(str2,'(i2.2)') it
	114	CALL histdef(nid_tra, "INV"//str2, "INV"//str2, "CM-3",
	115	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	116	. "ave(X)", zsto,zout)
	117	ENDDO
	118
	119	#else
	120	CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
	121	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	122	. "ave(X)", zsto,zout)
	123	CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
	124	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	125	. "ave(X)", zsto,zout)
	126	CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
	127	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	128	. "ave(X)", zsto,zout)
	129	CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
	130	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	131	. "ave(X)", zsto,zout)
	132	CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
	133	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	134	. "ave(X)", zsto,zout)
	135	CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
	136	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	137	. "ave(X)", zsto,zout)
	138	#endif
	139
	140	DO it=1, grpcnt
	141	CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
	142	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	143	. "ave(X)", zsto,zout)
	144	ENDDO
	145	#endif
	146
	147	#ifdef INCA_AER
	148
	149	CALL histdef(nid_tra2, "scavcoef_st","scavcoef_st", "S-1",
	150	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	151	. "ave(X)", zsto,zout)
	152	CALL histdef(nid_tra2, "scavcoef_cv","scavcoef_cv", "S-1",
	153	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	154	. "ave(X)", zsto,zout)
	155	CALL histdef(nid_tra2, "AngstroemComp","AngstroemComp",
	156	. "angs comp", iim,jjm+1,nhori, 1,1,1, -99, 32,
	157	. "ave(X)", zsto,zout)
	158
[646]	159	CALL histdef(nid_tra2, "TOTAERH2O", "Total aerosol water",
	160	. "kg/m2", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
	161	. "ave(X)", zsto,zout)
[524]	162	#endif
	163	#endif
	164	DO it=1,nqmax
	165	C champ 2D
	166	#ifdef INCA
	167	IF ( prt_flag_ts(it) == 0 ) CYCLE
	168
	169	CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),
	170	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	171	. "ave(X)", zsto,zout)
	172	CALL histdef(nid_tra, "Dep_"//solsym(it), "Dep_"//solsym(it),
	173	. "cm/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	174	. "ave(X)", zsto,zout)
	175	#ifdef INCA_AER
[646]	176
	177	CALL histdef(nid_tra2, "Dry_"//solsym(it), "Dry_"//solsym(it),
	178	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	179	. "ave(X)", zsto,zout)
	180
[524]	181	IF ((it .ge. trmx) .and. (it .le. trnx)) then
	182	CALL histdef(nid_tra2, "Sed_"//solsym(it), "Sed_"//solsym(it),
	183	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	184	. "ave(X)", zsto,zout)
	185	CALL histdef(nid_tra2, "Wet_"//solsym(it), "Wet_"//solsym(it),
	186	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	187	. "ave(X)", zsto,zout)
	188	CALL histdef(nid_tra2, "WetST_"//solsym(it), "WetST_"//solsym(it),
	189	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	190	. "ave(X)", zsto,zout)
	191	CALL histdef(nid_tra2, "WetCV_"//solsym(it), "WetCV_"//solsym(it),
	192	. "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
	193	. "ave(X)", zsto,zout)
	194	CALL histdef(nid_tra2, "Emi_alt_"//solsym(it), "Emi_alt_"//solsym(it),
	195	. "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	196	. "ave(X)", zsto,zout)
	197	CALL histdef(nid_tra2, "Load_"//solsym(it), "Load_"//solsym(it),
	198	. "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,
	199	. "ave(X)", zsto,zout)
	200	CALL histdef(nid_tra2, "SConc_"//solsym(it), "SConc_"//solsym(it),
	201	. "kg/m3", iim,jjm+1,nhori, 1,1,1, -99, 32,
	202	. "ave(X)", zsto,zout)
	203	do la=1,las
[646]	204	CALL histdef(nid_tra2, "OD"//cla(la)//"_"//solsym(it), "OD"//cla(la)//"_"//solsym(it),
[524]	205	. "opt. depth", iim,jjm+1,nhori, 1,1,1, -99, 32,
	206	. "ave(X)", zsto,zout)
[646]	207
	208	CALL histdef(nid_tra2, "OD3D"//cla(la)//"_"//solsym(it), "OD3d_"//cla(la)//"_"//solsym(it),
	209	. "opt. depth 3D", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
	210	. "ave(X)", zsto,zout)
[524]	211	enddo
[646]	212
	213	CALL histdef(nid_tra2, "CLOAD_"//solsym(it), "3D LOAD_"//solsym(it),
	214	. "kg/m2 resp. #/m2", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
	215	. "ave(X)", zsto,zout)
	216
[524]	217	CALL histdef(nid_tra2, "MD_"//solsym(it), "MD_"//solsym(it),
[646]	218	. "median diameter [m]", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
[524]	219	. "ave(X)", zsto,zout)
[646]	220
	221	CALL histdef(nid_tra2, "MDW_"//solsym(it), "MDW_"//solsym(it),
	222	. "wet median diameter [m]", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
	223	. "ave(X)", zsto,zout)
	224
	225
	226
	227	CALL histdef(nid_tra3, "Inst_Load_"//solsym(it), "Inst_Load_"//solsym(it),
[524]	228	. "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,
	229	. "inst(X)", zout,zout)
[646]	230
[524]	231	endif
[646]	232
[524]	233	#endif
	234	CALL histdef(nid_tra, solsym(it), solsym(it), "VMR",
	235	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	236	. "ave(X)", zsto,zout)
	237	#else
	238	iq=it+2
	239	iiq=niadv(iq)
	240	CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga",
	241	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	242	. "ave(X)", zsto,zout)
	243	if (lessivage) THEN
	244	CALL histdef(nid_tra, "fl"//tnom(iq),"Flux "//ttext(iiq),
	245	. "U/m2/s",iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	246	. "ave(X)", zsto,zout)
	247	endif
[549]	248
	249	c---Ajout Olivia
	250	CALL histdef(nid_tra, "d_tr_th_"//tnom(iq),
	251	. "tendance thermique"// ttext(iiq), "?",
	252	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	253	. "ave(X)", zsto,zout)
	254	CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq),
	255	. "tendance convection"// ttext(iiq), "?",
	256	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	257	. "ave(X)", zsto,zout)
	258	CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq),
	259	. "tendance couche limite"// ttext(iiq), "?",
	260	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	261	. "ave(X)", zsto,zout)
	262	c---fin Olivia
	263
[524]	264	#endif
	265	ENDDO
[549]	266
[524]	267	#ifdef INCA
	268	CALL histdef(nid_tra, "O3_column", "O3_column",
	269	. "DU", iim,jjm+1,nhori, 1,1,1, -99, 32,
	270	. "ave(X)", zsto,zout)
	271	CALL histdef(nid_tra, "CO_column", "CO_column",
	272	. "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
	273	. "ave(X)", zsto,zout)
	274	CALL histdef(nid_tra, "CH4_column", "CH4_column",
	275	. "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
	276	. "ave(X)", zsto,zout)
	277	CALL histdef(nid_tra, "NO2_column", "NO2_column",
	278	. "10^15 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
	279	. "ave(X)", zsto,zout)
	280	CALL histdef(nid_tra, "O3_ste", "O3_ste",
	281	. "CM-2 S-1", iim,jjm+1,nhori, 1,1,1, -99, 32,
	282	. "ave(X)", zsto,zout)
	283	CALL histdef(nid_tra, "O3_prod", "O3_prod", "CM-3 S-1",
	284	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	285	. "ave(X)", zsto,zout)
	286	CALL histdef(nid_tra, "O3_loss", "O3_loss", "CM-3 S-1",
	287	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	288	. "ave(X)", zsto,zout)
	289
	290	! Special variables for daytime averaging
	291	! CALL histdef(nid_tra, "day_cnt", "day_cnt", "-",
	292	! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	293	! . "t_sum(X)", zsto,zout)
	294	! CALL histdef(nid_tra, "NO_day", "NO_day", "VMR",
	295	! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	296	! . "t_sum(X)", zsto,zout)
	297
[616]	298
[646]	299
[624]	300	!! Ajout Anne
	301	#ifdef INCA_AER
	302	! for sulfur cycle
	303	CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",
	304	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	305	. "ave(X)", zsto,zout)
	306	CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",
	307	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	308	. "ave(X)", zsto,zout)
	309	CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",
	310	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	311	. "ave(X)", zsto,zout)
	312	CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",
	313	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	314	. "ave(X)", zsto,zout)
	315	CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",
	316	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	317	. "ave(X)", zsto,zout)
	318	CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",
	319	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	320	. "ave(X)", zsto,zout)
	321	CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",
	322	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	323	. "ave(X)", zsto,zout)
	324	CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",
	325	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	326	. "ave(X)", zsto,zout)
	327	CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",
	328	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	329	. "ave(X)", zsto,zout)
	330
	331	CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",
	332	. "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	333	. "ave(X)", zsto,zout)
	334	CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",
	335	. "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	336	. "ave(X)", zsto,zout)
	337	CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",
	338	. "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	339	. "ave(X)", zsto,zout)
	340	CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",
	341	. "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	342	. "ave(X)", zsto,zout)
	343
	344	CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",
	345	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	346	. "ave(X)", zsto,zout)
	347	#endif
	348
	349	#ifdef INCA_NMHC
	350	CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
	351	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	352	. "ave(X)", zsto,zout)
	353	CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",
	354	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	355	. "ave(X)", zsto,zout)
	356	CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",
	357	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	358	. "ave(X)", zsto,zout)
	359	CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",
	360	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	361	. "ave(X)", zsto,zout)
	362	CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",
	363	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	364	. "ave(X)", zsto,zout)
	365	CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",
	366	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	367	. "ave(X)", zsto,zout)
	368	CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",
	369	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	370	. "ave(X)", zsto,zout)
	371	#endif
	372	! Fin ajout Anne
[524]	373	#else
	374	CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
	375	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	376	. "inst(X)", zsto,zout)
	377
	378	CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
	379	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	380	. "inst(X)", zsto,zout)
	381	CALL histdef(nid_tra, "psrf1", "nature sol", "-",
	382	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	383	. "inst(X)", zsto,zout)
	384	CALL histdef(nid_tra, "psrf2", "nature sol", "-",
	385	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	386	. "inst(X)", zsto,zout)
	387	CALL histdef(nid_tra, "psrf3", "nature sol", "-",
	388	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	389	. "inst(X)", zsto,zout)
	390	CALL histdef(nid_tra, "psrf4", "nature sol", "-",
	391	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	392	. "inst(X)", zsto,zout)
	393	CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
	394	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	395	. "inst(X)", zsto,zout)
	396	CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
	397	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	398	. "inst(X)", zsto,zout)
	399	CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
	400	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	401	. "inst", zsto,zout)
	402	CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
	403	. iim,jjm+1,nhori, 1,1,1, -99, 32,
	404	. "inst(X)", zsto,zout)
	405	CALL histdef(nid_tra, "pplay", "flux u mont","-",
	406	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	407	. "inst(X)", zsto,zout)
	408	CALL histdef(nid_tra, "t", "flux u mont","-",
	409	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	410	. "inst(X)", zsto,zout)
	411	CALL histdef(nid_tra, "mfu", "flux u mont","-",
	412	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	413	. "ave(X)", zsto,zout)
	414	CALL histdef(nid_tra, "mfd", "flux u decen","-",
	415	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	416	. "ave(X)", zsto,zout)
	417	CALL histdef(nid_tra, "en_u", "flux u mont","-",
	418	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	419	. "ave(X)", zsto,zout)
	420	CALL histdef(nid_tra, "en_d", "flux u mont","-",
	421	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	422	. "ave(X)", zsto,zout)
	423	CALL histdef(nid_tra, "de_d", "flux u mont","-",
	424	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	425	. "ave(X)", zsto,zout)
	426	CALL histdef(nid_tra, "de_u", "flux u decen","-",
	427	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	428	. "ave(X)", zsto,zout)
	429	CALL histdef(nid_tra, "coefh", "turbulent coef","-",
	430	. iim,jjm+1,nhori, klev,1,klev,nvert, 32,
	431	. "ave(X)", zsto,zout)
	432
	433	#endif
	434	c
	435	CALL histend(nid_tra)
	436	#ifdef INCA_AER
	437	CALL histend(nid_tra2)
	438	CALL histend(nid_tra3)
	439	#endif
	440	ndex2d = 0
	441	ndex3d = 0
	442	ndex = 0

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

Download in other formats: