chaaged the plot

Jean-Philippe Bernard
1 parent e62a03b3
Showing 1 changed file with 50 additions and 1 deletions Show diff stats
src/idl/dustem_brute_force_fit.pro
 FUNCTION dustem_brute_force_fit,sed,table_name,filters,normalize=normalize,fact=fact,chi2=chi2,rchi2=rchi2,show_sed=show_sed
  
+;+
+; NAME:
+;    dustem_brute_force_fit
+; PURPOSE:
+;    does brut force fit from a pre-computed model grid
+; CATEGORY:
+;    DustEM
+; CALLING SEQUENCE:
+;    res=dustem_brute_force_fit(sed,table_name,filters,normalize=normalize,fact=fact,chi2=chi2,rchi2=rchi2,show_sed=show_sed)
+; INPUTS:
+;    None
+; OPTIONAL INPUT PARAMETERS:
+;    key  = input parameter number
+;           First parameter: E(B-V) for extinction applied to the template, as used in Phangs
+;           Following 13*6 parameters are weights by which to multiply the MILES templates [Msun/pc^2]
+;    val  = corresponding input parameter value
+; OUTPUTS:
+;    stellar_brightness = stellar spectrum spectrum resampled on dustem wavelengths [MJy/sr]
+; OPTIONAL OUTPUT PARAMETERS:
+;    scope = scope of the plugin
+;    paramdefault = default values of parameters
+;    paramtag = plugin parameter names as strings
+; ACCEPTED KEY-WORDS:
+;    help                  = if set, print this help
+;    method                = SSP used. can be EMILES or CB19. default='EMILES'
+; COMMON BLOCKS:
+;    None
+; SIDE EFFECTS:
+;    None
+; RESTRICTIONS:
+;    The DustEMWrap IDL code must be installed
+; PROCEDURE:
+;    This is a DustEMWrap plugin for phangs
+;    It differs from dustem_plugin_emiles_stellar_continuum.pro only by the way the extinction is applied to the output spectrum
+; EXAMPLES
+;    dustem_init
+;    vec=dustem_plugin_phangs_stellar_continuum(scope=scope)
+; MODIFICATION HISTORY:
+;    Written by JPB June 2023
+;    Evolution details on the DustEMWrap gitlab.
+;    See http://dustemwrap.irap.omp.eu/ for FAQ and help.  
+;-
+
+IF keyword_set(help) THEN BEGIN
+  doc_library,'dustem_plugin_phangs_stellar_continuum'
+  output=0.
+  goto,the_end
+ENDIF
+
 ;stop
 ;sed must be normalized to NH=1e20 H/cm2
  
@@ -85,7 +134,7 @@ IF keyword_set(show_sed) THEN BEGIN
 	xtit='Wavelength [mic]'
 	ytit='SED (model=red,data=blue)'
 	;should also plot uncertainties
-	cgplot,sed.wave,sed.stokesI,psym='Filled Circle',color='blue',/ylog,yrange=[0.01,10],/ysty,xtit=xtit,ytit=ytit,/xlog
+	cgplot,sed.wave,sed.stokesI,psym='Filled Circle',color='blue',/ylog,/ysty,xtit=xtit,ytit=ytit,/xlog
 	cgoplot,sed.wave,grid_seds[ind[0],*]*fact,psym='Filled Square',color='red'
 ENDIF