Administrace komentářů

Bry vecer, snazim se pri vykreslovani grafu v gnuplotu (4.2) dostat body nad (do popredi) krivky tak aby ji prekryvaly. Jde o to ze mam treba body tvorene "prazdnyma" ctvereckama, a ty jsou protnuty errorbarem a krivkou (pokud na ni lezi). Viz obrazek. Nevypada to moc dobre, kdyz jsou tam ty ctvercky, kundicky, a pentagony preskrtnute errorbarem. Bohuzel nemuzu pouzit jen "vyplnene" symboly pro body...je tam moc dat, ktere je treba rozlisit. Nevite, prosim vas, jak je dostat nad tu krivku? Diky, Milan.

Obrazek je generovan nasledujicimi dvema skripty. Upraveni vzheldu obrazku (viewgama.gp):

set terminal wxt enhanced
set timestamp "generated at %H:%M:%S on date %d-%m-%Y by `whoami`"
###########################################################
# set common properties for all output terminals          #
###########################################################
set style line 1 lt 1 lw 2
set style line 2 lt -1 lw 2
set border 15 ls 2
set bar 1.75
set size 1,1
set pointsize 2
Tabs=273.15
Tmax=385.15
Tmin=273.15
Tincr=20
fxaxis(x)=1/x

# xtics
t1=Tmin
t2=t1+Tincr
t3=t2+Tincr
t4=t3+Tincr
t5=t4+Tincr
t6=t5+Tincr
t7=t6+Tincr

# minor xtics
tm1=Tmin + Tincr/2
tm2=t1+Tincr/2
tm3=t2+Tincr/2
tm4=t3+Tincr/2
tm5=t4+Tincr/2
tm6=t5+Tincr/2

###########################################################
# set primary graph  (x1axis in 1/T)                      #
###########################################################
set grid
set xrange [fxaxis(Tmax):fxaxis(Tmin)]
set xtics 0.0026,0.0002,0.00365 nomirror
set mxtics 2
set xlabel "{/=22 1000.T ^{ -1} / K ^{ -1}}"
set format x "%.1s"
set mytics 2
set ylabel "{/=22 ln {/Symbol=22 \g@^{/=22 \245}_{1}}}"
set key outside noreverse
call 'params.dat'
############################################################
# set secondary graph  (x2axis in C^{o})                   #
############################################################
set x2tics ('0' fxaxis(t1),"20" fxaxis(t2),"40" \
fxaxis(t3),"60" fxaxis(t4),"80" fxaxis(t5),"100"\
fxaxis(t6))
set x2label "{/=22 t / ^{o}C}"
replot
############################################################
# write modified (pretty  output to *.eps file          #
############################################################
set notimestamp
set terminal postscript enhanced monochrome "Times-Romain" 22
set output 'gama.eps'
replot
set output 'gama-nokey.eps'
set key default
set key off
replot
pause -1

A pro samotne vykresleni params.dat(tento skrip je volan predchozim skriptem viewgama.gp):

plot "work.tmp" index 0:0 u ($$5/$$1):2:3 notitle w yerrorbars 1,\
 "work.tmp" index 0:0 u ($$5/$$1):2 title "92HWA/OLS" w p ps 1.75 pt 4

replot "work.tmp" index 1:1 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 1:1 u ($$5/$$1):2 title "99BEN/DOH" w p ps 1.75 pt 6

replot "work.tmp" index 2:2 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 2:2 u ($$5/$$1):2 title "16DAV" w p ps 1.75 pt 8

replot "work.tmp" index 3:3 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 3:3 u ($$5/$$1):2 title "31GRO/SAY" w p ps 1.75 pt 10

replot "work.tmp" index 4:4 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 4:4 u ($$5/$$1):2 title "41VER/WER" w p ps 1.75 pt 12

replot "work.tmp" index 5:5 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 5:5 u ($$5/$$1):2 title "80BAN/YAL" w p ps 1.75 pt 14

replot "work.tmp" index 6:6 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 6:6 u ($$5/$$1):2 title "03BER/DOH" w p ps 1.75 pt 5

replot "work.tmp" index 7:7 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 7:7 u ($$5/$$1):2 title "04BER/DOH" w p ps 1.75 pt 7

replot (36.391525069038+(-28.248913125533)/(1/(x*298.15))+(-29.889623983045)*log(1/(x*298.15))+(0)*(1/(x*298.15)))\
title "FIT with 3 params" with lines ls 1
set pointsize 2
set xrange [0.0026:0.0036]