m=936;
n=1062;
k=90;
suf=['_' int2str(m) 'x' int2str(n)];

% read grid infomration files
pn='~dmenemen/llc_4320/regions/Boxes/Box12/grid/';
for fn={'Depth','RAC','XC','YC','XG','YG'}
    eval([fn{1} '=readbin([pn fn{1} suf],[m n]);']);
end
hFacC=readbin([pn 'hFacC' suf 'x' int2str(k)],[m n k]);