PRO make_sed_phangs_tables,test=test,show_seds=show_seds,grid_type=grid_type,isrf_class=isrf_class

;make_sed_phangs_tables,/show   ;That's a test

;make_sed_phangs_tables,grid_type=1,/show_seds,/test   ;That's a test

;make_sed_phangs_tables,/test,grid_type=2,/show_seds,isrf_class=15

;make_sed_phangs_tables,/test,grid_type=3,/show_seds

;=== This is needed for Phangs work

;defsysv,'!phangs_data_dir','/Volumes/PILOT_FLIGHT1/PHANGS/'
;defsysv,'!phangs_data_dir','/data/projects/phangs/'

;=== This is where the grids will be written

;dir=!dustem_wrap_soft_dir+'/Grids/'

dir=!phangs_data_dir+'/ISRF/GRIDS/'

dustem_define_la_common

use_grid_type=0

IF keyword_set(grid_type) THEN BEGIN
  use_grid_type=grid_type
ENDIF

test_str=''
IF keyword_set(test) THEN BEGIN
  test_str='TEST_'
ENDIF

CASE use_grid_type OF

   0: BEGIN ;======= This is for test only
      !quiet=1
      model='DBP90'
      ;===== This is to do a small test_table
      pd = [ $
           '(*!dustem_params).G0' $  ;G0
            ]
      iv_min =   [0.1]
      iv_max =   [10]
      iv_Nvalues=[2]

      plog=[1]

      fpd=['(*!dustem_params).grains(0).mdust_o_mh',$          ;PAH0 mass fraction
           '(*!dustem_params).grains(1).mdust_o_mh']
      fiv=[1.e-3,1.e-2]
      table_name=dir+'TEST_'+model+'_G0.fits'
      dustem_init,model=model

      ;=== select filters to be used for the grid (default=ALL)
      ;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,plog=plog,show_seds=show_seds,/print_params,filters=filters

   END

   1: BEGIN ;======= This is to produce a grid with the DBP90 model + Mathis field for G0 and PAH abundance

      !quiet=1
      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] 

      plog=[1,1,1]

      fpd=[]
      fiv=[]

      IF keyword_set(test) THEN BEGIN

         ;iv_Nvalues=[5,3,3]
         iv_Nvalues=[2,2,2]

         table_name=dir+'TEST_'+model+'_JWST_G0_YPAH_YVSG_4Phangs.fits'

      ENDIF ELSE BEGIN
      	iv_Nvalues=[50,20,20]

      	table_name=dir+model+'_JWST_G0_YPAH_YVSG_4Phangs.fits'

      ENDELSE

      dustem_init,model=model
      ;=== select filters to be used for the grid
      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,plog=plog,show_seds=show_seds

   END

   2: BEGIN          ;This is DPB90 model with PHANGS ISRF classes

      !quiet=1

      model='DBP90'   ;This is to use the DBP90 dust model

      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2

      pd = [ $
           'dustem_plugin_phangs_class_isrf_2', $             ;G0 factor for Phangs ISRF classes
           '(*!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] 

      plog=[1,1,1]

      Nclass=31L

      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4']             ;ISRF class to be used

      dustem_init,model=model

      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number

      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too

      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)

      filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names, $

               (*!dustem_filters).iras.filter_names,(*!dustem_filters).pacs.filter_names,(*!dustem_filters).spire.filter_names]

      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table

         use_isrf_class=isrf_class

      ENDIF ELSE BEGIN

         use_isrf_class=0  ;This is for Mathis field

      ENDELSE

      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. alos removes ionising photons.

      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN

         ;iv_Nvalues=[2,2,2]
         iv_Nvalues=[4,4,6,10]

         table_name=dir+'TEST_'+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs_noionis'+isrf_class_str+'.fits'

      ENDIF ELSE BEGIN
         iv_Nvalues=[50,20,20]

         table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs_noionis'+isrf_class_str+'.fits'

      ENDELSE

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

   END

   3: BEGIN          ;This is DPB90 model with DL07 ISRF prescription (a la Chastanet)
      !quiet=1
      model='DBP90'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $

           'dustem_plugin_dl07_isrf_model_3', $             ;Umin factor [0.1-100] 6 values
           'dustem_plugin_dl07_isrf_model_1', $             ;gamma factor [0.0,0.02] 5 values
           '(*!dustem_params).grains(0).mdust_o_mh',$       ;PAH0 mass fraction [1.e-4,1.e-1] 20 values
           '(*!dustem_params).grains(1).mdust_o_mh' $       ;VSG mass fraction [1.e-4,1.e-1] 20 values

            ]

      ;iv_min =   [0.1,0.  ,1.e-4,1.e-4]      ;does not run for values =0., because of use of pmin (division by initial value)
      iv_min =   [0.1,1.e-10  ,1.e-4,1.e-4]

      iv_max =   [100.,0.02,1.e-1,1.e-1]

      plog=[1,0,1,1]

      fpd=['dustem_plugin_dl07_isrf_model_2', $ ;alpha=1

           'dustem_plugin_dl07_isrf_model_4' $  ;Umax=1.e7
           ]

      fiv=[2.,1.e7]     ;This sets the fixed parameters

      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number

      !dustem_verbose=1
      !quiet=0
      !dustem_nocatch=1
      ;!dustem_verbose=0
      ;(*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too

      ;=== select filters to be used for the grid (no definition -> all known filters used)
      ;filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names, $
      ;         (*!dustem_filters).iras.filter_names,(*!dustem_filters).pacs.filter_names,(*!dustem_filters).spire.filter_names]

      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN

         ;iv_Nvalues=[2,2,3,3]

         iv_Nvalues=[2,2,3,3]
         ;iv_Nvalues=[4,4,6,10]

         table_name=dir+'TEST_'+model+'_DL07ISRF_YPAH_YVSG_Umin_gamma_4Phangs'+'.fits'

      ENDIF ELSE BEGIN
         iv_Nvalues=[6,5,20,20]

         table_name=dir+model+'_DL07ISRF_YPAH_YVSG_Umin_gamma_4Phangs'+'.fits'

      ENDELSE

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

   END

   4: BEGIN          ;This is DL07 model with PHANGS ISRF classes
      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2', $             ;G0 factor for Phangs ISRF classes
           '(*!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] 
      iv_min =   [0.1,5.e-5,5.e-5]
      iv_max =   [1.e6,5.e-1,5.e-1]
      plog=[1,1,1]
      Nclass=31L
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4']             ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)
      filters=dustem_get_all_filter_names()
      ;filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names, $
      ;         (*!dustem_filters).iras.filter_names,(*!dustem_filters).pacs.filter_names,(*!dustem_filters).spire.filter_names]
      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. alos removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         ;iv_Nvalues=[2,2,2]
         iv_Nvalues=[2,2,3,5]
         table_name=dir+'TEST_'+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs_noionis'+isrf_class_str+'.fits'
      ENDIF ELSE BEGIN

         iv_Nvalues=[100,20,20]

         table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs_noionis'+isrf_class_str+'.fits'
      ENDELSE
      dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,plog=plog,show_seds=show_seds
   END

   5: BEGIN          ;This is DL07 model with PHANGS ISRF classes, with only G0 as free parameter. This takes way too long

      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2' $             ;G0 factor for Phangs ISRF classes
           ]          
      iv_min =   [0.1]
      iv_max =   [1.e6]
      plog=[1]

      IF keyword_set(test) THEN BEGIN
         iv_min =   [1.]
         iv_max =   [3.]
         plog=[1]
      ENDIF

      ;Nclass=31L
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4'] ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)

      ;filters=dustem_get_all_filter_names()
      filters=get_phangs_filter_names(/order_by_wavelengths)

      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. also removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         iv_Nvalues=[2]

         ;table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'

      ENDIF ELSE BEGIN
         iv_Nvalues=[100]

      ENDELSE

      table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'

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

      ;goto,one_by_one

      ;adding linear parameters abundance
      message,'Adding Dust abundances to grid',/continue
      pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[1].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))
      ;iv_min=[default_value/10.] & iv_max=[default_value*100.]
      iv_min=[default_value/10.]

      iv_max=[default_value*10.]
      iv_Nvalues=[20,20,20]

      plog=[1,1,1]
      IF keyword_set(test) THEN BEGIN
         iv_min=[default_value/2.] & iv_max=[default_value*2.]
         iv_Nvalues=[3,3,3]
         plog=[0,0,0]
      ENDIF
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypah1added_Ypah2added_Yvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'

      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds

      goto,the_end

      one_by_one:

      ;adding linear parameters abundance

      message,'Adding PAH0 abundances to grid',/continue
      pd=['(*!dustem_params).grains[0].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))[0]
      iv_min=[default_value/10.] & iv_max=[default_value*100.] & iv_Nvalues=[100] & plog=[1]
      IF keyword_set(test) THEN iv_Nvalues=[3]
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypahadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds

      ;adding linear parameter PAH1 abundance

      message,'Adding PAH1 abundances to grid',/continue
      table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypahadded_4Phangs_noionis'+isrf_class_str+'.fits'
      pd=['(*!dustem_params).grains[1].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))[0]
      iv_min=[default_value/10.] & iv_max=[default_value*100.] & iv_Nvalues=[100] & plog=[1]
      IF keyword_set(test) THEN iv_Nvalues=[3]
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypah1added_Ypah2added_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds

      ;adding linear parameter VSG abundance

      message,'Adding VSG abundances to grid',/continue
      table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypah1added_Ypah2added_4Phangs_noionis'+isrf_class_str+'.fits'
      pd=['(*!dustem_params).grains[3].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))[0]
      iv_min=[default_value/10.] & iv_max=[default_value*100.] & iv_Nvalues=[100] & plog=[1]
      IF keyword_set(test) THEN iv_Nvalues=[3]
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypah1added_Ypah2added_Yvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds

   END

   22: BEGIN          ;This is DL07 model with PHANGS ISRF classes, with only G0 as free parameter. Same as 2, but for DL07 model
      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2' $             ;G0 factor for Phangs ISRF classes
           ]          
      iv_min =   [0.1]
      iv_max =   [100.]
      plog=[1]
      IF keyword_set(test) THEN BEGIN
         iv_min =   [.1]
         iv_max =   [100.]
         plog=[1]
      ENDIF
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4'] ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)
      ;filters=dustem_get_all_filter_names()
      filters=get_phangs_filter_names(/order_by_wavelengths)
      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. also removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         iv_Nvalues=[4]
         ;table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'
      ENDIF ELSE BEGIN

         ;iv_Nvalues=[19]
         iv_Nvalues=[22]

      ENDELSE
      table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,plog=plog,show_seds=show_seds
      ;adding linear parameters (small grains abundances)
      message,'Adding Dust abundances to grid',/continue
      pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[1].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))
      iv_min=[default_value/10.]
      iv_max=[default_value*10.]

      ;iv_Nvalues=[11,11,11]
      iv_Nvalues=[13,13,13]

      plog=[1,1,1]
      IF keyword_set(test) THEN BEGIN
         iv_min=[default_value/2.]
         iv_max=[default_value*2.]
         iv_Nvalues=[3,3,3]
         plog=[1,1,1]
      ENDIF
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_G0_Ypah1added_Ypah2added_Yvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds
      goto,the_end
   END

   221: BEGIN          ;This is same as 22, but adapted for observed range in NGC0628
      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2' $             ;G0 factor for Phangs ISRF classes
           ]          
      iv_min =   [0.1]
      iv_max =   [10.]
      plog=[1]
      IF keyword_set(test) THEN BEGIN
         iv_min =   [.1]
         iv_max =   [10.]
         plog=[1]
      ENDIF
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4'] ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)
      ;filters=dustem_get_all_filter_names()
      filters=get_phangs_filter_names(/order_by_wavelengths)
      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. also removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         iv_Nvalues=[4]
         ;table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'
      ENDIF ELSE BEGIN
         ;iv_Nvalues=[19]
         iv_Nvalues=[12]
      ENDELSE
      table_name=dir+test_str+model+'_MuseISRF_JWST_adaptedG0_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,plog=plog,show_seds=show_seds
      ;adding linear parameters (small grains abundances)
      message,'Adding Dust abundances to grid',/continue
      pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[1].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))
      iv_min=[default_value/10.]
      iv_max=[default_value[0:1]*10.,default_value[2]*100.]
      ;iv_Nvalues=[11,11,11]
      iv_Nvalues=[13,13,13]
      plog=[1,1,1]
      IF keyword_set(test) THEN BEGIN
         iv_min=[default_value/2.]
         iv_max=[default_value*2.]
         iv_Nvalues=[3,3,3]
         plog=[1,1,1]
      ENDIF
      output_table_name=dir+test_str+model+'_MuseISRF_JWST_adaptedG0_Ypah1added_Ypah2added_adaptedYvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds
      goto,the_end
   END

   222: BEGIN          ;This is same as 221, but using the standard DL07 GRAIN.DAT, ie, Ypah+ set by Fortran and PAH keywords set to mix-logn
                        ;Note: This does not really work, because PAH+ contribution does not get multiplied linearly by abundance ...
      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2' $             ;G0 factor for Phangs ISRF classes
           ]          
      iv_min =   [0.1]
      iv_max =   [10.]
      plog=[1]
      IF keyword_set(test) THEN BEGIN
         iv_min =   [.1]
         iv_max =   [10.]
         plog=[1]
      ENDIF
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4'] ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;stop
      (*!dustem_params).grains[0].TYPE_KEYWORDS='logn-mix'
      (*!dustem_params).grains[1].TYPE_KEYWORDS='logn-mix'
      ;(*!dustem_params).grains[0].TYPE_KEYWORDS='logn'
      ;(*!dustem_params).grains[1].TYPE_KEYWORDS='logn'
      grain_keywords=[(*!dustem_params).grains.TYPE_KEYWORDS]   ;This is to inform dustem_make_sed_table, which actually runs the models
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)
      ;filters=dustem_get_all_filter_names()
      filters=get_phangs_filter_names(/order_by_wavelengths)
      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. also removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         iv_Nvalues=[4]
         ;table_name=dir+test_str+model+'_MuseISRF_JWST_G0_4Phangs_noionis'+isrf_class_str+'.fits'
      ENDIF ELSE BEGIN
         ;iv_Nvalues=[19]
         iv_Nvalues=[12]
      ENDELSE
      table_name=dir+test_str+model+'_orig_MuseISRF_JWST_adaptedG0_4Phangs_noionis'+isrf_class_str+'.fits'
      ;stop
      dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,plog=plog,show_seds=show_seds, $
                           grain_keywords=grain_keywords
      ;adding linear parameters (small grains abundances)
      message,'Adding Dust abundances to grid',/continue
      ;pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[1].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))
      iv_min=[default_value/10.]
      ;iv_max=[default_value[0:1]*10.,default_value[2]*100.]
      iv_max=[default_value[0]*10.,default_value[1]*100.]
      ;iv_Nvalues=[11,11,11]
      ;iv_Nvalues=[13,13,13]
      iv_Nvalues=[13,13]
      ;plog=[1,1,1]
      plog=[1,1]
      IF keyword_set(test) THEN BEGIN
         iv_min=[default_value/2.]
         iv_max=[default_value*2.]
         ;iv_Nvalues=[3,3,3]
         ;plog=[1,1,1]
         iv_Nvalues=[3,3]
         plog=[1,1]
      ENDIF
      output_table_name=dir+test_str+model+'_orig_MuseISRF_JWST_adaptedG0_Ypah1added_Yvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds
      goto,the_end
   END
   223: BEGIN          ;This is same as 221, but using the standard DL07 GRAIN.DAT, ie, YPAH0 and Ypah1 keywords set to logn
                        ;Note: This assumes a variable Xion which is different from DL07 and no size depednence of Xion
      !quiet=1
      model='DL07'   ;This is to use the DBP90 dust model
      ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
      ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
      pd = [ $
           'dustem_plugin_phangs_class_isrf_2' $             ;G0 factor for Phangs ISRF classes
           ]          
      iv_min =   [0.1]
      iv_max =   [10.]
      plog=[1]
      IF keyword_set(test) THEN BEGIN
         iv_min =   [.1]
         iv_max =   [10.]
         plog=[1]
      ENDIF
      fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1','dustem_plugin_phangs_class_isrf_4'] ;ISRF class to be used
      dustem_init,model=model
      fortran_user=dustem_set_up_fortran(/random_name)   ;use a random fortran number
      !dustem_verbose=0
      (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS    ;This makes Fortran be quiet too
      ;stop
      (*!dustem_params).grains[0].TYPE_KEYWORDS='logn'
      (*!dustem_params).grains[1].TYPE_KEYWORDS='logn'
      grain_keywords=[(*!dustem_params).grains.TYPE_KEYWORDS]   ;This is to inform dustem_make_sed_table, which actually runs the models
      ;=== select filters to be used for the grid (HOOPS, we need MUSE also, and maybe all filters ...)
      ;filters=dustem_get_all_filter_names()
      filters=get_phangs_filter_names(/order_by_wavelengths)
      IF keyword_set(isrf_class) THEN BEGIN  ;do just one table
         use_isrf_class=isrf_class
      ENDIF ELSE BEGIN
         use_isrf_class=0  ;This is for Mathis field
      ENDELSE
      fiv=[-1.,1.,use_isrf_class,0.]     ;This sets the ISRF class to the requested value. also removes ionising photons.
      isrf_class_str='_isrfclass'+strtrim(use_isrf_class,2)
      ;define the number of free parameters in the grid and the grid fits table name
      IF keyword_set(test) THEN BEGIN
         iv_Nvalues=[4]
      ENDIF ELSE BEGIN
         iv_Nvalues=[12]
      ENDELSE
      table_name=dir+test_str+model+'_logn_MuseISRF_JWST_adaptedG0_4Phangs_noionis'+isrf_class_str+'.fits'
      ;stop
      dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,plog=plog,show_seds=show_seds, $
                           grain_keywords=grain_keywords
      ;adding linear parameters (small grains abundances)
      message,'Adding Dust abundances to grid',/continue
      pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[1].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      ;pd=['(*!dustem_params).grains[0].MDUST_O_MH','(*!dustem_params).grains[2].MDUST_O_MH']
      default_value=(dustem_get_param_values_default(pd))
      iv_min=[default_value/10.]
      iv_max=[default_value[0:1]*10.,default_value[2]*100.]
      iv_Nvalues=[13,13,13]
      plog=[1,1,1]
      IF keyword_set(test) THEN BEGIN
         iv_min=[default_value/2.]
         iv_max=[default_value*2.]
         iv_Nvalues=[3,3,3]
         plog=[1,1,1]
      ENDIF
      output_table_name=dir+test_str+model+'_logn_MuseISRF_JWST_adaptedG0_Ypah1added_Yvsgadded_4Phangs_noionis'+isrf_class_str+'.fits'
      dustem_add_linear_params2grid,table_name,pd,iv_min,iv_max,iv_Nvalues,out_filename=output_table_name,plog=plog,show_seds=show_seds
      goto,the_end
   END

ENDCASE

the_end:

END