nx=20; ny=16; nt=14;
fin='tair.labsea1979';
fld=readbin(fin,[nx ny nt]);
for i=1:20
  tmp=fld;
  tmp(10,10,1)=tmp(10,10,1)+.01*i;
  fout=['tair' myint2str(i) '.labsea1979'];
  writebin(fout,tmp)
end