Context Navigation

source: trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h @ 1543

Last change on this file since 1543 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).

EM

File size: 13.5 KB

Rev	Line
[3]	1	!
	2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histins.h,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
	3	!
	4	IF (ok_instan) THEN
[892]	5
[902]	6	zsto = dtime * ecrit_ins
	7	zout = dtime * ecrit_ins
[892]	8
[3]	9	idayref = day_ref
[97]	10	CALL ymds2ju(annee_ref, 1, idayref, zero, zjulian)
[892]	11
	12	call histbeg_phy("histins.nc",itau_phy,
[902]	13	. zjulian,dtime,nhori,nid_ins)
[892]	14
	15	!$OMP MASTER
	16	CALL histvert(nid_ins, "presnivs", "Vertical levels", "Pa",
[97]	17	. klev, presnivs, nvert)
[3]	18
	19	c-------------------------------------------------------
[1305]	20	IF(lev_histins.GE.1) THEN
[888]	21	c
	22	ccccccccccccc 2D fields, basics
	23	c
	24	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
[1530]	25	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	26	. "once", zsto,zout)
	27	c
	28	CALL histdef(nid_ins, "aire", "Grid area", "-",
[1530]	29	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	30	. "once", zsto,zout)
	31	c
	32	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
[1530]	33	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	34	. "ins(X)", zsto,zout)
	35	c
	36	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
[1530]	37	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	38	. "ins(X)", zsto,zout)
	39	c
	40	c CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
[1530]	41	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	42	c . "ins(X)", zsto,zout)
	43	c
	44	c CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
[1530]	45	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	46	c . "ins(X)", zsto,zout)
	47	c
[902]	48	c CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-",
[1530]	49	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[902]	50	c . "ins(X)", zsto,zout)
[888]	51	c
[902]	52	c CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-",
[1530]	53	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[902]	54	c . "ins(X)", zsto,zout)
[888]	55	c
[1305]	56	ENDIF !lev_histins.GE.1
[888]	57	c
[3]	58	c-------------------------------------------------------
[1305]	59	IF(lev_histins.GE.2) THEN
[888]	60	c
	61	ccccccccccccc 3D fields, basics
	62	c
	63	CALL histdef(nid_ins, "temp", "Air temperature", "K",
[1530]	64	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	65	. "ins(X)", zsto,zout)
	66	c
	67	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
[1530]	68	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	69	. "ins(X)", zsto,zout)
	70	c
	71	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
[1530]	72	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	73	. "ins(X)", zsto,zout)
	74	c
	75	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
[1530]	76	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	77	. "ins(X)", zsto,zout)
	78	c
	79	CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
[1530]	80	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	81	. "ins(X)", zsto,zout)
	82	c
	83	CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
[1530]	84	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	85	. "ins(X)", zsto,zout)
	86	c
	87	CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
[1530]	88	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	89	. "ins(X)", zsto,zout)
	90	c
	91	CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
[1530]	92	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	93	. "ins(X)", zsto,zout)
	94	c
	95	c CALL histdef(nid_ins, "mang", "Angular momentum", "kg m2/s",
[1530]	96	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	97	c . "ins(X)", zsto,zout)
	98	c
[902]	99	c CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s",
[1530]	100	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[902]	101	c . "ins(X)", zsto,zout)
[888]	102	c
[1453]	103	CALL histdef(nid_ins, "mmean", "Mean molecular mass", "g/mol",
[1530]	104	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[1468]	105	. "ins(X)", zsto,zout)
[1452]	106
[1453]	107	CALL histdef(nid_ins, "rho", "Air density [mass/Vol]", "kg/m3",
[1530]	108	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[1468]	109	. "ins(X)", zsto,zout)
[1452]	110
	111
	112	c plusieurs traceurs !!!outputs in [vmr]
	113	IF (iflag_trac.eq.1) THEN
[888]	114	DO iq=1,nqmax
	115	IF (iq.LE.99) THEN
	116	WRITE(str2,'(i2.2)') iq
[1442]	117	CALL histdef(nid_ins, tname(iq), ttext(iq), "mol/mol",
[1530]	118	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	119	. "ins(X)", zsto,zout)
	120	ELSE
	121	PRINT*, "Trop de traceurs"
	122	CALL abort
	123	ENDIF
	124	ENDDO
