hi everyone,

          i am using the following ascii file conversion script,

################################################################

files=findfile('*.hdf')
n = n_elements(files)
print,'No. of files to process:',n

for i=0,n-1 do begin            & $
   print,files(i)                     & $
   ofile=files(i)  + '.txt'          & $
   load,files(i),ftype='mapped',prod_name=['new band']     & $
   display, band_no=i+1,fbuf=i+1                                               & $
      out_ascii,ofile,fbuf=i+1,vars=[3,4,2], format='(F20.10,F20.10,F20.10)', region=[1917,1377,2271,1694]  & $
endfor
clear_up

###########################################################

          i am getting geo data with corresponding lat & lon from this script....

I want to include a line which will give the output as,

                       1. least number in the series(geo data) excludes zero's and negatives.

                       2. largest number in the series(geo data).

                       3. Average number of all geo data's in the series.
         
How to get these three things in that script ? which terms will give this as output ?

Cheers

Suresh M