1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
openu,1,unit
a=assoc(1,bytarr(dim,dim,/nozero))
dens=a(0)
openu,2,"scale-idl"
b=assoc(2,bytarr(30,256,/nozero))
sca=b(0)
d=size(dens)
s=size(sca)
if ( color eq 0 ) then begin
for i=0, dim-1 do begin
for j=0, dim-1 do begin
if (dens(i,j) eq 0) then dens(i,j)=1
endfor
endfor
max = 255
expo = 2
r = findgen(max+1)
g = findgen(max+1)
for i=0, max do begin
g(i) = max * ( 1 - r(i)/max ) ^ expo
endfor
tvlct, g, g, g
black = 255
white = 0
ccol_1 = black
ccol_2 = black
ccol_3 = black
endif else begin
for i=0, dim-1 do begin
for j=0, dim-1 do begin
if (dens(i,j) eq 0) then dens(i,j)=1
endfor
endfor
loadct,color
max = 255
r = findgen(max+1)
g = findgen(max+1)
b = findgen(max+1)
tvlct, r, g, b, /get
r(0) = 255
g(0) = 255
b(0) = 255
tvlct, r, g, b
black = 1
white = 0
ccol_1 = black
ccol_2 = black
ccol_3 = black
endelse
set_plot,"ps"
device, /encapsulated, file=idl2eps, /color, bits=8
xoff = !d.x_px_cm*2
yoff = 0
s1x = !d.x_px_cm*10
s1y = !d.y_px_cm*10
; s2x = !d.x_px_cm
s2x = !d.x_px_cm*0.5
; s2y = !d.y_px_cm*10
s2y = !d.y_px_cm*8
; dist = !d.x_px_cm*4
dist = !d.x_px_cm
tv, dens, xsize=s1x, ysize=s1y, xoff, yoff
x = xmin + (xmax-xmin)/d(1)*findgen(d(1))
y = ymin + (ymax-ymin)/d(2)*findgen(d(2))
zcrit = 256/zmax
contour, dens, /noerase, xstyle=4, ystyle=4, $
level = [level1*zcrit, level2*zcrit, level3*zcrit], $
c_color = [ccol_1, ccol_2, ccol_3], $
position=[xoff,yoff,xoff+s1x,yoff+s1y],/device
plot, x, y, xstyle=1, ystyle=1, xrange=[xmin,xmax], yrange=[ymin,ymax],$
xtitle=xname, ytitle=yname,$
charsize=1.75, color=black, /nodata, /noerase, /device,$
position=[xoff,yoff,xoff+s1x,yoff+s1y]
; xoff = xoff + s1x + dist
; yoff = yoff
xoff = xoff + !d.x_px_cm * 0.5
yoff = yoff + !d.y_px_cm
tv, sca, xsize=s2x, ysize=s2y, xoff, yoff
x = 1.0/s(1)*findgen(s(1))
y = zmax/s(2)*findgen(s(2))
; contour, sca, /noerase, xstyle=4, ystyle=1, yrange=[zmin,zmax], charsize = 2, $
contour, sca, /noerase, xstyle=4, ystyle=4, yrange=[zmin,zmax], $
levels = [level1, level2, level3], $
c_color = [ccol_1, ccol_2, ccol_3], $
color = black, $
position = [xoff,yoff,xoff+s2x,yoff+s2y],/device
; plot, x, y, xstyle = 4, ystyle = 1, yrange = [0,zmax],$
; ytitle = yname, charsize = 2,$
; color = black, /nodata, /noerase, /device,$
; position = [xoff,yoff,xoff+s2x,yoff+s2y]
axis, xaxis=0, color=black, xticks=1, xtickname=replicate(' ',2), xminor=1
axis, xaxis=1, color=black, xticks=1, xtickname=replicate(' ',2), xminor=1
axis, yaxis=0, color=black, yticks=1, ytickname=replicate(' ',2)
axis, yaxis=1, color=black, yrange=[0,zmax], ystyle=1, ytitle=zname, ycharsize=1.5
device, /close
close,1
close,2
end
|