Context Navigation

source: trunk/libf/phytitan/ini_histmth.h @ 119

Last change on this file since 119 was 106, checked in by slebonnois, 14 years ago
SL : mise a jour de phytitan pour etre conforme aux sources actuelles utilisees sur gnome.
File size: 10.5 KB

Rev	Line
[3]	1	IF (ok_mensuel) THEN
	2	c
	3	zsto = dtime
	4	zout = dtime * FLOAT(ecrit_mth)
	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("histmth", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	18	. 1,iim,1,jjmp1, itau_phy, zjulian, zsto,
[3]	19	. nhori, nid_mth)
	20	write(,)'Journee ', itau_phy, zjulian
	21	CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb",
[97]	22	. klev, presnivs, nvert)
[3]	23
	24	c-------------------------------------------------------
	25	IF(lev_histmth.GE.1) THEN
	26	c
	27	ccccccccccccc 2D fields, invariables
	28	c
	29	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
	30	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	31	. "once", zsto,zout)
	32	c
	33	CALL histdef(nid_mth, "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_mth, "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_mth, "tsol", "Surface Temperature", "K",
	45	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	46	. "ave(X)", zsto,zout)
	47	c
	48	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
	49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	50	. "ave(X)", zsto,zout)
	51	c
[106]	52	c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
	53	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	54	c . "ave(X)", zsto,zout)
[3]	55	c
[106]	56	c CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
	57	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	58	c . "ave(X)", zsto,zout)
[3]	59	c
[106]	60	c CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
	61	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	62	c . "ave(X)", zsto,zout)
	63	c
	64	c CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
	65	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	66	c . "ave(X)", zsto,zout)
	67	c
[3]	68	ENDIF !lev_histmth.GE.1
	69	c
	70	c-------------------------------------------------------
	71	IF(lev_histmth.GE.2) THEN
	72	c
	73	ccccccccccccc 3D fields, basics
	74	c
	75	CALL histdef(nid_mth, "temp", "Air temperature", "K",
	76	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	77	. "ave(X)", zsto,zout)
	78	c
	79	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
	80	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	81	. "ave(X)", zsto,zout)
	82	c
	83	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
	84	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	85	. "ave(X)", zsto,zout)
	86	c
	87	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
	88	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	89	. "ave(X)", zsto,zout)
	90	c
	91	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
	92	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	93	. "ave(X)", zsto,zout)
	94	c
	95	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
	96	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	97	. "ave(X)", zsto,zout)
	98	c
[106]	99	c CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
	100	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	101	c . "ave(X)", zsto,zout)
	102	c
[3]	103	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
	104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	105	. "ave(X)", zsto1,zout)
	106	c
	107	if (iflag_trac.eq.1) THEN
[97]	108	c if (microfi.eq.1) then
	109	c DO iq=1,nmicro
	110	c CALL histdef(nid_mth, tname(iq), ttext(iq), "n/m2",
	111	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	112	c . "ave(X)", zsto,zout)
	113	c ENDDO
	114	c endif
	115	if (nmicro.lt.nqmax) then
[3]	116	DO iq=nmicro+1,nqmax
	117	CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
	118	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	119	. "ave(X)", zsto,zout)
	120	ENDDO
	121	c Condensation:
	122	c DO iq=nmicro+1,nqmax
	123	c CALL histdef(nid_mth, "c_"//tname(iq), "c_"//ttext(iq),
	124	c . "ppm/s",iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	125	c . "ave(X)", zsto,zout)
	126	c ENDDO
[97]	127	endif
[3]	128	endif
	129	c
	130	ENDIF !lev_histmth.GE.2
	131	c
	132	c-------------------------------------------------------
	133	IF(lev_histmth.GE.3) THEN
	134	c
	135	cccccccccccccccccc Radiative transfer
	136	c
	137	c 2D
	138	c
	139	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
	140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	141	. "ave(X)", zsto1,zout)
	142	c
	143	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
	144	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	145	. "ave(X)", zsto1,zout)
	146	c
	147	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
	148	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	149	. "ave(X)", zsto1,zout)
	150	c
	151	c 3D
	152	c
	153	CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
	154	. iim,jjmp1,nhori, klev,1,klev,nvert,
	155	. 32, "ave(X)", zsto1,zout)
	156	c
	157	CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
	158	. iim,jjmp1,nhori, klev,1,klev,nvert,
	159	. 32, "ave(X)", zsto1,zout)
	160	c
[106]	161	DO k=7,NSPECV,10
[97]	162	write(str1,'(i2.2)') k
	163	CALL histdef(nid_mth,"thv"//str1,"Haze Opa Vis",
	164	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	165	. "ave(X)",zsto1,zout)
	166	ENDDO
	167	c
[106]	168	DO k=7,NSPECV,10
[97]	169	write(str1,'(i2.2)') k
	170	CALL histdef(nid_mth,"khv"//str1,"Haze ext Vis ",
	171	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	172	. "ave(X)",zsto1,zout)
	173	ENDDO
	174	c
[106]	175	DO k=7,NSPECV,10
	176	write(str1,'(i2.2)') k
	177	CALL histdef(nid_mth,"tgv"//str1,"Haze Opa Vis",
	178	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	179	. "ave(X)",zsto1,zout)
	180	ENDDO
	181	c
	182	DO k=7,NSPECV,10
	183	write(str1,'(i2.2)') k
	184	CALL histdef(nid_mth,"kgv"//str1,"Haze ext Vis ",
	185	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	186	. "ave(X)",zsto1,zout)
	187	ENDDO
	188	c
	189	DO k=8,NSPECI,10
	190	write(str1,'(i2.2)') k
	191	CALL histdef(nid_mth,"thi"//str1,"Haze Opa IR",
	192	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	193	. "ave(X)",zsto1,zout)
	194	ENDDO
	195	c
	196	DO k=8,NSPECI,10
	197	write(str1,'(i2.2)') k
	198	CALL histdef(nid_mth,"khi"//str1,"Haze ext IR ",
	199	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	200	. "ave(X)",zsto1,zout)
	201	ENDDO
	202	c
	203	DO k=8,NSPECI,10
	204	write(str1,'(i2.2)') k
	205	CALL histdef(nid_mth,"tgi"//str1,"Haze Opa IR",
	206	. "--",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	207	. "ave(X)",zsto1,zout)
	208	ENDDO
	209	c
	210	DO k=8,NSPECI,10
	211	write(str1,'(i2.2)') k
	212	CALL histdef(nid_mth,"kgi"//str1,"Haze ext IR ",
	213	. "m-1",iim,jjmp1,nhori,klev,1,klev,nvert,32,
	214	. "ave(X)",zsto1,zout)
	215	ENDDO
	216	c
[3]	217	ENDIF !lev_histmth.GE.3
	218	c
	219	c-------------------------------------------------------
	220	IF(lev_histmth.GE.4) THEN
	221	c
	222	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
	223	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	224	. "ave(X)", zsto,zout)
	225	c
	226	CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
	227	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	228	. "ave(X)", zsto,zout)
	229	c
	230	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
	231	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	232	. "ave(X)", zsto,zout)
	233	c
	234	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
	235	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	236	. "ave(X)", zsto,zout)
	237	c
	238	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
	239	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	240	. "ave(X)", zsto,zout)
	241	c
	242	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
	243	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	244	. "ave(X)", zsto,zout)
	245	c
	246	c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
	247	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	248	c . "ave(X)", zsto,zout)
	249	c
	250	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
	251	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	252	. "ave(X)", zsto,zout)
	253	c
	254	CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
	255	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	256	. "ave(X)", zsto,zout)
	257	c
	258	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
	259	c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	260	c . "ave(X)", zsto,zout)
	261	c
	262	ENDIF !lev_histmth.GE.4
	263	c
	264	c-------------------------------------------------------
	265	IF(lev_histmth.GE.5) THEN
	266	c
	267	c
	268	c call histdef(nid_mth, "taux",
	269	c $ "Zonal wind stress", "Pa",
	270	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	271	c $ "ave(X)", zsto,zout)
	272	c
	273	c call histdef(nid_mth, "tauy",
	274	c $ "Meridional xind stress", "Pa",
	275	c $ iim,jjmp1,nhori, 1,1,1, -99, 32,
	276	c $ "ave(X)", zsto,zout)
	277	c
	278	c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
	279	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	280	c . "ave(X)", zsto,zout)
	281	c
	282	c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
	283	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	284	c . "ave(X)", zsto,zout)
	285	c
	286	ENDIF !lev_histmth.GE.5
	287	c-------------------------------------------------------
	288	c
	289	CALL histend(nid_mth)
	290	c
	291	ndex2d = 0
	292	ndex3d = 0
	293	c
	294	ENDIF ! fin de test sur ok_journe

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

Download in other formats: