Some upper and lower limit santiy checking added

Annie Hughes
1 parent a8bf7c1c
Showing 1 changed file with 31 additions and 0 deletions Show diff stats
src/idl/dustem_init_params.pro
@@ -93,6 +93,16 @@ if keyword_set(ulims) then begin
       message, 'Number of upper limits must equal number of parameters (use ulimed to specify ON/OFF for each parameter)',/info
       stop
    end
+   ; compare ulims and ivs
+   on=where(ulimed eq 1,ct)
+   if ct gt 0 then begin
+      test=total(ulims[on] lt iv[on])
+      if test ne 0 then begin
+         message, 'Upper limits are less than requested initial values',/info
+         stop
+      end
+      
+   end
 end
  
 if keyword_set(llimed) then begin
@@ -109,6 +119,27 @@ if keyword_set(llims) then begin
       message, 'Number of lower limits must equal number of parameters (use llimed to specify ON/OFF for each parameter)',/info
       stop
    end
+   ; compare llims and ivs
+   on=where(llimed eq 1,ct)
+   if ct gt 0 then begin
+      test=total(llims[on] gt iv[on])
+      if test ne 0 then begin
+         message, 'Lower limits are greater than requested initial values',/info
+         stop
+      end
+      
+   end
+end
+
+if keyword_set(llims) and keyword_set(ulims) then begin
+   on=where(llimed eq 1 and ulimed eq 1,ct)
+   if ct gt 0 then begin
+      test=total(llims[on] gt ulims[on])
+      if test ne 0 then begin
+         message, 'Requested lower limits are greater than upper limits',/info
+         stop
+      end
+   end
 end