[5570] | 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 | |
---|