galps
 device psfile gal              
 gal
 hardcopy
 device psfile galsim
 galsim
 hardcopy
 device psfile galpalfa
 galpalfa
 hardcopy
 device psfile galsimpalfa
 galsimpalfa
 hardcopy
 device psfile gallines 
 gallines
 hardcopy
 device psfile galsimlines 
 galsimlines     
 hardcopy
 device psfile galpalfalines
 galpalfalines   
 hardcopy
 device psfile galsimpalfalines 
 galsimpalfalines 
 hardcopy
 device x11
 !convert gal.ps gal.png
 !convert galsim.ps galsim.png
 !convert galpalfa.ps galpalfa.png
 !convert galsimpalfa.ps galsimpalfa.png

gal
 galinit
 plotpsr 0 0 0 0
 galfinish

gallines    
  gal 
  galaolimits
galsimlines 
  galsim 
  galaolimits
galsimpalfalines
  galsimpalfa
  galaolimits
galpalfalines
  galpalfa
  galaolimits

galsim
 galinit
 plotpsr 0 0 0 0
 plotsim
 galfinish

galpalfa
 galinit
 plotpsr 0 0 0 0
 plotpalfa
 galfinish

galsimpalfa
 galinit
 plotpsr 0 0 0 0
 plotsim
 plotpalfa
 galfinish


galinit
 define exp_def 1.5
 define exp_psr 0.8
 define exp_sim 0.8
 define exp_palfa 4.0
 define color_def black
 define color_arm black
 define color_psr red 
 define color_sim blue
 define color_palfa black
 define rad_gc 0.4
 define rad_sun_1 0.1
 define rad_sun_2 0.4

 define r0 8.5

 define TeX_strings 1
 
 erase
 ctype $color_def
 expand $exp_def
 limits -13 13 -13 13
 box
 xlabel x (kpc)
 ylabel y (kpc)

galfinish
 plotarm
 symbolgc
 symbolsun

plotarm
 data log_arms.out
 lines 2 99999
 read {n 1 x 3 y 4}
 ctype $color_arm
 do i = 1, 5 { 
   connect x y if (n==$i) 
 }

plotpsr 4
 data psrxyz.dat
 read {flag 1 x 2 y 3 ar1 8 ar2 9 ar3 10 ar4 11}
 ctype $color_psr
 ptype 20 3
 expand $exp_psr
 points x y if (flag==0 && ar1>=$1 && ar2>=$2 && ar3>=$3 && ar4>=$4)

plotsim
 data survey.det
 lines 3 9999
 read {x 9 y 10}
 ctype $color_sim
 ptype 20 3
 expand $exp_sim
 points x y

plotpalfa
 data palfa.xyz
 read {flag 1 x 2 y 3}
 ctype $color_palfa
 ptype 5 2
 do i = 0, 10 {
   expand $!(($i/10.)*$exp_palfa)
   points x y if (flag==0)
 }

symbolgc
 lweight 3
 ctype $color_def
 set theta = 0, 3.1416/2., .01
 set x = -$rad_gc * cos(theta)
 set y = $rad_gc * sin(theta)
 connect x y
 set x = -x
 set y = -y
 connect x y
 set x = x concat {0}
 set y = y concat {0}
 set x = -x
 shade 0 x y
 set x = -x
 set y = -y
 shade 0 x y
 lweight 1

symbolsun
 lweight 3
 ctype $color_def
 set theta = 0, 3.1416*2., .01
 set x = $rad_sun_2 * cos(theta)
 set y = $rad_sun_2 * sin(theta) + $r0
 connect x y
 set x = $rad_sun_1 * cos(theta)
 set y = $rad_sun_1 * sin(theta) + $r0
 shade 1 x y
 lweight 1

galaolimits
 lline 32
 lline 77
 lline 168
 lline 214

lline 1
 relocate 0 $r0
 define _tmp_lline_1 (50.*sin($1*3.14159265/180.))
 define _tmp_lline_2 ($r0 - 50.*cos($1*3.14159265/180.))
 draw $_tmp_lline_1 $_tmp_lline_2