pro st_dev, x, mean, sigma
;+
;
;  Name : st_dev
;  
;  Purpose :
;       Compute the standard deviation
;
;  Syntax : st_dev x, mean, num, sigma
;
;  Input :
;      x : data
;      mean : mean value
;
;  Output :
;      sigma : standard deviation
;
;  Keywords :
;
;  Required routines :
;
;  History
;       2010 Oct, Dong-uk Song
;
;-

;Input value ********************************
;x=[2,2,2,2,2,2,2,2,2]
;mean=mean(x)
;********************************************

num = n_elements(x)
mx=fltarr(num)
for i = 0, num-1 do begin
mx[i]=(x[i]-mean)^2
sumx=total(mx)
xx=sumx/num
sigma=sqrt(xx)                        ; the value of standard deviation
endfor

end