Context Navigation

← Previous Revision
Next Revision →
Blame
Revision Log

coare_validation_net_1992_93_all_v10.pro

Last change on this file was 203, checked in by pinsard, 10 years ago
fix thanks to coding rules
Property svn:keywords set to `Id URL`
File size: 20.1 KB

Line
1	;+
2	; .. _coare_validation_net_1992_93_all_v10.pro:
3	;
4	; ========================================
5	; coare_validation_net_1992_93_all_v10.pro
6	; ========================================
7	;
8	; DESCRIPTION
9	; ===========
10	;
11	; .. graphviz::
12	;
13	; digraph coare_validation_net_1992_93_all_v10 {
14	;
15	; met_data [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met_data_as_expt_199_p.txt"];
16	; swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_swr_19910101_20091231_oafluxgrid.nc"];
17	; lwr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lwr_19910101_20091231_oafluxgrid.nc"];
18	; lhf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lhf_19910101_20091231_oafluxgrid.nc"];
19	; shf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_shf_19910101_20091231_oafluxgrid.nc"];
20	; tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"];
21	; swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"];
22	; swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_oafluxgrid_1985_2007.nc"];
23	; lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"];
24	; lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2007.nc"];
25	; shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2007.nc"];
26	; swr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"];
27	; lwr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"];
28	; lhf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"];
29	; shf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"];
30	;
31	; ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ncep1_flux_19890101_20091231.nc"];
32	;
33	; figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/coare_validation_net_1992_93_all_v11.ps"];
34	;
35	; coare_validation_net_1992_93_all_v10 [shape=box,
36	; fontname=Courier,
37	; color=blue,
38	; URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/coare_validation_net_1992_93_all_v10.pro",
39	; label="${TROPFLUX}/src/paper01/fig14/coare_validation_net_1992_93_all_v10.pro"];
40	;
41	; {met_data swr_erai lwr_erai lhf_erai shf_erai tropflux_general swr_tropflux_nrt swr_oafluxgrid lwr_oafluxgrid lhf_oafluxgrid shf_oafluxgrid swr_ncep2 lwr_ncep2 lhf_ncep2 shf_ncep2 ncep1} -> {coare_validation_net_1992_93_all_v10} -> {figure}
42	; }
43	;
44	; SEE ALSO
45	; ========
46	;
47	; :ref:`project_profile.sh`
48	; :ref:`project_init.pro`
49	; :ref:`cm_project.pro`
50	;
51	; :ref:`statistics.pro`
52	;
53	; EXAMPLES
54	; ========
55	;
56	; .. code-block:: idl
57	;
58	; coare_validation_net_1992_93_all_v10
59	;
60	; TODO
61	; ====
62	;
63	; make it work on cratos : missing data
64	;
65	; coding rules
66	;
67	; complete description
68	;
69	; handle IO error
70	;
71	; EVOLUTIONS
72	; ==========
73	;
74	; $Id$
75	;
76	; $URL$
77	;
78	; - fplod 20110420T120449Z aedon.locean-ipsl.upmc.fr (Darwin)
79	;
80	; * remove hard coding path
81	; * add graphviz
82	;
83	; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
84	;
85	; * minimal header
86	;
87	;-
88	pro coare_validation_net_1992_93_all_v10
89	@cm_general
90	@cm_project
91	reinitplt, /z,/invert
92	key_portrait = 1
93	coefpalit=.9
94	;
95	openps, FILENAME = project_od_env+"coare_validation_net_1992_93_all_v11.ps"
96	; partie a changer
97	;
98	box=[155.5, 156.5, -2, -1]
99	;
100	; first leg
101	; start data - 11/11/92 and end date 03/12/92
102	fi=project_id_env+"met_data_coare_expt_1992_p1.txt"
103	res=read_ascii(fi,data_start=1)
104	ff=res.field1
105	date=reform(ff(0,*))
106	sw92_1=reform(ff(7,*))/0.93
107	lw92_1=reform(ff(8,*))
108	lh92_1=reform(ff(5,*))
109	sh92_1=reform(ff(6,*))
110	;
111	net1=sw92_1+lw92_1+lh92_1+sh92_1
112	;
113	st1=921111
114	en1=921203
115	;
116	; second leg
117	;
118	fi=project_id_env+"met_data_coare_expt_1992_p2.txt"
119	res=read_ascii(fi,data_start=1)
120	ff=res.field1
121	sw92_2=reform(ff(7,*))/0.93
122	lw92_2=reform(ff(8,*))
123	lh92_2=reform(ff(5,*))
124	sh92_2=reform(ff(6,*))
125	;
126	net2=sw92_2+lw92_2+lh92_2+sh92_2
127	;
128	; start data - 17/12/92 and end date 31/12/92
129	st2=19921217
130	en2=19921231
131	;
132	; third leg
133	;
134	fi=project_id_env+"met_data_coare_expt_1993_p1.txt"
135	res=read_ascii(fi,data_start=1)
136	ff=res.field1
137	sw93_1=reform(ff(7,*))/0.93
138	lw93_1=reform(ff(8,*))
139	lh93_1=reform(ff(5,*))
140	sh93_1=reform(ff(6,*))
141	;
142	net3=sw93_1+lw93_1+lh93_1+sh93_1
143	;
144	; start data - 01/01/93 and end date 11/01/93
145	st3=19930101
146	en3=19930111
147	;
148	; fourth leg
149	;
150	fi=project_id_env+"met_data_coare_expt_1993_p2.txt"
151	res=read_ascii(fi,data_start=1)
152	ff=res.field1
153	sw93_2=reform(ff(7,*))
154	lw93_2=reform(ff(8,*))
155	lh93_2=reform(ff(5,*))
156	sh93_2=reform(ff(6,*))
157	;
158	net4=sw93_2+lw93_2+lh93_2+sh93_2
159	;
160	; start data - 28/01/93 and end date 16/02/93
161	st4=19930128
162	en4=19930216
163	;
164	net=[net1, net2, net3, net4]
165	help, net
166	;
167	; reading other data sets
168	;
169	st1=19921111
170	en1=19921203
171	st2=19921217
172	en2=19921231
173	st3=19930101
174	en3=19930111
175	st4=19930128
176	en4=19930216
177	;
178	file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc"
179	initncdf, file
180	swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)
181	swr1=grossemoyenne(swr1, "xy", box=box,/nan)
182	swr1=reform(swr1)
183	swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)
184	swr2=grossemoyenne(swr2, "xy", box=box,/nan)
185	swr2=reform(swr2)
186	swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box)
187	swr3=grossemoyenne(swr3, "xy", box=box,/nan)
188	swr3=reform(swr3)
189	swr4=read_ncdf('swr',st4,en4,file=file,/nostr, box=box)
190	swr4=grossemoyenne(swr4, "xy", box=box,/nan)
191	swr4=reform(swr4)
192	;
193	file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc'
194	initncdf, file
195	lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box)
196	lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)
197	lwr1=reform(lwr1)
198	lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box)
199	lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)
200	lwr2=reform(lwr2)
201	lwr3=read_ncdf('lwr',st3-1,en3,file=file,/nostr, box=box)
202	lwr3=grossemoyenne(lwr3, "xy", box=box,/nan)
203	lwr3=reform(lwr3)
204	lwr4=read_ncdf('lwr',st4-1,en4,file=file,/nostr, box=box)
205	lwr4=grossemoyenne(lwr4, "xy", box=box,/nan)
206	lwr4=reform(lwr4)
207	;
208	file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc'
209	initncdf, file
210	lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)
211	lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)
212	lhf1=reform(lhf1)
213	lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)
214	lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)
215	lhf2=reform(lhf2)
216	lhf3=-1*read_ncdf('lhf',st3,en3,file=file,/nostr, box=box)
217	lhf3=grossemoyenne(lhf3, "xy", box=box,/nan)
218	lhf3=reform(lhf3)
219	lhf4=-1*read_ncdf('lhf',st4,en4,file=file,/nostr, box=box)
220	lhf4=grossemoyenne(lhf4, "xy", box=box,/nan)
221	lhf4=reform(lhf4)
222	;
223	file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc'
224	initncdf, file
225	shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box)
226	shf1=grossemoyenne(shf1, "xy", box=box,/nan)
227	shf1=reform(shf1)
228	shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box)
229	shf2=grossemoyenne(shf2, "xy", box=box,/nan)
230	shf2=reform(shf2)
231	shf3=-1*read_ncdf('shf',st3,en3,file=file,/nostr, box=box)
232	shf3=grossemoyenne(shf3, "xy", box=box,/nan)
233	shf3=reform(shf3)
234	shf4=-1*read_ncdf('shf',st4,en4,file=file,/nostr, box=box)
235	shf4=grossemoyenne(shf4, "xy", box=box,/nan)
236	shf4=reform(shf4)
237	;
238	net1=swr1+lwr1+lhf1+shf1
239	net2=swr2+lwr2+lhf2+shf2
240	net3=swr3+lwr3+lhf3+shf3
241	net4=swr4+lwr4+lhf4+shf4
242	era=[net1,net2,net3,net4]
243	;
244	help, era
245	swr1=0
246	swr2=0
247	swr3=0
248	swr4=0
249	lwr1=0
250	lwr2=0
251	lwr3=0
252	lwr4=0
253	lhf1=0
254	lhf2=0
255	lhf3=0
256	lhf4=0
257	shf1=0
258	shf2=0
259	shf3=0
260	shf4=0
261	net1=0
262	net2=0
263	net3=0
264	net4=0
265	;
266	file=project_id_env+"TropFlux_19890101_20091231.nc"
267	initncdf, file
268	swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)
269	swr1=grossemoyenne(swr1, "xy",/nan, box=box)
270	swr1=reform(swr1)
271	swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)
272	swr2=grossemoyenne(swr2, "xy",/nan, box=box)
273	swr2=reform(swr2)
274	swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box)
275	swr3=grossemoyenne(swr3, "xy",/nan, box=box)
276	swr3=reform(swr3)
277	swr4=read_ncdf('swr',st4,en4,file=file,/nostr, box=box)
278	swr4=grossemoyenne(swr4, "xy",/nan, box=box)
279	swr4=reform(swr4)
280	;
281	lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box)
282	lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
283	lwr1=reform(lwr1)
284	lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box)
285	lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
286	lwr2=reform(lwr2)
287	lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box)
288	lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)
289	lwr3=reform(lwr3)
290	lwr4=read_ncdf('lwr',st4,en4,file=file,/nostr, box=box)
291	lwr4=grossemoyenne(lwr4, "xy",/nan, box=box)
292	lwr4=reform(lwr4)
293	;
294	lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)
295	lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
296	lhf1=reform(lhf1)
297	lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)
298	lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
299	lhf2=reform(lhf2)
300	lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box)
301	lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)
302	lhf3=reform(lhf3)
303	lhf4=read_ncdf('lhf',st4,en4,file=file,/nostr, box=box)
304	lhf4=grossemoyenne(lhf4, "xy",/nan, box=box)
305	lhf4=reform(lhf4)
306	;
307	shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box)
308	shf1=grossemoyenne(shf1, "xy",/nan, box=box)
309	shf1=reform(shf1)
310	shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box)
311	shf2=grossemoyenne(shf2, "xy",/nan, box=box)
312	shf2=reform(shf2)
313	shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box)
314	shf3=grossemoyenne(shf3, "xy",/nan, box=box)
315	shf3=reform(shf3)
316	shf4=read_ncdf('shf',st4,en4,file=file,/nostr, box=box)
317	shf4=grossemoyenne(shf4, "xy",/nan, box=box)
318	shf4=reform(shf4)
319	;
320	net1=swr1+lwr1+lhf1+shf1
321	net2=swr2+lwr2+lhf2+shf2
322	net3=swr3+lwr3+lhf3+shf3
323	net4=swr4+lwr4+lhf4+shf4
324	;
325	file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc"
326	initncdf, file
327	olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box)
328	olr1=grossemoyenne(olr1, "xy",/nan, box=box)
329	olr1=reform(olr1)
330	olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box)
331	olr2=grossemoyenne(olr2, "xy",/nan, box=box)
332	olr2=reform(olr2)
333	olr3=read_ncdf("sw", st3-1, en3, file=file,/nostr, box=box)
334	olr3=grossemoyenne(olr3, "xy",/nan, box=box)
335	olr3=reform(olr3)
336	olr4=read_ncdf("sw", st4-1, en4, file=file,/nostr, box=box)
337	olr4=grossemoyenne(olr4, "xy",/nan, box=box)
338	olr4=reform(olr4)
339	;
340	neto_1=olr1+lwr1+lhf1+shf1
341	neto_2=olr2+lwr2+lhf2+shf2
342	neto_3=olr3+lwr3+lhf3+shf3
343	neto_4=olr4+lwr4+lhf4+shf4
344	;
345	trop=[net1,net2,net3,net4]
346	olr=[neto_1,neto_2,neto_3,neto_4]
347	help, trop, olr
348	;
349	swr1=0
350	swr2=0
351	swr3=0
352	swr4=0
353	lwr1=0
354	lwr2=0
355	lwr3=0
356	lwr4=0
357	lhf1=0
358	lhf2=0
359	lhf3=0
360	lhf4=0
361	shf1=0
362	shf2=0
363	shf3=0
364	shf4=0
365	net1=0
366	net2=0
367	net3=0
368	net4=0
369	neto_1=0
370	neto_2=0
371	neto_3=0
372	neto_4=0
373	;
374	file=project_id_env+'swr_oafluxgrid_1985_2007.nc'
375	initncdf, file
376	swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box)
377	swr1=grossemoyenne(swr1, "xy",/nan, box=box)
378	swr1=reform(swr1)
379	swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box)
380	swr2=grossemoyenne(swr2, "xy",/nan, box=box)
381	swr2=reform(swr2)
382	swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box)
383	swr3=grossemoyenne(swr3, "xy",/nan, box=box)
384	swr3=reform(swr3)
385	swr4=read_ncdf("swr", st4, en4, file=file,/nostr, box=box)
386	swr4=grossemoyenne(swr4, "xy",/nan, box=box)
387	swr4=reform(swr4)
388	;
389	file=project_id_env+'lwr_oafluxgrid_1985_2007.nc'
390	initncdf, file
391	lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box)
392	lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
393	lwr1=reform(lwr1)
394	lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box)
395	lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
396	lwr2=reform(lwr2)
397	lwr3=-1*read_ncdf("lwr", st3, en3, file=file,/nostr, box=box)
398	lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)
399	lwr3=reform(lwr3)
400	lwr4=-1*read_ncdf("lwr", st4, en4, file=file,/nostr, box=box)
401	lwr4=grossemoyenne(lwr4, "xy",/nan, box=box)
402	lwr4=reform(lwr4)
403	;
404	file=project_id_env+'lhf_oafluxgrid_1985_2009.nc'
405	initncdf, file
406	lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box)
407	lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
408	lhf1=reform(lhf1)
409	lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box)
410	lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
411	lhf2=reform(lhf2)
412	lhf3=-1*read_ncdf("lhf", st3, en3, file=file,/nostr, box=box)
413	lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)
414	lhf3=reform(lhf3)
415	lhf4=-1*read_ncdf("lhf", st4, en4, file=file,/nostr, box=box)
416	lhf4=grossemoyenne(lhf4, "xy",/nan, box=box)
417	lhf4=reform(lhf4)
418	;
419	file=project_id_env+'shf_oafluxgrid_1985_2009.nc'
420	initncdf, file
421	shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box)
422	shf1=grossemoyenne(shf1, "xy",/nan, box=box)
423	shf1=reform(shf1)
424	shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box)
425	shf2=grossemoyenne(shf2, "xy",/nan, box=box)
426	shf2=reform(shf2)
427	shf3=-1*read_ncdf("shf", st3, en3, file=file,/nostr, box=box)
428	shf3=grossemoyenne(shf3, "xy",/nan, box=box)
429	shf3=reform(shf3)
430	shf4=-1*read_ncdf("shf", st4, en4, file=file,/nostr, box=box)
431	shf4=grossemoyenne(shf4, "xy",/nan, box=box)
432	shf4=reform(shf4)
433	;
434	net1=swr1+lwr1+lhf1+shf1
435	net2=swr2+lwr2+lhf2+shf2
436	net3=swr3+lwr3+lhf3+shf3
437	net4=swr4+lwr4+lhf4+shf4
438	;
439	oaf=[net1,net2,net3,net4]
440	;
441	help, oaf
442	swr1=0
443	swr2=0
444	swr3=0
445	swr4=0
446	lwr1=0
447	lwr2=0
448	lwr3=0
449	lwr4=0
450	lhf1=0
451	lhf2=0
452	lhf3=0
453	lhf4=0
454	shf1=0
455	shf2=0
456	shf3=0
457	shf4=0
458	net1=0
459	net2=0
460	net3=0
461	net4=0
462	;
463	fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc'
464	initncdf, fi
465	swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box)
466	swr1=grossemoyenne(swr1, "xy",/nan, box=box)
467	swr1=reform(swr1)
468	swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box)
469	swr2=grossemoyenne(swr2, "xy",/nan, box=box)
470	swr2=reform(swr2)
471	swr3=read_ncdf("swr", st3, en3, file=fi,/nostr, box=box)
472	swr3=grossemoyenne(swr3, "xy",/nan, box=box)
473	swr3=reform(swr3)
474	swr4=read_ncdf("swr", st4, en4, file=fi,/nostr, box=box)
475	swr4=grossemoyenne(swr4, "xy",/nan, box=box)
476	swr4=reform(swr4)
477	;
478	fi=project_id_env+"NCEP2_flux_19890101_20090729.nc"
479	initncdf, fi
480	lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box)
481	lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
482	lwr1=reform(lwr1)
483	lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box)
484	lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
485	lwr2=reform(lwr2)
486	lwr3=-1*read_ncdf("lwr", st3, en3, file=fi,/nostr, box=box)
487	lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)
488	lwr3=reform(lwr3)
489	lwr4=-1*read_ncdf("lwr", st4, en4, file=fi,/nostr, box=box)
490	lwr4=grossemoyenne(lwr4, "xy",/nan, box=box)
491	lwr4=reform(lwr4)
492	;
493	fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc'
494	initncdf, fi
495	lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box)
496	lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
497	lhf1=reform(lhf1)
498	lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box)
499	lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
500	lhf2=reform(lhf2)
501	lhf3=-1*read_ncdf("lhf", st3, en3, file=fi,/nostr, box=box)
502	lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)
503	lhf3=reform(lhf3)
504	lhf4=-1*read_ncdf("lhf", st4, en4, file=fi,/nostr, box=box)
505	lhf4=grossemoyenne(lhf4, "xy",/nan, box=box)
506	lhf4=reform(lhf4)
507	;
508	fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc'
509	initncdf, fi
510	shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box)
511	shf1=grossemoyenne(shf1, "xy",/nan, box=box)
512	shf1=reform(shf1)
513	shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box)
514	shf2=grossemoyenne(shf2, "xy",/nan, box=box)
515	shf2=reform(shf2)
516	shf3=-1*read_ncdf("shf", st3, en3, file=fi,/nostr, box=box)
517	shf3=grossemoyenne(shf3, "xy",/nan, box=box)
518	shf3=reform(shf3)
519	shf4=-1*read_ncdf("shf", st4, en4, file=fi,/nostr, box=box)
520	shf4=grossemoyenne(shf4, "xy",/nan, box=box)
521	shf4=reform(shf4)
522	;
523	net1=swr1+lwr1+lhf1+shf1
524	net2=swr2+lwr2+lhf2+shf2
525	net3=swr3+lwr3+lhf3+shf3
526	net4=swr4+lwr4+lhf4+shf4
527	;
528	nce2=[net1,net2,net3,net4]
529	;
530	help, nce2
531	swr1=0
532	swr2=0
533	swr3=0
534	swr4=0
535	lwr1=0
536	lwr2=0
537	lwr3=0
538	lwr4=0
539	lhf1=0
540	lhf2=0
541	lhf3=0
542	lhf4=0
543	shf1=0
544	shf2=0
545	shf3=0
546	shf4=0
547	net1=0
548	net2=0
549	net3=0
550	net4=0
551	;
552	file=project_id_env+"ncep1_flux_19890101_20091231.nc"
553	initncdf, file
554	swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box)
555	swr1=grossemoyenne(swr1, "xy",/nan, box=box)
556	swr1=reform(swr1)
557	swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box)
558	swr2=grossemoyenne(swr2, "xy",/nan, box=box)
559	swr2=reform(swr2)
560	swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box)
561	swr3=grossemoyenne(swr3, "xy",/nan, box=box)
562	swr3=reform(swr3)
563	swr4=read_ncdf("swr", st4, en4, file=file,/nostr, box=box)
564	swr4=grossemoyenne(swr4, "xy",/nan, box=box)
565	swr4=reform(swr4)
566	;
567	lwr1=read_ncdf("lwr", st1, en1, file=file,/nostr, box=box)
568	lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
569	lwr1=reform(lwr1)
570	lwr2=read_ncdf("lwr", st2, en2, file=file,/nostr, box=box)
571	lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
572	lwr2=reform(lwr2)
573	lwr3=read_ncdf("lwr", st3, en3, file=file,/nostr, box=box)
574	lwr3=grossemoyenne(lwr3, "xy",/nan, box=box)
575	lwr3=reform(lwr3)
576	lwr4=read_ncdf("lwr", st4, en4, file=file,/nostr, box=box)
577	lwr4=grossemoyenne(lwr4, "xy",/nan, box=box)
578	lwr4=reform(lwr4)
579	;
580	lhf1=read_ncdf("lhf", st1, en1, file=file,/nostr, box=box)
581	lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
582	lhf1=reform(lhf1)
583	lhf2=read_ncdf("lhf", st2, en2, file=file,/nostr, box=box)
584	lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
585	lhf2=reform(lhf2)
586	lhf3=read_ncdf("lhf", st3, en3, file=file,/nostr, box=box)
587	lhf3=grossemoyenne(lhf3, "xy",/nan, box=box)
588	lhf3=reform(lhf3)
589	lhf4=read_ncdf("lhf", st4, en4, file=file,/nostr, box=box)
590	lhf4=grossemoyenne(lhf4, "xy",/nan, box=box)
591	lhf4=reform(lhf4)
592	;
593	shf1=read_ncdf("shf", st1, en1, file=file,/nostr, box=box)
594	shf1=grossemoyenne(shf1, "xy",/nan, box=box)
595	shf1=reform(shf1)
596	shf2=read_ncdf("shf", st2, en2, file=file,/nostr, box=box)
597	shf2=grossemoyenne(shf2, "xy",/nan, box=box)
598	shf2=reform(shf2)
599	shf3=read_ncdf("shf", st3, en3, file=file,/nostr, box=box)
600	shf3=grossemoyenne(shf3, "xy",/nan, box=box)
601	shf3=reform(shf3)
602	shf4=read_ncdf("shf", st4, en4, file=file,/nostr, box=box)
603	shf4=grossemoyenne(shf4, "xy",/nan, box=box)
604	shf4=reform(shf4)
605	;
606	net1=swr1+lwr1+lhf1+shf1
607	net2=swr2+lwr2+lhf2+shf2
608	net3=swr3+lwr3+lhf3+shf3
609	net4=swr4+lwr4+lhf4+shf4
610	;
611	nce1=[net1,net2,net3,net4]
612	nce1=nce1
613	;
614	help, nce1
615	swr1=0
616	swr2=0
617	swr3=0
618	swr4=0
619	lwr1=0
620	lwr2=0
621	lwr3=0
622	lwr4=0
623	lhf1=0
624	lhf2=0
625	lhf3=0
626	lhf4=0
627	shf1=0
628	shf2=0
629	shf3=0
630	shf4=0
631	net1=0
632	net2=0
633	net3=0
634	net4=0
635	;
636	ind=where(finite(net))
637	net=net(ind)
638	era=era(ind)
639	trop=trop(ind)
640	oaf=oaf(ind)
641	nce1=nce1(ind)
642	nce2=nce2(ind)
643	olr=olr(ind)
644	;
645	help, net, era, trop, oaf, nce1, nce2, olr
646	;
647	statistics, net, era, $
648	cor, bias, std, rmsd
649	print, cor, bias, std, rmsd
650	cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
651	;
652	statistics, net, trop, $
653	cor, bias, std, rmsd
654	print, cor, bias, std, rmsd
655	cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
656	;
657	statistics, net, oaf, $
658	cor, bias, std, rmsd
659	print, cor, bias, std, rmsd
660	cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
661	;
662	statistics, net, nce1, $
663	cor, bias, std, rmsd
664	print, cor, bias, std, rmsd
665	cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
666	;
667	statistics, net, nce2, $
668	cor, bias, std, rmsd
669	print, cor, bias, std, rmsd
670	cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
671	;
672	statistics, net, olr, $
673	cor, bias, std, rmsd
674	print, cor, bias, std, rmsd
675	cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
676	;
677	splot, net, era, charsize=1.1, title='Net Flux - COARE Vs ERAI', $
678	xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1
679	oplot, [-200,200], [-200,200]
680	xyouts, -150,-175, cstat_era, charsize=1.
681	;
682	splot, net, trop, charsize=1.1, title='Net Flux - COARE Vs TropFlux',/noer, $
683	xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1
684	oplot, [-200,200], [-200,200]
685	xyouts, -150,-175, cstat_trop, charsize=1.
686	;
687	splot, net, oaf, charsize=1.1, title='Net Flux - COARE Vs OAFlux',/noer, $
688	xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1
689	oplot, [-200,200], [-200,200]
690	xyouts, -150,-175, cstat_oaf, charsize=1.
691	;
692	splot, net, nce1, charsize=1.1, title='Net Flux - COARE Vs NCEP1',/noer, $
693	xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1
694	oplot, [-200,200], [-200,200]
695	xyouts, -150,-175, cstat_nce1, charsize=1.
696	;
697	splot, net, nce2, charsize=1.1, title='Net Flux - COARE Vs NCEP2',/noer, $
698	xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1
699	oplot, [-200,200], [-200,200]
700	xyouts, -150,-175, cstat_nce2, charsize=1.
701	;
702	splot, net, olr, charsize=1.1, title='Net Flux - COARE Vs TropFlux_NRT',/noer, $
703	xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1
704	oplot, [-200,200], [-200,200]
705	xyouts, -150,-175, cstat_olr, charsize=1.
706	;
707	closeps
708	;
709	end

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

Download in other formats:

Original Format