Context Navigation

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

Last change on this file since 418 was 418, checked in by (none), 22 years ago
This commit was manufactured by cvs2svn to create branch 'rel-LF'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 8.9 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
70	CALL histdef(nid_ins, "qsol", "Surface humidity", "mm",
71	. iim,jjmp1,nhori, 1,1,1, -99, 32,
72	. "inst(X)", zsto,zout)
73	c
74	CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
75	. iim,jjmp1,nhori, 1,1,1, -99, 32,
76	. "inst(X)", zsto,zout)
77	c
78	CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
79	. iim,jjmp1,nhori, 1,1,1, -99, 32,
80	. "inst(X)", zsto,zout)
81	c
82	CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
83	. "kg/s",
84	. iim,jjmp1,nhori, 1,1,1, -99, 32,
85	. "inst(X)", zsto,zout)
86	c
87	CALL histdef(nid_ins, "snow", "Snow fall", "kg/s",
88	. iim,jjmp1,nhori, 1,1,1, -99, 32,
89	. "inst(X)", zsto,zout)
90	c
91	CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
92	. iim,jjmp1,nhori, 1,1,1, -99, 32,
93	. "inst(X)", zsto,zout)
94	c
95	CALL histdef(nid_ins, "topl", "OLR", "W/m2",
96	. iim,jjmp1,nhori, 1,1,1, -99, 32,
97	. "inst(X)", zsto,zout)
98	c
99	CALL histdef(nid_ins, "evap", "Evaporation", "kg/s",
100	. iim,jjmp1,nhori, 1,1,1, -99, 32,
101	. "inst(X)", zsto,zout)
102	c
103	CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
104	. iim,jjmp1,nhori, 1,1,1, -99, 32,
105	. "inst(X)", zsto,zout)
106	c
107	CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
108	. iim,jjmp1,nhori, 1,1,1, -99, 32,
109	. "inst(X)", zsto,zout)
110	c
111	CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
112	. "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
113	. "inst(X)", zsto,zout)
114	c
115	CALL histdef(nid_ins, "bils", "Surf. total 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, "sens", "Sensible heat flux", "W/m2",
120	. iim,jjmp1,nhori, 1,1,1, -99, 32,
121	. "inst(X)", zsto,zout)
122	c
123	CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
124	. iim,jjmp1,nhori, 1,1,1, -99, 32,
125	. "inst(X)", zsto,zout)
126	c
127	CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
128	. iim,jjmp1,nhori, 1,1,1, -99, 32,
129	. "inst(X)", zsto,zout)
130	c
131	CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
132	. iim,jjmp1,nhori, 1,1,1, -99, 32,
133	. "inst(X)", zsto,zout)
134	c
135	CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
136	. iim,jjmp1,nhori, 1,1,1, -99, 32,
137	. "inst(X)", zsto,zout)
138	c
139	CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
140	. iim,jjmp1,nhori, 1,1,1, -99, 32,
141	. "inst(X)", zsto,zout)
142
143	DO nsrf = 1, nbsrf
144	C
145	call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
146	$ "Fraction"//clnsurf(nsrf), "W/m2",
147	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
148	$ "inst(X)", zsto,zout)
149
150	call histdef(nid_ins, "sens_"//clnsurf(nsrf),
151	$ "Sensible heat flux "//clnsurf(nsrf), "W/m2",
152	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
153	$ "inst(X)", zsto,zout)
154	c
155	call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
156	$ "Surface Temperature"//clnsurf(nsrf), "W/m2",
157	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
158	$ "inst(X)", zsto,zout)
159	c
160	call histdef(nid_ins, "lat_"//clnsurf(nsrf),
161	$ "Latent heat flux "//clnsurf(nsrf), "W/m2",
162	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
163	$ "inst(X)", zsto,zout)
164	C
165	call histdef(nid_ins, "taux_"//clnsurf(nsrf),
166	$ "Zonal wind stress"//clnsurf(nsrf),"Pa",
167	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
168	$ "inst(X)", zsto,zout)
169
170	call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
171	$ "Meridional xind stress "//clnsurf(nsrf), "Pa",
172	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
173	$ "inst(X)", zsto,zout)
174	c
175	call histdef(nid_ins, "albe_"//clnsurf(nsrf),
176	$ "Albedo "//clnsurf(nsrf), "-",
177	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
178	$ "inst(X)", zsto,zout)
179	c
180	call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
181	$ "rugosite "//clnsurf(nsrf), "-",
182	$ iim,jjmp1,nhori, 1,1,1, -99, 32,
183	$ "inst(X)", zsto,zout)
184	CXXX
185	END DO
186	CALL histdef(nid_ins, "rugs", "rugosity", "-",
187	. iim,jjmp1,nhori, 1,1,1, -99, 32,
188	. "inst(X)", zsto,zout)
189
190	c
191	CALL histdef(nid_ins, "albs", "Surface albedo", "-",
192	. iim,jjmp1,nhori, 1,1,1, -99, 32,
193	. "inst(X)", zsto,zout)
194	CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
195	. iim,jjmp1,nhori, 1,1,1, -99, 32,
196	. "inst(X)", zsto,zout)
197	c
198	c
199	c Champs 3D:
200	c
201	CALL histdef(nid_ins, "temp", "Temperature", "K",
202	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
203	. "inst(X)", zsto,zout)
204	c
205	CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
206	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
207	. "inst(X)", zsto,zout)
208	c
209	CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
210	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
211	. "inst(X)", zsto,zout)
212	c
213	CALL histdef(nid_ins, "geop", "Geopotential height", "m",
214	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
215	. "inst(X)", zsto,zout)
216	c
217	CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
218	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
219	. "inst(X)", zsto,zout)
220	c
221	CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
222	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
223	. "inst(X)", zsto,zout)
224	c
225	CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
226	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
227	. "inst(X)", zsto,zout)
228	c
229
230	CALL histend(nid_ins)
231	c
232	ndex2d = 0
233	ndex3d = 0
234	c
235	ENDIF

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

Download in other formats: