Context Navigation

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

Last change on this file since 201 was 175, checked in by slebonnois, 14 years ago

S.LEBONNOIS:

Revision majeure de la physique Titan => ajout des nuages version 10 bins (Jeremie Burgalat) Cette version reste a tester mais avec clouds=0, on reste sur l'ancienne.
Quelques ajouts dans la doc.

File size: 13.2 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
	13	zx_lon(i,1) = rlond(i+1)
	14	zx_lon(i,jjmp1) = rlond(i+1)
	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,:),
	18	. 1,iim,1,jjmp1, itau_phy, zjulian, zsto,
[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)
	115	CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
	116	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	117	. "ave(X)", zsto,zout)
	118	CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2",
	119	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	120	. "ave(X)", zsto,zout)
	121	CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2",
	122	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	123	. "ave(X)", zsto,zout)
	124	CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2",
	125	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	126	. "ave(X)", zsto,zout)
	127	c--------------
	128	c ----- SATURATION ESP NUAGES
	129	if (clouds.eq.1) then
	130	CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
[3]	131	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	132	. "ave(X)", zsto,zout)
[175]	133	CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--",
	134	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	135	. "ave(X)", zsto,zout)
	136	CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--",
	137	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	138	. "ave(X)", zsto,zout)
	139	c --------------
	140	c ----- RESERVOIR DE SURFACE
	141	CALL histdef(nid_day, "reserv", "Reservoir surface","m",
	142	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	143	. "ave(X)", zsto,zout)
	144	c --------------
	145	c ----- PRECIPITATIONS (precipitations cumulatives)
	146	CALL histdef(nid_day,"prech4","Precip CH4","m",
	147	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	148	. "ave(X)", zsto,zout)
	149	CALL histdef(nid_day,"prec2h6","Precip C2H6",
	150	. "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
	151	. "ave(X)", zsto,zout)
	152	CALL histdef(nid_day,"prec2h2","Precip C2H2",
	153	. "m",iim,jjmp1,nhori, 1,1,1, -99, 32,
	154	. "ave(X)", zsto,zout)
	155	c --------------
	156	c ----- FLUX GLACE
	157	CALL histdef(nid_day,"flxgl1", "flux gl CH4",
	158	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	159	. "ave(X)", zsto,zout)
	160	CALL histdef(nid_day,"flxgl2", "flux gl C2H6",
	161	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	162	. "ave(X)", zsto,zout)
	163	CALL histdef(nid_day,"flxgl3", "flux gl C2H2",
	164	. "kg/m2/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	165	. "ave(X)", zsto,zout)
	166	c --------------
	167	c ----- RAYON DES GOUTTES
	168	CALL histdef(nid_day,"rcldbar", "rayon moyen goutte",
	169	. "m",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	170	. "ave(X)", zsto,zout)
	171	endif
[3]	172	endif
[175]	173	c --------------
	174	c ----- TRACEURS CHIMIQUES
[3]	175	if (nmicro.lt.nqmax) then
	176	DO iq=nmicro+1,nqmax
	177	CALL histdef(nid_day, tname(iq), ttext(iq), "ppm",
	178	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	179	. "ave(X)", zsto,zout)
	180	ENDDO
	181	endif
	182	endif
	183	c
[175]	184	ENDIF !lev_histday.GE.2
	185	c
	186	c-------------------------------------------------------
	187	IF(lev_histday.GE.3) THEN
	188	c
[3]	189	cccccccccccccccccc Radiative transfer
	190	c
	191	c 2D
	192	c
	193	CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
	194	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	195	. "ave(X)", zsto1,zout)
	196	c
	197	CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2",
	198	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	199	. "ave(X)", zsto1,zout)
	200	c
	201	CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2",
	202	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	203	. "ave(X)", zsto1,zout)
	204	c
	205	c 3D
	206	c
	207	CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2",
	208	. iim,jjmp1,nhori, klev,1,klev,nvert,
	209	. 32, "ave(X)", zsto1,zout)
	210	c
	211	CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2",
	212	. iim,jjmp1,nhori, klev,1,klev,nvert,
	213	. 32, "ave(X)", zsto1,zout)
	214	c
[175]	215	c --------------
	216	c ----- OPACITE BRUME
[106]	217	DO k=7,NSPECV,10
	218	write(str1,'(i2.2)') k
	219	CALL histdef(nid_day,"thv"//str1,"Haze Opa Vis",
	220	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	221	. "ave(X)",zsto1,zout)
	222	ENDDO
	223	c
