! device x11 -g 800x800-100+10 -bg RoyalBlue -fg white
 device postencap Pkevolution.eps
macro prf 5 { lines $1 $2
              read k $4
		        read pk $5
	           ltype 0
	           ltype $3
              connect (LG(k)) (LG(pk)) if (pk > 0)
              ltype 0
	      }
location 5000 32000 3000 32000
lweight 2
expand 1.1
	      limits -1 2. -4.60 3.
	      ticksize -1 10 -1 10
	      box
 	      xlabel k(h/Mpc)
		   ylabel P(k) (Mpc/h)^{3}
data Sp128x256.init.DAT
prf 61 282 0 2 4
data Sp128x256.052.DAT
prf 61 282 0 2 4
data Sp128x256.092.DAT
prf 61 282 0 2 4
data Sp128x256.192.DAT
prf 61 282 0 2 4
data Sp128x256.a100.DAT
prf 61 282 0 2 4

data model.dat
prf 17 126 2 1 2
prf 129 238 2 1 2
data model.052.dat
prf 129 238 2 1 2
data model.092.dat
prf 129 238 2 1 2
data model.192.dat
prf 129 238 2 1 2

expand 0.9
relocate 0.3 2.8 label \Lambda CDM model: h=0.7 \Omega_0=0.3
relocate 0.3 2.6 label PM run with 128^3 particles
relocate 0.7 2.4 label 256^3 mesh
relocate 0.7 2.2 label 20h^{-1}Mpc box
		ltype 1
relocate 1.303 -3  draw 1.303 2.2
		ltype 0
		relocate 1.34 1.85 label Nyquist
		relocate 1.34 1.6 label Frequency
relocate -0.6 -0.3 label z=30
relocate  -0.50 0.87 label 18
relocate  -0.50 1.34 label 9.9
relocate  -0.50 2.02 label 4.2
relocate -0.6 2.6		label z=0
cursor

location 8300 16000 4500 15400
expand 0.801
	      limits 0.9 1.302 -4.0 -2.5
	      ticksize -1 2 -1 10
	      box
relocate 0.94 -3.52 label step=0
relocate 0.91 -2.85 label step=10
relocate 1.14 -2.59 label step=20		
data		Sp128x256.init.DAT
prf 61 282 0 2 4
data Sp128x256.052.DAT
prf 61 282 0 2 4
data Sp128x256.072.DAT
prf 61 282 0 2 4

data model.dat
prf 129 238 2 1 2
data model.052.dat
prf 129 238 2 1 2
data model.072.dat
prf 129 238 2 1 2
 hardcopy
 end