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