kill all

E0=-1272.800


read xny "e"
with g1.gr1.s1
y=(y-E0)*1000
legend "bcc"

title on
  title "Fe in GGA"
  title position 1.17716 , -0.0648753
  title restriction none

yaxislabel "Energy [mHa]"
xaxislabel "a [a$_B$.]"
autoscale
x auto tic off
x major tic spacing 0.4
vx0=0.152519
vy0=0.166
vw=0.29084
vh=0.66584
view vx0,vy0,vw,vh


read xny "s" into g2
with g2.gr1.s1
legend "bcc"
title off
yaxislabel "moment [$~m$_B$.]"
xaxislabel "a [a$_B$.]"
autoscale
x auto tic off
x major tic spacing 0.4
vx0=0.56626
view vx0,vy0,vw,vh
