go page_new -i 1 1 0 100 0 100

let count=($01)
let flag_ref=($03)

let norm_max=1.7


set window 1

go taylor_frame.jnl norm_max `norm_max/6,p=1` " " "Normalized Standard Deviation" 2


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
repeat/i=1:`count` ( \ 
let norm_ref=vec_norm_ref[i=`i`]; \
let norm_eval=vec_norm_eval[i=`i`]; \
let correl_ref=vec_correl_ref[i=`i`]; \
let correl_eval=vec_correl_eval[i=`i`]; \
let id_site=vec_id_site[i=`i`] ; \
let num_site=`i` ; \
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" ; \
polygon/over/nolabel/noaxis/nokey/pal="white"/line plt_x_arrow, plt_y_arrow, 1 ; \
go taylor_polymark `norm_ref` `correl_ref` /line/fill/nokey/pal="white"/over 1 circle 0.5 ; \
let mon_lab_x = `norm_ref`*cos(acos(`correl_ref`))+`norm_max/50` ; \
let mon_lab_y = `norm_ref`*sin(acos(`correl_ref`))+`norm_max/50` ; \
label `mon_lab_x`,`mon_lab_y`,0,0,0.1 "`num_site`"  ; \
let xpos=10*(0.7+INT((`num_site`-1)/10)*1.5) ; let ypos=85-mod(`num_site`-1,10)*3 ; \
let xpos_corr=if(xpos gt 100) then (-1) else xpos ; \
go text_put xpos_corr ypos "`num_site` `id_site`" -1 0.2 )

go text_put 50 95 "($02)-to-($02) ($04) variability" 0 0.6

label 0, `norm_max/3*-1,p=1`,1,0,0.2 "@ASVersion ($06)"
label 0, `norm_max/2*-1,p=1`,1,0,0.2 "For ($07) sites"
repeat/range=1:`flag_ref`/name=ref ( \ 
go world_arrow_put 0 `norm_max/3*-1,p=1` `norm_max/3,p=1` `norm_max/3*-1,p=1`;\
label `norm_max/3,p=1`, `norm_max/3*-1,p=1`,-1,0,0.2 "@AS($05)")