chaaged the plot

Jean-Philippe Bernard
1 parent b4dbd43c
Showing 1 changed file with 6 additions and 3 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,nostop=nostop
+FUNCTION dustem_brute_force_fit,sed,table_name,filters,normalize=normalize,fact=fact,chi2=chi2,rchi2=rchi2,best_sed=best_sed,show_sed=show_sed,nostop=nostop
  
 ;+
 ; NAME:
@@ -129,13 +129,16 @@ params=grid_param_values[ind[0],*]
 fact=facts[ind[0]]
 Nfreedom=Nseds-Nparams-1
 rchi2=chi2/Nfreedom
+best_sed=grid_seds[ind[0],*]*fact
  
 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,/ysty,xtit=xtit,ytit=ytit,/xlog
-	cgoplot,sed.wave,grid_seds[ind[0],*]*fact,psym='Filled Square',color='red'
+	minv=la_min([sed.stokesI,best_sed])
+	maxv=la_max([sed.stokesI,best_sed])
+	cgplot,sed.wave,sed.stokesI,psym='Filled Circle',color='blue',/ylog,/ysty,xtit=xtit,ytit=ytit,/xlog,yrange=[minv,maxv]
+	cgoplot,sed.wave,best_sed,psym='Filled Square',color='red'
 	if not keyword_set(nostop) THEN stop
 ENDIF