PRO Read_JPL, in_path, $
              thisDATE, $   
              level, $  ; L1a / L1b / L2
              output_file, $   
              return_code=return_code

return_code=0

path_array = strsplit(in_path, '/', /extract)
mission = path_array[-1]
thisYEAR=strmid(thisDATE, 0, 4)

if            mission eq 'cosmic2' then begin
  file_patten=in_path+'/c2f?/v2qp*0d_ncep_10Kp1N.20km/y'+thisYEAR+'/'+thisDATE+'/'+level+'/txt/*.txt.gz'

endif else if mission eq 'cosmic' then begin
  file_patten=in_path+'/cosmic?/postproc/y'+thisYEAR+'/'+thisDATE+'/'+level+'/txt/*.txt.gz'

endif else begin
  file_patten=in_path+'/postproc/y'+thisYEAR+'/'+thisDATE+'/'+level+'/txt/*.txt.gz'

endelse

file_list=file_search(file_patten, count=file_count)

print,file_patten, string(10B)
print,'file_count= ',file_count
if file_count eq 0 then begin
  return_code=1
  goto, jump_end
endif  

; ------------------------------------------------
; PRN_ID is the commen id (01-32) used by UCAR
; GPS_ID is the uniq id (01-73)
Extract_GPS_almanacs_sem, GPS_ID, PRN_ID, thisDATE

; ------------------------------------------------
; reading
data=[]
for file_loop = 0, file_count-1 do begin

  print,file_list[file_loop]

  if      level eq 'L1a' then status = Extract_JPLatmPhs(file_list[file_loop], data_this, thisDATE, mission, GPS_ID, PRN_ID) $
  else if level eq 'L1b' then status = Extract_JPLdry(file_list[file_loop], data_this, thisDATE, mission, GPS_ID, PRN_ID) $
  else if level eq 'L2'  then status = Extract_JPLwet(file_list[file_loop], data_this, thisDATE, mission, GPS_ID, PRN_ID)

  print, status
  if status eq 0 then continue
  ;print, data_this

  if n_elements(data_this) gt 0 then data =[data,  data_this]

endfor ; file_loop  

; no ref
if n_elements(data) gt 0 then save, data, filename=output_file else print,'No data input'

print,'end at : ',systime() 
return_code=2  

jump_end:
stop
END