Context Navigation

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

Last change on this file since 1243 was 1056, checked in by slebonnois, 11 years ago
SL: Titan runs ! see DOC/chantiers/commit_importants.log
File size: 13.4 KB

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

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

Download in other formats: