ms
	dev x11 -fg black -bg white

exit
	quit

bh0easy
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define y_gutter 0 
 window 1 2 1 1 limits -10 10 -5 5 #limits -4 4 -2 2     
 #---axis
        SET dimen(vlab)=6.s\n SET vlab[1]='1'\n SET vlab[2]='2'\n SET vlab[3]='6'\n SET vlab[4]='9'\n SET vlab[5]='10'
	SET vlab[0]='0'
        SET vbig = {0 1 2 6 9 10}
        SET vsma = {0 1 2 3 4 5 6 7 8 9 10 100 1E3 1E4 1E5}
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 #---BH
 	SET x=-2,2,0.01 
 	SET y=sqrt(4-x**2) ctype ltgrey SHADE 0 x y
 	SET y=-sqrt(4-x**2) ctype ltgrey SHADE 0 x y
 ctype black
 relocate 0 1.5\n putlabel 5 J=0
 relocate 0 3.5\n putlabel 5 black\n relocate 0 3.0\n putlabel 5 hole
 #---disk ISCOs (Bardeen et al. 1972 ApJ 178, 347, Eq. 2.21: 
 set y=-1,1,0.01 set x=6+20*y**2. ctype ltgrey SHADE 0 x y set x=-x SHADE 0 x y
 angle 90
 ctype ltgrey relocate 6.00 -2.0\n putlabel 5 "J=0"
 angle 0 ctype grey
 relocate 6.0 3.5 putlabel 5 accretion disk
 relocate 6.0 3.0 putlabel 5 with
 relocate 6.0 2.5 putlabel 5 innermost stable
 relocate 6.0 2.0 putlabel 5 circular orbits
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 limits -100 100 -100 100 #AU
 SET dimen(vlab)=9.s\n SET vlab[0]='Eris'\n SET vlab[1]='P'\n SET vlab[2]='N'\n SET vlab[3]='U'\n SET vlab[4]='S'\n SET vlab[5]='J'\n SET vlab[6]='0'\n SET vlab[7]='10'\n SET vlab[8]='100'
 SET vbig = {-67.67 -40 -30 -20 -10 -5 0 10 100}
 SET vsma = {0 10 20 30 40 50 60 70 80 90 100 }
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy2 $($gx2-$gx1) 1 1	# top
 ctype black xlabel \raise+16000 radius in units of AU for M_{BH}=10^9 M_{Sun}
 ctype black

bh0
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 #define x_gutter 0 
 define y_gutter 0 
 window 1 2 1 1 limits -10 10 -5 5 #limits -4 4 -2 2     
 #window 1 1 1 1 limits -10 10 -10 10
 #---axis
 #ticksize 10 10 10 10
 #box 0 4 4 4
        SET dimen(vlab)=6.s\n SET vlab[1]='1'\n SET vlab[2]='2'\n SET vlab[3]='6'\n SET vlab[4]='9'\n SET vlab[5]='10'
	SET vlab[0]='0'
        SET vbig = {0 1 2 6 9 10}
        SET vsma = {0 1 2 3 4 5 6 7 8 9 10 100 1E3 1E4 1E5}
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 #---BH
 #SET a=0.835 SET theta=0,2*pi,0.01*pi # ---max ergosphere
 #SET arr=(1+sqrt(1-((a)**2.)))
 #SET x=arr*(1+a**2)*sin(theta) SET y=arr*cos(theta) ctype grey connect x y
 	SET x=-2,2,0.01 
 	SET y=sqrt(4-x**2) ctype ltgrey SHADE 0 x y
 	SET y=-sqrt(4-x**2) ctype ltgrey SHADE 0 x y
 SET a=1.000 SET theta=0,2*pi,0.01*pi 
 SET arr=(1+sqrt(1-((a)**2.)))
 SET x=arr*(1+a**2)*sin(theta) SET y=arr*cos(theta) ctype grey SHADE 0 x y
 ##echo $(arr*(1+a**2))
 SET x=-arr,arr,0.005*arr
 SET y=sqrt(arr**2.-x**2) ctype black SHADE 0 x y
 SET y=-sqrt(arr**2.-x**2) ctype black SHADE 0 x y
 ##ctype black SET x=-2,2,0.01 SET y=sqrt(4-x**2) connect x y
 relocate 0 0\n ctype white\n putlabel 5 J=1\n ctype black
 relocate 0 1.5\n putlabel 5 J=0
 #relocate 2 1.5\n ctype grey\n putlabel 6 J=0.835\n ctype black
 relocate 0 3.5\n putlabel 5 black\n relocate 0 3.0\n putlabel 5 hole
 #---disk ISCOs (Bardeen et al. 1972 ApJ 178, 347, Eq. 2.21: 
 set y=-1,1,0.01 set x=1+10*y**2. ctype dkgrey SHADE 0 x y set x=-x SHADE 0 x y
 set y=-1,1,0.01 set x=6+20*y**2. ctype ltgrey SHADE 0 x y set x=-x SHADE 0 x y
 set y=-1,1,0.01 set x=9+40*y**2. ctype grey SHADE 0 x y set x=-x SHADE 0 x y
 angle 90
 ctype dkgrey relocate 0.5 -5\n label "  J=1"
 ctype dkgrey relocate 1.0 -5\n label "  co-"
 ctype dkgrey relocate 1.5 -5\n label "  rotating"
 ctype ltgrey relocate 6.00 -2.0\n putlabel 5 "J=0"
 ctype grey relocate 8.75 -5\n label "  J=1"
 ctype grey relocate 9.25 -5\n label "  counter-"
 ctype grey relocate 9.75 -5\n label "  rotating"
 angle 0 ctype grey
 relocate 6.0 3.5 putlabel 5 accretion disk
 relocate 6.0 3.0 putlabel 5 with
 relocate 6.0 2.5 putlabel 5 innermost stable
 relocate 6.0 2.0 putlabel 5 circular orbits
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 limits -100 100 -100 100 #AU
 SET dimen(vlab)=9.s\n SET vlab[0]='Eris'\n SET vlab[1]='P'\n SET vlab[2]='N'\n SET vlab[3]='U'\n SET vlab[4]='S'\n SET vlab[5]='J'\n SET vlab[6]='0'\n SET vlab[7]='10'\n SET vlab[8]='100'
 SET vbig = {-67.67 -40 -30 -20 -10 -5 0 10 100}
 SET vsma = {0 10 20 30 40 50 60 70 80 90 100 }
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy2 $($gx2-$gx1) 1 1	# top
 ctype black xlabel \raise+16000 radius in units of AU for M_{BH}=10^9 M_{Sun}
 #limits -840 840 -840 840 #light-minutes
 #AXIS $fx1 $fx2 0 0 $gx1 $gy2 $($gx2-$gx1) 1 1	# top
 #ctype black xlabel \raise+16000 radius in units of light-minutes for M_{BH}=10^9 M_{Sun}
 ctype black

bh1
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 10 0 10
        SET dimen(vlab)=6.s\n SET vlab[1]='1'\n SET vlab[2]='2'\n SET vlab[3]='6'\n SET vlab[4]='9'\n SET vlab[5]='10'
	SET vlab[0]='0'
        SET vbig = {0 1 2 6 9 10}
        SET vsma = {0 1 2 3 4 5 6 7 8 9 10 100 1E3 1E4 1E5}
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 #---BH
 	SET x=-2,2,0.01 
 	SET y=sqrt(4-x**2) ctype ltgrey SHADE 0 x y
 SET a=1.000 SET theta=0,2*pi,0.01*pi 
 SET arr=(1+sqrt(1-((a)**2.)))
 SET x=arr*(1+a**2)*sin(theta) SET y=arr*cos(theta) ctype grey SHADE 0 x y
 SET x=-arr,arr,0.005*arr
 SET y=sqrt(arr**2.-x**2) ctype black SHADE 0 x y
 relocate 0.5 0.375\n ctype white\n putlabel 5 "J=1"\n ctype black
 relocate 0 1.65\n putlabel 6 " BH"
 relocate 0 1.25\n putlabel 6 " J=0"
 #---disk ISCOs (Bardeen et al. 1972 ApJ 178, 347, Eq. 2.21: 
 set y=-1,1,0.01 set x=1+10*y**2. ctype dkgrey SHADE 0 x y set x=-x SHADE 0 x y
 set y=-1,1,0.01 set x=6+20*y**2. ctype ltgrey SHADE 0 x y set x=-x SHADE 0 x y
 set y=-1,1,0.01 set x=9+40*y**2. ctype grey SHADE 0 x y set x=-x SHADE 0 x y
 angle 90
 angle 0 ctype grey
 relocate 5.5 2.5 putlabel 5 X-ray emission region
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 angle 90 AXIS $fy1 $fy2 0 0 $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 limits 0 10 0 10 #---orbital
 SET vborb = 0 concat 1 concat 6 concat 10 SET vsorb=vborb
 SET dimen(vlorb)=4.s\n SET vlorb[1] = '0.36'\n SET vlorb[2] = '5.25'\n SET vlorb[3] = '11.5'
 AXIS $fx1 $fx2 vsorb vborb vlorb $gx1 $($gy1/2.) $($gx2-$gx1) 1 0 # basement
 ctype black xlabel orbital time in units of days for M_{BH}=10^9 M_{Sun}
 limits 0 100 0 100
 angle 90 AXIS $fy1 $fy2 0 0 $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of AU for M_{BH}=10^9 M_{Sun}

bhlog1
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 1 0 1
 #---BH
 SET lr=lg(2) SET theta=0,1.99*pi,0.01*pi 
 set rx=lr*cos(theta) set ry=lr*sin(theta) ctype ltgrey SHADE 0 rx ry 
 SET a=1.000 SET theta=0.01,0.50*pi,0.01*pi 
 SET arr=(1+sqrt(1-((a)**2.)))
 SET ex=arr*(1+a**2)*sin(theta) SET ey=arr*cos(theta) ctype grey 
 set theta=atan(ey/ex) set arr=lg(sqrt(ex**2.+ey**2.))
 set ex=arr*cos(theta) set ey=arr*sin(theta) 
 SHADE histogram 0 ex ey
 #set lx=lg(x) set ly=lg(y) SHADE 0 lx ly
 ctype black relocate $(lg(1.35)) $(lg(1.5)) putlabel 5 BH
 relocate $(lg(1.35)) $(lg(1.35)) putlabel 5 (J=0)
 #---disk ISCOs (Bardeen et al. 1972 ApJ 178, 347, Eq. 2.21: 
 set y=0,1,0.001 set x=100*y**2. ctype dkgrey SHADE histogram 0 x y 
 set y=0,1,0.001 set x=lg(6)+75*y**2. ctype ltgrey SHADE histogram 0 x y 
 set y=0,1,0.001 set x=lg(9)+50*y**2. ctype grey SHADE histogram 0 x y
 angle 0 ctype grey
 relocate 0.78 0.25 putlabel 5 X-ray
 relocate 0.78 0.2 putlabel 5 emission
 relocate 0.78 0.15 putlabel 5 region
 SET dimen(vlab)=5.s\n SET vlab[1]='2'\n SET vlab[2]='6'\n SET vlab[3]='9'\n SET vlab[4]='10'
 SET vlab[0]='1'
 SET vbig = 0 concat $(lg(2)) concat $(lg(6)) concat $(lg(9)) concat $(lg(10)) 
 SET vsma = 0 concat $(lg(2)) concat $(lg(3)) concat $(lg(4)) concat $(lg(5)) concat $(lg(6)) concat $(lg(7)) concat $(lg(8)) concat $(lg(9)) concat $(lg(10))
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0   # bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vbig = 0 concat $(lg(10)) 
 SET dimen(vlab)=2.s\n SET vlab[0]='1'\n SET vlab[1]='10'
 angle 90 AXIS $fy1 $fy2 vsma vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 limits $(lg(10.)) $(lg(100.)) $(lg(10.)) $(lg(100.))
 SET vaub = $(lg(10)) concat $(lg(100))
 SET vaus = 0 concat $(lg(20)) concat $(lg(30)) concat $(lg(40)) concat $(lg(50)) concat $(lg(60)) concat $(lg(70)) concat $(lg(80)) concat $(lg(90)) concat $(lg(100))
 SET dimen(vaul)=2.s\n SET vaul[0] = '10'\n SET vaul[1] = '100'
 angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of AU for M_{BH}=10^9 M_{Sun}

bh2
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 100 0 100
 #---BH
 SET lr=2 SET theta=0,2*pi,0.01*pi set x=lr*cos(theta) set y=lr*sin(theta) 
 ctype ltgrey SHADE 0 x y 
 SET lr=1 SET theta=0,2*pi,0.01*pi set x=lr*cos(theta) set y=lr*sin(theta) 
 ctype black SHADE 0 x y 
 #---disk: 
 set x=0,100,0.01 set y=0.066*x ctype dkgrey shade histogram 0 x y
 set x=(6),100,0.01 set y=0.1*(x-(6)) ctype ltgrey shade histogram 0 x y
 set x=(9),100,0.01 set y=0.11*(x-(9)) ctype grey shade histogram 0 x y
 relocate $((55)) $((20.0)) putlabel 5 UV/optical
 relocate $((55)) $((16.0)) putlabel 5 continuum
 relocate $((55)) $((12.0)) putlabel 5 emission
 relocate $((55)) $((8.0)) putlabel 5 region
 ctype black angle 0
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.40,2,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x) 
 set r=(sqrt(x**2.+y**2.)) set theta=atan(y/x) 
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=px/3 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=2*px connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 14 51\n putlabel 6 (possible)
 relocate 14 46\n putlabel 6 jet launching and
 relocate 14 42\n putlabel 6 magnetohydrodynamic
 relocate 14 38\n putlabel 6 collimation region
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]='10'\n SET vlab[2]='100'\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 10 100 1E3 1E4 1E5}
 SET vsma = 0 concat $((10)) concat $((20)) concat $((30)) concat $((40)) concat $((50)) concat $((60)) concat $((70)) concat $((80)) concat $((90)) concat $((100)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 limits 0 1 0 1 #---orbital
 SET vborb = 0 concat 0.1 concat 0.5 concat 1 SET vsorb=vborb
 SET dimen(vlorb)=4.s\n SET vlorb[1] = '1.64'\n SET vlorb[2] = '18.1'\n SET vlorb[3] = '52'
 AXIS $fx1 $fx2 vsorb vborb vlorb $gx1 $($gy1/2.) $($gx2-$gx1) 1 0 # basement
 ctype black xlabel orbital time in units of weeks for M_{BH}=10^9 M_{Sun}
 limits $((0)) $((140.)) $((0)) $((140.))
 SET vlmb = $((0)) concat $((14.)) concat $((140.))
 SET vlms = vlmb
 SET dimen(vlml)=3.s\n SET vlml[0] = '0'\n SET vlml[1] = '14'\n SET vlml[2] = '140'
 angle 90 AXIS $fy1 $fy2 vlms vlmb vlml $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-hours for M_{BH}=10^9 M_{Sun}
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 ptype 20 3\n limits $((0.0)) $((5.8)) $((0.0)) $((5.8))
 ctype grey\n set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=(xs+3.25) \n set ys0=4.25,61.25,9.5
 set lys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 points lxs lys
 limits 0 100 0 100
 relocate 80 80 label stars\n ctype black
 relocate 81 10 angle 45 label shielding\n relocate 87 10 label gas\n angle 90
 relocate 100 12 label He II BELR
 angle 0

bhlog2
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 2 0 2
 ##---axis
 #ticksize 10 10 10 10 box 0 4 4 4
 #---BH
 SET lr=lg(2) SET theta=0,2*pi,0.01*pi 
 set x=lr*cos(theta) set y=lr*sin(theta) ctype ltgrey SHADE 0 x y ctype black
 relocate $(lg(1.35)) $(lg(1.5)) putlabel 5 BH
 relocate $(lg(1.35)) $(lg(1.25)) putlabel 5 (J=0)
 #---disk: 
 set x=0,2,0.01 set y=0.066*x ctype dkgrey shade histogram 0 x y
 set x=lg(6),2,0.01 set y=0.1*(x-lg(6)) ctype ltgrey shade histogram 0 x y
 set x=lg(9),2,0.01 set y=0.11*(x-lg(9)) ctype grey shade histogram 0 x y
 relocate $(lg(30)) $(lg(2.83)) putlabel 5 UV/optical
 relocate $(lg(30)) $(lg(2.38)) putlabel 5 continuum
 relocate $(lg(30)) $(lg(2.00)) putlabel 5 emission
 relocate $(lg(30)) $(lg(1.68)) putlabel 5 region
 ctype dkgrey 
 relocate $(lg(4)) $(lg(2.0)) putlabel 5 X-ray 
 relocate $(lg(4)) $(lg(1.68)) putlabel 5 emission
 relocate $(lg(4)) $(lg(1.41)) putlabel 5 region
 ctype black angle 0
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.50,2,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x) 
 set r=lg(sqrt(x**2.+y**2.)) set theta=atan(y/x) 
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=px/3 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=2*px connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 0.5 1.3\n putlabel 5 jet launching and
 relocate 0.5 1.2\n putlabel 5 magnetohydrodynamic
 relocate 0.5 1.1\n putlabel 5 collimation region
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]='10'\n SET vlab[2]='100'\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 1 2 3 4 5}
 SET vsma = 0 concat $(lg(2)) concat $(lg(6)) concat $(lg(9)) concat $(lg(10)) concat $(lg(20)) concat $(lg(30)) concat $(lg(40)) concat $(lg(50)) concat $(lg(60)) concat $(lg(70)) concat $(lg(80)) concat $(lg(90)) concat $(lg(100)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 #AXIS $fx1 $fx2 vsma vbig vlab $gx1 $(0.5*$gy1) $($gx2-$gx1) 1 0 # basement
 #limits 14.47 16.47 14.47 16.47
 #angle 90 AXIS $fy1 $fy2 0 0 $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 #ctype black ylabel \raise+750 log(radius) in units of cm for M_{BH}=10^9 M_{Sun}
 #limits 9.47 11.47 9.47 11.47
 #angle 90 AXIS $fy1 $fy2 0 0 $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 #ctype black ylabel \raise+750 log(radius) in units of km for M_{BH}=10^9 M_{Sun}
 #limits $(lg(10.)) $(lg(1000.)) $(lg(10.)) $(lg(1000.))
 #SET vaub = $(lg(10)) concat $(lg(100)) concat $(lg(1000))
 #SET vaus = vaub
 #SET dimen(vaul)=3.s\n SET vaul[0] = '10'\n SET vaul[1] = '100'\n SET vaul[2] = '10^3'
 #angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 #ctype black ylabel \raise+750 radius in units of AU for M_{BH}=10^9 M_{Sun}
 limits $(lg(1.4)) $(lg(140.)) $(lg(1.4)) $(lg(140.))
 SET vlmb = $(lg(1.4)) concat $(lg(14.)) concat $(lg(140.))
 SET vlms = vlmb
 SET dimen(vlml)=3.s\n SET vlml[0] = '1.4'\n SET vlml[1] = '14'\n SET vlml[2] = '140'
 angle 90 AXIS $fy1 $fy2 vlms vlmb vlml $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-hours for M_{BH}=10^9 M_{Sun}
 limits 0 2 0 2

bh3 # blr
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 
 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 1000 0 1000
 #---BH
 SET theta=0,2*pi,0.01*pi 
 set x=cos(theta) set y=sin(theta) ctype black SHADE 0 x y
 #---disk: 
 set x=(6),1E3,0.01 set y=0.1*(x-(6)) ctype ltgrey shade histogram 0 x y
 ctype black angle 0
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.40,3,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x)
 set r=(sqrt(x**2.+y**2.)) set theta=atan(y/x)
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=px/2 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 46 850\n putlabel 6 (possible)
 relocate 46 800\n putlabel 6 jet 
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]=''\n SET vlab[2]='100'\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 0 100 1E3 1E4 1E5}
 SET vsma = 0 concat $((100)) concat $((200)) concat $((300)) concat $((400)) concat $((500)) concat $((600)) concat $((700)) concat $((800)) concat $((900)) concat $((1000)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 limits 0 31 0 31 #---orbital
 SET vborb = 0 concat 3.1 concat 15.5 concat 31 SET vsorb=vborb
 SET dimen(vlorb)=4.s\n SET vlorb[1] = '1'\n SET vlorb[2] = '11'\n SET vlorb[3] = '31'
 AXIS $fx1 $fx2 vsorb vborb vlorb $gx1 $($gy1/2.) $($gx2-$gx1) 1 0 # basement
 ctype black xlabel orbital time in units of years for M_{BH}=10^9 M_{Sun}
 limits $((0.0)) $((58.)) $((0.0)) $((58.))
 SET vldb = $((0.0)) concat $((0.0)) concat $((5.8)) concat $((58.))
 SET vlds = vldb
 SET dimen(vldl)=4.s\n SET vldl[0] = '0'\n SET vldl[1] = ''\n SET vldl[2] = '5.8'\n SET vldl[3] = '58.'
 angle 90 AXIS $fy1 $fy2 vlds vldb vldl $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-days for M_{BH}=10^9 M_{Sun}
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 ptype 20 3\n limits $((0.0)) $((58.)) $((0.0)) $((58.))
 ctype grey\n set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=(xs+3.25) \n set ys0=4.25,61.25,9.5
 set lys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 expand 0.6 points lxs lys expand 1.2
 limits 0 1000 0 1000
 relocate 760 900 label stars\n ctype black
 #---BLR
 angle 63
 ctype grey set x=50,1000,1 set y=0.48*(x-50)+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 set z=0.1*x+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 set yz=y concat reverse(z) set xx = x concat reverse(x)
 shade 899 xx yz
 ctype black angle 90
 relocate $(10**2) $(10**1.50) putlabel 6 He II
 relocate $(10**2.48) $(10**1.66) putlabel 6 C IV, Si IV
 relocate $(10**2.95) $(10**2.66) putlabel 6 H\beta, Mg II, C III]
 angle 45\n relocate 500 110 putlabel 6 Broad Line Region
 ctype black angle 0
 #---torus:
 ctype dkgrey
 set x1=550,1000,0.01 set y1=1.000*(x1-550.)+0.1*(x1-6)
 set x2=750,1000,0.01 set y2=1.000*(x2-750.)+0.1*(x2-6)
 set x12 = x1 concat reverse(x2)\n set y12 = y1 concat reverse(y2)
 shade 0 x12 y12
 ctype black
 angle 45\n relocate 675 110 ctype white putlabel 6 Obscuring Torus
 ctype black angle 0

