Context Navigation

source: trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h @ 1055

Last change on this file since 1055 was 902, checked in by slebonnois, 12 years ago
SL: small modifs on Venus outputs
File size: 13.5 KB

Rev	Line
[3]	1	IF (ok_journe) THEN
	2	c
	3	zsto = dtime
	4	zout = dtime * FLOAT(ecrit_day)
	5	c zsto1: pour des flux radiatifs calcules tous les radpas appels physiq
	6	zsto1= dtime * FLOAT(radpas)
	7	c
	8	idayref = day_ref
[97]	9	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
[3]	10	c
	11	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
	12	DO i = 1, iim
[888]	13	zx_lon(i,1) = rlond(i+jjmp1-jjm)
	14	zx_lon(i,jjmp1) = rlond(i+jjmp1-jjm)
[3]	15	ENDDO
	16	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)
[97]	17	CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
[902]	18	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
[3]	19	. nhori, nid_day)
	20	write(,)'Journee ', itau_phy, zjulian
	21	CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
[97]	22	. klev, presnivs, nvert)
[3]	23
	24	c-------------------------------------------------------
	25	IF(lev_histday.GE.1) THEN
	26	c
	27	ccccccccccccc 2D fields, invariables
	28	c
	29	CALL histdef(nid_day, "phis", "Surface geop. height", "-",
	30	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	31	. "once", zsto,zout)
	32	c
	33	CALL histdef(nid_day, "aire", "Grid area", "-",
	34	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	35	. "once", zsto,zout)
	36	c
	37	ccccccc axe Ls
	38	CALL histdef(nid_day, "ls", "Solar longitude", "degrees",
	39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	40	. "ave(X)", zsto,zout)
	41	c
	42	ccccccccccccc 2D fields, variables
	43	c
	44	CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
	45	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	46	. "ave(X)", zsto,zout)
	47	c
	48	CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
	49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	50	. "ave(X)", zsto,zout)
	51	c
	52	c CALL histdef(nid_day, "ue", "Zonal energy transport", "-",
	53	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	54	c . "ave(X)", zsto,zout)
	55	c
	56	c CALL histdef(nid_day, "ve", "Merid energy transport", "-",
	57	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	58	c . "ave(X)", zsto,zout)
	59	c
	60	ENDIF !lev_histday.GE.1
	61	c
	62	c-------------------------------------------------------
	63	IF(lev_histday.GE.2) THEN
	64	c
	65	ccccccccccccc 3D fields, basics
	66	c
	67	CALL histdef(nid_day, "temp", "Air temperature", "K",
	68	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	69	. "ave(X)", zsto,zout)
	70	c
	71	CALL histdef(nid_day, "pres", "Air pressure", "Pa",
	72	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	73	. "ave(X)", zsto,zout)
	74	c
	75	CALL histdef(nid_day, "geop", "Geopotential height", "m",
	76	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	77	. "ave(X)", zsto,zout)
	78	c
	79	CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
	80	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	81	. "ave(X)", zsto,zout)
	82	c
	83	CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
	84	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	85	. "ave(X)", zsto,zout)
	86	c
	87	CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
	88	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	89	. "ave(X)", zsto,zout)
	90	c
	91	CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	93	. "ave(X)", zsto1,zout)
	94	c
	95	CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2",
	96	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	97	. "ave(X)", zsto,zout)
	98	c
	99	CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2",
	100	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	101	. "ave(X)", zsto,zout)
	102	c
[106]	103	cccccccccccccccccc Tracers
	104	c
[3]	105	if (iflag_trac.eq.1) THEN
[175]	106	if (microfi.ge.1) then
	107	c DO iq=1,nmicro
	108	c CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2",
	109	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	110	c . "ave(X)", zsto,zout)
	111	c ENDDO
	112	CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2",
	113	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	114	. "ave(X)", zsto,zout)
[306]	115
[474]	116	if (clouds.eq.1) then
[175]	117	CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
	118	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	119	. "ave(X)", zsto,zout)
	120	CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
	121	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	122	. "ave(X)", zsto,zout)
	123	CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
	124	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	125	. "ave(X)", zsto,zout)
	126	CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
	127	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	128	. "ave(X)", zsto,zout)
	129	c--------------
	130	c ----- SATURATION ESP NUAGES
	131	CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
[3]	132	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	133	. "ave(X)", zsto,zout)
[175]	134	CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
	135	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	136	. "ave(X)", zsto,zout)
	137	CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
	138	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	139	. "ave(X)", zsto,zout)
	140	c --------------
	141	c ----- RESERVOIR DE SURFACE
	142	CALL histdef(nid_day, "reserv", "Reservoir surface","m",
	143	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	144	. "ave(X)", zsto,zout)
	145	c --------------
[474]	146	c ----- ECHANGE GAZ SURF/ATM (evaporation)
	147	CALL histdef(nid_day, "evapch4", "Evaporation CH4","m",
	148	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	149	. "ave(X)", zsto,zout)
	150	c --------------
[175]	151	c ----- PRECIPITATIONS (precipitations cumulatives)
	152	CALL histdef(nid_day,"prech4","Precip CH4","m",
	153	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	154	. "ave(X)", zsto,zout)
	155	CALL histdef(nid_day,"prec2h6","Precip C2H6",
	156	. "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	. "ave(X)", zsto,zout)
	158	CALL histdef(nid_day,"prec2h2","Precip C2H2",
	159	. "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
	160	. "ave(X)", zsto,zout)
	161	c --------------
	162	c ----- FLUX GLACE
	163	CALL histdef(nid_day,"flxgl1", "flux gl CH4",
	164	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	165	. "ave(X)", zsto,zout)
	166	CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
	167	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	168	. "ave(X)", zsto,zout)
	169	CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
	170	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	171	. "ave(X)", zsto,zout)
	172	c --------------
	173	c ----- RAYON DES GOUTTES
	174	CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
	175	. "m",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	176	. "ave(X)", zsto,zout)
[474]	177	endif
[3]	178	endif
[175]	179	c --------------
	180	c ----- TRACEURS CHIMIQUES
[3]	181	if (nmicro.lt.nqmax) then
	182	DO iq=nmicro+1,nqmax
	183	CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
	184	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	185	. "ave(X)", zsto,zout)
	186	ENDDO
	187	endif
	188	endif
	189	c
[175]	190	ENDIF !lev_histday.GE.2
	191	c
	192	c-------------------------------------------------------
	193	IF(lev_histday.GE.3) THEN
	194	c
[3]	195	cccccccccccccccccc Radiative transfer
	196	c
	197	c 2D
	198	c
	199	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
	200	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	201	. "ave(X)", zsto1,zout)
	202	c
	203	CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
	204	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	205	. "ave(X)", zsto1,zout)
	206	c
	207	CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
	208	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	209	. "ave(X)", zsto1,zout)
	210	c
	211	c 3D
	212	c
	213	CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
	214	. iim,jjmp1,nhori, klev,1,klev,nvert,
	215	. 32, "ave(X)", zsto1,zout)
	216	c
	217	CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
	218	. iim,jjmp1,nhori, klev,1,klev,nvert,
	219	. 32, "ave(X)", zsto1,zout)
	220	c
[175]	221	c --------------
	222	c ----- OPACITE BRUME
[106]	223	DO k=7,NSPECV,10
	224	write(str1,'(i2.2)') k
	225	CALL histdef(nid_day,"thv"//str1,"Haze Opa Vis",
	226	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	227	. "ave(X)",zsto1,zout)
	228	ENDDO
	229	c
[175]	230	DO k=8,NSPECI,10
[106]	231	write(str1,'(i2.2)') k
[175]	232	CALL histdef(nid_day,"thi"//str1,"Haze Opa IR",
	233	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	234	. "ave(X)",zsto1,zout)
	235	ENDDO
	236	c
[175]	237	c --------------
	238	c ----- EXTINCTION BRUME
[106]	239	DO k=7,NSPECV,10
	240	write(str1,'(i2.2)') k
[175]	241	CALL histdef(nid_day,"khv"//str1,"Haze ext Vis ",
	242	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	243	. "ave(X)",zsto1,zout)
	244	ENDDO
	245	c
[175]	246	DO k=8,NSPECI,10
[106]	247	write(str1,'(i2.2)') k
[175]	248	CALL histdef(nid_day,"khi"//str1,"Haze ext IR ",
[106]	249	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	250	. "ave(X)",zsto1,zout)
	251	ENDDO
	252	c
[175]	253	c --------------
	254	c ----- OPACITE GAZ
	255	DO k=7,NSPECV,10
[106]	256	write(str1,'(i2.2)') k
[175]	257	CALL histdef(nid_day,"tgv"//str1,"Gas Opa Vis",
[106]	258	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	259	. "ave(X)",zsto1,zout)
	260	ENDDO
	261	c
	262	DO k=8,NSPECI,10
	263	write(str1,'(i2.2)') k
[175]	264	CALL histdef(nid_day,"tgi"//str1,"Gas Opa IR",
	265	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	266	. "ave(X)",zsto1,zout)
	267	ENDDO
	268	c
[175]	269	c --------------
	270	c ----- EXTINCTION GAZ
	271	DO k=7,NSPECV,10
[106]	272	write(str1,'(i2.2)') k
[175]	273	CALL histdef(nid_day,"kgv"//str1,"Gas ext Vis ",
	274	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	275	. "ave(X)",zsto1,zout)
	276	ENDDO
	277	c
	278	DO k=8,NSPECI,10
	279	write(str1,'(i2.2)') k
[175]	280	CALL histdef(nid_day,"kgi"//str1,"Gas ext IR ",
[106]	281	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	282	. "ave(X)",zsto1,zout)
	283	ENDDO
	284	c
[175]	285	c --------------
	286	c ----- OPACITE NUAGES
	287	if (clouds.eq.1) then
	288	CALL histdef(nid_day,"tcld","Cld Opa proxy",
	289	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	290	. "ave(X)",zsto,zout)
	291	c
	292	c --------------
	293	c ----- EXTINCTION NUAGES
	294	CALL histdef(nid_day,"kcld","Cld Ext proxy",
	295	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	296	. "ave(X)",zsto,zout)
	297	endif
	298	c
[3]	299	ENDIF !lev_histday.GE.3
	300	c
	301	c-------------------------------------------------------
	302	IF(lev_histday.GE.4) THEN
	303	c
	304	CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
	305	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	306	. "ave(X)", zsto,zout)
	307	c
	308	CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
	309	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	310	. "ave(X)", zsto,zout)
	311	c
	312	CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
	313	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	314	. "ave(X)", zsto,zout)
	315	c
	316	CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
	317	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	318	. "ave(X)", zsto,zout)
	319	c
	320	CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
	321	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	322	. "ave(X)", zsto,zout)
	323	c
	324	CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
	325	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	326	. "ave(X)", zsto,zout)
	327	c
[106]	328	c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
	329	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	330	c . "ave(X)", zsto,zout)
[3]	331	c
	332	ENDIF !lev_histday.GE.4
	333	c
	334	c-------------------------------------------------------
	335	IF(lev_histday.GE.5) THEN
	336	c
	337	c
	338	c call histdef(nid_day, "taux",
	339	c $ "Zonal wind stress", "Pa",
	340	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	341	c $ "ave(X)", zsto,zout)
	342	c
	343	c call histdef(nid_day, "tauy",
	344	c $ "Meridional xind stress", "Pa",
	345	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	346	c $ "ave(X)", zsto,zout)
	347	c
	348	c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
	349	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	350	c . "ave(X)", zsto,zout)
	351	c
	352	c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
	353	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	354	c . "ave(X)", zsto,zout)
	355	c
	356	ENDIF !lev_histday.GE.5
	357	c-------------------------------------------------------
	358	c
	359	CALL histend(nid_day)
	360	c
	361	ndex2d = 0
	362	ndex3d = 0
	363	c
	364	ENDIF ! fin de test sur ok_journe

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

Download in other formats: