Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~palfa/plots/gal/gal.sm Дата изменения: Sun Jul 16 05:52:39 2006 Дата индексирования: Mon Apr 11 00:59:43 2016 Кодировка: Поисковые слова: solar eclipse

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