Dear Gnuplot users,
I'm trying to use a pm3d map to plot the isotopic distribution as a
function of atomic number (Z) and mass number (A).
I have three main doubts:
First.- I would like to assign the colours directly onto the integer
values (centred in the values: Z,A), extending to half a unit up and
down in both directions.
Currently, by using corners2colors {c1|c2|c3|c4} I can obtain the value
of the crossing points in one of the intermediate regions, but never
centred around the crossing point (which would be useful to read the
value as a function of Z and A).
Second.- I would like to expand the x (A) axis in order to obtain square
pixels (not the rectangular ones I get now). That is to set a one-to-one
scale between both axis.
Three.- To force Gnuplot to create an output file (eps or png) which
holds the whole picture! (currently, my eps is cutting the right boundary).
I also include the macro file I use.
Thank you very much in advance!
------------------------------------------------------------------------
# Execute this macro file with: load "fiename"
#set terminal svg
#set term png color transparent nocrop enhanced font arial 8 size 420,320
#set output 'residuals.png'
set term postscript landscape eps color
set output 'residuals.eps'
#set output 'test.svg'
# Limits
MIN=1E+3 # Minimum value
MAX=1E+11 # Maximum value
CPD=4 # Colors per decade
# Labels & Titles
#set title 'FISSION PRODUCT RATES' font ",15"
set ylabel "Z (atomic number)"
set xlabel "A (mass number)"
set cblabel "Isotopic yields (isotope/cm3/s/MW of beam)"
#set label "Isotopic yields (isotope/cm3/s/MW of beam)" rotate by -90 at 0,0 right
# Settings
set ytics 5
set xtics 10
set mytics 5
set mxtics 5
set grid xtics ytics # mxtics mytics
set ticslevel 0
#set style line 100 lt 5 lw 0.01
unset key
set colorbox vertical user ori 0.865,0.2 size 0.03,0.6
#set pm3d at bs
set pm3d map explicit corners2color c4 # hidden3d 100
set logscale cb
COLORS=(log10(MAX)-log10(MIN))*CPD+1
set palette defined ( 1 "white", 2 "magenta", 3 "blue", \
4 "green", 5 "yellow", 6 "red", 7 "black") maxcolors COLORS
set cbrange [MIN:MAX]
set yrange [25:65]
set xrange [60:170]
splot 'residuals.dat' us 2:1:3, 'stable.dat' us 3:2:(0) with points lt -1 pt 5 ps 1 #lines lt -1 lw 0
# pause -1 "Hit return to continue"
--
-----------------------------------------------------------------------------------------------------------------
Adonai Herrera Martinez
-----------------------------------------------------------------------------------------------------------------
AB Dept. ATB/EET Group
CERN, European Organization for Nuclear Research
CH-1211 Geneva
Switzerland
-----------------------------------------------------------------------------------------------------------------
tel: +41 22 767 81 65
fax: +41 22 766 95 28
adonai.herrera.martinez@...
http://www.cern.ch/adonai
|