% extract December 2011 SST
nx=4320;
region_name='SST';
pin='~dmenemen/llc_4320/MITgcm/run/';
pout=['~dmenemen/llc_4320/regions/' region_name '/'];
eval(['cd ' pout])
for ts=283536:144:390528
    fin=[pin 'Theta.' myint2str(ts,10) '.data'];
    dy=ts2dte(ts,25,2011,9,10,30);
    fout=['SST.' dy];
    fld=quikread_llc(fin,nx,1,'real*4');
    writebin(fout,fld);
end
