;;;plot the # of passes; # of counts; # of counts/passes in each bin ;;;with equal size (MLAT*MLT). start_end_year = indgen(10)+2001 ;abs_cut1 = 0.1e-14 ;abs_cut2 = 3.0e-14;1.0e-14 percent_cut = 30. ctime = fltarr(80000L)*!Values.F_NaN cmaxden = fltarr(80000L)*!Values.F_NaN cdelta_den = fltarr(80000L)*!Values.F_NaN cglat = fltarr(80000L)*!Values.F_NaN cglon = fltarr(80000L)*!Values.F_NaN cglt = fltarr(80000L)*!Values.F_NaN cmlat = fltarr(80000L)*!Values.F_NaN cmlon = fltarr(80000L)*!Values.F_NaN cmlt = fltarr(80000L)*!Values.F_NaN cby = fltarr(80000L)*!Values.F_NaN cbz = fltarr(80000L)*!Values.F_NaN cpress = fltarr(80000L)*!Values.F_NaN csymh = fltarr(80000L)*!Values.F_NaN ;cpass_number = fltarr(80000L)*!Values.F_NaN npts0 = 0 for jj = 0, n_elements(start_end_year)-1 do begin ;begin major analysis loop, 1 file/loop. print, start_end_year(jj) ;;;;;read in ---------------------------------------------- champ = read_champ_maxden('./CHAMP_360km_AACGM_peaks_30%_'+strtrim(start_end_year(jj),2)+'.dat') time = champ.time maxden = champ.maxden delta_den = champ.delta_den percent = champ.percent glat = champ.glat glon = champ.glon glt = champ.glt mlat = champ.mlat mlon = champ.mlon mlt = champ.mlt by = champ.by bz = champ.bz pre =champ.press symh = champ.symh ;pass_number = champ.number_pass npts = n_elements(time) ctime(npts0:npts0+npts-1) = time cmaxden(npts0:npts0+npts-1) = maxden cdelta_den(npts0:npts0+npts-1) = delta_den cglat(npts0:npts0+npts-1) = glat cglon(npts0:npts0+npts-1) = glon cglt(npts0:npts0+npts-1) = glt cmlat(npts0:npts0+npts-1) = mlat cmlon(npts0:npts0+npts-1) = mlon cmlt(npts0:npts0+npts-1) = mlt cby(npts0:npts0+npts-1) = by cbz(npts0:npts0+npts-1) = bz cpress(npts0:npts0+npts-1) = pre csymh(npts0:npts0+npts-1) = symh ;cpass_number(npts0:npts0+npts-1) = pass_number npts0 = npts0+npts endfor ;end major analysis loop. ;;;read in original pass coordinates cpass_mlat = fltarr(15000000L)*!Values.F_NaN cpass_mlt = fltarr(15000000L)*!Values.F_NaN npts0 = 0 for jj = 0, n_elements(start_end_year)-1 do begin print, start_end_year(jj) filename = file_search('../CHAMP_MSIS_400_360_320km_'+strtrim(start_end_year(jj),2)+'.dat') data = read_champ_msis(filename) pass_mlat = data.mlat pass_mlt = data.mlt npts = n_elements(pass_mlat) cpass_mlat(npts0:npts0+npts-1) = pass_mlat cpass_mlt(npts0:npts0+npts-1) = pass_mlt npts0 = npts0+npts endfor ;;;;data binning ;Define bin size mlat_window = 1.0 ;5.0 binsz = 180./mlat_window ;2.0 MLat(90) mlt_window = 1.0 mlat_bin = (findgen(binsz+1)-(binsz/2.))*mlat_window mlt_bin = findgen(24./mlt_window+1) ;binave_peak = fltarr(24,binsz)*!Values.F_NaN bincount_peak = fltarr(24,binsz) ;*!Values.F_NaN bincount_pass = fltarr(24,binsz) bincount = fltarr(24,binsz) ;normalized peak counts max_cb = 100. ;125 min_cb = 0. ;100 max_cb2 = 10000. min_cb2 = 0. max_cb3 = .04 min_cb3 = 0. ;---------------------Figure(mp)---------------------------------------------------------------- set_plot,'ps' ;device,/inches,xsize=10.90,ysize=5.347,yoffset=10.9,xoffset=1.0,/color,/landscape device,filename='CHAMP_2001-2010_maxden_equalbin.ps',/portrait device,/color loadct, 33 ;xyouts,-49.0,44.0,'Poynting Flux(mW/m!U2!N)',charsize=1. plot, mlt_bin, mlt_bin*0.,/nodata, xr = [0, 24], xstyle = 1, yr=[-90,90], ystyle = 1, $ xticks = 24, yticks = 18,xtitle='MLT',ytitle='MLAT',title='# of peak counts (CHAMP)' for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 fp_bin = where((cmlat ge mlat_bin(imlat)) and (cmlat le mlat_bin(imlat+1)) $ and (cmlt ge mlt_bin(imlt)) and (cmlt le mlt_bin(imlt+1)), count_bin) if (count_bin ne 0) then begin ;temp = cpass_number(fp_bin) ;;print, 'pass_number = ', n_elements(temp(uniq(temp,sort(temp)))) ;;print, 'peak counts = ', count_bin ;binave(imlt,imlat) = n_elements(temp(uniq(temp,sort(temp)))) ;number of pass in a bin ;;binave(imlt,imlat) = 100.0*count_bin/n_elements(temp(uniq(temp,sort(temp)))) ;normalized counts bincount_peak(imlt,imlat) = count_bin bin_color = (bincount_peak(imlt,imlat)-min_cb)/(max_cb-min_cb)*256 ;bin_color = (binave(imlt,imlat)-min_cb)/(max_cb-min_cb)*256 bin_xpos1 = mlt_bin(imlt) bin_xpos2 = mlt_bin(imlt+1) bin_xpos3 = mlt_bin(imlt+1) bin_xpos4 = mlt_bin(imlt) bin_ypos1 = mlat_bin(imlat) bin_ypos2 = mlat_bin(imlat) bin_ypos3 = mlat_bin(imlat+1) bin_ypos4 = mlat_bin(imlat+1) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor ;;colorbar number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb ; max(alog10(temps),/nan) minData = min_cb ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;polar plot polar_dial_n, 1, 1 xyouts,-49.0,48.0,'peak counts',charsize=1.1 xyouts,-49.0,44.0,'Northern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 ;fp_bin = where((cmlat ge mlat_bin(imlat)) and (cmlat le mlat_bin(imlat+1)) $ ; and (cmlt ge mlt_bin(imlt)) and (cmlt le mlt_bin(imlt+1)), count_bin) if ((bincount_peak(imlt,imlat) ne 0) and (mlat_bin(imlat) ge 50.)) then begin ;temp = cpass_number(fp_bin) ;binave(imlt,imlat) = n_elements(temp(uniq(temp,sort(temp)))) ;number of pass in a bin ;;binave(imlt,imlat) = 100.0*count_bin/n_elements(temp(uniq(temp,sort(temp)))) ;normalized counts ;bincount_peak(imlt,imlat) = count_bin bin_color = (bincount_peak(imlt,imlat)-min_cb)/(max_cb-min_cb)*256 ;bin_color = (binave(imlt,imlat)-min_cb)/(max_cb-min_cb)*256 bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(90-i),2) + '!9%!5' ;northern hemisphere if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb ; max(alog10(temps),/nan) minData = min_cb ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 polar_dial_s, 1, 2 xyouts,-49.0,48.0,'peak counts',charsize=1.1 xyouts,-49.0,44.0,'Southern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 ;fp_bin = where((cmlat ge mlat_bin(imlat)) and (cmlat le mlat_bin(imlat+1)) $ ; and (cmlt ge mlt_bin(imlt)) and (cmlt le mlt_bin(imlt+1)), count_bin) if ((bincount_peak(imlt,imlat) ne 0) and (mlat_bin(imlat+1) le -50.)) then begin ;bincount_peak(imlt,imlat) = count_bin bin_color = (bincount_peak(imlt,imlat)-min_cb)/(max_cb-min_cb)*256 bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(i-90),2) + '!9%!5' if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb ; max(alog10(temps),/nan) minData = min_cb ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;histogram loadct, 39 !p.multi = [0,2,3] xbin_mlat = (findgen(binsz)-(binsz/2.))*mlat_window+mlat_window/2. mlat_bincount_peak = fltarr(binsz)*0.0 for imlat = 0, binsz-1 do begin mlat_bincount_peak(imlat) = total(bincount_peak(*,imlat)) endfor indexn = where(xbin_mlat ge 50.0) plot, xbin_mlat, mlat_bincount_peak, xtitle = 'mlat(NH)',ytitle='counts',charsize = 1.5,xr=[50, 90],xstyle=1,yr=[0,.2e4] temp = total(xbin_mlat(indexn)*mlat_bincount_peak(indexn))/total(mlat_bincount_peak(indexn)) xyouts,55,1000,'mean='+string(temp,format='(f6.2)'),charsize = 1. indexs = where(xbin_mlat le -50.0) plot, xbin_mlat, mlat_bincount_peak, xtitle = 'mlat(SH)',ytitle='counts',charsize = 1.5,xr=[-90,-50],xstyle=1,yr=[0,.2e4] temp = total(xbin_mlat(indexs)*mlat_bincount_peak(indexs))/total(mlat_bincount_peak(indexs)) xyouts,-70,1000,'mean='+string(temp,format='(f6.2)'),charsize = 1. xbin_mlt = findgen(24)*mlt_window+mlt_window/2. mlt_bincount_peak_nh = fltarr(24)*0.0 mlt_bincount_peak_sh = fltarr(24)*0.0 for imlt = 0, 23 do begin ;indexn = where(xbin_mlat ge 50.0) mlt_bincount_peak_nh(imlt) = total(bincount_peak(imlt,indexn)) ;indexs = where(xbin_mlat le -50.0) mlt_bincount_peak_sh(imlt) = total(bincount_peak(imlt,indexs)) endfor plot, xbin_mlt, mlt_bincount_peak_nh, xtitle = 'mlt(NH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1,yr=[0,1.e3],xticks = 4 temp = total(xbin_mlt*mlt_bincount_peak_nh)/total(mlt_bincount_peak_nh) xyouts,3,800,'mean='+string(temp,format='(f6.2)'),charsize = 1. plot, xbin_mlt, mlt_bincount_peak_sh, xtitle = 'mlt(SH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1,yr=[0,1.e3],xticks = 4 temp = total(xbin_mlt*mlt_bincount_peak_sh)/total(mlt_bincount_peak_sh) xyouts,12,800,'mean='+string(temp,format='(f6.2)'),charsize = 1. !p.multi = 0 ;;====================================pass counts=========================== loadct, 33 plot, mlt_bin, mlt_bin*0.,/nodata, xr = [0, 24], xstyle = 1, yr=[-90,90], ystyle = 1, $ xticks = 24, yticks = 18,xtitle='MLT',ytitle='MLAT',title='# of pass counts (CHAMP)' for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 fp_bin = where((cpass_mlat ge mlat_bin(imlat)) and (cpass_mlat le mlat_bin(imlat+1)) $ and (cpass_mlt ge mlt_bin(imlt)) and (cpass_mlt le mlt_bin(imlt+1)), count_bin) if (count_bin ne 0) then begin bincount_pass(imlt,imlat) = count_bin bin_color = (bincount_pass(imlt,imlat)-min_cb2)/(max_cb2-min_cb2)*256 bin_xpos1 = mlt_bin(imlt) bin_xpos2 = mlt_bin(imlt+1) bin_xpos3 = mlt_bin(imlt+1) bin_xpos4 = mlt_bin(imlt) bin_ypos1 = mlat_bin(imlat) bin_ypos2 = mlat_bin(imlat) bin_ypos3 = mlat_bin(imlat+1) bin_ypos4 = mlat_bin(imlat+1) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor ;;colorbar number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb2 ; max(alog10(temps),/nan) minData = min_cb2 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;polar plot polar_dial_n, 1, 1 xyouts,-49.0,48.0,'pass counts',charsize=1.1 xyouts,-49.0,44.0,'Northern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 if ((bincount_pass(imlt,imlat) ne 0) and (mlat_bin(imlat) ge 50.)) then begin ;bincount_pass(imlt,imlat) = count_bin bin_color = (bincount_pass(imlt,imlat)-min_cb2)/(max_cb2-min_cb2)*256 bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(90-i),2) + '!9%!5' ;northern hemisphere if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb2 ; max(alog10(temps),/nan) minData = min_cb2 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 polar_dial_s, 1, 2 xyouts,-49.0,48.0,'pass counts',charsize=1.1 xyouts,-49.0,44.0,'Southern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 if ((bincount_pass(imlt,imlat) ne 0) and (mlat_bin(imlat+1) le -50.)) then begin bin_color = (bincount_pass(imlt,imlat)-min_cb2)/(max_cb2-min_cb2)*256 bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(i-90),2) + '!9%!5' if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb2 ; max(alog10(temps),/nan) minData = min_cb2 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;histogram loadct, 39 !p.multi = [0,2,3] xbin_mlat = (findgen(binsz)-(binsz/2.))*mlat_window+mlat_window/2. mlat_bincount_pass = fltarr(binsz)*0.0 for imlat = 0, binsz-1 do begin mlat_bincount_pass(imlat) = total(bincount_pass(*,imlat)) endfor indexn = where(xbin_mlat ge 50.0) plot, xbin_mlat, mlat_bincount_pass, xtitle = 'mlat(NH)',ytitle='counts',charsize = 1.5,xr=[50, 90],xstyle=1,yr=[0,1.2e5] temp = total(xbin_mlat(indexn)*mlat_bincount_pass(indexn))/total(mlat_bincount_pass(indexn)) xyouts,55,50000,'mean='+string(temp,format='(f6.2)'),charsize = 1. indexs = where(xbin_mlat le -50.0) plot, xbin_mlat, mlat_bincount_pass, xtitle = 'mlat(SH)',ytitle='counts',charsize = 1.5,xr=[-90,-50],xstyle=1,yr=[0,1.2e5] temp = total(xbin_mlat(indexs)*mlat_bincount_pass(indexs))/total(mlat_bincount_pass(indexs)) xyouts,-70,50000,'mean='+string(temp,format='(f6.2)'),charsize = 1. xbin_mlt = findgen(24)*mlt_window+mlt_window/2. mlt_bincount_pass_nh = fltarr(24)*0.0 mlt_bincount_pass_sh = fltarr(24)*0.0 for imlt = 0, 23 do begin ;indexn = where(xbin_mlat ge 50.0) mlt_bincount_pass_nh(imlt) = total(bincount_pass(imlt,indexn)) ;indexs = where(xbin_mlat le -50.0) mlt_bincount_pass_sh(imlt) = total(bincount_pass(imlt,indexs)) endfor plot, xbin_mlt, mlt_bincount_pass_nh, xtitle = 'mlt(NH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1, $ yr=[5.6e4,6.2e4],ystyle = 1, xticks = 4 temp = total(xbin_mlt*mlt_bincount_pass_nh)/total(mlt_bincount_pass_nh) xyouts,3,6.e4,'mean='+string(temp,format='(f6.2)'),charsize = 1. plot, xbin_mlt, mlt_bincount_pass_sh, xtitle = 'mlt(SH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1, $ yr=[5.6e4,6.2e4],ystyle = 1, xticks = 4 temp = total(xbin_mlt*mlt_bincount_pass_sh)/total(mlt_bincount_pass_sh) xyouts,12,6.e4,'mean='+string(temp,format='(f6.2)'),charsize = 1. !p.multi = 0 ;;;;========================normalized peak;;;;counts==================================== loadct, 33 plot, mlt_bin, mlt_bin*0.,/nodata, xr = [0, 24], xstyle = 1, yr=[-90,90], ystyle = 1, $ xticks = 24, yticks = 18,xtitle='MLT',ytitle='MLAT',title='normalized peak counts (CHAMP)' for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 if ((bincount_peak(imlt,imlat) ne 0) and (bincount_pass(imlt,imlat) ne 0)) then begin bincount(imlt,imlat) = 1.0*bincount_peak(imlt,imlat)/bincount_pass(imlt,imlat) bin_color = (bincount(imlt,imlat)-min_cb3)/(max_cb3-min_cb3)*256 bin_xpos1 = mlt_bin(imlt) bin_xpos2 = mlt_bin(imlt+1) bin_xpos3 = mlt_bin(imlt+1) bin_xpos4 = mlt_bin(imlt) bin_ypos1 = mlat_bin(imlat) bin_ypos2 = mlat_bin(imlat) bin_ypos3 = mlat_bin(imlat+1) bin_ypos4 = mlat_bin(imlat+1) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor ;;colorbar number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb3 ; max(alog10(temps),/nan) minData = min_cb3 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;polar plot polar_dial_n, 1, 1 xyouts,-49.0,48.0,'normalized peak counts',charsize=1.1 xyouts,-49.0,44.0,'Northern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 if ((bincount_peak(imlt,imlat) ne 0) and (bincount_pass(imlt,imlat) ne 0) and (mlat_bin(imlat) ge 50.)) then begin bin_color = (bincount(imlt,imlat)-min_cb3)/(max_cb3-min_cb3)*256 bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.-mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(90-i),2) + '!9%!5' ;northern hemisphere if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb3 ; max(alog10(temps),/nan) minData = min_cb3 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 polar_dial_s, 1, 2 xyouts,-49.0,48.0,'normalized peak counts',charsize=1.1 xyouts,-49.0,44.0,'Southern Hemisphere',charsize=1.1 for imlt = 0, n_elements(mlt_bin)-2 do begin ;24 for imlat = 0, n_elements(mlat_bin)-2 do begin ;90 if ((bincount_peak(imlt,imlat) ne 0) and (bincount_pass(imlt,imlat) ne 0) and (mlat_bin(imlat+1) le -50.)) then begin bin_color = (bincount(imlt,imlat)-min_cb3)/(max_cb3-min_cb3)*256 bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos1 = bin_radius*cos(bin_angle) bin_ypos1 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt)-6.)*!pi/12. bin_xpos2 = bin_radius*cos(bin_angle) bin_ypos2 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat+1)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos3 = bin_radius*cos(bin_angle) bin_ypos3 = bin_radius*sin(bin_angle) bin_radius = (90.+mlat_bin(imlat)) bin_angle = (mlt_bin(imlt+1)-6.)*!pi/12. bin_xpos4 = bin_radius*cos(bin_angle) bin_ypos4 = bin_radius*sin(bin_angle) xx = [bin_xpos1,bin_xpos2,bin_xpos3,bin_xpos4] yy = [bin_ypos1,bin_ypos2,bin_ypos3,bin_ypos4] polyfill,xx,yy,color = bin_color endif endfor endfor loadct, 39 ;additional parts pi = 3.141592654 angle = 2.*pi*indgen(721)/720. for i=10.,40.,10. do begin radius = replicate(i,721) oplot,/polar,radius,angle,linestyle =2 endfor oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] for i=10,40,10 do begin tang = -45.*pi/180. txt = ' '+strtrim(string(i-90),2) + '!9%!5' if (i eq 40) then txt = txt + ' MLAT' xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize endfor xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5 xyouts,46.,0.,'6',charsize=csize,alignment=0. xyouts,0.,46.,'12',charsize=csize,alignment=0.5 xyouts,-46.,0.,'18',charsize=csize,alignment=1. ;;colorbar loadct, 33 number_levels = 33 levels_contour = dblarr(number_levels) colors_contour = dblarr(number_levels) color_table_offset =.0 maxData = max_cb3 ; max(alog10(temps),/nan) minData = min_cb3 ; min(alog10(temps),/nan) for i = 0, number_levels-1 do begin levels_contour(i) = (maxData-minData)*i/(number_levels-1.0)+minData colors_contour(i) = (255.0-color_table_offset)*i/(number_levels-1.0)+color_table_offset endfor contour, [[levels_contour],[levels_contour]],levels_contour,indgen(2),charsize=0.8,charthick=1, $ levels=levels_contour,c_color=colors_contour,/fill,position=[0.67,0.005,0.97,0.035],/noerase,ystyle=4, $ xr=[minData,maxData],xstyle=1,xticks = 5 ;;;histogram loadct, 39 !p.multi = [0,2,3] xbin = (findgen(binsz)-(binsz/2.))*mlat_window+mlat_window/2. mlat_bincount = fltarr(binsz)*0.0 for imlat = 0, binsz-1 do begin mlat_bincount(imlat) = total(bincount(*,imlat)) endfor indexn = where(xbin ge 50.0) plot, xbin, mlat_bincount, xtitle = 'mlat(NH)',ytitle='counts',charsize = 1.5,xr=[50, 90],xstyle=1,yr=[0,0.5] temp = total(xbin(indexn)*mlat_bincount(indexn))/total(mlat_bincount(indexn)) xyouts,55,0.2,'mean='+string(temp,format='(f6.2)'),charsize = 1. indexs = where(xbin le -50.0) plot, xbin, mlat_bincount, xtitle = 'mlat(SH)',ytitle='counts',charsize = 1.5,xr=[-90,-50],xstyle=1,yr=[0,0.5] temp = total(xbin(indexs)*mlat_bincount(indexs))/total(mlat_bincount(indexs)) xyouts,-70,0.2,'mean='+string(temp,format='(f6.2)'),charsize = 1. xbin_mlt = findgen(24)*mlt_window+mlt_window/2. mlt_bincount_nh = fltarr(24)*0.0 mlt_bincount_sh = fltarr(24)*0.0 for imlt = 0, 23 do begin ;indexn = where(xbin_mlat ge 50.0) mlt_bincount_nh(imlt) = total(bincount(imlt,indexn)) ;indexs = where(xbin_mlat le -50.0) mlt_bincount_sh(imlt) = total(bincount(imlt,indexs)) endfor plot, xbin_mlt, mlt_bincount_nh, xtitle = 'mlt(NH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1,yr=[0.1,.6], xticks = 4 temp = total(xbin_mlt*mlt_bincount_nh)/total(mlt_bincount_nh) xyouts,3,.4,'mean='+string(temp,format='(f6.2)'),charsize = 1. plot, xbin_mlt, mlt_bincount_sh, xtitle = 'mlt(SH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1,yr=[0.1,.6], xticks = 4 temp = total(xbin_mlt*mlt_bincount_sh)/total(mlt_bincount_sh) xyouts,12,.4,'mean='+string(temp,format='(f6.2)'),charsize = 1. !p.multi = 0 device,/close set_plot,'x' end