Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histins.h @ 451

Last change on this file since 451 was 451, checked in by lmdzadmin, 21 years ago
L'unite de Precipitation,Evaporation, Snow Fall passe de kg/s a kg/(s*m2) JP/NT IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 8.8 KB

Line
1	IF (ok_instan) THEN
2	c
3	idayref = day_ref
4	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
5	c
6	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
7	DO i = 1, iim
8	zx_lon(i,1) = rlon(i+1)
9	zx_lon(i,jjmp1) = rlon(i+1)
10	ENDDO
11	DO ll=1,klev
12	znivsig(ll)=float(ll)
13	ENDDO
14	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
15	CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
16	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
17	. nhori, nid_ins)
18	write(,)'Inst ', itau_phy, zjulian
19	CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb",
20	. klev, presnivs, nvert)
21	c call histvert(nid_ins, 'sig_s', 'Niveaux sigma','-',
22	c . klev, znivsig, nvert)
23	c
24	c
25	zsto = dtime * ecrit_ins
26	zout = dtime * ecrit_ins
27	C
28	CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
29	. iim,jjmp1,nhori, 1,1,1, -99, 32,
30	. "once", zsto,zout)
31	c
32	CALL histdef(nid_ins, "aire", "Grid area", "-",
33	. iim,jjmp1,nhori, 1,1,1, -99, 32,
34	. "once", zsto,zout)
35	c
36	c Champs 2D:
37	c
38	CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
39	. iim,jjmp1,nhori, 1,1,1, -99, 32,
40	. "inst(X)", zsto,zout)
41	c
42	CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
43	. iim,jjmp1,nhori, 1,1,1, -99, 32,
44	. "inst(X)", zsto,zout)
45	c
46	CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
47	. iim,jjmp1,nhori, 1,1,1, -99, 32,
48	. "inst(X)", zsto,zout)
49	c
50	CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
51	. iim,jjmp1,nhori, 1,1,1, -99, 32,
52	. "inst(X)", zsto,zout)
53	c
54	CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
55	. iim,jjmp1,nhori, 1,1,1, -99, 32,
56	. "inst(X)", zsto,zout)
57	c
58	CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
59	. iim,jjmp1,nhori, 1,1,1, -99, 32,
60	. "inst(X)", zsto,zout)
61	c
62	CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
63	. iim,jjmp1,nhori, 1,1,1, -99, 32,
64	. "inst(X)", zsto,zout)
65	c
66	CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
67	. iim,jjmp1,nhori, 1,1,1, -99, 32,
68	. "inst(X)", zsto,zout)
69	c
70	CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
72	. "inst(X)", zsto,zout)
73	c
74	CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
76	. "inst(X)", zsto,zout)
77	c
78	CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
79	. "kg/(s*m2)",
80	. iim,jjmp1,nhori, 1,1,1, -99, 32,
81	. "inst(X)", zsto,zout)
82	c
83	CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)",
84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
85	. "inst(X)", zsto,zout)
86	c
87	CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
89	. "inst(X)", zsto,zout)
90	c
91	CALL histdef(nid_ins, "topl", "OLR", "W/m2",
92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
93	. "inst(X)", zsto,zout)
94	c
95	CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)",
96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
97	. "inst(X)", zsto,zout)
98	c
99	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
101	. "inst(X)", zsto,zout)
102	c
103	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
105	. "inst(X)", zsto,zout)
106	c
107	CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
108	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
109	. "inst(X)", zsto,zout)
110	c
111	CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
112	. iim,jjmp1,nhori, 1,1,1, -99, 32,
113	. "inst(X)", zsto,zout)
114	c
115	CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2",
116	. iim,jjmp1,nhori, 1,1,1, -99, 32,
117	. "inst(X)", zsto,zout)
118	c
119	CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
121	. "inst(X)", zsto,zout)
122	c
123	CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
124	. iim,jjmp1,nhori, 1,1,1, -99, 32,
125	. "inst(X)", zsto,zout)
126	c
127	CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
129	. "inst(X)", zsto,zout)
130	c
131	CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
133	. "inst(X)", zsto,zout)
134	c
135	CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
137	. "inst(X)", zsto,zout)
138
139	DO nsrf = 1, nbsrf
140	C
141	call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
142	$ "Fraction"//clnsurf(nsrf), "W/m2",
143	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
144	$ "inst(X)", zsto,zout)
145
146	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
147	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
148	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
149	$ "inst(X)", zsto,zout)
150	c
151	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
152	$ "Surface Temperature"//clnsurf(nsrf), "W/m2",
153	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
154	$ "inst(X)", zsto,zout)
155	c
156	call histdef(nid_ins, "lat_"//clnsurf(nsrf),
157	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
158	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
159	$ "inst(X)", zsto,zout)
160	C
161	call histdef(nid_ins, "taux_"//clnsurf(nsrf),
162	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
163	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
164	$ "inst(X)", zsto,zout)
165
166	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
167	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
168	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
169	$ "inst(X)", zsto,zout)
170	c
171	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
172	$ "Albedo "//clnsurf(nsrf), "-",
173	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
174	$ "inst(X)", zsto,zout)
175	c
176	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
177	$ "rugosite "//clnsurf(nsrf), "-",
178	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
179	$ "inst(X)", zsto,zout)
180	CXXX
181	END DO
182	CALL histdef(nid_ins, "rugs", "rugosity", "-",
183	. iim,jjmp1,nhori, 1,1,1, -99, 32,
184	. "inst(X)", zsto,zout)
185
186	c
187	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
188	. iim,jjmp1,nhori, 1,1,1, -99, 32,
189	. "inst(X)", zsto,zout)
190	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
191	. iim,jjmp1,nhori, 1,1,1, -99, 32,
192	. "inst(X)", zsto,zout)
193	c
194	c
195	c Champs 3D:
196	c
197	CALL histdef(nid_ins, "temp", "Temperature", "K",
198	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
199	. "inst(X)", zsto,zout)
200	c
201	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
202	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
203	. "inst(X)", zsto,zout)
204	c
205	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
206	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
207	. "inst(X)", zsto,zout)
208	c
209	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
210	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
211	. "inst(X)", zsto,zout)
212	c
213	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
214	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
215	. "inst(X)", zsto,zout)
216	c
217	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
218	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
219	. "inst(X)", zsto,zout)
220	c
221	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
222	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
223	. "inst(X)", zsto,zout)
224	c
225
226	CALL histend(nid_ins)
227	c
228	ndex2d = 0
229	ndex3d = 0
230	c
231	ENDIF

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

Download in other formats: