[Gnuplot-info] pm3d map, splotting doubts

[Adonai H. M. <ado...@ce...>]

  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

ado...@ce...

http://www.cern.ch/adonai