atm_gray/inp_ll/gendata.m


fname='northEquatMask.bin';
nx=128; ny=64;

msk=zeros(nx,ny);
msk(:,33:36)=1.;
msk(41:nx,33:36)=2.;

%figure(1);clf;
%imagesc(msk'); set(gca,'YDir','normal');
%colorbar

%fid=fopen(fname,'w','b'); fwrite(fid,msk,'real*8'); fclose(fid);
%fprintf(['write file',fname,'\n']);

%-- cos shape SST:
yy=[1:ny]; yy=yy-mean(yy);
yy=yy*pi*2/ny;
sst=ones(nx,1)*cos(yy);
sst=273.15+9+19*sst;

var=rand([nx,ny]); var=var-mean(var(:));
var=var.*(1+ones(nx,1)*cos(yy));
figure(1);clf;
imagesc(var'); set(gca,'YDir','normal');
colorbar

noise=0.;
%noise=1.e-3;
sst=sst+noise*var;

%rDir='res_n26/';
%namf='SST.check';
%sst=rdmds([rDir,namf]);

figure(3);clf;
var=sst;
zav=mean(var); var=var-ones(nx,1)*zav;
subplot(211);
imagesc(var'); set(gca,'YDir','normal');
colorbar
subplot(212);
plot([1:ny],zav'-273.15);

%fname='SST_sym.bin';
fname='SST_cos0.bin';
%fname='SST_cos1.bin'; % with: noise=1.e-3;
fid=fopen(fname,'w','b'); fwrite(fid,sst,'real*8'); fclose(fid);
fprintf(['write file',fname,'\n']);
1
2	fname='northEquatMask.bin';
3	nx=128; ny=64;
4
5	msk=zeros(nx,ny);
6	msk(:,33:36)=1.;
7	msk(41:nx,33:36)=2.;
8
9	%figure(1);clf;
10	%imagesc(msk'); set(gca,'YDir','normal');
11	%colorbar
12
13	%fid=fopen(fname,'w','b'); fwrite(fid,msk,'real*8'); fclose(fid);
14	%fprintf(['write file',fname,'\n']);
15
16	%-- cos shape SST:
17	yy=[1:ny]; yy=yy-mean(yy);
18	yy=yypi2/ny;
19	sst=ones(nx,1)*cos(yy);
20	sst=273.15+9+19*sst;
21
22	var=rand([nx,ny]); var=var-mean(var(:));
23	var=var.(1+ones(nx,1)cos(yy));
24	figure(1);clf;
25	imagesc(var'); set(gca,'YDir','normal');
26	colorbar
27
28	noise=0.;
29	%noise=1.e-3;
30	sst=sst+noise*var;
31
32	%rDir='res_n26/';
33	%namf='SST.check';
34	%sst=rdmds([rDir,namf]);
35
36	figure(3);clf;
37	var=sst;
38	zav=mean(var); var=var-ones(nx,1)*zav;
39	subplot(211);
40	imagesc(var'); set(gca,'YDir','normal');
41	colorbar
42	subplot(212);
43	plot([1:ny],zav'-273.15);
44
45	%fname='SST_sym.bin';
46	fname='SST_cos0.bin';
47	%fname='SST_cos1.bin'; % with: noise=1.e-3;
48	fid=fopen(fname,'w','b'); fwrite(fid,sst,'real*8'); fclose(fid);
49	fprintf(['write file',fname,'\n']);