Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h @ 894

Last change on this file since 894 was 892, checked in by slebonnois, 12 years ago
SL: Important commit ! Adaptation of Venus physics to parallel computation / template for arch on the LMD servers using ifort / documentation for 1D column physics and for parallel computations
File size: 9.2 KB

Rev	Line
[3]	1	!
	2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histmth.h,v 1.4 2005/03/09 12:30:16 fairhead Exp $
	3	!
	4	IF (ok_mensuel) THEN
[892]	5
[3]	6	zsto = dtime
	7	zout = dtime * ecrit_mth
[892]	8
[3]	9	c zsto1: pour des valeurs "instantanees" mensuelles
	10	zsto1 = dtime * ecrit_mth
	11	c zsto2: pour des flux radiatifs calcules avec un autre pas de temps
	12	zsto2 = dtime * radpas
	13	PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout
[892]	14
[3]	15	idayref = day_ref
[97]	16	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
[892]	17
	18	call histbeg_phy("histmth.nc",itau_phy,
	19	. zjulian,zsto,nhori,nid_mth)
	20
	21	!$OMP MASTER
	22	CALL histvert(nid_mth, "presnivs", "Vertical levels", "Pa",
[97]	23	. klev, presnivs, nvert)
[3]	24
	25	c-------------------------------------------------------
	26	IF(lev_histmth.GE.1) THEN
	27	c
	28	ccccccccccccc 2D fields, basics
	29	c
	30	CALL histdef(nid_mth, "phis", "Surface geop. height", "-",
[892]	31	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	32	. "once", zsto,zout)
	33	c
	34	CALL histdef(nid_mth, "aire", "Grid area", "-",
[892]	35	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	36	. "once", zsto,zout)
	37	c
	38	CALL histdef(nid_mth, "tsol", "Surface Temperature", "K",
[892]	39	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	40	. "ave(X)", zsto,zout)
	41	c
	42	CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa",
[892]	43	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	44	. "ave(X)", zsto,zout)
	45	c
	46	c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-",
[892]	47	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	48	c . "ave(X)", zsto,zout)
	49	c
	50	c CALL histdef(nid_mth, "ve", "Merid energy transport", "-",
[892]	51	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	52	c . "ave(X)", zsto,zout)
	53	c
	54	CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-",
[892]	55	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	56	. "ave(X)", zsto,zout)
	57	c
	58	CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-",
[892]	59	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	60	. "ave(X)", zsto,zout)
	61	c
	62	ENDIF !lev_histmth.GE.1
	63	c
	64	c-------------------------------------------------------
	65	IF(lev_histmth.GE.2) THEN
	66	c
	67	ccccccccccccc 3D fields, basics
	68	c
	69	CALL histdef(nid_mth, "temp", "Air temperature", "K",
[892]	70	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	71	. "ave(X)", zsto,zout)
	72	c
	73	CALL histdef(nid_mth, "pres", "Air pressure", "Pa",
[892]	74	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	75	. "ave(X)", zsto,zout)
	76	c
	77	CALL histdef(nid_mth, "geop", "Geopotential height", "m",
[892]	78	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	79	. "ave(X)", zsto,zout)
	80	c
	81	CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",
[892]	82	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	83	. "ave(X)", zsto,zout)
	84	c
	85	CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",
[892]	86	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	87	. "ave(X)", zsto,zout)
	88	c
	89	CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s",
[892]	90	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	91	. "ave(X)", zsto,zout)
	92	c
	93	CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2",
[892]	94	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	95	. "ave(X)", zsto,zout)
	96	c
	97	CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2",
[892]	98	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	99	. "ave(X)", zsto,zout)
	100	c
	101	c CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s",
[892]	102	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	103	c . "ave(X)", zsto,zout)
	104	c
	105	CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s",
[892]	106	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	107	. "ave(X)", zsto,zout)
	108	c
	109	c plusieurs traceurs
	110	if (iflag_trac.eq.1) THEN
	111	DO iq=1,nqmax
	112	IF (iq.LE.99) THEN
	113	WRITE(str2,'(i2.2)') iq
