% Recreate Depth.data from hFacC.data and thk90.mat
clear, close all
cd ~dmenemen/llc_4320/grid
nx=4320; nz=90;
load thk90
tmp=quikread_llc('hFacC.data',4320);
Depth=tmp*thk90(1);
for k=2:nz, disp(k)
    tmp=quikread_llc('hFacC.data',4320,k);
    Depth=Depth+tmp*thk90(k);
end
writebin('Depth.data',Depth);