bhlog3 # blr
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 
 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 3 0 3
 #---BH
 SET lr=lg(2) SET theta=0,2*pi,0.01*pi 
 set x=lr*cos(theta) set y=lr*sin(theta) ctype ltgrey SHADE 0 x y ctype black
 relocate $(lg(1.3)) $(lg(1.3)) putlabel 5 BH
 #---BLR
 angle 45
 relocate 2.5 1.1 putlabel 5 Broad Line Region
 angle 63
 ctype grey set x=1.5,4,0.01 set y=x-1.5 shade histogram 899 x y
 ctype black angle 90
 relocate 2 0.16 putlabel 6 He II
 relocate 2.48 0.248 putlabel 6 C IV, Si IV
 relocate 2.95 0.330 putlabel 6 H\beta, Mg II, C III]
 ctype black angle 0
 #---disk: 
 #set x=0,3,0.01 set y=0.066*x ctype dkgrey shade histogram 0 x y
 set x=lg(6),3,0.01 set y=0.1*(x-lg(6)) ctype ltgrey shade histogram 0 x y
 #set x=lg(9),3,0.01 set y=0.11*(x-lg(9)) ctype grey shade histogram 0 x y
 relocate $(lg(30)) $(lg(2.83)) putlabel 5 UV/optical
 relocate $(lg(30)) $(lg(2.38)) putlabel 5 continuum
 relocate $(lg(30)) $(lg(2.00)) putlabel 5 emission
 relocate $(lg(30)) $(lg(1.68)) putlabel 5 region
 ctype dkgrey 
 relocate $(lg(4)) $(lg(2.0)) putlabel 5 X-ray 
 relocate $(lg(4)) $(lg(1.68)) putlabel 5 emission
 relocate $(lg(4)) $(lg(1.41)) putlabel 5 region
 ctype black angle 0
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.50,2,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x) 
 set ly=ly concat 3\n set y=y concat 1000
 set lx=lx concat 1.785\n set x=x concat 61
 set r=lg(sqrt(x**2.+y**2.)) set theta=atan(y/x) 
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=px/3 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=2*px connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 0.5 1.4\n putlabel 5 "    (possible)"
 relocate 0.5 1.3\n putlabel 5 "    jet launching and"
 relocate 0.5 1.2\n putlabel 5 "    collimation region"
 relocate 0.4 2.9\n putlabel 5 jet
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]='10'\n SET vlab[2]='100'\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 1 2 3 4 5}
 SET vsma = 0 concat $(lg(2)) concat $(lg(6)) concat $(lg(9)) concat $(lg(10)) concat $(lg(20)) concat $(lg(30)) concat $(lg(40)) concat $(lg(50)) concat $(lg(60)) concat $(lg(70)) concat $(lg(80)) concat $(lg(90)) concat $(lg(100)) concat $(lg(200)) concat $(lg(300)) concat $(lg(400)) concat $(lg(500)) concat $(lg(600)) concat $(lg(700)) concat $(lg(800)) concat $(lg(900)) concat $(lg(1000)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 #limits $(lg(10.)) $(lg(10000.)) $(lg(10.)) $(lg(10000.))
 #SET vaub = $(lg(10)) concat $(lg(100)) concat $(lg(1000)) concat $(lg(10000))
 #SET vaus = vaub
 #SET dimen(vaul)=4.s\n SET vaul[0] = '10'\n SET vaul[1] = '100'\n SET vaul[2] = '10^3'\n SET vaul[3]='10^4'
 #angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 #ctype black ylabel \raise+750 radius in units of AU for M_{BH}=10^9 M_{Sun}
 limits $(lg(0.058)) $(lg(58.)) $(lg(0.0579)) $(lg(58.))
 SET vldb = $(lg(0.058)) concat $(lg(0.58)) concat $(lg(5.8)) concat $(lg(58.))
 SET vlds = vldb
 SET dimen(vldl)=4.s\n SET vldl[0] = '0.058'\n SET vldl[1] = '0.58'\n SET vldl[2] = '5.8'\n SET vldl[3] = '58.'
 angle 90 AXIS $fy1 $fy2 vlds vldb vldl $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-days for M_{BH}=10^9 M_{Sun}
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 ptype 20 3\n expand 0.26
 ctype grey\n set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=lg(xs+3.25) \n set ys0=4.25,61.25,9.5
 set ys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 set lys=lg(ys) \n points lxs lys
 limits 0 3 0 3 expand 1.2
 relocate 1.5 2.9 label stars\n ctype black

bh4 # torus
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 
 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 1E4 0 1E4
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 ptype 20 3\n expand 0.26
 limits $((0.0)) $((580.)) $((0.0)) $((580.)) \n ctype ltgrey
 set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=lg(xs) \n set ys0=4.25,61.25,9.5
 set ys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 set i=0,8 foreach var i {
  set lys=(ys+$var*65.5)\n points lxs lys
  set j=0,8 foreach rav j {
   set lxs=(xs+$rav*65.5)\n points lxs lys
  }
 }
 limits 0 1E4 0 1E4
 expand 1.2\n relocate 8600 9600\n ctype grey\n putlabel 6 stars\n ctype black
 relocate 9800 6100\n putlabel 4 inner
 relocate 9800 5750\n putlabel 4 edge
 relocate 9800 5400\n putlabel 4 of
 relocate 9800 5050\n putlabel 4 narrow
 relocate 9800 4700\n putlabel 4 line
 relocate 9800 4350\n putlabel 4 region
 #---disk:
 set x=(6),7E3,0.01 set y=0.1*(x-(6)) ctype ltgrey shade histogram 0 x y
 ctype black angle 0
 relocate $(10**3.93) 1000 putlabel 5 Clumpy
 relocate $(10**3.93)  500 putlabel 5 disk?
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.40,4,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x)
 set r=(sqrt(x**2.+y**2.)) set theta=atan(y/x)
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 460 9700\n putlabel 6 (possible)
 relocate 460 9300\n putlabel 6 jet
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]=''\n SET vlab[2]=''\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 0 0 1E3 1E4 1E5}
 SET vsma = 0 concat $((1000)) concat $((2000)) concat $((3000)) concat $((4000)) concat $((5000)) concat $((6000)) concat $((7000)) concat $((8000)) concat $((9000)) concat $((10000)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 limits 0 981 0 981 #---orbital
 SET vborb = 0 concat 98.1 concat 490.5 concat 981 SET vsorb=vborb
 SET dimen(vlorb)=4.s\n SET vlorb[1] = '31'\n SET vlorb[2] = '347'\n SET vlorb[3] = '981'
 AXIS $fx1 $fx2 vsorb vborb vlorb $gx1 $($gy1/2.) $($gx2-$gx1) 1 0 # basement
 ctype black xlabel orbital time in units of years for M_{BH}=10^9 M_{Sun}
 limits $((0.0)) $((30.)) $((0.0)) $((3E1))
 SET vaub = $((0.0)) concat $((0.0)) concat $((0.0)) concat $((3.)) concat $((29.999))
 SET vaus = vaub
 SET dimen(vaul)=5.s\n SET vaul[0] = '0'\n SET vaul[1] = ''\n SET vaul[2] = ''\n SET vaul[3]='0.16'\n SET vaul[4]='1.6'
 angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-years for M_{BH}=10^9 M_{Sun}
 limits 2 5.5 1.5 5.0
 ctype grey circle 2.6 2 0.09 circle 3.1 2.1 0.06 circle 3.9 3.1 0.05
 ellipse 3.1 2.8 0.7 0.1 53\n circle 3.4 2.4 0.07 circle 3.55 2.75 0.08
 relocate 3.4 2.9 angle 53 putlabel 5 BAL clouds
 limits 0.5 4.5 0.0 4.0
 ctype grey circle 2.6 2 0.07 circle 3.1 2.1 0.05 circle 3.9 3.1 0.03 
 ellipse 3.1 2.8 0.35 0.1 53\n circle 3.4 2.4 0.06 circle 3.55 2.75 0.04 
 limits 0.5 4.5 0.0 4.0
 ctype grey circle 2.6 2 0.07 circle 3.1 2.1 0.05 circle 3.9 3.1 0.03 
 limits 0.5 4.5 0.8 4.8\n circle 3.4 2.4 0.02 circle 3.55 2.75 0.06
 ctype grey circle 2.6 2 0.04 circle 3.1 2.1 0.05 circle 3.9 3.1 0.03
 angle 0 ctype black
 limits 0 1E4 0 1E4
 #---BLR
 angle 63
 ctype grey set x=50,1000,1 set y=0.48*(x-50)+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 set z=0.1*x+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 set yz=y concat reverse(z) set xx = x concat reverse(x)
 shade 899 xx yz
 ctype black angle 90
 angle 30\n relocate 550 400 putlabel 6 BELR
 ctype black angle 0
 #---torus:
 ctype dkgrey
 set x1=550,7500,0.01 set y1=0.400*(x1-550.)+0.1*(x1-6)
 set x2=750,7500,0.01 set y2=0.400*(x2-750.)+0.1*(x2-6)
 set x12 = x1 concat reverse(x2)\n set y12 = y1 concat reverse(y2)
 shade 0 x12 y12 
 ctype black
 angle 27\n relocate 4000 1350 ctype black putlabel 6 Obscuring Torus
 ctype black angle 0

bh4b 3 # torus NLRlabely BLRangle BALyoff
 # bh4b 4500 12 -0.2
 # bh4b 6000 22 0
 # bh4b 7000 33 0.3
 # bh4b 9000 45 0.66
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 
 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 1E4 0 1E4
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 expand 0.26
 limits $((0.0)) $((580.)) $((0.0)) $((580.)) \n ctype ltgrey
 set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=lg(xs) \n set ys0=4.25,61.25,9.5
 set ys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 set i=0,8 foreach var i {
  set lys=(ys+$var*65.5)\n points lxs lys
  set j=0,8 foreach rav j {
   set lxs=(xs+$rav*65.5)\n points lxs lys
  }
 }
 limits 0 1E4 0 1E4
 expand 1.2\n relocate 8600 9600\n ctype grey\n putlabel 6 stars\n ctype black
 relocate 9800 $1\n putlabel 4 inner
 relocate 9800 $($1-350)\n putlabel 4 edge
 relocate 9800 $($1-700)\n putlabel 4 of
 relocate 9800 $($1-1050)\n putlabel 4 narrow
 relocate 9800 $($1-1400)\n putlabel 4 line
 relocate 9800 $($1-1750)\n putlabel 4 region
 #---disk:
 set x=(6),7E3,0.01 set y=0.1*(x-(6)) ctype ltgrey shade histogram 0 x y
 ctype black angle 0
 relocate $(10**3.93) 1000 putlabel 5 Clumpy
 relocate $(10**3.93)  500 putlabel 5 disk?
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.40,4,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x)
 set r=(sqrt(x**2.+y**2.)) set theta=atan(y/x)
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 460 9700\n putlabel 6 (possible)
 relocate 460 9300\n putlabel 6 jet
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]=''\n SET vlab[2]=''\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 0 0 1E3 1E4 1E5}
 SET vsma = 0 concat $((1000)) concat $((2000)) concat $((3000)) concat $((4000)) concat $((5000)) concat $((6000)) concat $((7000)) concat $((8000)) concat $((9000)) concat $((10000)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 limits 0 981 0 981 #---orbital
 SET vborb = 0 concat 98.1 concat 490.5 concat 981 SET vsorb=vborb
 SET dimen(vlorb)=4.s\n SET vlorb[1] = '31'\n SET vlorb[2] = '347'\n SET vlorb[3] = '981'
 AXIS $fx1 $fx2 vsorb vborb vlorb $gx1 $($gy1/2.) $($gx2-$gx1) 1 0 # basement
 ctype black xlabel orbital time in units of years for M_{BH}=10^9 M_{Sun}
 limits $((0.0)) $((30.)) $((0.0)) $((3E1))
 SET vaub = $((0.0)) concat $((0.0)) concat $((0.0)) concat $((3.)) concat $((29.999))
 SET vaus = vaub
 SET dimen(vaul)=5.s\n SET vaul[0] = '0'\n SET vaul[1] = ''\n SET vaul[2] = ''\n SET vaul[3]='0.16'\n SET vaul[4]='1.6'
 angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-years for M_{BH}=10^9 M_{Sun}
 limits 2 5.5 1.5 5.0
 ctype grey circle 2.6 $($3+2) 0.09 circle 3.1 $($3+2.1) 0.06 circle 3.9 $($3+3.1) 0.05
 ellipse 3.1 $($3+2.8) 0.7 0.1 53\n circle 3.4 $($3+2.4) 0.07 circle 3.55 $($3+2.75) 0.08
 relocate 3.4 $($3+2.9) angle 53 putlabel 5 BAL clouds
 limits 0.5 4.5 0.0 4.0
 ctype grey circle 2.6 $($3*1.5+2) 0.07 circle 3.1 $($3*1.5+2.1) 0.05 circle 3.9 $($3*1.5+3.1) 0.03 
 ellipse 3.1 $($3*1.5+2.8) 0.35 0.1 53\n circle 3.4 $($3*1.5+2.4) 0.06 circle 3.55 $($3*1.5+2.75) 0.04 
 ctype grey circle 2.6 $($3*1.5+2) 0.07 circle 3.1 $($3*1.5+2.1) 0.05 circle 3.9 $($3*1.5+3.1) 0.03 
 limits 0.5 4.5 0.8 4.8\n circle 3.4 $($3*2+2.4) 0.02 circle 3.55 $($3*2+2.75) 0.06
 ctype grey circle 2.6 $($3*2+2) 0.04 circle 3.1 $($3*2+2.1) 0.05 circle 3.9 $($3*2+3.1) 0.03
 angle 0 ctype black
 limits 0 1E4 0 1E4
 #---BLR
 #angle 63
 #ctype grey set x=50,1000,1 set y=0.48*(x-50)+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 #set z=0.1*x+0.048*(x-500)*10**(abs(x-500.5)/(x-500.5))
 #set yz=y concat reverse(z) set xx = x concat reverse(x)
 #shade 899 xx yz
 ctype black angle 90
 angle $(atand(0.1+tand($2)))\n relocate 350 450 putlabel 6 BELR
 ctype black angle 0
 #---torus:
 ctype dkgrey
 set xlim=550.+9500./$(tand($(atand(0.1+tand($2)))))
 set uplim=7000.\n if ($(uplim) < $(xlim)) {set xlim=$(uplim)}
 set x1=550,$(xlim),0.01 set y1=tand($2)*(x1-550.)+0.1*(x1-6)
 set x2=750,$(xlim),0.01 set y2=tand($2)*(x2-750.)+0.1*(x2-6)
 set x12 = x1 concat reverse(x2)\n set y12 = y1 concat reverse(y2)
 shade 0 x12 y12 
 ctype black
 angle $(atand(0.1+tand($2)))
 set xpos=550.+3450./$(tand($(atand(0.1+tand($2)))))
 set yr=reverse(y2)\n set ypos=$((yr[0])*0.9)
 relocate $(xlim) $(ypos)\n putlabel 1 Obscuring Torus
 ctype black angle 0

bhlog4 # torus
 erase ltype 0 lweight 3 caray expand 1.2 ctype black
 define x_gutter 0 
 define y_gutter 0 
 window 10 10 3:10 3:10 limits 0 4 0 4
 #---BH
 SET lr=lg(2) SET theta=0,2*pi,0.01*pi 
 set x=lr*cos(theta) set y=lr*sin(theta) ctype ltgrey SHADE 0 x y ctype black
 relocate $(lg(1.3)) $(lg(1.3)) putlabel 5 BH
 #---stars, based on MW: Schodel+ adsabs.harvard.edu/abs/2007A%26A...469..125S
 expand 0.26
 ctype ltgrey\n limits $(lg(0.058)) $(lg(580.)) $(lg(0.0579)) $(lg(580.))
 set xs={1 1 1 1 1 1 1 20 20 20 20 20 20 20 39 39 39 39 39 39 39 58 58 58 58 58 58 58 10.5 10.5 10.5 10.5 10.5 10.5 10.5 29.5 29.5 29.5 29.5 29.5 29.5 29.5 48.5 48.5 48.5 48.5 48.5 48.5 48.5} \n set lxs=lg(xs) \n set ys0=4.25,61.25,9.5
 set ys=ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0 concat ys0
 set i=0,8 foreach var i {
  set lys=lg(ys+$var*65.5)\n points lxs lys
  set j=0,8 foreach rav j {
   set lxs=lg(xs+$rav*65.5)\n points lxs lys
  }
 }
 limits 0 4 0 4 expand 1.2
 relocate 1.5 3.9 label stars\n ctype black
 #---BLR
 angle 63
 relocate 2.68 1.9 putlabel 5 "Broad Line Region               "
 ctype grey 
 set x0=1.5,3.90,0.01 set y0=1.25*(x0-1.5) 
 set x3=2.5,3.90,0.01 set y3=2*(x3-2.5) 
 set x03 = x0 concat reverse(x3)\n set y03 = y0 concat reverse(y3)
 shade histogram 899 x03 y03
 ctype black angle 90
 relocate 2 0.16 putlabel 6 He II
 relocate 2.48 0.248 putlabel 6 C IV, Si IV
 #relocate 2.95 0.330 putlabel 6 H\beta, Mg II
 ctype black angle 0
 #---torus:
 ctype dkgrey 
 set x1=2.45,4,0.01 set y1=2.000*(x1-2.45) 
 set x2=2.95,4,0.01 set y2=2.000*(x2-2.95) 
 set x12 = x1 concat reverse(x2)\n set y12 = y1 concat reverse(y2)
 shade histogram 0 x12 y12
 ctype black
 angle 90 relocate 2.95 1.200 putlabel 6 H\beta, Mg II\n angle 0
 angle 63 relocate 3.2 1 ctype white label Obscuring `Torus'
 ctype black angle 0
 #---disk: 
 #set x=0,4,0.01 set y=0.066*x ctype dkgrey shade histogram 0 x y
 set x=lg(6),4,0.01 set y=0.1*(x-lg(6)) ctype ltgrey shade histogram 0 x y
 #set x=lg(9),4,0.01 set y=0.11*(x-lg(9)) ctype grey shade histogram 0 x y
 relocate $(lg(20)) $(lg(3.90)) putlabel 5 UV/optical
 relocate $(lg(20)) $(lg(2.90)) putlabel 5 continuum
 relocate $(lg(20)) $(lg(2.15)) putlabel 5 emission
 relocate $(lg(20)) $(lg(1.50)) putlabel 5 region
 ctype dkgrey 
 relocate $(lg(4)) $(lg(2.0)) putlabel 5 X-ray 
 relocate $(lg(4)) $(lg(1.68)) putlabel 5 emission
 relocate $(lg(4)) $(lg(1.31)) putlabel 5 region
 ctype black angle 0
 relocate 3.68 0.2 putlabel 5 Clumpy
 relocate 3.68 0.1 putlabel 5 disk?
 #---jet: launched at 28-46 R_g in M87; Hada et al. 2011
 set ly=1.50,2,0.01 set y=10**ly set x=y/sqrt(y-23) set lx=lg(x) 
 set ly=ly concat 4\n set y=y concat 10000
 set lx=lx concat 2.785\n set x=x concat 610
 set r=lg(sqrt(x**2.+y**2.)) set theta=atan(y/x) 
 set px=r*cos(theta) set py=r*sin(theta)
 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=px/3 connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 set px=2*px connect px py set rx = reverse(px) set ry = reverse(py)
 draw_arrow $(rx[1]) $(ry[1]) $(rx[0]) $(ry[0])
 relocate 0.5 1.3\n putlabel 5 "        jet launching and"
 relocate 0.5 1.2\n putlabel 5 "        collimation region"
 relocate 0.5 3.9\n putlabel 5 jet
 #---axes
 SET dimen(vlab)=6.s\n SET vlab[1]='10'\n SET vlab[2]='100'\n SET vlab[3]='10^3'\n SET vlab[4]='10^4'\n SET vlab[5]='10^5'
 #SET vlab[0]='0'
 SET vbig = {0 1 2 3 4 5}
 SET vsma = 0 concat $(lg(2)) concat $(lg(6)) concat $(lg(9)) concat $(lg(10)) concat $(lg(20)) concat $(lg(30)) concat $(lg(40)) concat $(lg(50)) concat $(lg(60)) concat $(lg(70)) concat $(lg(80)) concat $(lg(90)) concat $(lg(100)) concat $(lg(200)) concat $(lg(300)) concat $(lg(400)) concat $(lg(500)) concat $(lg(600)) concat $(lg(700)) concat $(lg(800)) concat $(lg(900)) concat $(lg(1000)) 
 ctype black AXIS $fx1 $fx2 vsma vbig vlab $gx1 $gy1 $($gx2-$gx1) 1 0	# bottom
 ctype black xlabel radius in units of r_g=GM_{BH}/c^2
 SET vsm = vsma
 angle 90 AXIS $fy1 $fy2 vsm vbig vlab $gy1 $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of r_g=GM_{BH}/c^2
 #limits $(lg(10.)) $(lg(100000.)) $(lg(10.)) $(lg(100000.))
 #SET vaub = $(lg(10)) concat $(lg(100)) concat $(lg(1000)) concat $(lg(10000)) concat $(lg(100000))
 #SET vaus = vaub
 #SET dimen(vaul)=5.s\n SET vaul[0] = '10'\n SET vaul[1] = '100'\n SET vaul[2] = '10^3'\n SET vaul[3]='10^4'\n SET vaul[4]='10^5'
 #angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 #ctype black ylabel \raise+750 radius in units of AU for M_{BH}=10^9 M_{Sun}
 limits $(lg(0.003)) $(lg(30.)) $(lg(0.003)) $(lg(3E1))
 SET vaub = $(lg(0.003)) concat $(lg(0.03)) concat $(lg(0.3)) concat $(lg(3.)) concat $(lg(29.999))
 SET vaus = vaub
 SET dimen(vaul)=5.s\n SET vaul[0] = '1.6E-4'\n SET vaul[1] = '1.6E-3'\n SET vaul[2] = '0.016'\n SET vaul[3]='0.16'\n SET vaul[4]='1.6'
 angle 90 AXIS $fy1 $fy2 vaus vaub vaul $(0.5*$gy1) $gx1 $($gy2-$gy1) 1 1 angle 0 # left
 ctype black ylabel \raise+750 radius in units of light-years for M_{BH}=10^9 M_{Sun}
 limits 0 4 0 4
 ctype grey circle 2.9 3 0.09 circle 3.1 2.1 0.06 circle 3.9 3.1 0.05
 ellipse 3.1 2.8 0.7 0.1 53\n circle 3.4 2.4 0.07 circle 3.55 2.75 0.08
 relocate 3.4 2.9 angle 53 putlabel 5 BAL clouds
 angle 0 ctype black
 limits 0 4 0 4

caray
 add_ctype a 255 0   0   add_ctype b 255 80  0   add_ctype d 255 128 0
 add_ctype f 255 208 0   add_ctype g 255 255 0   add_ctype h 215 255 0
 add_ctype j 128 255 0   add_ctype l 0   255 0   add_ctype n 0   255 176
 add_ctype p 0   255 255 add_ctype q 0   192 255 add_ctype r 0   128 255
 add_ctype s 0   64  255 add_ctype t 0   0   255 add_ctype u 80  0   255
 add_ctype v 176 0   255 add_ctype x 255 0   255 add_ctype y 255 0   212
 add_ctype z 255 0   128 add_ctype omega 0 0 0	 add_ctype brown 255 99 99
 set carray = { a b d f g h j l n p q r s t u v x y z omega } # 0 to 19 (20 entries)
 add_ctype ltgrey 192 192 192 \n add_ctype grey 128 128 128 \n add_ctype dkgrey 64 64 64


wlq11b # GN11B WLQ fig 3
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.0
 define file_type fits \n define table_type bintable
 window 1 3 1 1:3 \n xlabel Rest-frame Wavelength in \AA
 ylabel Normalized F_\lambda
 #ylabel F_\lambda\ (10^{-17} erg s^{-1} cm^{-2} \AA^{-1})

 window 1 3 1 3 \n limits 912 1575 0 2.75 #14
 image "/Users/phall/CURRENT/WILD/SKYSUB/0895/spSpec-52581-0895-271_skysubDR.fit"
 define ZHW image\n define NAXIS1 image\n define COEFF0 image\n define COEFF1 image
 set zed = $ZHW \n set _i = 0, $($NAXIS1 - 6) \n set wl = $COEFF0 + $COEFF1*_i
 set wl = 10.**(wl)/(1+zed) \n set fl = image[_i,0] \n set er = image[_i,1]
 set fl = fl/5
 box 0 2 0 0 \n smooth fl fl7 7 \n connect wl fl7
 relocate 1275 0.45 \n putlabel 6 J083122 #1.5
 data "/Users/phall/IRAF/PHIST/sdss_qso_median.tab1"
 ltype 0 \n lweight 1 \n read { wl 1 fl 2 } \n set fl = fl / 6.5 \n connect wl fl \n lweight 2 \n ltype 0
 relocate 1225 2.45 \n putlabel 6 Ly\alpha+NV
 relocate 1535 1.75 \n putlabel 5 CIV

 window 1 3 1 2 \n limits 3750 9250 0 20 \n limits 912 1575 0 1.75 #3.8
 image "/Users/phall/CURRENT/WILD/SKYSUB/0740/spSpec-52263-0740-627_skysubDR.fit"
 define ZHW image\n define NAXIS1 image\n define COEFF0 image\n define COEFF1 image
 set zed = $ZHW \n set _i = 0, $($NAXIS1 - 6) \n set wl = $COEFF0 + $COEFF1*_i
 set wl = 10.**(wl)/(1+zed) \n set fl = image[_i,0] \n set er = image[_i,1]
 set fl = fl/2
 box 0 2 0 0 \n smooth fl fl13 13 \n connect wl fl13
 relocate 1275 0.30 \n putlabel 6 J225246 #0.75

 window 1 3 1 1 \n limits 3750 9250 0 20 \n limits 912 1575 0 1.75 #7
 image "/Users/phall/CURRENT/WILD/SKYSUB/0655/spSpec-52162-0655-482_skysubDR.fit"
 define ZHW image\n define NAXIS1 image\n define COEFF0 image\n define COEFF1 image
 set zed = $ZHW \n set _i = 0, $($NAXIS1 - 6) \n set wl = $COEFF0 + $COEFF1*_i
 set wl = 10.**(wl)/(1+zed) \n set fl = image[_i,0] \n set er = image[_i,1]
 set fl=fl/4
 box 1 2 0 0 \n smooth fl fl7 7 \n connect wl fl7
 relocate 1275 0.30 \n putlabel 6 J004054 #1.25

 #VDB way - check it
        #define CRVAL1 image
        #define CD1_1 image
        #set s_log_l_min  = sprintf('%f',$CRVAL1)
        #set s_log_dl     = sprintf('%f',$CD1_1)
        #set fl = image[*,0]  # 'SPECTRUM'/ units of (10^-17 erg/s/cm^2/A)
        #set s_log_l = fl        do i=0, dimen(fl)-1 {
          #set s_log_l[$i] = $(s_log_l_min) + $i * $(s_log_dl)
        #}
        #set wl       = 10**(s_log_l)
        #set er = image[*,1]


leighlify # PHL 1811 analog & WLQ unification, Leighly shout-out
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 
 add_ctype dgrey 192 192 192 \n add_ctype bgrey 128 128 128 \n add_ctype lgrey 64 64 64
 #add_ctype lgrey 190 190 190 \n add_ctype dgrey 140 140 140 \n add_ctype bgrey 65 65 65

 do i = 2, 1, -1 {
  window 1 2 1 $i \n limits -1.00 1.00 -0.66 1.14

 # regions
  if ($i == 2) {	# top
   xlabel \raise+20700 \bf Leighly variable-height X-ray source model #Normal Quasars
   #set xs = {0.075 0.1 0.150 0.150} \n set ys = {0 0.20 0.300 0} \n set zs = -1*ys
   set xx = {-0.025 -0.05 0.05 0.025} \n set yx = {0.15 0.70 0.70 0.15} \n set zx = -1*yx
  }
  if ($i == 1) {	# bottom
   xlabel \raise+11500 \bf Unified WLQs / PHL 1811 Analogs	# or \raise+2000
   #set xs = {0.075 0.1125 0.150 0.150} \n set ys = {0 0.25 0.350 0} \n set zs = -1*ys
   set xx = {-0.025 -0.05 0.05 0.025} \n set yx = {0.05 0.10 0.10 0.05} \n set zx = -1*yx
  }
  ctype lgrey \n shade 0 xx yx \n shade 0 xx zx \n ctype default # X-ray source
  set xs = {0.075 0.11 0.165 0.165} \n set ys = {0 0.20 0.275 0} \n set zs = -1*ys
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  set xe = {0.220 0.242 0.440 0.440} \n set ye = {0.0 0.330 0.600 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe \n set xa = -1*xa
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  if ($i == 2) { # top
   expand 1.2 \n ctype white \n angle 90 \n relocate 0 -0.40 \n putlabel 5 X-ray
   relocate 0 0.40 \n putlabel 5 source\n angle 0 \n ctype default \n expand 1.5
   ltype 3
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.3125 -0.225 0.375 -0.300
   draw_arrow 0.3125  0.225 0.375  0.300
   expand 1.5
  }
  if ($i == 1) { # bottom
   expand 1.2 \n ctype black \n angle 90 \n relocate 0 -0.40 \n putlabel 5 X-ray
   relocate 0 0.40 \n putlabel 5 source\n angle 0 \n ctype default \n expand 1.5
   ltype 3
   draw_arrow -0.05 0.2 -0.2 0.9 #\n draw_arrow 0.0625 0.2 0.25 0.9 
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.0 1.0 \n putlabel 5 Weak-Line Quasars
   draw_arrow 0.05 0.055 0.65 0.715
   draw_arrow 0.05 0.075 0.6 0.9
   relocate 0.625 0.925 \n putlabel 6 PHL 1811 \n relocate 0.675 0.725 \n putlabel 6 analogs
   ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.250 -0.225 0.420 -0.400
   draw_arrow 0.250  0.200 0.420  0.375
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.440 0 \n draw 0.440 0	# disk
  set xt = { -0.44 -0.88 -0.88 -0.44 -0.44 } \n set yt = {  0.33  0.66 -0.66 -0.33  0.33 }
  set zt = -1*xt \n ctype default \n shade 0 xt yt \n shade 0 zt yt	# torus
  relocate 0 0 \n ptype 32 3 \n expand 1.1 \n dot \n expand 1.5	# X-ray source
  relocate -0.325 0.075 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  ctype white \n relocate -0.53 0 \n expand 1.1 \n putlabel 5 torus \n expand 1.5
  ctype default
  relocate 0 -0.80 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.75 -0.1 -0.15 \n draw_arrow 0.075 -0.75 0.1 -0.15 \n expand 1.5
 }


progafy # PHL 1811 analog & WLQ unification, Proga shout-out
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 
 add_ctype dgrey 192 192 192 \n add_ctype bgrey 128 128 128 \n add_ctype lgrey 64 64 64

 do i = 2, 1, -1 {
  window 1 2 1 $i 
 # regions
  if ($i == 2) {	# top
   limits -1.00 1.00 -0.95 0.85
   xlabel \raise+20700 \bf Normal Quasars a la Proga  # Leave comment symbol as is (centers text)
   set xs = {0.075 0.1 0.120 0.125} \n set ys = {0 0.025 0.040 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 normal\n expand 0.3
   set xd = {0.075 0.1 0.150 0.155} \n set yd = {0 0.080 0.120 0} \n set zd = -1*yd \n set wd = -1*xd
   ltype 1 \n connect xd yd \n connect xd zd \n connect wd yd \n connect wd zd \n ltype 0	# dotted outline
  }
  if ($i == 1) {	# bottom
   limits -1.00 1.00 -0.66 1.14
   xlabel \raise+11500 \bf Unified WLQs / PHL 1811 Analogs    # or \raise+2000
   set xs = {0.075 0.1125 0.150 0.165} \n set ys = {0 0.225 0.325 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 high-f_C\n expand 0.3
  }
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  #set xe = {0.200 0.220 0.440 0.440} \n set ye = {0.0 0.300 0.600 0.0} \n set ze = -1*ye
  set xe = {0.200 0.220 0.48 0.400} \n set ye = {0.0 0.300 0.64 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe 
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  relocate 0 -0.600 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.375 -0.1 -0.175 \n draw_arrow 0.075 -0.375 0.1 -0.175 \n expand 1.5
  if ($i == 2) { # top
   ltype 3
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.3125 -0.225 0.375 -0.300
   draw_arrow 0.3125  0.225 0.375  0.300
   expand 1.5
  }
  if ($i == 1) { # bottom
   ltype 3
   draw_arrow -0.05 0.2 -0.2 0.9 #\n draw_arrow 0.0625 0.2 0.25 0.9 
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate -0.2 1.0 \n putlabel 5 Weak-Line Quasars
   draw_arrow 0.05 0.055 0.65 0.715
   draw_arrow 0.05 0.075 0.6 0.9
   relocate 0.555 1.0 \n putlabel 6 PHL 1811 \n relocate 0.605 0.8 \n putlabel 6 analogs
   ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.250 -0.225 0.420 -0.400
   draw_arrow 0.250  0.200 0.420  0.375
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.440 0 \n draw 0.440 0	# disk
  #draw 0.44 0.33 \n draw 0.65 0.33 \n draw 0.65 -0.33 \n draw 0.44 -0.33 \n draw 0.44 0
  #relocate -0.44 0.33 \n draw -0.65 0.33 \n draw -0.65 -0.33 \n draw -0.44 -0.33 \n draw -0.44 0.33
  #set xt = { -0.44 -0.88 -0.88 -0.44 -0.44 } \n set yt = {  0.33  0.66 -0.66 -0.33  0.33 }
  set xt = {-0.4 -0.44 -0.88 -0.88 -0.44 -0.4} \n set yt = {0 0.33 0.66 -0.66 -0.33 0}
  	set zt = -1*xt \n ctype lgrey \n shade 0 xt yt \n shade 0 zt yt \n ctype default	# TORUS
  relocate 0 0 \n ptype 32 3 \n expand 2 \n dot \n expand 1.5	# X-ray source
  relocate -0.3000 0.125 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  #relocate -0.55 0.075 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
  	ctype white \n relocate -0.63 0 \n expand 1.1 \n putlabel 5 TORUS \n expand 1.5
  ctype default
 }


unify4 # PHL 1811 analog & WLQ unification
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 
 add_ctype dgrey 192 192 192 \n add_ctype bgrey 128 128 128 \n add_ctype lgrey 64 64 64

 do i = 2, 1, -1 {
  window 1 2 1 $i 
 # regions
  if ($i == 2) {	# top
   limits -1.00 1.00 -0.95 0.85
   xlabel \raise+20700 \bf Normal Quasars  # Leave comment symbol as is (centers text)
   set xs = {0.075 0.1 0.120 0.125} \n set ys = {0 0.025 0.040 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 normal\n expand 0.3
   set xd = {0.075 0.1 0.150 0.155} \n set yd = {0 0.080 0.120 0} \n set zd = -1*yd \n set wd = -1*xd
   ltype 1 \n connect xd yd \n connect xd zd \n connect wd yd \n connect wd zd \n ltype 0	# dotted outline
  }
  if ($i == 1) {	# bottom
   limits -1.00 1.00 -0.66 1.14
   xlabel \raise+10750 \bf Unified WLQs / PHL 1811 Analogs    # or \raise+2000
   set xs = {0.075 0.1125 0.150 0.165} \n set ys = {0 0.225 0.325 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 high-f_C\n expand 0.3
  }
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  #set xe = {0.200 0.220 0.440 0.440} \n set ye = {0.0 0.300 0.600 0.0} \n set ze = -1*ye
  set xe = {0.200 0.220 0.48 0.400} \n set ye = {0.0 0.300 0.64 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe 
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  relocate 0 -0.600 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.375 -0.1 -0.175 \n draw_arrow 0.075 -0.375 0.1 -0.175 \n expand 1.5
  if ($i == 2) { # top
   ltype 3
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.3125 -0.225 0.375 -0.300
   draw_arrow 0.3125  0.225 0.375  0.300
   expand 1.5
  }
  if ($i == 1) { # bottom
   ltype 3
   draw_arrow -0.05 0.2 -0.2 0.9 #\n draw_arrow 0.0625 0.2 0.25 0.9 
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate -0.2 1.0 \n putlabel 5 Weak-Line Quasars
   draw_arrow 0.05 0.055 0.65 0.715
   draw_arrow 0.05 0.075 0.6 0.9
   relocate 0.555 1.0 \n putlabel 6 PHL 1811 \n relocate 0.605 0.8 \n putlabel 6 analogs
   ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.250 -0.225 0.420 -0.400
   draw_arrow 0.250  0.200 0.420  0.375
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.440 0 \n draw 0.440 0	# disk
  #draw 0.44 0.33 \n draw 0.65 0.33 \n draw 0.65 -0.33 \n draw 0.44 -0.33 \n draw 0.44 0
  #relocate -0.44 0.33 \n draw -0.65 0.33 \n draw -0.65 -0.33 \n draw -0.44 -0.33 \n draw -0.44 0.33
  #set xt = { -0.44 -0.88 -0.88 -0.44 -0.44 } \n set yt = {  0.33  0.66 -0.66 -0.33  0.33 }
  set xt = {-0.4 -0.44 -0.88 -0.88 -0.44 -0.4} \n set yt = {0 0.33 0.66 -0.66 -0.33 0}
  	set zt = -1*xt \n ctype lgrey \n shade 0 xt yt \n shade 0 zt yt \n ctype default	# TORUS
  relocate 0 0 \n ptype 32 3 \n expand 2 \n dot \n expand 1.5	# X-ray source
  relocate -0.3000 0.125 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  #relocate -0.55 0.075 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
  	ctype white \n relocate -0.63 0 \n expand 1.1 \n putlabel 5 TORUS \n expand 1.5
  ctype default
 }


unify3 # PHL 1811 analog & WLQ unification
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 
 add_ctype dgrey 192 192 192 \n add_ctype bgrey 128 128 128 \n add_ctype lgrey 64 64 64

 do i = 2, 1, -1 {
  window 1 2 1 $i 
 # regions
  if ($i == 2) {	# top
   limits -1.00 1.00 -0.95 0.85
   xlabel \raise+20700 \bf Normal Quasars  # Leave comment symbol as is (centers text)
   set xs = {0.075 0.1 0.150 0.155} \n set ys = {0 0.05 0.075 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 normal\n expand 0.3
  }
  if ($i == 1) {	# bottom
   limits -1.00 1.00 -0.66 1.14
   xlabel \raise+11500 \bf Unified WLQs / PHL 1811 Analogs    # or \raise+2000
   set xs = {0.075 0.1125 0.150 0.165} \n set ys = {0 0.225 0.325 0} \n set zs = -1*ys
   relocate 0 -0.450 \n expand 1.2 \n putlabel 5 high-f_C\n expand 0.3
  }
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  set xe = {0.200 0.220 0.440 0.440} \n set ye = {0.0 0.300 0.600 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe \n set xa = -1*xa
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  relocate 0 -0.600 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.375 -0.1 -0.175 \n draw_arrow 0.075 -0.375 0.1 -0.175 \n expand 1.5
  if ($i == 2) { # top
   ltype 3
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.3125 -0.225 0.375 -0.300
   draw_arrow 0.3125  0.225 0.375  0.300
   expand 1.5
  }
  if ($i == 1) { # bottom
   ltype 3
   draw_arrow -0.05 0.2 -0.2 0.9 #\n draw_arrow 0.0625 0.2 0.25 0.9 
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.0 1.0 \n putlabel 5 Weak-Line Quasars
   draw_arrow 0.05 0.055 0.65 0.715
   draw_arrow 0.05 0.075 0.6 0.9
   relocate 0.625 0.925 \n putlabel 6 PHL 1811 \n relocate 0.675 0.725 \n putlabel 6 analogs
   ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.250 -0.225 0.420 -0.400
   draw_arrow 0.250  0.200 0.420  0.375
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.440 0 \n draw 0.440 0	# disk
  #draw 0.44 0.33 \n draw 0.65 0.33 \n draw 0.65 -0.33 \n draw 0.44 -0.33 \n draw 0.44 0
  #relocate -0.44 0.33 \n draw -0.65 0.33 \n draw -0.65 -0.33 \n draw -0.44 -0.33 \n draw -0.44 0.33
  set xt = { -0.44 -0.88 -0.88 -0.44 -0.44 } \n set yt = {  0.33  0.66 -0.66 -0.33  0.33 }
  set zt = -1*xt \n ctype lgrey \n shade 0 xt yt \n shade 0 zt yt \n ctype default	# TORUS
  relocate 0 0 \n ptype 32 3 \n expand 2 \n dot \n expand 1.5	# X-ray source
  relocate -0.3125 0.125 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  #relocate -0.55 0.075 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
  ctype white \n relocate -0.63 0 \n expand 1.1 \n putlabel 5 TORUS \n expand 1.5
  ctype default
 }


unify2 # PHL 1811 analog & WLQ unification
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 \n add_ctype lgrey 190 190 190
 add_ctype dgrey 192 192 192 \n add_ctype bgrey 128 128 128
 #add_ctype dgrey 140 140 140 \n add_ctype bgrey 65 65 65 

 do i = 2, 1, -1 {
  window 1 2 1 $i \n limits -0.65 1.35 -0.65 1.15

 # regions
  if ($i == 2) {	# top
   xlabel \raise+20700 \bf Normal Quasars
   set xs = {0.075 0.1 0.150 0.150} \n set ys = {0 0.05 0.075 0} \n set zs = -1*ys
  # BALR clouds
   #set xa = {0.5 0.585 0.55} \n set ya = {0.100 0.05 -0.175} 
   #set wa = -1*xa \n set za = -1*ya \n ctype lgrey \n expand 3 
   #ptype 7 3 \n points xa ya \n points wa za \n expand 1.5 \n ctype default
  }
  if ($i == 1) {	# bottom
   xlabel \raise+11500 \bf Unified WLQs / PHL 1811 Analogs	# or \raise+2000
   set xs = {0.075 0.1125 0.150 0.150} \n set ys = {0 0.25 0.350 0} \n set zs = -1*ys
  # BALR clouds
   #set xa = {0.5 0.55 0.585 0.5 0.55} \n set ya = {0.15 0.45 0.05 -0.45 -0.200} 
   #set wa = -1*xa \n set za = -1*ya \n ctype lgrey \n expand 3 
   #ptype 7 3 \n points xa ya \n points wa za \n expand 1.5 \n ctype default
  }
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  set xe = {0.220 0.220 0.440 0.440} \n set ye = {0.0 0.300 0.600 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe \n set xa = -1*xa
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  relocate 0 -0.45 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.35 -0.1 -0.15 \n draw_arrow 0.075 -0.35 0.1 -0.15 \n expand 1.5
  if ($i == 2) { # top
   ltype 3
   #draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.100 1.0 \n putlabel 5 Normal Quasar, Low Blueshift
   #expand 1.3 
   #draw_arrow 0.067 0.05  0.8 0.6 \n relocate 0.8 0.6 \n putlabel 6 Normal Quasar,
   #relocate 0.8 0.45 \n putlabel 6 Medium Blueshift 
   #draw_arrow 0.06 -0.024 0.8 -0.325 \n relocate 0.8 -0.325 \n putlabel 6 BAL Quasar
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   #draw_arrow 0.300 -0.125 0.375 -0.150
   #draw_arrow 0.300  0.125 0.375  0.150
   draw_arrow 0.3125 -0.225 0.375 -0.300
   draw_arrow 0.3125  0.225 0.375  0.300
   expand 1.5
  }
  if ($i == 1) { # bottom
   ltype 3
   #draw_arrow -0.25 0.4 -0.20 0.9 \n draw_arrow 0.35 0.4  0.40 0.9
   draw_arrow -0.05 0.2 -0.2 0.9 #\n draw_arrow 0.0625 0.2 0.25 0.9 
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.0 1.0 \n putlabel 5 Weak-Line Quasars
   draw_arrow 0.05 0.055 0.65 0.715
   draw_arrow 0.05 0.075 0.6 0.9
   relocate 0.625 0.925 \n putlabel 6 PHL 1811 \n relocate 0.675 0.725 \n putlabel 6 analogs
   #draw_arrow 0.05 -0.025  0.8 -0.325
   #expand 1.25 \n relocate 0.8 -0.25 \n putlabel 6 PHL 1811 \n relocate 0.8 -0.4 \n putlabel 6 analog \n relocate 0.8 -0.55 \n putlabel 6 BAL Quasar \n expand 1.5
   ltype 0
  # outflow vectors
   expand 0.5 
   #draw_arrow 0.250 -0.0625 0.420 -0.1000
   #draw_arrow 0.250  0.0625 0.420  0.1000
   draw_arrow 0.250 -0.225 0.420 -0.400
   draw_arrow 0.250  0.200 0.420  0.375
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.440 0 \n draw 0.440 0	# disk
  #draw 0.44 0.33 \n draw 0.65 0.33 \n draw 0.65 -0.33 \n draw 0.44 -0.33 \n draw 0.44 0
  #relocate -0.44 0.33 \n draw -0.65 0.33 \n draw -0.65 -0.33 \n draw -0.44 -0.33 \n draw -0.44 0.33
  set xt = { -0.44 -0.65 -0.65 -0.44 -0.44 } \n set yt = {  0.33  0.33 -0.33 -0.33  0.33 }
  set zt = -1*xt \n ctype default \n shade 0 xt yt \n shade 0 zt yt
  relocate 0 0 \n ptype 32 3 \n expand 2 \n dot \n expand 1.5	# X-ray source
  relocate -0.325 0.075 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  #relocate -0.55 0.075 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
  ctype white \n relocate -0.53 0 \n expand 1.1 \n putlabel 5 torus \n expand 1.5
  ctype default
 }


unify1 # PHL 1811 analog & WLQ unification
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 \n add_ctype bgrey 65 65 65 
 add_ctype dgrey 140 140 140 \n add_ctype lgrey 190 190 190

 do i = 2, 1, -1 {
  window 1 2 1 $i \n limits -0.65 1.35 -0.65 1.15

 # regions
  if ($i == 2) {	# top
   xlabel \raise+20700 \bf Normal Quasars
   set xs = {0.075 0.100 0.150 0.150} \n set ys = {0 0.05 0.0625 0} \n set zs = -1*ys
  # BALR clouds
   set xa = {0.5 0.585 0.55} \n set ya = {0.100 0.05 -0.175} 
   set wa = -1*xa \n set za = -1*ya \n ctype lgrey \n expand 3 
   ptype 7 3 \n points xa ya \n points wa za \n expand 1.5 \n ctype default
  }
  if ($i == 1) {	# bottom
   xlabel \raise+11500 \bf Unified WLQs / PHL 1811 Analogs	# or \raise+2000
   set xs = {0.075 0.125 0.150 0.150} \n set ys = {0 0.175 0.225 0} \n set zs = -1*ys
  # BALR clouds
   set xa = {0.5 0.55 0.585 0.5 0.55} \n set ya = {0.15 0.45 0.05 -0.45 -0.200} 
   set wa = -1*xa \n set za = -1*ya \n ctype lgrey \n expand 3 
   ptype 7 3 \n points xa ya \n points wa za \n expand 1.5 \n ctype default
  }
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  #set xe = {0.225 0.225 0.450 0.450} \n set ye = {-0.1 0.225 0.450 -0.1} \n set ze = -1*ye
  set xe = {0.225 0.225 0.450 0.450} \n set ye = {0.0 0.225 0.450 0.0} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xs = -1*xs \n set xe = -1*xe \n set xa = -1*xa
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR

 # labels and arrows 
  relocate 0 -0.45 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.3
  draw_arrow -0.075 -0.35 -0.1 -0.15 \n draw_arrow 0.075 -0.35 0.1 -0.15 \n expand 1.5
  if ($i == 2) { # top
   ltype 3
   #draw_arrow -0.25 0.4 -0.20 0.9 \n draw_arrow 0.35 0.4  0.40 0.9
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.100 1.0 \n putlabel 5 Normal Quasar, Low Blueshift
   expand 1.3 
   draw_arrow 0.067 0.05  0.8 0.6 \n relocate 0.8 0.6 \n putlabel 6 Normal Quasar,
   relocate 0.8 0.45 \n putlabel 6 Medium Blueshift 
   draw_arrow 0.06 -0.024 0.8 -0.325 \n relocate 0.8 -0.325 \n putlabel 6 BAL Quasar
   expand 1.5 \n ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.300 -0.075 0.375 -0.100
   draw_arrow 0.300  0.075 0.375  0.100
   draw_arrow 0.300 -0.175 0.375 -0.225
   draw_arrow 0.300  0.175 0.375  0.225
   expand 1.5
  }
  if ($i == 1) { # bottom
   ltype 3
   #draw_arrow -0.25 0.4 -0.20 0.9 \n draw_arrow 0.35 0.4  0.40 0.9
   draw_arrow 0.025 0.2  0.10 0.9 \n relocate 0.100 1.0 \n putlabel 5 Weak-Line Quasar
   draw_arrow 0.05 0.025  0.8 0.4
   relocate 0.8 0.4 \n putlabel 6 PHL 1811 \n relocate 0.8 0.2 \n putlabel 6 analog
   draw_arrow 0.05 -0.025  0.8 -0.325
   expand 1.25 \n relocate 0.8 -0.25 \n putlabel 6 PHL 1811 \n relocate 0.8 -0.4 \n putlabel 6 analog \n relocate 0.8 -0.55 \n putlabel 6 BAL Quasar \n expand 1.5
   ltype 0
  # outflow vectors
   expand 0.5 
   draw_arrow 0.275 -0.0625 0.400 -0.1000
   draw_arrow 0.275  0.0625 0.400  0.1000
   draw_arrow 0.275 -0.175 0.400 -0.275
   draw_arrow 0.275  0.175 0.400  0.275
   expand 1.5
  }
 
  box 4 4 3 4
  relocate -0.725 0 \n draw 0.725 0	# disk
  relocate 0 0 \n ptype 32 3 \n expand 2 \n dot \n expand 1.5	# X-ray source
  relocate -0.325 0.075 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  relocate -0.55 0.075 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
 }


unify0 # PHL 1811 analog & WLQ unification FIRST DRAFT
 erase \n ltype 0 \n lweight 2 \n ctype default \n expand 1.5
 define TeXstrings 1 \n define y_gutter 0.60 \n add_ctype bgrey 65 65 65 
 add_ctype dgrey 140 140 140 \n add_ctype lgrey 190 190 190

 do i = 2, 1, -1 {
  window 1 2 1 $i \n limits -0.65 1.35 -0.65 1.15

  if ($i == 2) { # top
   xlabel \raise+20700 \bf Normal Quasar
   set xs = {0.075 0.125 0.175 0.175} \n set ys = {0 0.05 0.05 0} \n set zs = -1*ys
  }
  if ($i == 1) { # bottom
   xlabel \raise+11500 \bf Unified WLQ / PHL 1811 Analog # or \raise+2000
   set xs = {0.075 0.125 0.175 0.175} \n set ys = {0 0.175 0.200 0} \n set zs = -1*ys
  }

  relocate 0 -0.45 \n expand 1.2 \n putlabel 5 shielding gas \n expand 0.5
  draw_arrow -0.075 -0.35 -0.1 -0.15 \n draw_arrow 0.075 -0.35 0.1 -0.15 \n expand 1.5

  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  set xe = {0.225 0.225 0.450 0.450} \n set ye = {-0.1 0.225 0.450 -0.1} \n set ze = -1*ye
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  set xa = {0.45 0.45 0.60 0.60} \n set ya = {0 0.45 0.60 0} \n set za = -1*ya 
  #ctype lgrey \n shade 0 xa ya \n shade 0 xa za \n ctype default	# BALR
  set xs = -1*xs \n set xe = -1*xe \n set xa = -1*xa
  ctype bgrey \n shade 0 xs ys \n shade 0 xs zs \n ctype default	# shield
  ctype dgrey \n shade 0 xe ye \n shade 0 xe ze \n ctype default	# BELR
  #ctype lgrey \n shade 0 xa ya \n shade 0 xa za \n ctype default	# BALR
 # BALR clouds
  set xb = {0.5 0.55 0.585 0.5 0.55} \n set yb = {0.15 0.45 0 -0.45 -0.275} 
  set wb = -1*xb \n set zb = -1*yb \n ctype lgrey \n expand 3 
  ptype 7 3 \n points xb yb \n points wb zb \n expand 1.5 \n ctype default

 # labels and arrows 
  if ($i == 2) { # top
   draw_arrow 0.1 0.05  0.8 0.4
   draw_arrow 0.1 -0.05 0.8 -0.4
   relocate 0.8 -0.4 \n putlabel 6 BAL Quasar
  }
  if ($i == 1) { # bottom
   draw_arrow -0.25 0.4 -0.20 0.9
   draw_arrow 0.35 0.4  0.40 0.9
   draw_arrow 0.025 0.2  0.10 0.9
   relocate 0.100 1.0 \n putlabel 5 Weak-Line Quasar
   draw_arrow 0.05 0.025  0.8 0.4
   relocate 0.8 0.4 \n putlabel 6 PHL 1811 \n relocate 0.8 0.2 \n putlabel 6 analog
   draw_arrow 0.05 -0.025  0.8 -0.4
   relocate 0.8 -0.2 \n putlabel 6 PHL 1811 \n relocate 0.8 -0.4 \n putlabel 6 analog \n relocate 0.8 -0.6 \n putlabel 6 BAL Quasar
  }
 
  box 4 4 3 4
  relocate -0.225 0 \n draw 0.225 0	# disk
  relocate 0 0 \n ptype 32 3 \n expand 4 \n dot \n expand 1.5	# X-ray source
  relocate 0.35 0 \n expand 1.2 \n putlabel 5 BELR \n expand 1.5
  relocate 0.55 -0.125 \n expand 1.1 \n putlabel 5 BALR \n expand 1.5
 }


sidehist 2 #Draw a histogram on the side of a figure
            #sidehist x y
   define _size (dimen($1) - 2)
   define _startx ($1[0])
   define _starty ($2[0])
   relocate $_starty $_startx
   do _i = 1, $_size{
    define _nextx ($1[$_i])
    define _nexty ($2[$_i])
    draw $_nexty $_nextx
    define _nextx ($1[$($_i + 1)])
    draw $_nexty $_nextx
   }
   foreach i (_size _startx _starty _i _nextx _nexty) {
   define $i delete}

zM      #plot z vs. M_g (sidehist demonstration)
        
        era \n lw 3 \n expand 1.8 \n add_ctype grey 180 180 180

        da dr3qso.targfinal.i.lfinput.dat \n read {zem 1 i 2 Mi 3}

        window -16 -16 3:16 3:16 \n limits 0 5.3 -22 -29.99
        box 0 0 0 0 \n ptype 1 1
        ctype default \n points zem Mi \n ptype 4 1

        set y1 = <-22.5 -29.7>
        set u1 = <0.3 5.0>

        ctype grey

        define vv -29.7
        while {$vv<=-22.5} {
          set v1 = <$vv $vv>
          define vv $($vv+0.3)
           connect u1 v1
        }

        define vv 0.3
        while {$vv<=2.2} {
          set v1 = <$vv $vv>
          define vv $($vv+0.38)
           connect v1 y1
        }

        define vv 2.2
        while {$vv<=3.0} {
          set v1 = <$vv $vv>
          define vv $($vv+0.4)
           connect v1 y1
        }
        
        define vv 3.0
        while {$vv<=5.0} {
          set v1 = <$vv $vv>
          define vv $($vv+0.5)
           connect v1 y1
        }

        
        set zz = 0,6,0.01
          
        define HoN 70  #HoN stands for New Hubble constant.
        define alpha -0.50
        
        set K = -2.5*(1.0+$alpha)*LG(1.0+zz) +
2.5*(1.0+$alpha)*LG(1.0+2.0)

        define s $((0.7/0.3)**(1/3))
        define s2 $($s*$s)
        define s3 $(0.7/0.3)
        define s4 $($s2*$s2)
        define lead $(2.0*sqrt($s3+1))
        define last $(0.066941*$s4)

        set n1 =
$lead*(1-0.1540*$s+0.4304*$s2+0.19097*$s3+$last)**(-0.125)
        set n2 =
$lead*((1+zz)**4-0.1540*$s*(1+zz)**3+0.4304*$s2*(1+zz)**2+0.19097*$s3*(1+zz)+$last)**(-0.125)

        set Dl = ((3e11*(1+zz))/$HoN)*(n1-n2)

        set Mi150 = 15.0 + 5 - 5*LG(Dl) - K
        set Mi191 = 19.1 + 5 - 5*LG(Dl) - K
        set Mi202 = 20.2 + 5 - 5*LG(Dl) - K

        ltype 2
        connect zz Mi150
        connect zz Mi191
        connect zz Mi202
        ltype 0
        
        ctype default
        
        da dps.kcorr.gfix.dat
        read {z 1 corrg 3 corri 5}

        set kcorrcont = -2.5*(1.0-0.436)*lg(1.0+z)
        set corri2 = corri-kcorrcont
        set Mi191k = Mi191 - corri2 if (zz < 3.0)
        set zz1 = zz if (zz < 3.0)
        
        set Mi202k = Mi202 - corri2 if (zz >= 3.0 && zz < 5.5)
        set zz2 = zz if (zz >= 3.0 && zz < 5.5)

        ctype red
        connect zz1 Mi191k
        connect zz2 Mi202k
        ctype default

        window -16 -16 1:2 3:16
        limits 0 2000 -22 -29.99
        box 0 1 1 0
        set mbin = -30.0,-18.0,0.25
        set mhist = histogram(Mi:mbin)
        sidehist mbin mhist
        xlabel N
        ylabel \raise+500M_i

        window -16 -16 3:16 1:2
        limits 0 5.3 0 500
        box 1 0 0 1
        set zbin2 = 0,5,0.05
        set zhist = histogram(zem:zbin2)
        histogram zbin2 zhist
        xlabel Redshift
        
        window 1 1 1 1
        lw 1
        expand 1.0001