36 |
% Temperature profile |
% Temperature profile |
37 |
[sprintf('Tref =') sprintf(' %8.6g,',Tz*z-mean(Tz*z))] |
[sprintf('Tref =') sprintf(' %8.6g,',Tz*z-mean(Tz*z))] |
38 |
|
|
39 |
|
ieee='b'; |
40 |
|
accuracy='real*8'; |
41 |
|
|
42 |
% Gaussian bump |
% Gaussian bump |
43 |
h=-H+dh*exp( -(X.^2+Y.^2)/(2*(L^2)) ); |
h=-H+dh*exp( -(X.^2+Y.^2)/(2*(L^2)) ); |
44 |
f77write('topog.bump',h,'real*8','b'); |
fid=fopen('topog.bump','w',ieee); fwrite(fid,h,accuracy); fclose(fid); |
45 |
|
|
46 |
% Side walls + bump |
% Side walls + bump |
47 |
h(:,1)=0; |
h(:,1)=0; |
48 |
h(:,ny)=0; |
h(:,ny)=0; |
49 |
f77write('topog.bumpchannel',h,'real*8','b'); |
fid=fopen('topog.bumpchannel','w',ieee); fwrite(fid,h,accuracy); fclose(fid); |
50 |
|
|
51 |
% Simple channel |
% Simple channel |
52 |
h(:,1)=0; |
h(:,1)=0; |
53 |
h(:,2:ny-1)=-H; |
h(:,2:ny-1)=-H; |
54 |
h(:,ny)=0; |
h(:,ny)=0; |
55 |
f77write('topog.channel',h,'real*8','b'); |
fid=fopen('topog.channel','w',ieee); fwrite(fid,h,accuracy); fclose(fid); |