% base directory
pn='/data8/llc/llc_2160/regions/latlon/';
eval(['cd ' pn 'matlab'])

% file dimensions
nx=8640;
ny=4320;
nz=1;
suf2d=['_' int2str(nx) 'x' int2str(ny)]; % 2d suffix
suf3d=[suf2d 'x' int2str(nz)];       % 3d suffix

% plotting parameters
dt='22-Apr-2013 6:00:00';     % date to plot
cx=[-1 1]/2;

% U/V figure
fn=[pn 'U/U' suf3d '.' datestr(datenum(dt),30)];
U=readbin(fn,[nx ny]);
fn=[pn 'V/V' suf3d '.' datestr(datenum(dt),30)];
V=readbin(fn,[nx ny]);
clf reset
subplot(211)
mypcolor(U')
caxis(cx)
thincolorbar
subplot(212)
mypcolor(V')
caxis(cx)
thincolorbar
print -dpng ../figs/UV