[1452]	125	ENDIF
	126
	127	IF (callthermos .and. ok_chem) THEN
[1453]	128	CALL histdef(nid_ins, "d_qmoldif CO2", "Dif molec" , "kg/kg",
[1530]	129	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[1468]	130	. "ins(X)", zsto,zout)
[1453]	131	CALL histdef(nid_ins, "d_qmoldif O3p", "Dif molec" , "kg/kg",
[1530]	132	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[1468]	133	. "ins(X)", zsto,zout)
[1453]	134	CALL histdef(nid_ins, "d_qmoldif N2", "Dif molec" , "kg/kg",
[1530]	135	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[1468]	136	. "ins(X)", zsto,zout)
[1452]	137	ENDIF
[888]	138	c
	139	CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
[1530]	140	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	141	. "ins(X)", zsto,zout)
	142	c
[1442]	143	if (ok_cloud) THEN
	144
	145	if (nb_mode.GE.1) THEN
	146
	147	c
	148	CALL histdef(nid_ins, "NBRTOTm1", "Nbr total droplet",
[1530]	149	. "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
[1305]	150	. "ins(X)", zsto,zout)
[1442]	151	c
	152
	153	c
	154	c CALL histdef(nid_ins, "R_MEDIANm1", "Median radius
	155	c . for log normal distribution" ,
	156	c . "fraction",
[1530]	157	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	158	c . "ins(X)", zsto,zout)
	159	c
	160
	161	c
	162	c CALL histdef(nid_ins, "STDDEVm1", "Std Deviation
	163	c . for log normal distribution",
	164	c . "fraction",
[1530]	165	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	166	c . "ins(X)", zsto,zout)
	167	c
	168
	169	if (nb_mode.GE.2) THEN
	170
	171	c
	172	CALL histdef(nid_ins, "NBRTOTm2", "Nbr total droplet",
[1530]	173	. "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
[1305]	174	. "ins(X)", zsto,zout)
[1442]	175	c
	176
	177	c
	178	c CALL histdef(nid_ins, "R_MEDIANm2", "Median radius
	179	c . for log normal distribution" ,
	180	c . "fraction",
[1530]	181	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	182	c . "ins(X)", zsto,zout)
	183	c
	184
	185	c
	186	c CALL histdef(nid_ins, "STDDEVm2", "Std Deviation
	187	c . for log normal distribution",
	188	c . "fraction",
[1530]	189	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	190	c . "ins(X)", zsto,zout)
	191	c
	192
	193	if (nb_mode.GE.3) THEN
	194
	195	c
	196	CALL histdef(nid_ins, "NBRTOTm3", "Nbr total droplet", "#/cm3",
[1530]	197	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	198	. "ins(X)", zsto,zout)
[1442]	199	c
	200
	201	c
	202	c CALL histdef(nid_ins, "R_MEDIANm3", "Median radius
	203	c . for log normal distribution" ,
	204	c . "fraction",
[1530]	205	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	206	c . "ins(X)", zsto,zout)
	207	c
	208
	209	c
	210	c CALL histdef(nid_ins, "STDDEVm3", "Std Deviation
	211	c . for log normal distribution",
	212	c . "fraction",
[1530]	213	c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1442]	214	c . "ins(X)", zsto,zout)
	215	c
	216
	217	ENDIF
	218	ENDIF
	219	ENDIF
	220
	221	c
	222	CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4",
	223	. "fraction",
[1530]	224	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	225	. "ins(X)", zsto,zout)
[1442]	226	c
	227
	228	c
	229	CALL histdef(nid_ins, "rho_droplet", "density cloud droplet",
	230	. "kg.m-3",
[1530]	231	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	232	. "ins(X)", zsto,zout)
[1442]	233	c
[1305]	234
[1442]	235	ENDIF
	236
	237	if (ok_sedim) THEN
	238	c
	239	CALL histdef(nid_ins, "d_tr_sed_H2SO4", "var mmr from sedim",
	240	. "kg/kg",
[1530]	241	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	242	. "ins(X)", zsto,zout)
[1442]	243	c
	244
	245	c
	246	CALL histdef(nid_ins, "d_tr_sed_H2O", "var mmr from sedim",
	247	. "kg/kg",
[1530]	248	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	249	. "ins(X)", zsto,zout)
[1442]	250	c
	251
	252	c
	253	CALL histdef(nid_ins, "F_sedim", "tendency from sedim",
	254	. "kg.m-2.s-1",
[1530]	255	. nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32,
[1305]	256	. "ins(X)", zsto,zout)
[1442]	257	c
	258	ENDIF
