Context Navigation

source: LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h @ 704

Last change on this file since 704 was 704, checked in by Laurent Fairhead, 18 years ago
Inclusion des modifs de Y. Meurdesoif pour la version V3 LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 21.4 KB

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

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

Download in other formats: