set   redc = {  255 0  255    0    0  200 255  255 200 255   0   0 100  90 150 200 200 150   0 255}
set   grnc = {  255 0   0    255   0  200 255  255 255   0 255   0 100 100 150 200 220   0 150 255}
set   bluc  ={  255 0   0    0   255    0   0  150 255 255 255  55 255 100 150 200 255   0   0 255}
		ctype = redc + 256*(grnc +256*bluc)
		ctype = { white black red green blue dyellow yellow byellow bbyellow magenta cyan 
		               darkgrey grey lightgrey steel
		               lblue lightblue
		               darkred darkgreen wwhite}
!dev x11 -g 1000x1000 -bg black
!dev postfile halosxy.ps
dev postscript_bg  ":SY=echo File is \$F:OF=/europa1/aklypin/PROJECTS/PM/TEST32x128/points.ps:"
erase
macro circ 4 {# circles of different radius
	expand 0.2
        points $1 $2 if( (abs((LG($3))-10.) < 0.75 ))
	expand 0.5           
        points $1 $2 if( (abs((LG($3))-11.0) < 0.25 ))
	expand 0.9           
        points $1 $2 if( (abs((LG($3))-11.5) < 0.25 ))
	expand 1.3           
        points $1 $2 if( (abs((LG($3))-12.0) < 0.25 ))
	expand 1.9           
        points $1 $2 if( (abs((LG($3))-12.5) < 0.25 ))
	expand 2.8           
        points $1 $2 if( (abs((LG($3))-14.0) < 1.25 ))
             }
macro pnts 3 {# points of different color
		ctype grey        points $1 $2 if( (v < 1.0) )
		ctype blue        points $1 $2 if( (v> 1.0)   &&  (v< 1.250) )
		ctype darkgreen   points $1 $2 if( (v> 1.250) &&  (v< 1.50) )
		ctype darkred     points $1 $2 if( (v> 1.50)  &&  (v< 1.750) )
		ctype red         points $1 $2 if( (v> 1.7500)&&  (v< 2.0) )
		ctype magenta     points $1 $2 if( (v> 2.00)  &&  (v< 2.250) )
 		ctype green       points $1 $2 if( (v> 2.250) &&  (v< 2.50) )
		ctype yellow      points $1 $2 if( (v> 2.50)  &&  (v< 2.6) )
 		ctype byellow     points $1 $2 if( (v> 2.6)   &&  (v< 2.70) )
 		ctype bbyellow    points $1 $2 if( (v> 2.7)   &&  (v< 2.8) )
 		ctype lightblue   points $1 $2 if( (v> 2.8)   &&  (v< 3.0) )
 		ctype wwhite      points $1 $2 if( (v> 3.0) )
               }
location 2000	29000 2000 29000
data	coordinates.ascii.a100
lines          6 50000
read		{x 2 y 3  z 4 vx 5 vy 6 vz 7}
set v = (LG(sqrt(vx**2+vy**2+vz**2)))

expand 1.0001
limits    0. 20. 0. 20.
ctype		red
box
expand	1.001
		relocate 16. -0.6 label Y(Mpc/h)
		expand 0.8 ctype cyan
		relocate 0.1 20.3
		label LCDM h=0.7 \Omega_0=0.3 Z=0 N=32^3x128^3
		relocate 9.1 20.3
		ctype lightblue label Velocity(km/s):
		ctype red    relocate 13. 20.3 label 100
		ctype magenta relocate 14. 20.3 label 150
		ctype green  relocate 15.0 20.3 label 200
		ctype yellow relocate 16.0 20.3 label 350
		ctype byellow relocate 17. 20.3 label 450
		ctype bbyellow relocate 18. 20.3 label 550
		ctype wwhite  relocate 19.0 20.3 label 1000

ptype	0	0
pnts	 y z v
expand 1.000
 hardcopy
 end