function mean, data, missing if n_elements(missing) eq 0 then missing = 2.0*max(data) n_e = n_elements(data) if n_e gt 0 then begin n = 0 mean = 0.0 for i=0,n_e-1 do if data(i) ne missing then begin n = n+1 mean = mean + data(i) endif if n gt 0 then mean=mean/float(n) endif else mean = 1.0e+32 return, mean end