Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf.h @ 520

Last change on this file since 520 was 511, checked in by lmdzadmin, 21 years ago
Enleve region "1d" des fichiers ini*.h IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 6.3 KB

Rev	Line
[416]	1
	2	IF (ok_hf) THEN
	3	c
[506]	4	zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
	5	zsto = zout
	6	zsto1 = dtime
	7	PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf
	8	c
[453]	9	idayref = day_ref
	10	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
	11
[416]	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
	19	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
	20
	21	cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
	22	CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
[511]	23	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
[416]	24	. nhori, nid_hf)
	25
	26	CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
	27	. klev, presnivs, nvert)
[506]	28	c
[416]	29
[506]	30	IF(lev_histhf.GE.1) THEN
[416]	31	c
	32	c CALL histdef(nid_hf, "phis", "Surface geop. height", "-",
	33	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	34	c . "once", zsto,zout)
	35	c
	36	c CALL histdef(nid_hf, "aire", "Grid area", "-",
	37	c . iim,jjmp1,nhori, 1,1,1, -99, 32,
	38	c . "once", zsto,zout)
	39	c
[506]	40	CALL histdef(nid_hf, "aireTER","Grid area CONT","-",
	41	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	42	. "once", zsto,zout)
	43	c
	44	CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-",
	45	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	46	. "once", zsto,zout)
	47	c
	48	CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-",
	49	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	50	. "once", zsto,zout)
	51	c
[416]	52	c Champs 2D:
	53	c
[506]	54	CALL histdef(nid_hf, "t2m", "Temperature 2m", "K",
[416]	55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	56	. "inst(X)", zsto,zout)
	57	c
[506]	58	CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg",
	59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	60	. "inst(X)", zsto,zout)
	61	c
[416]	62	CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
	63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	64	. "inst(X)", zsto,zout)
	65	c
	66	print*,'ATTENTION METTRE AVE(X) POUR LES PRECIPS'
	67
[506]	68	CALL histdef(nid_hf, "rain", "Precipitation", "kg/m^2s",
[416]	69	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[506]	70	. "ave(X)", zsto1,zout)
	71	c
	72	cIM ENSEMBLES BEG
	73	c
	74	CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
	75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	76	. "inst(X)", zsto,zout)
	77	c
[506]	78	CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa",
	79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	80	. "inst(X)", zsto,zout)
[506]	81	c
	82	c
	83	CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s",
	84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	85	. "inst(X)", zsto,zout)
	86	c
[506]	87	CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s",
	88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	89	. "inst(X)", zsto,zout)
[506]	90	c
	91	CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s",
	92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	93	. "inst(X)", zsto,zout)
[506]	94	c
	95	DO k=1, nlevENS
	96	IF(clev(k).EQ."500") THEN
	97	CALL histdef(nid_hf, "phi"//clev(k),
	98	. "Geopotential"//clev(k)//"mb", "m2/s2",
	99	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	100	. "inst(X)", zsto,zout)
[506]	101	ENDIF !clev(k).EQ."500"
	102	ENDDO !k=1, nlevENS
	103	c
	104	ENDIF !lev_histhf.GE.1
	105	c
	106	IF(lev_histhf.GE.2) THEN
	107	c
	108	CALL histdef(nid_hf, "cldt", "Total cloudiness", "-",
	109	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	110	. "ave(X)", zsto1,zout)
	111	c
	112	c
	113	CALL histdef(nid_hf, "SWdownOR",
	114	. "Sfce incident SW radiation OR", "W/m^2",
	115	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	116	. "ave(X)", zsto1,zout)
	117	c
	118	CALL histdef(nid_hf, "LWdownOR",
	119	. "Sfce incident LW radiation OR", "W/m^2",
	120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	121	. "ave(X)", zsto1,zout)
	122	ENDIF !lev_histhf.GE.2
	123	c
	124	cIM ENSEMBLES END
	125	c
	126	IF(lev_histhf.GE.3) THEN
	127	c
	128	DO k=1, nlevENS
	129	c
	130	CALL histdef(nid_hf, "t"//clev(k),
	131	. "Temperature"//clev(k)//"mb","K",
	132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	133	. "inst(X)", zsto,zout)
[506]	134	c
	135	IF(clev(k).NE."500") THEN
	136	CALL histdef(nid_hf, "phi"//clev(k),
	137	. "Geopotential"//clev(k)//"mb", "m2/s2",
	138	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[416]	139	. "inst(X)", zsto,zout)
[506]	140	ENDIF
	141	c
	142	CALL histdef(nid_hf, "q"//clev(k),
	143	. "Specific humidity"//clev(k)//"mb","kg/kg",
[467]	144	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	145	. "inst(X)", zsto,zout)
[506]	146	c
	147	IF(1.EQ.0) THEN
	148	CALL histdef(nid_hf, "rh"//clev(k),
	149	. "Relative humidity"//clev(k)//"mb", "%",
	150	. iim,jjmp1,nhori, 1,1,1, -99, 32,
[467]	151	. "inst(X)", zsto,zout)
[506]	152	ENDIF
	153	c
	154	CALL histdef(nid_hf, "u"//clev(k),
	155	. "Zonal wind"//clev(k)//"mb","m/s",
	156	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	157	. "inst(X)", zsto,zout)
	158	c
	159	CALL histdef(nid_hf, "v"//clev(k),
	160	. "Meridional wind"//clev(k)//"mb","m/s",
	161	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	162	. "inst(X)", zsto,zout)
	163	c
	164	ENDDO !nlevENS
	165	IF(1.EQ.0) THEN
[467]	166	CALL histdef(nid_hf, "cdrm", " Momentum drag coef.", "-",
	167	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	168	. "inst(X)", zsto,zout)
[506]	169	c
[467]	170	CALL histdef(nid_hf, "cdrh", "Heat drag coef.", "-",
	171	. iim,jjmp1,nhori, 1,1,1, -99, 32,
	172	. "inst(X)", zsto,zout)
[506]	173	ENDIF !(1.EQ.0) THEN
[416]	174	c
[506]	175	ENDIF !lev_histhf.GE.3
	176	c
	177	IF(lev_histhf.GE.4) THEN
	178	#define histhf3d
	179	#ifdef histhf3d
	180	#include "ini_histhf3d.h"
	181	#endif
	182	ENDIF !lev_histhf.GE.4
	183	c
	184	c#define histhf3d
	185	c#ifdef histhf3d
	186	c#include "ini_histhf3d.h"
	187	c#endif
	188	c
[416]	189	CALL histend(nid_hf)
[506]	190	c
[416]	191	endif ! ok_hf

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

Download in other formats: