source: trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro @ 50

Last change on this file since 50 was 50, checked in by pinsard, 13 years ago

replace TROPFLUX by PROJECT
File size: 13.7 KB
Line 
1;+
2; .. _subdctn_ses_validation_net_1991_93_all_v10.pro:
3;
4; ==============================================
5; subdctn_ses_validation_net_1991_93_all_v10.pro
6; ==============================================
7;
8; DESCRIPTION
9; ===========
10;
11; .. graphviz::
12;
13;    digraph subdctn_ses_validation_net_1991_93_all_v10 {
14;       graph [
15;       rankdir="LR",
16;       ]
17;       met_data [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met_data_subduct_ses_9*.txt"];
18;       swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_swr_19910101_20091231_oafluxgrid.nc"];
19;       lwr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lwr_19910101_20091231_oafluxgrid.nc"];
20;       lhf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lhf_19910101_20091231_oafluxgrid.nc"];
21;       shf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_shf_19910101_20091231_oafluxgrid.nc"];
22;       tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"];
23;       swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"];
24;       swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_oafluxgrid_1985_2007.nc"];
25;       lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"];
26;       lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2007.nc"];
27;       shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2007.nc"];
28;       swr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"];
29;       lwr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"];
30;       lhf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"];
31;       shf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"];
32;
33;       ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ncep1_flux_19890101_20091231.nc"];
34;
35;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/subdctn_ses_validation_net_1991_93_all_v10.ps"];
36;
37;       subdctn_ses_validation_net_1991_93_all_v10 [shape=box,
38;       fontname=Courier,
39;       color=blue,
40;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro",
41;       label="${TROPFLUX}/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro"];
42;
43;       {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} -> {subdctn_ses_validation_net_1991_93_all_v10} -> {figure}
44;    }
45;
46; SEE ALSO
47; ========
48;
49; :ref:`project_profile.sh`
50; :ref:`project_init.pro`
51; :ref:`cm_project.pro`
52;
53; :ref:`statistics.pro`
54;
55; EXAMPLES
56; ========
57;
58; ::
59;
60;  IDL> subdctn_ses_validation_net_1991_93_all_v10
61;
62; TODO
63; ====
64;
65; make it work on cratos : missing data
66;
67; coding rules
68;
69; get rid of spwan
70;
71; complete description
72;
73; handle IO error
74;
75; EVOLUTIONS
76; ==========
77;
78; $Id$
79;
80; - fplod 20110420T121359Z aedon.locean-ipsl.upmc.fr (Darwin)
81;
82;   * remove hard coding path
83;   * add graphviz
84;
85; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
86;
87;   * minimal header
88;
89;-
90pro subdctn_ses_validation_net_1991_93_all_v10
91@cm_general
92@cm_project
93reinitplt, /z,/invert
94key_portrait = 1
95coefpalit=.9
96
97openps, FILENAME = 'idl.ps'
98; partie a changer
99
100;; subduction experiment in the atlantic in 1991-93
101box=[337.5, 338.5, 17.5, 18.5]
102
103;; first leg
104
105;; second leg
106
107fi=project_id_env+"met_data_subduct_ses_92.txt"
108res=read_ascii(fi,data_start=1)
109ff=res.field1
110swr_92=reform(ff(7,*))
111lhf_92=reform(ff(5,*))
112shf_92=reform(ff(6,*))
113lwr_92=reform(ff(8,*))
114net_92=[swr_92+lwr_92+shf_92+lhf_92]
115help, net_92
116
117;; second leg
118
119fi=project_id_env+"met_data_subduct_ses_93.txt"
120res=read_ascii(fi,data_start=1)
121ff=res.field1
122swr_93=reform(ff(7,*))
123lhf_93=reform(ff(5,*))
124shf_93=reform(ff(6,*))
125lwr_93=reform(ff(8,*))
126net_93=[swr_93+lwr_93+shf_93+lhf_93]
127help, net_93
128
129net=[net_92, net_93]
130help, net
131
132;; reading other data sets
133
134st1=19920209 & en1=19920911
135st2=19920921 & en2=19930618
136
137file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc"
138initncdf, file
139swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1)
140swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2)
141
142file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc'
143initncdf, file
144lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1)
145lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2)
146
147file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc'
148initncdf, file
149shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1)
150shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2)
151
152file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc'
153initncdf, file
154lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1)
155lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2)
156
157net1=swr1+(lwr1+shf1+lhf1) & net2=swr2+(lwr2+shf2+lhf2)
158era=[net1,net2]
159help, era
160
161lhf1=0. & lhf2=0. & shf1=0. & shf2=0.
162lwr1=0. & lwr2=0. & swr1=0. & swr2=0.
163net1=0. & net2=0.
164
165file=project_id_env+"TropFlux_19890101_20091231.nc"
166initncdf, file
167lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1)
168lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2)
169shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1)
170shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2)
171lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1)
172lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2)
173swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1)
174swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2)
175
176file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc"
177initncdf, file
178olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) & olr1=grossemoyenne(olr1, "xy",/nan, box=box) & olr1=reform(olr1)
179olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box) & olr2=grossemoyenne(olr2, "xy",/nan, box=box) & olr2=reform(olr2)
180
181net1=swr1+lwr1+(shf1+lhf1) & net2=swr2+lwr2+(shf2+lhf2)
182neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2
183trop=[net1,net2] & olr=[neto_1,neto_2]
184help, trop, olr
185
186file=project_id_env+'swr_oafluxgrid_1985_2007.nc'
187initncdf, file
188swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy",/nan, box=box) & swr1=reform(swr1)
189swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy",/nan, box=box) & swr2=reform(swr2)
190
191file=project_id_env+'lwr_oafluxgrid_1985_2007.nc'
192initncdf, file
193lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) & lwr1=reform(lwr1)
194lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) & lwr2=reform(lwr2)
195
196file=project_id_env+'lhf_oafluxgrid_1985_2009.nc'
197initncdf, file
198lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) & lhf1=reform(lhf1)
199lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) & lhf2=reform(lhf2)
200
201file=project_id_env+'shf_oafluxgrid_1985_2009.nc'
202initncdf, file
203shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy",/nan, box=box) & shf1=reform(shf1)
204shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy",/nan, box=box) & shf2=reform(shf2)
205
206net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2
207oaf=[net1,net2]
208help, oaf
209
210fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc'
211initncdf, fi
212swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy",/nan, box=box) & swr1=reform(swr1)
213swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy",/nan, box=box) & swr2=reform(swr2)
214
215fi=project_id_env+"NCEP2_flux_19890101_20090729.nc"
216initncdf, fi
217lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) & lwr1=reform(lwr1)
218lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) & lwr2=reform(lwr2)
219
220fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc'
221initncdf, fi
222lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) & lhf1=reform(lhf1)
223lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) & lhf2=reform(lhf2)
224
225fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc'
226initncdf, fi
227shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy",/nan, box=box) & shf1=reform(shf1)
228shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy",/nan, box=box) & shf2=reform(shf2)
229
230net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2
231nce2=[net1,net2]
232help, nce2
233
234file=project_id_env+"ncep1_flux_19890101_20091231.nc"
235initncdf, file
236lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1)
237lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2)
238shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1)
239shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2)
240lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1)
241lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2)
242swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1)
243swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2)
244net1=swr1+(lwr1+shf1+lhf1) & net2=swr2+(lwr2+shf2+lhf2)
245nce1=[net1,net2]
246
247help, nce1
248
249ind=where(finite(net))
250net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind)
251olr=olr(ind)
252
253statistics, net, era, $
254                 cor, bias, std, rmsd
255print, cor, bias, std,  rmsd
256cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
257
258statistics, net, trop, $
259                 cor, bias, std, rmsd
260print, cor, bias, std,  rmsd
261cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
262
263statistics, net, oaf, $
264                 cor, bias, std, rmsd
265print, cor, bias, std,  rmsd
266cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
267
268statistics, net, nce1, $
269                 cor, bias, std, rmsd
270print, cor, bias, std,  rmsd
271cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
272
273statistics, net, nce2, $
274                 cor, bias, std, rmsd
275print, cor, bias, std,  rmsd
276cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
277
278statistics, net, olr, $
279                 cor, bias, std, rmsd
280print, cor, bias, std,  rmsd
281cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
282
283
284splot, net, era, charsize=1.1, title='Net Flux - Subdn. ses Vs ERAI', $
285     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1
286oplot, [-200,200], [-200,200]
287xyouts, 0,-175, cstat_era, charsize=1.
288
289splot, net, trop, charsize=1.1, title='Net Flux - Subdn. ses Vs TropFlux',/noer, $
290     xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1
291oplot, [-200,200], [-200,200]
292xyouts, 0,-175, cstat_trop, charsize=1.
293
294splot, net, oaf, charsize=1.1, title='Net Flux - Subdn. ses Vs OAFlux',/noer, $
295     xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1
296oplot, [-200,200], [-200,200]
297xyouts, 0,-175, cstat_oaf, charsize=1.
298
299splot, net, nce1, charsize=1.1, title='Net Flux - Subdn. ses Vs NCEP1',/noer, $
300     xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1
301oplot, [-200,200], [-200,200]
302xyouts, 0,-175, cstat_nce1, charsize=1.
303
304splot, net, nce2, charsize=1.1, title='Net Flux - Subdn. ses Vs NCEP2',/noer, $
305     xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1
306oplot, [-200,200], [-200,200]
307xyouts, 0,-175, cstat_nce2, charsize=1.
308
309splot, net, olr, charsize=1.1, title='Net Flux - Subdn. ses Vs TropFlux_NRT',/noer, $
310     xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1
311oplot, [-200,200], [-200,200]
312xyouts, 0,-175, cstat_olr, charsize=1.
313
314closeps
315fig=project_od_env+"subdctn_ses_validation_net_1991_93_all_v10.ps"
316spawn, 'mv '+psdir+'idl.ps '+fig
317spawn, 'gv '+fig
318end
Note: See TracBrowser for help on using the repository browser.