fixed a bug

Jean-Philippe Bernard
1 parent 49c5df35
Showing 1 changed file with 7 additions and 7 deletions Show diff stats
src/idl/dustem_brute_force_fit.pro
@@ -126,13 +126,6 @@ ENDELSE
 chi2=min(chi2s)
 ind=where(chi2s EQ chi2,count)
  
-;==== should check bumping to table edges there
-params_hit=intarr(Nparams)
-FOR i=0L,Nparams-1 DO BEGIN
-  IF params[i] EQ minvalues[i] THEN params_hit[i]=-1
-  IF params[i] EQ maxvalues[i] THEN params_hit[i]=1
-ENDFOR
-
 params=grid_param_values[ind[0],*]  ;these are the best fit parameters
 fact=facts[ind[0]]   ; This is the corresponding scaling factor
 Nfreedom=Nseds-Nparams-1   ;This is the degree of freedom
@@ -140,6 +133,13 @@ rchi2=chi2/Nfreedom        ;This is the reduced chi2
 best_grid_sed=reform(grid_seds[ind[0],*])   ;This is the best SED fround in the grid.
 best_sed=best_grid_sed*fact  ;This is the best scaled grid SED matching the input SED.
  
+;==== should check bumping to table edges there
+params_hit=intarr(Nparams)
+FOR i=0L,Nparams-1 DO BEGIN
+  IF params[i] EQ minvalues[i] THEN params_hit[i]=-1
+  IF params[i] EQ maxvalues[i] THEN params_hit[i]=1
+ENDFOR
+
 ;==== should compute parameter uncertainties there
 conf_level=90./100.
 dchi2=delta_chi2(Nfreedom,conf_level)