Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_bilKP_ins.h @ 682

Last change on this file since 682 was 644, checked in by Laurent Fairhead, 19 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 10.4 KB

Rev	Line
[644]	1	c
	2	c $Header$
	3	c
	4	IF (ok_journe) THEN
	5	c
	6	zsto = dtime
	7	zout = dtime
	8	typeval=tinst
	9	c
	10	idayref = day_ref
	11	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	12	c
	13	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
	14	DO i = 1, iim
	15	zx_lon(i,1) = rlon(i+1)
	16	zx_lon(i,jjmp1) = rlon(i+1)
	17	ENDDO
	18	DO ll=1,klev
	19	znivsig(ll)=float(ll)
	20	ENDDO
	21	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	22	write(,)'zx_lon = ',zx_lon(:,1)
	23	write(,)'zx_lat = ',zx_lat(1,:)
	24	c
	25	cIM 280405 BEG
	26	c
	27	cIM cf. AM 081204 BEG region
	28	imin_ins=1
	29	imax_ins=iim
	30	jmin_ins=1
	31	jmax_ins=jjmp1
	32	do i=1,iim-1
	33	if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
	34	if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
	35	enddo
	36	do j=1,jjmp1-1
	37	if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
	38	if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
	39	enddo
	40	c
	41	print*,'On stoke le fichier bilKP instantanne sur ',
	42	s imin_ins,imax_ins,jmin_ins,jmax_ins
	43	print*,'On stoke le fichier bilKP instantanne sur ',
	44	s zx_lon(imin_ins,1),zx_lon(imax_ins,1),
	45	s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins)
	46	cIM cf. AM 081204 END region
	47	c
	48	cIM 280405 END
	49	c
	50	IF(1.EQ.0) THEN
	51	CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
	52	. zx_lat(1,:),
	53	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
	54	. nhori, nid_bilKPins)
	55	ENDIF
	56	c
	57	cIM 280405 BEG
	58	c
	59	cIM cf. AM 081204 BEG region
	60	CALL histbeg("histbilKP_ins", iim,zx_lon(:,1),
	61	. jjmp1,zx_lat(1,:),
	62	. imin_ins,imax_ins-imin_ins+1,
	63	. jmin_ins,jmax_ins-jmin_ins+1,
	64	. itau_phy, zjulian, dtime,
	65	. nhori, nid_bilKPins)
	66	cIM 081204 END
	67	c
	68	cIM 280405 END
	69	c
	70	write(,)'Journee ', itau_phy, zjulian
	71	CALL histvert(nid_bilKPins, "presnivs",
	72	. "Vertical levels","mb",
	73	. klev, presnivs/100., nvert)
	74	c
	75	c Champs 3D:
	76	c
	77	CALL histdef(nid_bilKPins,"ue",
	78	. "Zonal energy transport","-",
	79	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	80	. typeval, zsto,zout)
	81	c
	82	CALL histdef(nid_bilKPins,"ve",
	83	. "Merid energy transport","-",
	84	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	85	. typeval, zsto,zout)
	86	c
	87	CALL histdef(nid_bilKPins,"uq",
	88	. "Zonal humidity transport","-",
	89	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	90	. typeval, zsto,zout)
	91	c
	92	CALL histdef(nid_bilKPins,"vq",
	93	. "Merid humidity transport","-",
	94	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	95	. typeval, zsto,zout)
	96	c
	97	c Champs 3D:
	98	c
	99	CALL histdef(nid_bilKPins, "temp",
	100	. "Air temperature", "K",
	101	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	102	. typeval, zsto,zout)
	103	c
	104	CALL histdef(nid_bilKPins,"ovap",
	105	. "Specific humidity","Kg/Kg",
	106	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	107	. typeval, zsto,zout)
	108	c
	109	CALL histdef(nid_bilKPins,"geop",
	110	. "Geopotential height", "m",
	111	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	112	. typeval, zsto,zout)
	113	c
	114	CALL histdef(nid_bilKPins,"vitu",
	115	. "Zonal wind", "m/s",
	116	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	117	. typeval, zsto,zout)
	118	c
	119	CALL histdef(nid_bilKPins,"vitv",
	120	. "Meridional wind", "m/s",
	121	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	122	. typeval, zsto,zout)
	123	c
	124	CALL histdef(nid_bilKPins, "vitw",
	125	. "Vertical wind", "m/s",
	126	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	127	. typeval, zsto,zout)
	128	c
	129	CALL histdef(nid_bilKPins, "pres",
	130	. "Inter-Layer Air pressure",
	131	. "Pa",
	132	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	133	. typeval, zsto,zout)
	134	c
	135	CALL histdef(nid_bilKPins, "play",
	136	. "Mean-Layer Air pressure",
	137	. "Pa",
	138	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	139	. typeval, zsto,zout)
	140	c
	141	CALL histdef(nid_bilKPins, "oliq",
	142	. "Liquid water content",
	143	. "kg/kg",
	144	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	145	. typeval, zsto,zout)
	146	c
	147	CALL histdef(nid_bilKPins, "dtdyn",
	148	. "Dynamics dT", "K/s",
	149	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	150	. typeval, zsto,zout)
	151	c
	152	CALL histdef(nid_bilKPins, "dqdyn",
	153	. "Dynamics dQ", "Kg/Kg/s",
	154	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	155	. typeval, zsto,zout)
	156	c
	157	CALL histdef(nid_bilKPins, "dtcon",
	158	. "Convection dT", "K/s",
	159	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	160	. typeval, zsto,zout)
	161	c
	162	CALL histdef(nid_bilKPins, "ducon",
	163	. "Convection du", "m/s2",
	164	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	165	. typeval, zsto,zout)
	166	c
	167	CALL histdef(nid_bilKPins, "dvcon",
	168	. "Convection dv", "m/s2",
	169	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	170	. typeval, zsto,zout)
	171	c
	172	CALL histdef(nid_bilKPins,"dqcon",
	173	. "Convection dQ","Kg/Kg/s",
	174	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	175	. typeval, zsto,zout)
	176	c
	177	CALL histdef(nid_bilKPins, "dtlsc",
	178	. "Condensation dT", "K/s",
	179	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	180	. typeval, zsto,zout)
	181	c
	182	CALL histdef(nid_bilKPins,"dqlsc",
	183	. "Condensation dQ","Kg/Kg/s",
	184	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	185	. typeval, zsto,zout)
	186	c
	187	CALL histdef(nid_bilKPins,"dtvdf",
	188	. "Boundary-layer dT","K/s",
	189	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	190	. typeval, zsto,zout)
	191	c
	192	CALL histdef(nid_bilKPins, "dqvdf",
	193	. "Boundary-layer dQ",
	194	. "Kg/Kg/s",
	195	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	196	. typeval, zsto,zout)
	197	c
	198	CALL histdef(nid_bilKPins,"dtajs",
	199	. "Ajustement sec dT","K/s",
	200	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	201	. typeval, zsto,zout)
	202	c
	203	CALL histdef(nid_bilKPins,"dqajs",
	204	. "Ajustement sec dQ",
	205	. "Kg/Kg/s",
	206	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	207	. typeval, zsto,zout)
	208	c
	209	CALL histdef(nid_bilKPins,"dteva",
	210	. "Reevaporation dT","K/s",
	211	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	212	. typeval, zsto,zout)
	213	c
	214	CALL histdef(nid_bilKPins,"dqeva",
	215	. "Reevaporation dQ",
	216	. "Kg/Kg/s",
	217	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	218	. typeval, zsto,zout)
	219
	220	c
	221	CALL histdef(nid_bilKPins, "dtswr",
	222	. "SW radiation dT", "K/s",
	223	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	224	. typeval, zsto,zout)
	225	c
	226	CALL histdef(nid_bilKPins, "dtsw0",
	227	. "SW radiation dT", "K/s",
	228	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	229	. typeval, zsto,zout)
	230	c
	231	CALL histdef(nid_bilKPins, "dtlwr",
	232	. "LW radiation dT", "K/s",
	233	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	234	. typeval, zsto,zout)
	235	c
	236	CALL histdef(nid_bilKPins, "dtlw0",
	237	. "LW radiation dT", "K/s",
	238	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	239	. typeval, zsto,zout)
	240	c
	241	CALL histdef(nid_bilKPins,"duvdf",
	242	. "Boundary-layer dU","m/s2",
	243	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	244	. typeval, zsto,zout)
	245	c
	246	CALL histdef(nid_bilKPins,"dvvdf",
	247	. "Boundary-layer dV","m/s2",
	248	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	249	. typeval, zsto,zout)
	250	c
	251	IF (ok_orodr) THEN
	252	IF (ok_orolf) THEN
	253	CALL histdef(nid_bilKPins, "duoli",
	254	. "Orography dU", "m/s2",
	255	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	256	. typeval, zsto,zout)
	257	c
	258	CALL histdef(nid_bilKPins, "dvoli",
	259	. "Orography dV", "m/s2",
	260	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	261	. typeval, zsto,zout)
	262	c
	263	ENDIF
	264	ENDIF
	265	C
	266	CALL histdef(nid_bilKPins, "duphy",
	267	. "Physiq dU","-",
	268	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	269	. typeval, zsto,zout)
	270	C
	271	CALL histdef(nid_bilKPins, "dvphy",
	272	. "Physiq dV","-",
	273	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	274	. typeval, zsto,zout)
	275	C
	276	CALL histdef(nid_bilKPins, "dtphy",
	277	. "Physiq dT","-",
	278	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	279	. typeval, zsto,zout)
	280	C
	281	CALL histdef(nid_bilKPins, "dqphy",
	282	. "Physiq dQ","-",
	283	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	284	. typeval, zsto,zout)
	285	C
	286	CALL histdef(nid_bilKPins, "dqlphy",
	287	. "Physiq dQl","-",
	288	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
	289	. typeval, zsto,zout)
	290	cIM 280405 BEG
	291	c
	292	c Champs 2D:
	293	c
	294	c u850, v850
	295	c DO k=1, nlevSTD
	296	DO k=1, 12
	297	c
	298	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
	299	IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
	300	c
	301	IF(bb2.EQ."850") THEN
	302	c
	303	CALL histdef(nid_bilKPins, "u"//bb2,
	304	. "Zonal wind "//bb2//"mb","m/s",
	305	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	306	. typeval, zsto,zout)
	307	c
	308	CALL histdef(nid_bilKPins, "v"//bb2,
	309	. "Meridional wind "//bb2//"mb","m/s",
	310	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	311	. typeval, zsto,zout)
	312	c
	313	ENDIF !(bb2.EQ."850")
	314	c
	315	ENDDO !k=1, 12
	316	c
	317	cIM 280405 END
	318	c
	319	CALL histend(nid_bilKPins)
	320	c
	321	ndex2d = 0
	322	ndex3d = 0
	323	c
	324	ENDIF ! fin de test sur ok_journe

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

Download in other formats: