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 | go taylor_frame.jnl norm_max `norm_max/6,p=1` " " "Normalized Standard Deviation" 2 |
---|
13 | |
---|
14 | |
---|
15 | file/TYPE="TEXT,NUMERIC,NUMERIC,NUMERIC,NUMERIC"/variables="vec_id_site,vec_norm_eval,vec_correl_eval,vec_norm_ref,vec_correl_ref"/FORMAT=DELIMITERS/DELIM="\b" stat_($05)_($06)_($07)_($04)_($02).csv |
---|
16 | repeat/i=1:`count` ( \ |
---|
17 | let norm_ref=vec_norm_ref[i=`i`]; \ |
---|
18 | let norm_eval=vec_norm_eval[i=`i`]; \ |
---|
19 | let correl_ref=vec_correl_ref[i=`i`]; \ |
---|
20 | let correl_eval=vec_correl_eval[i=`i`]; \ |
---|
21 | let id_site=vec_id_site[i=`i`] ; \ |
---|
22 | let num_site=`i` ; \ |
---|
23 | 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" ; \ |
---|
24 | polygon/over/nolabel/noaxis/nokey/pal="white"/line plt_x_arrow, plt_y_arrow, 1 ; \ |
---|
25 | go taylor_polymark `norm_ref` `correl_ref` /line/fill/nokey/pal="white"/over 1 circle 0.5 ; \ |
---|
26 | let mon_lab_x = `norm_ref`*cos(acos(`correl_ref`))+`norm_max/50` ; \ |
---|
27 | let mon_lab_y = `norm_ref`*sin(acos(`correl_ref`))+`norm_max/50` ; \ |
---|
28 | label `mon_lab_x`,`mon_lab_y`,0,0,0.1 "`num_site`" ; \ |
---|
29 | let xpos=10*(0.7+INT((`num_site`-1)/10)*1.5) ; let ypos=85-mod(`num_site`-1,10)*3 ; \ |
---|
30 | let xpos_corr=if(xpos gt 100) then (-1) else xpos ; \ |
---|
31 | go text_put xpos_corr ypos "`num_site` `id_site`" -1 0.2 ) |
---|
32 | |
---|
33 | go text_put 50 95 "($02)-to-($02) ($04) variability" 0 0.6 |
---|
34 | |
---|
35 | label 0, `norm_max/3*-1,p=1`,1,0,0.2 "@ASVersion ($06)" |
---|
36 | label 0, `norm_max/2*-1,p=1`,1,0,0.2 "For ($07) sites" |
---|
37 | repeat/range=1:`flag_ref`/name=ref ( \ |
---|
38 | go world_arrow_put 0 `norm_max/3*-1,p=1` `norm_max/3,p=1` `norm_max/3*-1,p=1`;\ |
---|
39 | label `norm_max/3,p=1`, `norm_max/3*-1,p=1`,-1,0,0.2 "@AS($05)") |
---|