;;;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 ---------------------------------------------- filename = './CHAMP_360km_AACGM_peaks_30%_'+strtrim(start_end_year(jj),2)+'.dat' ;print, filename champ = read_champ_maxden(filename) 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. ;stop ;;;read in original pass coordinates ;cpass_time = fltarr(15000000L)*!Values.F_NaN cpass_mlat = fltarr(15000000L)*!Values.F_NaN cpass_mlt = fltarr(15000000L)*!Values.F_NaN cpass_number = intarr(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_time= data.doy pass_mlat = data.mlat pass_mlt = data.mlt npts = n_elements(pass_mlat) ;;;mark the pass number to each original data point. pass_number = intarr(npts)*!values.f_nan ;;seperate NH/SH lat_next = pass_time*!values.f_nan lat_next(0:npts-2) = pass_mlat(1:npts-1) multi = pass_mlat*lat_next fp_traj = where(multi le 0.0,ntraj) ; number of pass:ntraj+1 time_orbit = [pass_time(0), pass_time(fp_traj),pass_time(npts-1)] ; time for the change of NH/SH for ipoint = 0, npts-1 do begin pass_number(ipoint) = min(where(time_orbit ge pass_time(ipoint))) endfor ;cpass_time(npts0:npts0+npts-1) = pass_time cpass_mlat(npts0:npts0+npts-1) = pass_mlat cpass_mlt(npts0:npts0+npts-1) = pass_mlt cpass_number(npts0:npts0+npts-1) = pass_number npts0 = npts0+npts endfor ;stop ;;;;data binning ;Define bin size mlat_window = 1.0 ;2.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_equalbin2.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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount_peak(imlt,imlat) ge peak_counts) then begin peak_counts = bincount_peak(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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 ;stop loadct, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(i2)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount_peak(imlt,imlat) ge peak_counts) then begin peak_counts = bincount_peak(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(i2)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 xbin_mlat_array = fltarr(90000L)*!values.f_nan ;MLAT array with number of counts in each bin nn0 = 0. for imlat = 0, binsz-1 do begin mlat_bincount_peak(imlat) = total(bincount_peak(*,imlat)) nn = round(mlat_bincount_peak(imlat)) xbin_mlat_array(nn0:nn0+nn-1) = xbin_mlat(imlat) nn0 = nn0+nn endfor indexn = where(xbin_mlat_array 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. xyouts,55,1000,'mean='+string(mean(xbin_mlat_array(indexn),/nan),format='(f6.2)'),charsize = 1. xyouts,55,500,'median='+string(median(xbin_mlat_array(indexn)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount_peak(indexn),fp_max) ;xyouts,55,500,'max='+string(xbin_mlat(indexn(fp_max)),format='(f6.2)'),charsize = 1. indexs = where(xbin_mlat_array 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. xyouts,-70,1000,'mean='+string(mean(xbin_mlat_array(indexs),/nan),format='(f6.2)'),charsize = 1. xyouts,-70,500,'median='+string(median(xbin_mlat_array(indexs)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount_peak(indexs),fp_max) ;xyouts,-70,500,'max='+string(xbin_mlat(indexs(fp_max)),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 index1 = where(xbin_mlat ge 50.0) mlt_bincount_peak_nh(imlt) = total(bincount_peak(imlt,index1)) index2 = where(xbin_mlat le -50.0) mlt_bincount_peak_sh(imlt) = total(bincount_peak(imlt,index2)) 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. ;print, max(mlt_bincount_peak_nh,fp_max) ;xyouts,3,600,'max='+string(xbin_mlt(fp_max),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. ;print, max(mlt_bincount_peak_sh,fp_max) ;xyouts,12,600,'max='+string(xbin_mlt(fp_max),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) ;number of original data points in a bin if (count_bin ne 0) then begin pbin = cpass_number(fp_bin) bincount_pass(imlt,imlat) = n_elements(pbin(uniq(pbin,sort(pbin)))) 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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount_pass(imlt,imlat) ge peak_counts) then begin peak_counts = bincount_pass(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(i4)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount_pass(imlt,imlat) ge peak_counts) then begin peak_counts = bincount_pass(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(i4)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 xbin_mlat_array = fltarr(8000000L)*!values.f_nan ;MLAT array with number of counts in each bin nn0 = 0. for imlat = 0, binsz-1 do begin mlat_bincount_pass(imlat) = total(bincount_pass(*,imlat)) nn = round(mlat_bincount_pass(imlat)) xbin_mlat_array(nn0:nn0+nn-1) = xbin_mlat(imlat) nn0 = nn0+nn endfor indexn = where(xbin_mlat_array ge 50.0) plot, xbin_mlat, mlat_bincount_pass, xtitle = 'mlat(NH)',ytitle='counts',charsize = 1.5,xr=[50, 90],xstyle=1,yr=[0, 0.6e5] ;[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. xyouts,55,50000,'mean='+string(mean(xbin_mlat_array(indexn),/nan),format='(f6.2)'),charsize = 1. xyouts,55,20000,'median='+string(median(xbin_mlat_array(indexn)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount_pass(indexn),fp_max) ;xyouts,55,20000,'max='+string(xbin_mlat(indexn(fp_max)),format='(f6.2)'),charsize = 1. indexs = where(xbin_mlat_array le -50.0) plot, xbin_mlat, mlat_bincount_pass, xtitle = 'mlat(SH)',ytitle='counts',charsize = 1.5,xr=[-90,-50],xstyle=1,yr=[0, 0.6e5] ;[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. xyouts,-70,50000,'mean='+string(mean(xbin_mlat_array(indexs),/nan),format='(f6.2)'),charsize = 1. xyouts,-70,20000,'median='+string(median(xbin_mlat_array(indexs)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount_pass(indexs),fp_max) ;xyouts,-70,20000,'max='+string(xbin_mlat(indexs(fp_max)),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 index1 = where(xbin_mlat ge 50.0) mlt_bincount_pass_nh(imlt) = total(bincount_pass(imlt,index1)) index2 = where(xbin_mlat le -50.0) mlt_bincount_pass_sh(imlt) = total(bincount_pass(imlt,index2)) endfor plot, xbin_mlt, mlt_bincount_pass_nh, xtitle = 'mlt(NH)',ytitle='counts',charsize = 1.5,xr=[0, 24],xstyle=1, $ ystyle = 1, xticks = 4, yr=[5.2e4,5.8e4] ;[4.8e4,5.4e4] temp = total(xbin_mlt*mlt_bincount_pass_nh)/total(mlt_bincount_pass_nh) xyouts,3,5.65e4,'mean='+string(temp,format='(f6.2)'),charsize = 1. ;print, max(mlt_bincount_pass_nh,fp_max) ;xyouts,3,5.3e4,'max='+string(xbin_mlt(fp_max),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, $ ystyle = 1, xticks = 4, yr=[5.2e4,5.8e4] ;[4.8e4,5.4e4] temp = total(xbin_mlt*mlt_bincount_pass_sh)/total(mlt_bincount_pass_sh) xyouts,12,5.65e4,'mean='+string(temp,format='(f6.2)'),charsize = 1. ;print, max(mlt_bincount_pass_sh,fp_max) ;xyouts,12,5.3e4,'max='+string(xbin_mlt(fp_max),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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount(imlt,imlat) ge peak_counts) then begin peak_counts = bincount(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(f4.2)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 peak_counts = 0.0 & peak_mlat = 0.0 & peak_mlt = 0.0 loadct, 33 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 if (bincount(imlt,imlat) ge peak_counts) then begin peak_counts = bincount(imlt,imlat) peak_mlat = mlat_bin(imlat) peak_mlt = mlt_bin(imlt) endif 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, 10 ;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 =0, color = 200, thick = 6 endfor 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,charthick=4,color=200 endfor loadct, 39 oplot,[0.,0.],[-45.,45.] oplot,[-45.,45.],[0.,0.] 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. xyouts, -49, -44., 'peak counts:'+string(peak_counts,format='(f4.2)'),charsize=1.1 xyouts, -49, -48., 'peak MLAT:'+string(peak_mlat,format='(i3)'),charsize=1.1 xyouts, -49, -52, 'peak MLT:'+string(peak_mlt,format='(i2)'),charsize=1.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 xbin_array = fltarr(50000L)*!values.f_nan ;MLAT array with number of counts in each bin nn0 = 0 for imlat = 0, binsz-1 do begin mlat_bincount(imlat) = total(bincount(*,imlat)) nn = round(mlat_bincount(imlat)*1000.) xbin_array(nn0:nn0+nn-1L) = xbin(imlat) nn0 = nn0+nn endfor indexn = where(xbin_array 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.35,'mean='+string(temp,format='(f6.2)'),charsize = 1. xyouts,55,0.35,'mean='+string(mean(xbin_array(indexn),/nan),format='(f6.2)'),charsize = 1. xyouts,55,0.2,'median='+string(median(xbin_array(indexn)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount(indexn),fp_max) ;xyouts,55,0.2,'max='+string(xbin(indexn(fp_max)),format='(f6.2)'),charsize = 1. indexs = where(xbin_array 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.35,'mean='+string(temp,format='(f6.2)'),charsize = 1. xyouts,-70,0.35,'mean='+string(mean(xbin_array(indexs),/nan),format='(f6.2)'),charsize = 1. xyouts,-70,0.2,'median='+string(median(xbin_array(indexs)),format='(f6.2)'),charsize = 1. ;print, max(mlat_bincount(indexs),fp_max) ;xyouts,-70,0.2,'max='+string(xbin(indexs(fp_max)),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 index1 = where(xbin ge 50.0) mlt_bincount_nh(imlt) = total(bincount(imlt,index1)) index2 = where(xbin le -50.0) mlt_bincount_sh(imlt) = total(bincount(imlt,index2)) 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,0.4,'mean='+string(temp,format='(f6.2)'),charsize = 1. ;print, max(mlt_bincount_nh,fp_max) ;xyouts,3,0.2,'max='+string(xbin_mlt(fp_max),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,0.4,'mean='+string(temp,format='(f6.2)'),charsize = 1. ;print, max(mlt_bincount_sh,fp_max) ;xyouts,12,0.2,'max='+string(xbin_mlt(fp_max),format='(f6.2)'),charsize = 1. !p.multi = 0 ;;;MLAT vs. counts @ noon. temp = (bincount(11,*)+bincount(12,*))*0.5 ;;average of 11.5 and 12.5 MLT bins plot, xbin, temp,xr=[-90, 90], xstyle = 1, xticks = 12, $ xtitle ='MLAT', ytitle='Counts @ 12 MLT' print, max(temp(index1),fp_max1), max(temp(index2),fp_max2) xyouts, 30, 0.03, string(xbin(index1(fp_max1)),format='(f6.2)'),charsize = 1. xyouts, -60, 0.03, string(xbin(index2(fp_max2)),format='(f6.2)'),charsize = 1. device,/close set_plot,'x' ;;;;Write the counts @ noon into a file ;;;;-------------------- openw, 22,'CHAMP_360km_AACGM_peaks_30%_normalized@12MLT.dat',width = 200 ;printf, 22, 'bin_mlat counts' for imlat = 0, binsz-1 do begin printf, 22, xbin(imlat), temp(imlat) endfor close,22 ;;;;Write the normalized data into a file ;;;;-------------------- ; openw, 22,'CHAMP_360km_AACGM_peaks_30%_normalized_1MLAT.dat', width = 800 ; printf, 22, 'bin_mlat counts' ; for imlat = 0, binsz-1 do begin ; printf, 22, xbin(imlat), mlat_bincount(imlat) ; endfor ; close,22 end