exp1/input/gendata.m

ieee='b';
accuracy='real*8';

Ho=2000;
nx=60;
ny=60;

% Flat bottom at z=-Ho
h=-Ho*ones(nx,ny);
% Walls
h(end,:)=0;
h(:,end)=0;
fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);

% Wind-stress
tauMax=0.1;
x=((1:nx)-0.5)/(nx-1); % nx-1 accounts for a solid wall
y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall
[X,Y]=ndgrid(x,y);
tau=tauMax*sin(pi*Y);
fid=fopen('windx.sin_y','w',ieee); fwrite(fid,tau,accuracy); fclose(fid);
1	adcroft	1.2	ieee='b';
2			accuracy='real*8';
3
4	adcroft	1.1	Ho=2000;
5			nx=60;
6			ny=60;
7
8			% Flat bottom at z=-Ho
9			h=-Ho*ones(nx,ny);
10			% Walls
11			h(end,:)=0;
12			h(:,end)=0;
13	adcroft	1.2	fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
14	adcroft	1.1
15			% Wind-stress
16			tauMax=0.1;
17			x=((1:nx)-0.5)/(nx-1); % nx-1 accounts for a solid wall
18			y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall
19			[X,Y]=ndgrid(x,y);
20			tau=tauMaxsin(piY);
21	adcroft	1.2	fid=fopen('windx.sin_y','w',ieee); fwrite(fid,tau,accuracy); fclose(fid);