-
Notifications
You must be signed in to change notification settings - Fork 3
/
plot_qflux_histogram.pro
executable file
·103 lines (70 loc) · 2.75 KB
/
plot_qflux_histogram.pro
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
PRO plot_qflux_histogram,filename,qflux_final=qflux_final,ps=ps,chisq=chisq
restore,filename,/verbose
delvarx,te_map,em_map,sig_map,chi_map,tsig_errmap_symmetric,telog_errmap_symmetric,temperature_map,emission_map,sigma_map
big=size(aia_simul_map_cube.data)
qflux_final=fltarr(big[1],big[2],6)
!p.multi=[0,2,3]
wave=['131','171','193','211','335','94']
IF keyword_set(ps) THEN begin
set_plot,'ps'
device,encaps=1,bits_per_pixel=8,ysize=20,filename='qflux_histogram.ps'
ENDIF ELSE BEGIN
window,17,xsize=1024,ysize=1024
ENDELSE
for n=0,5 do begin
qflux_final(*,*,n)=float(aia_simul_map_cube(n).data(*,*))/float(aia_map_cube(n).data(*,*))
endfor
for i=0,big[1]-1 do begin
for j=0,big[2]-1 do begin
for k=0,5 do begin
IF (qflux_final(i,j,k) gt 10.) OR (qflux_final(i,j,k) le 0.0) THEN qflux_final(i,j,k)=10.
endfor
endfor
endfor
histo0=histogram(reform(qflux_final(*,*,0)),binsize=0.1,locations=locations0,max=3)
histo1=histogram(reform(qflux_final(*,*,1)),binsize=0.1,locations=locations1,max=3)
histo2=histogram(reform(qflux_final(*,*,2)),binsize=0.1,locations=locations2,max=3)
histo3=histogram(reform(qflux_final(*,*,3)),binsize=0.1,locations=locations3,max=3)
histo4=histogram(reform(qflux_final(*,*,4)),binsize=0.1,locations=locations4,max=3)
histo5=histogram(reform(qflux_final(*,*,5)),binsize=0.1,locations=locations5,max=3)
plot,locations0,histo0,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[0]+'A',/right,charsize=1.5
plot,locations1,histo1,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[1]+'A',/right,charsize=1.5
plot,locations2,histo2,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[2]+'A',/right,charsize=1.5
plot,locations3,histo3,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[3]+'A',/right,charsize=1.5
plot,locations4,histo4,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[4]+'A',/right,charsize=1.5
plot,locations5,histo5,psym=10,xrange=[0,3],charsize=1.5,ytitle='N',xtitle='qflux'
al_legend,wave[5]+'A',/right,charsize=1.5
IF keyword_set(ps) THEN BEGIN
device,/close
set_plot,'x'
ENDIF
IF keyword_set(chisq) THEN BEGIN
!p.multi=[0,1,1]
FOR i=0,big[1]-1 do begin
FOR j=0,big[2]-1 do begin
IF (chisq_map.data(i,j) gt 100.) THEN chisq_map.data(i,j)=100.
endfor
endfor
histo_chi=histogram((chisq_map.data(*,*)),binsize=0.02,locations=locations6,max=10)
IF keyword_set(ps) THEN BEGIN
set_plot,'ps'
device,encaps=1,bits_per_pixel=8,filename='chisq_histogram.ps'
ENDIF ELSE BEGIN
window,18,xsize=512,ysize=512
ENDELSE
plot,locations6,histo_chi,psym=10,xrange=[0,6],charsize=1.5,ytitle='N',xtitle='chisq'
IF keyword_set(ps) THEN BEGIN
device,/close
set_plot,'x'
ENDIF
ENDIF
;endfor
;device,/close
;set_plot,'x'
!p.multi=[0,1,1]
END