[175]	224	DO k=8,NSPECI,10
[106]	225	write(str1,'(i2.2)') k
[175]	226	CALL histdef(nid_day,"thi"//str1,"Haze Opa IR",
	227	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	228	. "ave(X)",zsto1,zout)
	229	ENDDO
	230	c
[175]	231	c --------------
	232	c ----- EXTINCTION BRUME
[106]	233	DO k=7,NSPECV,10
	234	write(str1,'(i2.2)') k
[175]	235	CALL histdef(nid_day,"khv"//str1,"Haze ext Vis ",
	236	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	237	. "ave(X)",zsto1,zout)
	238	ENDDO
	239	c
[175]	240	DO k=8,NSPECI,10
[106]	241	write(str1,'(i2.2)') k
[175]	242	CALL histdef(nid_day,"khi"//str1,"Haze ext IR ",
[106]	243	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	244	. "ave(X)",zsto1,zout)
	245	ENDDO
	246	c
[175]	247	c --------------
	248	c ----- OPACITE GAZ
	249	DO k=7,NSPECV,10
[106]	250	write(str1,'(i2.2)') k
[175]	251	CALL histdef(nid_day,"tgv"//str1,"Gas Opa Vis",
[106]	252	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	253	. "ave(X)",zsto1,zout)
	254	ENDDO
	255	c
	256	DO k=8,NSPECI,10
	257	write(str1,'(i2.2)') k
[175]	258	CALL histdef(nid_day,"tgi"//str1,"Gas Opa IR",
	259	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	260	. "ave(X)",zsto1,zout)
	261	ENDDO
	262	c
[175]	263	c --------------
	264	c ----- EXTINCTION GAZ
	265	DO k=7,NSPECV,10
[106]	266	write(str1,'(i2.2)') k
[175]	267	CALL histdef(nid_day,"kgv"//str1,"Gas ext Vis ",
	268	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
[106]	269	. "ave(X)",zsto1,zout)
	270	ENDDO
	271	c
	272	DO k=8,NSPECI,10
	273	write(str1,'(i2.2)') k
[175]	274	CALL histdef(nid_day,"kgi"//str1,"Gas ext IR ",
[106]	275	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	276	. "ave(X)",zsto1,zout)
	277	ENDDO
	278	c
[175]	279	c --------------
	280	c ----- OPACITE NUAGES
	281	if (clouds.eq.1) then
	282	CALL histdef(nid_day,"tcld","Cld Opa proxy",
	283	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	284	. "ave(X)",zsto,zout)
	285	c
	286	c --------------
	287	c ----- EXTINCTION NUAGES
	288	CALL histdef(nid_day,"kcld","Cld Ext proxy",
	289	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	290	. "ave(X)",zsto,zout)
	291	endif
	292	c
[3]	293	ENDIF !lev_histday.GE.3
	294	c
	295	c-------------------------------------------------------
	296	IF(lev_histday.GE.4) THEN
	297	c
	298	CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s",
	299	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	300	. "ave(X)", zsto,zout)
	301	c
	302	CALL histdef(nid_day, "dtphy", "Physics dT", "K/s",
	303	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	304	. "ave(X)", zsto,zout)
	305	c
	306	CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s",
	307	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	308	. "ave(X)", zsto,zout)
	309	c
	310	CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s",
	311	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	312	. "ave(X)", zsto,zout)
	313	c
	314	CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s",
	315	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	316	. "ave(X)", zsto,zout)
	317	c
	318	CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s",
	319	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	320	. "ave(X)", zsto,zout)
	321	c
[106]	322	c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s",
	323	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	324	c . "ave(X)", zsto,zout)
[3]	325	c
	326	ENDIF !lev_histday.GE.4
	327	c
	328	c-------------------------------------------------------
	329	IF(lev_histday.GE.5) THEN
	330	c
	331	c
	332	c call histdef(nid_day, "taux",
	333	c $ "Zonal wind stress", "Pa",
	334	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	335	c $ "ave(X)", zsto,zout)
	336	c
	337	c call histdef(nid_day, "tauy",
	338	c $ "Meridional xind stress", "Pa",
	339	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	340	c $ "ave(X)", zsto,zout)
	341	c
	342	c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-",
	343	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	344	c . "ave(X)", zsto,zout)
	345	c
	346	c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-",
	347	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	348	c . "ave(X)", zsto,zout)
	349	c
	350	ENDIF !lev_histday.GE.5
	351	c-------------------------------------------------------
	352	c
	353	CALL histend(nid_day)
	354	c
	355	ndex2d = 0
	356	ndex3d = 0
	357	c
	358	ENDIF ! fin de test sur ok_journe

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

Download in other formats: