% List of site names
site_names = {'Site16', 'Site17'};

% Initialize an empty struct to hold the data
site_info = struct();

% Loop through each site name
for i = 1:numel(site_names)
    site_name = site_names{i};
    mat_filename = [site_name '.mat'];
    
    % Load the MAT file
    load(mat_filename, 'sar_filenames', 'csv', 'ind');
    
    % Find the latitude and longitude boundaries
    [lat, lon] = FindLatLon(mat_filename);
    
    if ~exist('csv', 'var')
        csv = [];
    end
    if ~exist('ind', 'var')
        ind = [];
    end
    
    % Create a field in the struct for this site
    site_info.(site_name) = struct('sar_filenames', {sar_filenames}, ...
                                   'csv', csv, ...
                                   'ind', ind, ...
                                   'lat', lat, ...
                                   'lon', lon);
    clear sar_filenames csv ind lat lon
end