pro process_cips_days_l3, startDay, stopDay, preliminary=preliminary, pca_preliminary=pca_preliminary

  set_cips_production_vars
  
  set_production_flag, flag = (keyword_set(preliminary) ? 0 : 1)

  for index = startDay, stopDay do begin
    
    ;l3a = obj_new('level_3a_daisy_server', obj_new())
    
    start_day = index
    stop_day = start_day + 1
    print, "Begin processing day: "+start_day

    start_time = yd2usec(start_day)
    stop_time = yd2usec(stop_day)
    orbit_info_array = get_orbit_info(start_time, stop_time)
    ; let's check for valid orbits, ie are there any in the missing orbits file
    orbit_array = []
    foreach orbit, orbit_info_array do begin
      if is_orbit_valid( orbit.orbit_number ) then orbit_array = [orbit_array, orbit]
    endforeach

    ;result = l3a->get_data( orbit_info_array, /net_cdf )
    if keyword_set(pca_preliminary) then begin
       l3a_daisy_v52, orbit_array, /pca_preliminary
    endif else begin
       l3a_daisy_v52, orbit_array
    endelse
    ;str_free, result
    ;result=0
    ; clean up
    start_time = 0
    stop_time = 0
    str_free, orbit_info_array
    str_free, orbit_array
    orbit_info_array = 0
    orbit_array = 0
    heap_gc

  endfor
  
end