Context Navigation

ini_histday.h @ 3556

Last change on this file since 3556 was 1530, checked in by emillour, 9 years ago

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

File size: 9.2 KB

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

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

Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h @ 3556

Download in other formats: