;procedure to read Huntsville sonde .LE1 files
; by Wesley Cantrell
;wes.cantrell@nsstc.uah.edu




PRO   Read_sonde_LE1, numfiles, file, array_size, range, flight_num, date, time, station

  array_size=600.        ;all altitudes
  ;Data dummy variables
  o3mr=FLTARR(array_size)
  alt=o3mr  & temp=o3mr & press=o3mr & o3press=o3mr
  

 ;Header sizes
    HU_le10=15  & HU_le11=14 & l100=18  & HU_DAT=19     
  ;------------ 
  headertest=STRARR(HU_le10)
  
  
     OPENR,lun,file,/GET_LUN
     READF, lun, headertest
     SWITCH STRPOS(headertest[3],'Voemel') OF   
         55: header=STRARR(15)         
         55: date=8
         55: time=9
         55: station=5
         55: indx=7
         55: break
         -1: header=STRARR(14)
         -1: date=7
         -1: time=8
         -1: station=4
         -1: indx=6
         -1: break
     ENDSWITCH
     

     
     POINT_LUN, lun, 0

    READF, lun, header
    data_dum=DBLARR(13)
    station=STRTRIM(STRMID(header[station], 16),2)

    
    flight_num=STRTRIM(STRMID(header[indx],16),2)
    date=STRTRIM(STRMID(header[date],16),2)
    time=STRTRIM(STRMID(header[time],16),2)
    ;month=STRMID(STRTRIM(STRMID(STRTRIM(STRMID(header[date],16),2),2),2),13,3,/REVERSE_OFFSET)
    ;hour=STRMID(STRTRIM(STRMID(header[time],16),2),0,2)

    i=0
    WHILE NOT EOF(lun) DO BEGIN
      READF, lun, data_dum
      IF(data_dum[2] EQ 0. ) THEN BEGIN
        i++
        CONTINUE
      ENDIF
      o3mr[i]=data_dum[8]*1000  ;ppbv
      press[i]=data_dum[1]
      temp[i]=data_dum[4]
      o3press[i]=data_dum[7]
      alt[i]=data_dum[2]
      
      i++
    ENDWHILE
  FREE_LUN, lun


stop
END