PRO Read_RAOB_adpupa_bufr, decodefile, $
                         out_path, $
                         return_code=return_code

return_code=0
print,string(10B)+'input_file:'+decodefile

dirname = file_dirname(decodefile)
basename = file_basename(decodefile)    
dirname_array=strsplit(dirname, '/', /extract)
input_file_basename = dirname_array[-1]+'/'+basename

;--------------------------------------
; reading
status= Extract_RAOB_adpupa_bufr(decodefile, data_input, input_file_basename)
if status eq 0 then goto, jump_end

thisDATE = data_input.thisDATE
thisDATE_u = thisDATE[uniq(thisDATE, sort(thisDATE))]

for date_loop = 0, n_elements(thisDATE_u)-1 do begin

  data_today = data_input[where(thisDATE eq thisDATE_u[date_loop])]

    output_file = out_path+'/raob_'+thisDATE_u[date_loop]+'.sav'
    file_mkdir, out_path
    print,string(10B)+'output_file:',output_file

    data_previous=[]
    if file_test(output_file) then begin
      restore,output_file

      ; search the profiles which should not be replaced
      valid_id = where(data.filename ne input_file_basename, valid_num)
      if valid_num gt 0 then data_previous = data[valid_id]
    endif 

    data = [data_previous, data_today]
    save, data, filename=output_file
 
endfor ; date
 
;print,'end at : ',systime() 
return_code=2  

jump_end:
stop
END