Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
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')