1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69
|
{ # -*- ifm -*-
# template for the R portion of an Rk plot
# {$D->group} returns the ifeffit group name
# {$D->parameter} returns the value of that parameter
# {$P->parameter} returns the value of a Plot parameter
}
{
$this = $PT || $D;
$ylabel = $P->plot_rylabel;
$mag = $P->tempfile;
$this->points(file => $mag,
space => 'R',
suffix => 'chir_mag',
);
$magfit = $P->tempfile;
$this->points(file => $magfit,
part => 'fit',
space => 'R',
suffix => 'chir_mag',
);
$win = $P->tempfile;
$this->points(file => $win,
space => 'R',
scale => $C->get("window_size"),
suffix => 'rwin',
);
$re = $P->tempfile;
$this->points(file => $re,
space => 'R',
suffix => 'chir_re',
yoffset => $D->rmr_offset
);
$refit = $P->tempfile;
$this->points(file => $refit,
part => 'fit',
space => 'R',
suffix => 'chir_re',
yoffset => $D->rmr_offset
);
$thiskey = $D->name;
$thiskey =~ s{_}{\\\\_}g;
if (ref($PT) =~ m{Path}) { $thiskey = $PT->name || 'Path '.$PT->Index };
$legend = ($P->showlegend) ? '' : "\nunset key" ;
q{}
}
## gnuplot script
set multiplot layout 2,1
set size 1,0.55
set origin 0,0.01
{$legend}
unset title
set xlabel "\{/*1.2 Radial distance&\{aa\}(Å)\}"
set ylabel "\{/*1.2 {$ylabel}\}"
set xrange [ {$P->rmin} : {$P->rmax} ]
set key bottom right width 1 height 1 box
# {$this->name}: ifeffit group = {$D->group}
plot '{$mag}' using 1:2 with {$C->default("gnuplot", "datastyle")} ls {$P->increm+1} title "{$thiskey}",\
'{$magfit}' using 1:2 with {$C->default("gnuplot", "datastyle")} ls {$P->increm+2} title "fit",\
'{$win}' using 1:2 with {$C->default("gnuplot", "datastyle")} ls {$P->increm+3} title "window",\
'{$re}' using 1:2 with {$C->default("gnuplot", "datastyle")} ls {$P->increm+1} title "",\
'{$refit}' using 1:2 with {$C->default("gnuplot", "datastyle")} ls {$P->increm+2} title ""
|