[101]	114	CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
[892]	115	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	116	. "ave(X)", zsto,zout)
	117	ELSE
	118	PRINT*, "Trop de traceurs"
	119	CALL abort
	120	ENDIF
	121	ENDDO
	122	endif
	123	c
	124	CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
[892]	125	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	126	. "ave(X)", zsto2,zout)
	127	c
	128	ENDIF !lev_histmth.GE.2
	129	c
	130	c-------------------------------------------------------
	131	IF(lev_histmth.GE.3) THEN
	132	c
	133	cccccccccccccccccc Radiative transfer
	134	c
	135	c 2D
	136	c
	137	CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2",
[892]	138	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	139	. "ave(X)", zsto2,zout)
	140	c
	141	CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2",
[892]	142	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	143	. "ave(X)", zsto2,zout)
	144	c
	145	CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2",
[892]	146	. iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	147	. "ave(X)", zsto2,zout)
	148	c
	149	c 3D
	150	c
	151	CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2",
[892]	152	. iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	153	. 32, "ave(X)", zsto2,zout)
	154	c
	155	CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2",
[892]	156	. iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	157	. 32, "ave(X)", zsto2,zout)
	158	c
	159	CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2",
[892]	160	. iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	161	. 32, "ave(X)", zsto,zout)
	162	c
	163	CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2",
[892]	164	. iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	165	. 32, "ave(X)", zsto,zout)
	166	c
	167	CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2",
[892]	168	. iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	169	. 32, "ave(X)", zsto,zout)
	170	c
	171	c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2",
[892]	172	c . iim,jj_nb,nhori, klev,1,klev,nvert,
[3]	173	c . 32, "ave(X)", zsto,zout)
	174	c
	175	ENDIF !lev_histmth.GE.3
	176	c
	177	c-------------------------------------------------------
	178	IF(lev_histmth.GE.4) THEN
	179	c
	180	CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s",
[892]	181	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	182	. "ave(X)", zsto,zout)
	183	c
	184	c CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",
[892]	185	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	186	c . "ave(X)", zsto,zout)
	187	c
	188	CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s",
[892]	189	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	190	. "ave(X)", zsto,zout)
	191	c
	192	CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s",
[892]	193	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	194	. "ave(X)", zsto,zout)
	195	c
	196	CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s",
[892]	197	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	198	. "ave(X)", zsto,zout)
	199	c
	200	CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
[892]	201	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	202	. "ave(X)", zsto,zout)
	203	c
	204	c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
[892]	205	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	206	c . "ave(X)", zsto,zout)
	207	c
	208	CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2",
[892]	209	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	210	. "ave(X)", zsto,zout)
	211	c
	212	CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2",
[892]	213	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	214	. "ave(X)", zsto,zout)
	215	c
	216	CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2",
[892]	217	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	218	. "ave(X)", zsto,zout)
	219	c
	220	c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
[892]	221	c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[3]	222	c . "ave(X)", zsto,zout)
	223	c
	224	ENDIF !lev_histmth.GE.4
	225	c
	226	c-------------------------------------------------------
	227	IF(lev_histmth.GE.5) THEN
	228	c
	229	c call histdef(nid_mth, "taux",
	230	c $ "Zonal wind stress", "Pa",
[892]	231	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	232	c $ "ave(X)", zsto,zout)
	233	c
	234	c call histdef(nid_mth, "tauy",
	235	c $ "Meridional xind stress", "Pa",
[892]	236	c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	237	c $ "ave(X)", zsto,zout)
	238	c
	239	c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-",
[892]	240	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	241	c . "ave(X)", zsto,zout)
	242	c
	243	c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-",
[892]	244	c . iim,jj_nb,nhori, 1,1,1, nvert, 32,
[3]	245	c . "ave(X)", zsto,zout)
	246	c
	247	ENDIF !lev_histmth.GE.5
	248	c-------------------------------------------------------
[892]	249
[3]	250	CALL histend(nid_mth)
[892]	251	!$OMP END MASTER
	252
[3]	253	ENDIF ! fin de test sur ok_mensuel

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

Download in other formats: