PRO plot_phangs_ssps_against_uv_fast

source_name='ngc0628'
data_dir='/Volumes/PILOT_FLIGHT1/PHANGS-JWST/DR1/'

restore,data_dir+'ngc0628_astrosat_minus_prediction_fast.sav',/verb
restore,data_dir+'ngc0628_astrosat_voronoi_prediction_fast.sav',/verb


;stop

window,0,xsize=800,ysize=800
!p.multi=[0,1,2]

xtit='(data-model) [MJy/sr]'
IF do_normalize THEN xtit='(data-model) normalized to 1e21 H/cm2 [MJy/sr]'
NHthreshold=0. & tit=''
;NHthreshold=0.1 & tit='NH > 1e20 H/cm2'
;NHthreshold=0.5 & tit='NH > 5e20 H/cm2'
;NHthreshold=1. & tit='NH > 1e21 H/cm2'
;NHthreshold=5. & tit='NH > 5e21 H/cm2'
;NHthreshold=10. & tit='NH > 1e22 H/cm2'
;NHthreshold=50. & tit='NH > 5e22 H/cm2'
;vmin=-0.1 & vmax=0.1 & Nbins=50
vmin=-1. & vmax=1. & Nbins=50
;vmin=-5. & vmax=5. & Nbins=100
;vmin=-10. & vmax=10. & Nbins=100
;im=diff_maps[*,*,0]
ims=actual_seds-predicted_seds

im=ims[*,k]
ind=where(im NE la_undef())
xtit='(data-model) normalized to 1e21 H/cm2 [MJy/sr]',
res=histogram(im[ind],locations=x,min=vmin,max=vmax,Nbins=Nbins)
yrange=[max([min(res),1.]),1.e4]
colors=['violet','blue','yellow','brown','red']
cgplot,x,res,/ylog,xtit=xtit,ytit='N',psym=10,yrange=yrange,/ysty,color='violet',tit=tit,/no_data
FOR k=0L,4 DO BEGIN
	im=ims[*,k]
	ind=where(im NE la_undef())
	res=histogram(im[ind],locations=x,min=vmin,max=vmax,Nbins=Nbins)
   cgoplot,x1,res1,color=colors[k],psym=10
ENDFOR

stop

END