1 | |
---|
2 | go page_new -i 1 1 0 100 0 100 |
---|
3 | |
---|
4 | let count=($01) |
---|
5 | let flag_ref=($03) |
---|
6 | |
---|
7 | let norm_max=1.7 |
---|
8 | |
---|
9 | |
---|
10 | set window 1 |
---|
11 | |
---|
12 | DEFINE VIEWPORT/XLIM=0,1/YLIM=0.25,0.95 V1 |
---|
13 | DEFINE VIEWPORT/XLIM=0,1/YLIM=0,1 V2 |
---|
14 | DEFINE VIEWPORT/XLIM=0,1/YLIM=0.07,0.27 V3 |
---|
15 | |
---|
16 | SET VIEWPORT V1 |
---|
17 | go taylor_frame.jnl norm_max `norm_max/6,p=1` " " "Normalized Standard Deviation" 2 |
---|
18 | |
---|
19 | |
---|
20 | file/TYPE="TEXT,NUMERIC,NUMERIC,NUMERIC,NUMERIC,NUMERIC,NUMERIC"/variables="vec_id_site,vec_norm_eval,vec_correl_eval,vec_bias_eval,vec_norm_ref,vec_correl_ref,vec_bias_ref"/FORMAT=DELIMITERS/DELIM="\b" stat_($05)_($06)_($07)_($04)_($02).csv |
---|
21 | |
---|
22 | repeat/i=1:`count` ( \ |
---|
23 | let norm_ref=vec_norm_ref[i=`i`]; \ |
---|
24 | let norm_eval=vec_norm_eval[i=`i`]; \ |
---|
25 | let correl_ref=vec_correl_ref[i=`i`]; \ |
---|
26 | let correl_eval=vec_correl_eval[i=`i`]; \ |
---|
27 | let num_site=`i` ; \ |
---|
28 | go poly_vectors.jnl `norm_ref`*cos(acos(`correl_ref`)) `norm_ref`*sin(acos(`correl_ref`)) `norm_eval`*cos(acos(`correl_eval`))-`norm_ref`*cos(acos(`correl_ref`)) `norm_eval`*sin(acos(`correl_eval`))-`norm_ref`*sin(acos(`correl_ref`)) 0.2760 "arrow" ; \ |
---|
29 | polygon/over/nolabel/noaxis/nokey/pal="white"/line plt_x_arrow, plt_y_arrow, 1 ; \ |
---|
30 | go taylor_polymark `norm_ref` `correl_ref` /line/fill/nokey/pal="white"/over 1 circle 0.5 ; \ |
---|
31 | let mon_lab_x = `norm_ref`*cos(acos(`correl_ref`))+`norm_max/50` ; \ |
---|
32 | let mon_lab_y = `norm_ref`*sin(acos(`correl_ref`))+`norm_max/50` ; \ |
---|
33 | let mon_lab_x_corr =if(mon_lab_x gt norm_max*1.3) then (-1) else mon_lab_x ; \ |
---|
34 | label `mon_lab_x_corr`,`mon_lab_y`,0,0,0.1 "`num_site`" ; ) |
---|
35 | |
---|
36 | set VIEWPORT V2 |
---|
37 | repeat/i=1:`count` ( \ |
---|
38 | let id_site=vec_id_site[i=`i`] ; \ |
---|
39 | let num_site=`i` ; \ |
---|
40 | let xpos=10*(0.7+INT((`num_site`-1)/10)*1.5) ; let ypos=90-mod(`num_site`-1,10)*3 ; \ |
---|
41 | go text_put xpos ypos "`num_site` `id_site`" -1 0.2 ) |
---|
42 | |
---|
43 | go text_put 50 95 "($07) sites: ($02)-to-($02) ($04) variability" 0 0.5 |
---|
44 | |
---|
45 | set VIEWPORT V1 |
---|
46 | label 0, `norm_max/3*-1,p=1`,1,0,0.2 "@ASVersion ($06)" |
---|
47 | repeat/range=1:`flag_ref`/name=ref ( \ |
---|
48 | go world_arrow_put 0 `norm_max/3*-1,p=1` `norm_max/3,p=1` `norm_max/3*-1,p=1`;\ |
---|
49 | label `norm_max/3,p=1`, `norm_max/3*-1,p=1`,-1,0,0.2 "@AS($05)") |
---|
50 | |
---|
51 | |
---|
52 | set VIEWPORT V3 |
---|
53 | go margins_set 20 20 20 5 |
---|
54 | |
---|
55 | define axis/x=1:`count+1`:1/edges my_x |
---|
56 | |
---|
57 | let var_pft=if(j[j=1:2] le 1) then vec_bias_eval[gx=my_x@asn] else vec_bias_ref[gx=my_x@asn] |
---|
58 | shade/set_up/pal=inverse_greyscale/VGRAT=(line,col=black,THICK=3)/hlimits=1:`count+1`:1/AXES=1,1,0,0/nolabels/nokey var_pft |
---|
59 | |
---|
60 | |
---|
61 | PPL AXLABP 0,0 !turn off numeric labels |
---|
62 | PPL AXLSZE 0.95,0.25 |
---|
63 | PPL LABSET 0.3 0.3 0.3 0.3 |
---|
64 | ppl shakey 1, 0, 0.15, 5, 3, 6, ($PPL$XORG), `($PPL$XORG)+($PPL$XLEN)`,`($PPL$YORG)-1.3`,`($PPL$YORG)-0.9` |
---|
65 | ppl shade |
---|
66 | |
---|
67 | |
---|
68 | repeat/i=1:2 (\ |
---|
69 | let xlegend="@AS$`4+i`" ;\ |
---|
70 | let my_y=`i-0.5`*($PPL$YLEN)/2; \ |
---|
71 | let my_x=-0.5 ; \ |
---|
72 | label/nouser `my_x`,`my_y`,1,0,.15 `xlegend` ) |
---|
73 | |
---|
74 | repeat/i=1:`count` (\ |
---|
75 | let xlegend="@AS`i`" ;\ |
---|
76 | let my_x=`i-0.5`*($PPL$XLEN)/`count`; \ |
---|
77 | let my_y=-0.5 ; \ |
---|
78 | label/nouser `my_x`,`my_y`,0,0,.15 `xlegend` ) |
---|
79 | |
---|