Portál AbcLinuxu, 29. prosince 2025 00:16
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]
Nasel jsem to po tydnu hledani. V pointtype (zkaratka pt) je treba zvolit cislo z intervalu 70 az 75. Pak jsou symboly po exportu do eps (set term postscript enhanced monochrome) vyplnene, a tim padem nepreskrtnute errorbarem.
Tiskni
Sdílej:
ISSN 1214-1267, (c) 1999-2007 Stickfish s.r.o.