First commit

rmaris
1 parent 30881cdf
Showing 1 changed file with 32 additions and 0 deletions Show diff stats
LabTools/IRAP/RM/test_size_distribution.pro
@@ -0,0 +1,32 @@
+PRO test_size_distribution
+
+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
+model='DL07'
+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
+
+
+params = ['(*!dustem_params).grains(0).amax','(*!dustem_params).grains(1).amax']
+fiv=[-1., 1,0.,0]     ;This sets the ISRF class to the requested value. also removes ionising photons.
+; use_fit_filters_names=[dustem_instru2filters('NIRCAM'),dustem_instru2filters('MIRI'),dustem_instru2filters('PACS'),dustem_instru2filters('SPIRE')] 
+
+use_fit_filters_names = dustem_filter_names2filters(['F300M','F335M','F360M','F0770W','F1000W','F1130W','F2100W','PACS_BLUE','PACS_GREEN','PACS_RED','PSW','PMW','PLW'])
+path_save = 'tmp/sed_size.fits'
+
+
+dustem_make_sed_table,model,params, [1.2000E-07,1.2000E-07],[3.300E-08,3.300E-08],[5,5],fpd=fpd,fiv=fiv,plog=[0,0],/show_seds,filters=use_fit_filters_names,filename=path_save, $
+                     grain_keywords=grain_keywords
+
+
+
+
+the_end:
+message,'test_size_distribution',/info$
+
+END
 \ No newline at end of file