exp1/input/gendata.m

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

% Flat bottom at z=-Ho
h=-Ho*ones(nx,ny);
% Walls
h(end,:)=0;
h(:,end)=0;
f77write('topog.box',h,'real*8','b');

% 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);
f77write('windx.sin_y',tau,'real*8','b');
1	adcroft	1.1	Ho=2000;
2			nx=60;
3			ny=60;
4
5			% Flat bottom at z=-Ho
6			h=-Ho*ones(nx,ny);
7			% Walls
8			h(end,:)=0;
9			h(:,end)=0;
10			f77write('topog.box',h,'real*8','b');
11
12			% Wind-stress
13			tauMax=0.1;
14			x=((1:nx)-0.5)/(nx-1); % nx-1 accounts for a solid wall
15			y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall
16			[X,Y]=ndgrid(x,y);
17			tau=tauMaxsin(piY);
18			f77write('windx.sin_y',tau,'real*8','b');