Context Navigation

source: LMDZ4/trunk/libf/phytherm/ini_histins.h @ 1080

Last change on this file since 1080 was 814, checked in by Laurent Fairhead, 17 years ago
Rajout de la physique utilisant les thermiques FH LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 11.0 KB

Rev	Line
[814]	1	!
	2	! $Header$
	3	!
	4	IF (ok_instan) THEN
	5	c$OMP MASTER
	6	c
	7	zsto = ecrit_ins
	8	zout = ecrit_ins
	9	c
	10	idayref = day_ref
	11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	12	c
	13	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	14	cym DO i = 1, iim
	15	cym zx_lon(i,1) = rlon(i+1)
	16	cym zx_lon(i,jjmp1) = rlon(i+1)
	17	cym ENDDO
	18	DO ll=1,klev
	19	znivsig(ll)=float(ll)
	20	ENDDO
	21	cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	22	cym CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
	23	cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	24	cym . nhori, nid_ins)
	25	CALL histbeg_phy("histins", itau_phy, zjulian, dtime,
	26	. nhori, nid_ins)
	27	write(,)'Inst ', itau_phy, zjulian
	28	CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
	29	cIM . klev, presnivs, nvert)
	30	. klev, presnivs/100., nvert)
	31	c call histvert(nid_ins, 'sig_s', 'Niveaux sigma','-',
	32	c . klev, znivsig, nvert)
	33	c
	34	c
	35	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
	36	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	37	. "once", zsto,zout)
	38	c
	39	CALL histdef(nid_ins, "aire", "Grid area", "-",
	40	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	41	. "once", zsto,zout)
	42	c
	43	c Champs 2D:
	44	c
	45	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
	46	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	47	. "inst(X)", zsto,zout)
	48	c
	49	CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
	50	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	51	. "inst(X)", zsto,zout)
	52	c
	53	CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
	54	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	55	. "inst(X)", zsto,zout)
	56	c
	57	CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
	58	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	59	. "inst(X)", zsto,zout)
	60	c
	61	CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
	62	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	63	. "inst(X)", zsto,zout)
	64	c
	65	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
	66	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	67	. "inst(X)", zsto,zout)
	68	c
	69	CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
	70	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	71	. "inst(X)", zsto,zout)
	72	c
	73	CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
	74	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	75	. "inst(X)", zsto,zout)
	76	c
	77	CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
	78	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	79	. "inst(X)", zsto,zout)
	80	c
	81	CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
	82	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	83	. "inst(X)", zsto,zout)
	84	c
	85	CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
	86	. "kg/(s*m2)",
	87	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	88	. "inst(X)", zsto,zout)
	89	c
	90	CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)",
	91	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	92	. "inst(X)", zsto,zout)
	93	c
	94	c CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
	95	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	96	c . "inst(X)", zsto,zout)
	97	c
	98	CALL histdef(nid_ins, "topl", "OLR", "W/m2",
	99	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	100	. "inst(X)", zsto,zout)
	101	c
	102	CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)",
	103	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	104	. "inst(X)", zsto,zout)
	105	c
	106	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
	107	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	108	. "inst(X)", zsto,zout)
	109	c
	110	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
	111	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	112	. "inst(X)", zsto,zout)
	113	c
	114	CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
	115	. "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
	116	. "inst(X)", zsto,zout)
	117	c
	118	CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
	119	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	120	. "inst(X)", zsto,zout)
	121	c
	122	CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2",
	123	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	124	. "inst(X)", zsto,zout)
	125	c
	126	CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
	127	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	128	. "inst(X)", zsto,zout)
	129	c
	130	CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
	131	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	132	. "inst(X)", zsto,zout)
	133	c
	134	CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
	135	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	136	. "inst(X)", zsto,zout)
	137	c
	138	CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
	139	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	140	. "inst(X)", zsto,zout)
	141	c
	142	CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
	143	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	144	. "inst(X)", zsto,zout)
	145
	146	DO nsrf = 1, nbsrf
	147	C
	148	call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
	149	$ "% "//clnsurf(nsrf), "%",
	150	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	151	$ "inst(X)", zsto,zout)
	152
	153	call histdef(nid_ins, "fract_"//clnsurf(nsrf),
	154	$ "Fraction "//clnsurf(nsrf), "1",
	155	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	156	$ "inst(X)", zsto,zout)
	157
	158	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
	159	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
	160	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	161	$ "inst(X)", zsto,zout)
	162	c
	163	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
	164	$ "Surface Temperature"//clnsurf(nsrf), "W/m2",
	165	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	166	$ "inst(X)", zsto,zout)
	167	c
	168	call histdef(nid_ins, "lat_"//clnsurf(nsrf),
	169	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
	170	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	171	$ "inst(X)", zsto,zout)
	172	C
	173	call histdef(nid_ins, "taux_"//clnsurf(nsrf),
	174	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
	175	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	176	$ "inst(X)", zsto,zout)
	177
	178	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
	179	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
	180	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	181	$ "inst(X)", zsto,zout)
	182	c
	183	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
	184	$ "Albedo "//clnsurf(nsrf), "-",
	185	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	186	$ "inst(X)", zsto,zout)
	187	c
	188	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
	189	$ "rugosite "//clnsurf(nsrf), "-",
	190	$ iim,jj_nb,nhori, 1,1,1, -99, 32,
	191	$ "inst(X)", zsto,zout)
	192	CXXX
	193	END DO
	194	CALL histdef(nid_ins, "rugs", "rugosity", "-",
	195	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	196	. "inst(X)", zsto,zout)
	197
	198	c
	199	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
	200	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	201	. "inst(X)", zsto,zout)
	202	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
	203	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	204	. "inst(X)", zsto,zout)
	205	c
	206	c
	207	cIM cf. AM 081204 BEG
	208	c HBTM2
	209	CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m",
	210	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	211	. "inst(X)", zsto,zout)
	212	c
	213	CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",
	214	. "K",
	215	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	216	. "inst(X)", zsto,zout)
	217	c
	218	CALL histdef(nid_ins, "s_lcl", "Condensation level", "m",
	219	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	220	. "inst(X)", zsto,zout)
	221	c
	222	CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for
	223	. ABL", "J/m2",
	224	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	225	. "inst(X)", zsto,zout)
	226	c
	227	CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2",
	228	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	229	. "inst(X)", zsto,zout)
	230	c
	231	CALL histdef(nid_ins, "s_cteiCL", "Instability criteria
	232	. (ABL)", "K",
	233	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	234	. "inst(X)", zsto,zout)
	235	c
	236	CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K",
	237	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	238	. "inst(X)", zsto,zout)
	239	c
	240	CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2",
	241	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	242	. "inst(X)", zsto,zout)
	243	c
	244	CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2",
	245	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	246	. "inst(X)", zsto,zout)
	247	c
	248	CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m",
	249	. iim,jj_nb,nhori, 1,1,1, -99, 32,
	250	. "inst(X)", zsto,zout)
	251	c
	252	cIM cf. AM 081204 END
	253	c
	254	c Champs 3D:
	255	c
	256	CALL histdef(nid_ins, "temp", "Temperature", "K",
	257	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	258	. "inst(X)", zsto,zout)
	259	c
	260	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
	261	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	262	. "inst(X)", zsto,zout)
	263	c
	264	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
	265	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	266	. "inst(X)", zsto,zout)
	267	c
	268	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
	269	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	270	. "inst(X)", zsto,zout)
	271	c
	272	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
	273	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	274	. "inst(X)", zsto,zout)
	275	c
	276	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
	277	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	278	. "inst(X)", zsto,zout)
	279	c
	280	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
	281	. iim,jj_nb,nhori, klev,1,klev,nvert, 32,
	282	. "inst(X)", zsto,zout)
	283	c
	284
	285	CALL histend(nid_ins)
	286	c
	287	ndex2d = 0
	288	ndex3d = 0
	289	c
	290	c$OMP END MASTER
	291
	292	ENDIF

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

Download in other formats: