Context Navigation

ini_bilKP_ins.h @ 3632

Last change on this file since 3632 was 1910, checked in by Laurent Fairhead, 11 years ago
Merged trunk changes r1860:1909 into testing branch
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 10.7 KB

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

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

Context Navigation

source: LMDZ5/branches/testing/libf/phylmd/ini_bilKP_ins.h @ 3632

Download in other formats: