dsls=dir('netcdfs/*.nc')

for dd=1:length(dsls)
    ncload(dsls(dd).name)
end

%_t : total
%_m : mean
%_s : stationary eddy
%_e : transient eddy