function out = read_col_int(DDD, f_path)
%{
    Reads column integral data for MP3000 radiometer
        from parsed ColInt_Obs.txt file.
    Returns out array of JD,Int_vap_cm, Int_liq_mm, Cld_base_km
    Inputs: DDD - day-of-year
            f_path - string path to ColInt_Obs.txt file

    Input file format:
    Record ,Date_Time         ,Record_Type ,Int_Vapor_cm ,Int_Liquid_mm ,Cloud_Base_km
    33     ,02/10/15 00:01:43 ,301         ,5.193        ,0.000         ,7.750
    46     ,02/10/15 00:01:45 ,301         ,5.011        ,0.035         ,7.750
    47     ,02/10/15 00:01:45 ,301         ,5.003        ,0.015         ,7.750
    48     ,02/10/15 00:01:45 ,301         ,5.007        ,0.025         ,7.750
%}

%% initialize
disp(['Reading MP3000 data for day ',int2str(DDD)]);
fin = fopen(f_path);

temp = {};    % empty arrays for textscan
in = [];

%% read entire file
while feof(fin)==0   % read entire file into cell array
    temp = textscan(fin,'%*f %2f/%2f/%2f %2f:%2f:%2f %*f %f %f %f','delimiter',',','emptyvalue',NaN);
    in = [in cell2mat(temp)];  % N x 9 array: MM,DD,YY,hh,mm,ss,vap,liq,cldht
end
fclose(fin);
jd = md2yd_vn(in(:,3)+2000,in(:,1),in(:,2)) + in(:,4)/24 + in(:,5)/1440 + in(:,6)/86400;

ii = jd>DDD & jd<DDD+1;  % select data for desired day of year

out = [jd(ii),in(ii,7:9)];
    
end