PRO test_dustem_dl07_isrf_model

use_model='DL07'
dustem_init,model=use_model,polarization=use_polarization
!quiet=1

;set up some data (needed as will return an SED)
dir=!dustem_wrap_soft_dir+'/Data/EXAMPLE_OBSDATA/'
file=dir+'SED_DraineModel_U1.0.xcat'
spec=read_xcat(file,/silent)
dustem_set_data,m_fit=spec,m_show=spec

;key=[1]
;val=[0.2]
key=0 & val=0   ;default values only
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0,Umean=Umean)

yr=[1.e-26,5.e-24]
;yr=[1.e-13,1.e-09]
xr=[2,1.e3]
nu=3.e14/dustem_get_wavelengths()

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

xtit='Wav [mic]'
ytit='nuInu/<U> [whatever]'

;==== Umin=1

val=[0.0,2.,1.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0,Umean=Umean0)
cgplot,dustem_get_wavelengths(),dl_spectrum_0/Umean,/xlog,/ylog,yr=yr,/ysty,xr=xr,/xsty,title='Umin=1',xtit=xtit,ytit=ytit

;gamma,alpha,Umin,Umax
key=[1,2,3,4]
val=[0.02,2.,1.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p02,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p02/Umean,color='blue',linestyle=4

val=[0.01,2.,1.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p01,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p01/Umean,color='green',linestyle=4

val=[0.005,2.,1.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p005,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p005/Umean,color='red',linestyle=4

cgoplot,dustem_get_wavelengths(),dl_spectrum_0/Umean0,color='black'

;stop

;==== Umin=10
val=[0.0,2.,10.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0,Umean=Umean0)
cgplot,dustem_get_wavelengths(),dl_spectrum_0/Umean0,/xlog,/ylog,yr=yr,/ysty,xr=xr,/xsty,title='Umin=10',xtit=xtit


val=[0.04,2.,10.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p04,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p04/Umean,color='violet',linestyle=4

val=[0.02,2.,10.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p02,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p02/Umean,color='blue',linestyle=4

val=[0.01,2.,10.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p01,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p01/Umean,color='green',linestyle=4

val=[0.005,2.,10.,1.e6]
sed=dustem_plugin_dl07_isrf_model(key=key,val=val,spec=dl_spectrum_0p005,Umean=Umean)
cgoplot,dustem_get_wavelengths(),dl_spectrum_0p005/Umean,color='red',linestyle=4

cgoplot,dustem_get_wavelengths(),dl_spectrum_0/Umean0,color='black'

stop


END