function make_AoM2023_CL31_30s_v1_0(inname,outname)
% Script to make output ceilometer file from AoM23 expedition. Keep native
% resolution, set data period to be same as for WX, and change some of the
% missing data flags to be nan.
%
% JZP January 2019
% mod SM Feb 2024

%% load data
load(inname);

%% specify dates for which measurements could be made
% exclude periods when Oden in Helsingborg but include transit
startdate = datenum(2023,05,07,0,0,0);
enddate = datenum(2023,06,13,16,24,35);

ii = cl31.mday >= startdate & cl31.mday < enddate;

fields = fieldnames(cl31);
for n = 1:length(fields) - 1
    cl31.(fields{n}) = cl31.(fields{n})(ii,:);
end

cl31.time = (cl31.mday - datenum(1970,1,1,0,0,0))*86400; % = unix time


%% change cloud base heights to incorporate measurement height
cl31.base_ht = cl31.base_ht + 25;
cl31.sc_ht = cl31.sc_ht + 25;

%% Change missing data values in sky condition to nan. NOT NEEDED HERE
% ii = cl31.sc_frac == 99 | cl31.sc_frac < 0;
% cl31.sc_frac(ii) = nan; 
% cl31.sc_ht(ii) = nan;

%% save
save(outname,'cl31')