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
|
pro idlplot2ps
dim= 400
xmin = -1.0
xmax = 1.0
ymin = 0
ymax = 20
zmax = 0.02
level1 = 3
level2 = 4
level3 = 5
color = 1
openu,1,"dens.idl"
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)
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
if ( color eq 0 ) then begin
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
loadct,4
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, file="res.ps", /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
s2y = !d.y_px_cm*10
dist = !d.x_px_cm*4
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="!3x/!7k!3!D0!3", ytitle="!3t/!7s!3",$
charsize=2, color=black, /nodata, /noerase, /device,$
position=[xoff,yoff,xoff+s1x,yoff+s1y]
xoff = xoff + s1x + dist
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=[0,zmax],$
charsize = 2, $
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 = "!3E!N/E!Dr!3", charsize = 2,$
color = black, /nodata, /noerase, /device,$
position = [xoff,yoff,xoff+s2x,yoff+s2y]
device, /close
close,1
close,2
end
|