pro find_mag,value,mag
;=====================
; Finds what order of 
; magnitude a value is
;=====================

r=100.
mag =-15.
while (r ge 10.) do begin
	mag = mag + 1.
	r = abs(value)/(10^mag)
      ;  print,mag,r
endwhile
;print,'*',mag,r

return
end