if !d.name eq 'PS' then begin
  device,filename='idl.ps',xsize=18,ysize=14,yoffset=3,/color
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
!p.charsize=2
!x.margin=[10,3.5]
!y.margin=[3.2,.5]
!p.multi=[0,1,1]
;Define parameters
nbins=50 ;choose number of bins for histogram
onethird=1.0/3.0
microm=1e-6
;
;----Obtain the number density of the particles
pc_read_pdim,obj=pdim
npar=pdim.npar
pc_read_dim,obj=dim
nx=dim.nx
ny=dim.ny
nz=dim.nz
pc_read_param,obj=start
np_swarm0=start.np_swarm0
n0=npar*np_swarm0/(nx*ny*nz)
print,'n0=',n0
;calculate "norma" from the initial profile
restore,'../SW512condens0_coag_grav0_turb_shima_f2em2_ngp_a24_L025/pvar_file/pvar0.sav' 
pdfi_n,a,npswarm,xx,yy,n=nbins-1,/log
norma=total(yy)/n0
;number density------------------------------------------------------------
restore,'./pvar_file/pvar5.sav'
pdfi_n,a,npswarm,xx,yy,n=nbins-1,/log
radii_norm5=xx/microm
prob_norm5=yy/norma
save,file='./data/radii_norm5_c.sav',radii_norm5
save,file='./data/prob_norm5_c.sav',prob_norm5
END