OV-GSO-DC / dustem-wrapper_idl

Commit 42b6dd578aa10919d1adb25d45f2c010d537df93

Authored by Jean-Philippe Bernard
1 parent 0c48096d
Exists in master

first commit

Showing 1 changed file with 40 additions and 0 deletions   Show diff stats
Inline Side-by-side
LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro 0 → 100644
  Diff comments   View file @42b6dd5
... ... @@ -0,0 +1,40 @@
  1 +
  2 +PRO phangs_compare_seds_isrf
  3 +
  4 +dir=!phangs_data_dir+'/ISRF/GRIDS/'
  5 +isrf_class_str='_isrfclass15'
  6 +model='DBP90'
  7 +table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs'+isrf_class_str+'.fits'
  8 +st=mrdfits(table_name,1,h)
  9 +
  10 +filters=['NIRCAM1','MIRI2']
  11 +Nfilter=n_elements(filters)
  12 +
  13 +ii=100L
  14 +sed=fltarr(Nfilters)
  15 +tags=tagnames(st)
  16 +FOR i=0L,Nfilters-1 DO BEGIN
  17 + ind=where(tags EQ 'I'+filters[i],count)
  18 + sed[i]=st[ii].(ind[0])
  19 +ENDFOR
  20 +Nparams=3
  21 +true_params=fltarr(Nparams)
  22 +FOR i=0L,Nparams-1 DO BEGIN
  23 + true_params[i]=st[ii].(i)
  24 +ENDFOR
  25 +
  26 +
  27 +isrf_class_str='_isrfclass0'
  28 +table_name_2=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs'+isrf_class_str+'.fits'
  29 +
  30 +params=dustem_brute_force_fit(sed,table_name_2,filters,fact=fact,chi2=chi2,/normalize,rchi2=rchi2,show_sed=show_sed $
  31 + ,params_hit=params_hit,params_uncertainties=params_uncertainties,params_min=params_min,params_max=params_max $
  32 + ,fixed_parameters_description=fixed_parameters_description,fixed_parameters_values=fixed_parameters_values)
  33 +
  34 +print,true_params
  35 +print,params
  36 +
  37 +stop
  38 +
  39 +END
  40 +
... ...