MITgcm_contrib/lab_sea_test/mk_tair.m

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
1	nx=20; ny=16; nt=14;
2	fin='tair.labsea1979';
3	fld=readbin(fin,[nx ny nt]);
4	for i=1:20
5	tmp=fld;
6	tmp(10,10,1)=tmp(10,10,1)+.01*i;
7	fout=['tair' myint2str(i) '.labsea1979'];
8	writebin(fout,tmp)
9	end