[1305]	259	ENDIF !lev_histins.GE.2
[888]	260	c
	261	c-------------------------------------------------------
[1305]	262	IF(lev_histins.GE.3) THEN
[888]	263	c
	264	cccccccccccccccccc Radiative transfer
	265	c
	266	c 2D
	267	c
	268	CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
[1530]	269	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	270	. "ins(X)", zsto,zout)
	271	c
	272	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
[1530]	273	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	274	. "ins(X)", zsto,zout)
	275	c
	276	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
[1530]	277	. nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	278	. "ins(X)", zsto,zout)
	279	c
	280	c 3D
	281	c
	282	CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
[1530]	283	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	284	. 32, "ins(X)", zsto,zout)
	285	c
	286	CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
[1530]	287	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	288	. 32, "ins(X)", zsto,zout)
	289	c
	290	CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2",
[1530]	291	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	292	. 32, "ins(X)", zsto,zout)
	293	c
	294	CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2",
[1530]	295	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	296	. 32, "ins(X)", zsto,zout)
	297	c
	298	CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2",
[1530]	299	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	300	. 32, "ins(X)", zsto,zout)
	301	c
	302	c CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2",
[1530]	303	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
[888]	304	c . 32, "ins(X)", zsto,zout)
	305	c
[1305]	306	ENDIF !lev_histins.GE.3
[888]	307	c
	308	c-------------------------------------------------------
[1305]	309	IF(lev_histins.GE.4) THEN
[888]	310	c
	311	CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
[1530]	312	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	313	. "ins(X)", zsto,zout)
	314	c
	315	c CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
[1530]	316	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	317	c . "ins(X)", zsto,zout)
	318	c
	319	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
[1530]	320	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	321	. "ins(X)", zsto,zout)
	322	c
	323	CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
[1530]	324	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	325	. "ins(X)", zsto,zout)
	326	c
	327	CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
[1530]	328	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	329	. "ins(X)", zsto,zout)
	330	c
	331	CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
[1530]	332	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	333	. "ins(X)", zsto,zout)
	334	c
	335	c CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
[1530]	336	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	337	c . "ins(X)", zsto,zout)
	338	c
	339	CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2",
[1530]	340	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	341	. "ins(X)", zsto,zout)
	342	c
	343	CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2",
[1530]	344	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	345	. "ins(X)", zsto,zout)
	346	c
	347	CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2",
[1530]	348	. nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	349	. "ins(X)", zsto,zout)
	350	c
	351	c CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
[1530]	352	c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
[888]	353	c . "ins(X)", zsto,zout)
	354	c
[1305]	355	ENDIF !lev_histins.GE.4
[888]	356	c
	357	c-------------------------------------------------------
[1305]	358	IF(lev_histins.GE.5) THEN
[888]	359	c
	360	c call histdef(nid_ins, "taux",
	361	c $ "Zonal wind stress", "Pa",
[1530]	362	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	363	c $ "ins(X)", zsto,zout)
	364	c
	365	c call histdef(nid_ins, "tauy",
	366	c $ "Meridional xind stress", "Pa",
[1530]	367	c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	368	c $ "ins(X)", zsto,zout)
	369	c
	370	c CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
[1530]	371	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	372	c . "ins(X)", zsto,zout)
	373	c
	374	c CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
[1530]	375	c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
[888]	376	c . "ins(X)", zsto,zout)
	377	c
[1305]	378	ENDIF !lev_histins.GE.5
[888]	379	c-------------------------------------------------------
	380	c
[3]	381	CALL histend(nid_ins)
[892]	382	!$OMP END MASTER
	383
[3]	384	ENDIF

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

Download in other formats: