close all
clear all
clc


load FMT.mat;
load grid.mat;
[nxc nyc nzc] = size(hFacC);
nt = 3743;

%% old packing
time = rdslice(['TIME_CCS_2007_2017.bin'],[nt 1],1, fmt, Ieee);

for it = [1 nt]
    u =  rdslice(['U_MIT_CCS_2007_2017_FULL.bin'],[nxc nyc nzc],it, fmt, Ieee);
    v =  rdslice([ 'V_MIT_CCS_2007_2017_FULL.bin'],[nxc nyc nzc],it, fmt, Ieee);
end

