J13 grid

Georgina Sampson-olalde
1 parent 0ea960be
Showing 2 changed files with 76 additions and 2 deletions Show diff stats
LabTools/IRAP/GSO/make_grids_gso.pro
LabTools/IRAP/GSO/make_sed_tables_gso.pro
@@ -11,8 +11,8 @@ pro make_grids_gso
   t0 = systime(/sec)
  
   ; make_sed_tables_gso, grid_type = 'DBP90', /extend_trick
-  make_sed_tables_gso, grid_type = 'MC10', /extend_trick
-  ; make_sed_tables_gso, grid_type = 'DBP90', /extend_trick
+  ; make_sed_tables_gso, grid_type = 'MC10', /extend_trick
+  make_sed_tables_gso, grid_type = 'J13', /extend_trick
  
   !quiet = 0
  
@@ -216,6 +216,80 @@ pro make_sed_tables_gso, test = test, show_seds = show_seds, grid_type = grid_ty
         goto,the_end
       endelse
     end
+
+    'J13': begin ; ======= This is to produce a grid with the J13 model
+      message, 'Case J13', /info
+
+      !quiet = 1
+      model = 'J13'
+
+      ; Silence 
+      ; 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
+
+      dustem_init, model = model
+      if keyword_set(extend_trick) then begin
+        pd = [ $
+        '(*!dustem_params).G0' $ ; G0
+        ]
+        default_value=(dustem_get_param_values_default(pd))
+        iv_min=[default_value/10.] & iv_max=[default_value*10.]
+        plog = [1]
+
+        if keyword_set(test) then begin
+          iv_Nvalues = [3]
+        endif else begin
+          iv_Nvalues = [10]
+        endelse
+
+        fpd=['(*!dustem_params).grains(0).amin']
+        fiv=[4.0000E-08]
+        table_name = dir + test_str + model + '_G0.fits'
+      endif else begin
+        pd = [ $
+          '(*!dustem_params).G0', $ ; G0
+          '(*!dustem_params).grains(0).mdust_o_mh', $ ; CM20 mass fraction
+          '(*!dustem_params).grains(1).mdust_o_mh', $ ; CM20 mass fraction
+          '(*!dustem_params).grains(2).mdust_o_mh', $ ; aPyM5 mass fraction
+          '(*!dustem_params).grains(3).mdust_o_mh'] ; aOlM5 mass fraction
+        iv_min = [1.e-1, 1.70E-04, 6.30E-05, 2.55E-04, 2.55E-04]
+        iv_max = [1.e1, 1.70E-02, 6.30E-03, 2.55E-02, 2.55E-02]
+        plog = [1,1,1,1,1]
+
+        if keyword_set(test) then begin
+          iv_Nvalues = [3,3,3,3,3]
+        endif else begin
+          iv_Nvalues = [10,10,10,10,10]
+        endelse
+
+        fpd = [] & fiv = []
+        table_name = dir + test_str + model + '_G0_YCM20_YCM20_YaPyM5_YaOlM5.fits'
+      endelse
+
+
+      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
+
+
+      if keyword_set(extend_trick) then begin
+        ;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','(*!dustem_params).grains[3].MDUST_O_MH']
+        default_value=(dustem_get_param_values_default(pd))
+        iv_min=[default_value/10.] & iv_max=[default_value*10.]
+        iv_Nvalues=[10,10,10,10]
+        plog=[1,1,1,1]
+
+        IF keyword_set(test) THEN BEGIN
+          iv_Nvalues=[3,3,3,3]
+        ENDIF
+
+        output_table_name= dir + test_str + model + '_G0_YCM20added_YCM20added_YaPyM5added_YaOlM5added.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
+      endif else begin
+        goto,the_end
+      endelse
+    end
   endcase
  
   the_end: