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')