PRO make_sed_phangs_tables,test=test,show_seds=show_seds

;make_sed_phangs_tables

dustem_define_la_common
model='DBP90'

;===== This is to do a small test_table
pd = [ $
     '(*!dustem_params).G0', $                           ;G0
     '(*!dustem_params).grains(0).mdust_o_mh',$          ;PAH0 mass fraction
     '(*!dustem_params).grains(1).mdust_o_mh']          ;VSG mass fraction

iv_min =   [0.1,1.e-4,1.e-4]
iv_max =   [100,1.e-1,1.e-1] 
log=[1,1,1]

fpd=[]
fiv=[]

dir=!dustem_wrap_soft_dir+'/Grids/'

IF keyword_set(test) THEN BEGIN
   iv_Nvalues=[5,3,3]
   table_name=dir+'TEST_DBP90_JWST_G0_YPAH_YVSG_4Phangs.fits'
ENDIF ELSE BEGIN
	iv_Nvalues=[50,20,20]
	table_name=dir+'DBP90_JWST_G0_YPAH_YVSG_4Phangs.fits'
ENDELSE

dustem_init,model=model
filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names]

dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,log=log,show_seds=show_seds

stop

END