% animation of SST and W at 100 m
cd ~dmenemen/llc_4320/regions/Canary
clf reset, orient tall, wysiwyg
for ts=10368:144:279360;
    dy=ts2dte(ts,25,2011,9,10,30);
    fin=['Theta/Theta_240x577x85.' dy];
    theta=readbin(fin,[240 577]);
    fin=['W/W_240x577x85.' dy];
    w=readbin(fin,[240 577],1,'real*4',22);
    clf, subplot(211); mypcolor(theta');
    caxis([20 25]); thincolorbar
    title(['SST (deg C) on ' ts2dte(ts,25,2011,9,10)])
    subplot(212); mypcolor(1000*w');
    caxis([-1 1]*5); thincolorbar
    title('Vertical velocity (mm/s) at 100 m depth (positive=upwelling)')
    text(150,100,ts2dte(ts,25,2011,9,10))
    pause(.1)
end
