t=20.0*ones([20 20 20]); 
t(:,:,1)=t(:,:,1)+0.01*(rand([20 20 1])-0.5);
fid=fopen('tinit.bin','w','b');
fwrite(fid,t,'float64');