source: trunk/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro @ 182

Last change on this file since 182 was 182, checked in by pinsard, 12 years ago

fix some svn propset
  • Property svn:keywords set to Id URL
File size: 13.6 KB
Line 
1;+
2; .. _swr_statistics_map_2000_2009_v50.pro:
3;
4; ====================================
5; swr_statistics_map_2000_2009_v50.pro
6; ====================================
7;
8; DESCRIPTION
9; ===========
10;
11; .. graphviz::
12;
13;    digraph swr_statistics_map_2000_2009_v50 {
14;
15;       swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_erai_v50.txt"];
16;       swr_trop [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_trop_v50.txt"];
17;       swr_oaflx [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_oaflx_v50.txt"];
18;       swr_olr [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_olr_v50.txt"];
19;       mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/longwave_IO_mask.nc"];
20;
21;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/swr_statistics_map_2000_2009_v50.ps"];
22;
23;       swr_statistics_map_2000_2009_v50 [shape=box,
24;       fontname=Courier,
25;       color=blue,
26;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro",
27;       label="${TROPFLUX}/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro"];
28;
29;       {swr_erai swr_trop swr_oaflx swr_olr mask} -> {swr_statistics_map_2000_2009_v50} -> {figure}
30;    }
31;
32; SEE ALSO
33; ========
34;
35; :ref:`project_profile.sh`
36; :ref:`project_init.pro`
37; :ref:`cm_project.pro`
38;
39; EXAMPLES
40; ========
41;
42; ::
43;
44;  swr_statistics_map_2000_2009_v50
45;
46; TODO
47; ====
48;
49; make it work on cratos : missing data
50;
51; coding rules
52;
53; complete description
54;
55; handle IO error
56;
57; EVOLUTIONS
58; ==========
59;
60; $Id$
61;
62; $URL$
63;
64; - fplod 20110420T113846Z aedon.locean-ipsl.upmc.fr (Darwin)
65;
66;   * remove stop
67;   * remove hard coding path
68;   * add graphviz
69;
70; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
71;
72;   * minimal header
73;
74;-
75pro swr_statistics_map_2000_2009_v50
76@cm_general
77@cm_project
78reinitplt, /z,/invert
79key_portrait = 1
80coefpalit=.9
81;
82openps, FILENAME = project_od_env+'swr_statistics_map_2000_2009_v50.ps'
83; partie a changer
84bias_mi=-20
85bias_ma=20
86bias_int=2
87std_mi=0.7
88std_ma=1.31
89std_int=0.05
90rmsd_mi=10
91rmsd_ma=30
92rmsd_int=1.25
93cor_mi=0.5
94cor_ma=1.
95cor_int=0.02
96fi_swr_erai=project_id_env+'swr_2000_2009_erai_v50.txt'
97fi_swr_trop=project_id_env+'swr_2000_2009_trop_v50.txt'
98fi_swr_oaflx=project_id_env+'swr_2000_2009_oaflx_v50.txt'
99fi_swr_olr=project_id_env+'swr_2000_2009_olr_v50.txt'
100;
101res=read_ascii(fi_swr_erai,data_start=1)
102ff=res.field1
103lat=reform(ff(0,*))
104lon=reform(ff(1,*))
105cor_era=reform(ff(2,*))
106bias_era=reform(ff(3,*))
107std_era=reform(ff(4,*))
108rmsd_era=reform(ff(5,*))
109;
110ind=where(bias_era ge bias_ma)
111bias_era(ind)=bias_ma-0.5
112ind=where(rmsd_era ge rmsd_ma)
113rmsd_era(ind)=rmsd_ma-0.5
114;
115res=read_ascii(fi_swr_trop,data_start=1)
116ff=res.field1
117lat=reform(ff(0,*))
118lon=reform(ff(1,*))
119cor_trop=reform(ff(2,*))
120bias_trop=reform(ff(3,*))
121std_trop=reform(ff(4,*))
122rmsd_trop=reform(ff(5,*))
123ind=where(rmsd_trop ge rmsd_ma)
124rmsd_trop(ind)=rmsd_ma-0.5
125;
126res=read_ascii(fi_swr_oaflx,data_start=1)
127ff=res.field1
128lat=reform(ff(0,*))
129lon=reform(ff(1,*))
130cor_oaflx=reform(ff(2,*))
131bias_oaflx=reform(ff(3,*))
132std_oaflx=reform(ff(4,*))
133rmsd_oaflx=reform(ff(5,*))
134ind=where(rmsd_oaflx ge rmsd_ma)
135rmsd_oaflx(ind)=rmsd_ma-0.5
136;
137res=read_ascii(fi_swr_olr,data_start=1)
138ff=res.field1
139lat=reform(ff(0,*))
140lon=reform(ff(1,*))
141cor_olr=reform(ff(2,*))
142bias_olr=reform(ff(3,*))
143std_olr=reform(ff(4,*))
144rmsd_olr=reform(ff(5,*))
145ind=where(rmsd_olr ge rmsd_ma)
146rmsd_olr(ind)=rmsd_ma-0.5
147;
148;ind=where(std_olr ge std_ma)
149;std_olr(ind)=std_ma-0.01
150file=project_id_env+'longwave_IO_mask.nc'
151initncdf, file
152domdef, 30,390,-30,30
153msk=ncdf_lec(file,var='msk')
154;box=[0,290,-30,30]
155marge=[0,0,-0.8,0]
156marge1=[0,0,-5,0]
157plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=64,cor_mi, cor_ma,int=cor_int, $
158   title='1) SWR Correlation - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge
159;
160NN=n_elements(lat)
161;
162usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
163;
164for n=0,NN-1 do begin
165  x=lon(n)
166  y=lat(n)
167  c=cor_trop(n)
168  cmi=cor_mi
169  cma=cor_ma
170  dc=cma-cmi
171  col=((10+244*(c-cmi)/dc) > 10) < 254
172  plots, x,y,psym=8,symsize=1.5,color=col
173endfor
174;
175usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
176;
177for n=0,NN-1 do begin
178  x=lon(n)
179  y=lat(n)
180  c=cor_trop(n)
181;
182  plots, x,y,psym=8,symsize=1.5,color=0
183endfor
184plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $
185   title='2) SWR Correlation -  ERAI', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge
186;
187NN=n_elements(lat)
188;
189usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
190;
191for n=0,NN-1 do begin
192  x=lon(n)
193  y=lat(n)
194  c=cor_era(n)
195  cmi=cor_mi
196  cma=cor_ma
197  dc=cma-cmi
198  col=((10+244*(c-cmi)/dc) > 10) < 254
199  plots, x,y,psym=8,symsize=1.5,color=col
200endfor
201;
202usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
203;
204for n=0,NN-1 do begin
205  x=lon(n)
206  y=lat(n)
207  c=cor_era(n)
208;
209  plots, x,y,psym=8,symsize=1.5,color=0
210endfor
211plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $
212   title='3) SWR Correlation - OAFlux', subtitle='', small=[1,4,3],/noer,/rempl, marge=marge
213;
214NN=n_elements(lat)
215;
216usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
217;
218for n=0,NN-1 do begin
219  x=lon(n)
220  y=lat(n)
221  c=cor_oaflx(n)
222  cmi=cor_mi
223  cma=cor_ma
224  dc=cma-cmi
225  col=((10+244*(c-cmi)/dc) > 10) < 254
226  plots, x,y,psym=8,symsize=1.5,color=col
227endfor
228;
229usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
230;
231for n=0,NN-1 do begin
232  x=lon(n)
233  y=lat(n)
234  c=cor_oaflx(n)
235;
236  plots, x,y,psym=8,symsize=1.5,color=0
237endfor
238plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int, $
239   title='4) SWR Correlation - SWR_olr', subtitle='', small=[1,4,4],/noer,/rempl, marge=marge1
240;
241NN=n_elements(lat)
242;
243usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
244;
245for n=0,NN-1 do begin
246  x=lon(n)
247  y=lat(n)
248  c=cor_olr(n)
249  cmi=cor_mi
250  cma=cor_ma
251  dc=cma-cmi
252  col=((10+244*(c-cmi)/dc) > 10) < 254
253  plots, x,y,psym=8,symsize=1.5,color=col
254endfor
255;
256usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
257;
258for n=0,NN-1 do begin
259  x=lon(n)
260  y=lat(n)
261  c=cor_olr(n)
262;
263  plots, x,y,psym=8,symsize=1.5,color=0
264endfor
265;
266erase
267;
268plt, -5+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int, marge=marge,/nocolorb, $
269   title='1) SWR Mean bias - TropFlux', subtitle='', small=[1,4,1],/rempl,format='(i3)'
270;
271NN=n_elements(lat)
272;
273usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
274;
275for n=0,NN-1 do begin
276  x=lon(n)
277  y=lat(n)
278  c=bias_trop(n)
279  cmi=bias_mi
280  cma=bias_ma
281  dc=cma-cmi
282  col=((10+244*(c-cmi)/dc) > 10) < 254
283  plots, x,y,psym=8,symsize=1.5,color=col
284endfor
285;
286usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
287;
288for n=0,NN-1 do begin
289  x=lon(n)
290  y=lat(n)
291  c=bias_trop(n)
292;
293  plots, x,y,psym=8,symsize=1.5,color=0
294endfor
295plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $
296   title='2) SWR Mean bias - ERAI', subtitle='', small=[1,4,2],/rempl,format='(i3)',/nocolorb
297;
298NN=n_elements(lat)
299;
300usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
301;
302for n=0,NN-1 do begin
303  x=lon(n)
304  y=lat(n)
305  c=bias_era(n)
306  cmi=bias_mi
307  cma=bias_ma
308  dc=cma-cmi
309  col=((10+244*(c-cmi)/dc) > 10) < 254
310  plots, x,y,psym=8,symsize=1.5,color=col
311endfor
312;
313usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
314;
315for n=0,NN-1 do begin
316  x=lon(n)
317  y=lat(n)
318  c=bias_era(n)
319;
320  plots, x,y,psym=8,symsize=1.5,color=0
321endfor
322;
323plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $
324   title='3) SWR Mean bias - OAFlux', subtitle='', small=[1,4,3],/rempl,format='(i3)',/nocolorb
325;
326NN=n_elements(lat)
327;
328usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
329;
330for n=0,NN-1 do begin
331  x=lon(n)
332  y=lat(n)
333  c=bias_oaflx(n)
334  cmi=bias_mi
335  cma=bias_ma
336  dc=cma-cmi
337  col=((10+244*(c-cmi)/dc) > 10) < 254
338  plots, x,y,psym=8,symsize=1.5,color=col
339endfor
340;
341usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
342;
343for n=0,NN-1 do begin
344  x=lon(n)
345  y=lat(n)
346  c=bias_oaflx(n)
347;
348  plots, x,y,psym=8,symsize=1.5,color=0
349endfor
350;
351plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge1, $
352   title='4) SWR Mean bias - SWR_olr', subtitle='', small=[1,4,4],/rempl ;,format='(i4)'
353;
354NN=n_elements(lat)
355;
356usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
357;
358for n=0,NN-1 do begin
359  x=lon(n)
360  y=lat(n)
361  c=bias_olr(n)
362  cmi=bias_mi
363  cma=bias_ma
364  dc=cma-cmi
365  col=((10+244*(c-cmi)/dc) > 10) < 254
366  plots, x,y,psym=8,symsize=1.5,color=col
367endfor
368;
369usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
370;
371for n=0,NN-1 do begin
372  x=lon(n)
373  y=lat(n)
374  c=bias_olr(n)
375;
376  plots, x,y,psym=8,symsize=1.5,color=0
377endfor
378;
379erase
380;
381plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $
382   title='1) RMSD - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge
383;
384NN=n_elements(lat)
385;
386usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
387;
388for n=0,NN-1 do begin
389  x=lon(n)
390  y=lat(n)
391  c=rmsd_trop(n)
392  cmi=rmsd_mi
393  cma=rmsd_ma
394  dc=cma-cmi
395  col=((10+244*(c-cmi)/dc) > 10) < 254
396  plots, x,y,psym=8,symsize=1.5,color=col
397endfor
398;
399usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
400;
401for n=0,NN-1 do begin
402  x=lon(n)
403  y=lat(n)
404  c=rmsd_trop(n)
405;
406  plots, x,y,psym=8,symsize=1.5,color=0
407endfor
408plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $
409   title='2) RMSD - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge
410;
411NN=n_elements(lat)
412;
413usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
414;
415for n=0,NN-1 do begin
416  x=lon(n)
417  y=lat(n)
418  c=rmsd_era(n)
419  cmi=rmsd_mi
420  cma=rmsd_ma
421  dc=cma-cmi
422  col=((10+244*(c-cmi)/dc) > 10) < 254
423  plots, x,y,psym=8,symsize=1.5,color=col
424endfor
425;
426usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
427;
428for n=0,NN-1 do begin
429  x=lon(n)
430  y=lat(n)
431  c=rmsd_era(n)
432;
433  plots, x,y,psym=8,symsize=1.5,color=0
434endfor
435plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $
436   title='3) RMSD - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge
437;
438NN=n_elements(lat)
439;
440usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
441;
442for n=0,NN-1 do begin
443  x=lon(n)
444  y=lat(n)
445  c=rmsd_oaflx(n)
446  cmi=rmsd_mi
447  cma=rmsd_ma
448  dc=cma-cmi
449  col=((10+244*(c-cmi)/dc) > 10) < 254
450  plots, x,y,psym=8,symsize=1.5,color=col
451endfor
452;
453usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
454;
455for n=0,NN-1 do begin
456  x=lon(n)
457  y=lat(n)
458  c=rmsd_oaflx(n)
459;
460  plots, x,y,psym=8,symsize=1.5,color=0
461endfor
462;
463plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $
464   title='4) RMSD - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1
465;
466NN=n_elements(lat)
467;
468usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
469;
470for n=0,NN-1 do begin
471  x=lon(n)
472  y=lat(n)
473  c=rmsd_olr(n)
474  cmi=rmsd_mi
475  cma=rmsd_ma
476  dc=cma-cmi
477  col=((10+244*(c-cmi)/dc) > 10) < 254
478  plots, x,y,psym=8,symsize=1.5,color=col
479endfor
480;
481usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
482;
483for n=0,NN-1 do begin
484  x=lon(n)
485  y=lat(n)
486  c=rmsd_olr(n)
487;
488  plots, x,y,psym=8,symsize=1.5,color=0
489endfor
490erase
491plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $
492   title='1) STD ratio - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge
493;
494NN=n_elements(lat)
495;
496usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
497;
498for n=0,NN-1 do begin
499  x=lon(n)
500  y=lat(n)
501  c=std_trop(n)
502  cmi=std_mi
503  cma=std_ma
504  dc=cma-cmi
505  col=((10+244*(c-cmi)/dc) > 10) < 254
506  plots, x,y,psym=8,symsize=1.5,color=col
507endfor
508;
509usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
510;
511for n=0,NN-1 do begin
512  x=lon(n)
513  y=lat(n)
514  c=std_trop(n)
515;
516  plots, x,y,psym=8,symsize=1.5,color=0
517endfor
518plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $
519   title='2) STD ratop - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge
520;
521NN=n_elements(lat)
522;
523usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
524;
525for n=0,NN-1 do begin
526  x=lon(n)
527  y=lat(n)
528  c=std_era(n)
529  cmi=std_mi
530  cma=std_ma
531  dc=cma-cmi
532  col=((10+244*(c-cmi)/dc) > 10) < 254
533  plots, x,y,psym=8,symsize=1.5,color=col
534endfor
535;
536usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
537;
538for n=0,NN-1 do begin
539  x=lon(n)
540  y=lat(n)
541  c=std_era(n)
542;
543  plots, x,y,psym=8,symsize=1.5,color=0
544endfor
545;
546plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $
547   title='3) STD rato - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge
548;
549NN=n_elements(lat)
550;
551usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
552;
553for n=0,NN-1 do begin
554  x=lon(n)
555  y=lat(n)
556  c=std_oaflx(n)
557  cmi=std_mi
558  cma=std_ma
559  dc=cma-cmi
560  col=((10+244*(c-cmi)/dc) > 10) < 254
561  plots, x,y,psym=8,symsize=1.5,color=col
562endfor
563;
564usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
565;
566for n=0,NN-1 do begin
567  x=lon(n)
568  y=lat(n)
569  c=std_oaflx(n)
570;
571  plots, x,y,psym=8,symsize=1.5,color=0
572endfor
573;
574plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $
575   title='4) STD rato - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1
576;
577NN=n_elements(lat)
578;
579usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
580;
581for n=0,NN-1 do begin
582  x=lon(n)
583  y=lat(n)
584  c=std_olr(n)
585  cmi=std_mi
586  cma=std_ma
587  dc=cma-cmi
588  col=((10+244*(c-cmi)/dc) > 10) < 254
589  plots, x,y,psym=8,symsize=1.5,color=col
590endfor
591;
592usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
593;
594for n=0,NN-1 do begin
595  x=lon(n)
596  y=lat(n)
597  c=std_olr(n)
598;
599  plots, x,y,psym=8,symsize=1.5,color=0
600endfor
601;
602closeps
603;
604end
Note: See TracBrowser for help on using the repository browser.