killall

col=0xffffff

pw=950
ph=500
paper size pw,ph

h=0.7

w=h*ph/pw

vx0=0.1
vy0=0.12
view vx0,vy0,w,h
title "penetration depth 1 block"

g2 on
view vx0+0.5,vy0,w,h
title "penetration depth 10 blocks"

######################333
with g1
legend off

read akbl "slabifyres/+fs_sf_1block.spin1"

with gr1.s1
interpolation depth 1

min cutoff 1
null cutoff 0.12
max cutoff 1
min color 0x00
null color 0xff
max color 0xff5000
zpower 1

g1 irregular tic off
g1.gr1 irregular tic off
x tic side normal
x major tic line color col 
x minor tic line color col 
x tic side both
x tic label side normal
y major tic line color col 
y minor tic line color col 
y tic side both

xaxislabel on
xaxislabel "kx"
yaxislabel "ky"

#########################

with g2

legend off

read akbl "slabifyres/+fs_sf_10block.spin1"



with gr1.s1
interpolation depth 1

min cutoff 1
null cutoff 7.8
max cutoff 0.6
min color 0x00
null color 0xff
max color 0xff5000
zpower 1
g2 irregular tic off
g2.gr1 irregular tic off
x tic side normal
x major tic line color col 
x minor tic line color col 
x tic side both
x tic label side normal
y major tic line color col 
y minor tic line color col 
y tic side both

xaxislabel on
xaxislabel "kx"
yaxislabel "ky"


print to "fssemi.png" dpi 75