1 |
|
2 |
%- To allow to split the domain in several tiles, change the domain size |
3 |
% from 31x31 to 32x32 by adding 1 row (on Southern side) + 1 column (on |
4 |
% Western side) of land points. |
5 |
|
6 |
%- read old input files (31x31) and write new ones: |
7 |
nx1=31; ny1=31; |
8 |
nx2=1+nx1; |
9 |
ny2=1+ny1; |
10 |
prec='real*4'; |
11 |
|
12 |
fnam1='topog.box'; |
13 |
fnam2='topo_box.bin'; |
14 |
fprintf('converting %s to %s ...',fnam1,fnam2) |
15 |
vv1=rdda(fnam1,[nx1 ny1],1,prec,'b'); |
16 |
vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1; |
17 |
fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid); |
18 |
fprintf(' done\n'); |
19 |
%figure(1);clf; imagesc(vv2');set(gca,'YDir','normal');colorbar |
20 |
|
21 |
fnam1='windx.cos_y'; |
22 |
fnam2='taux_cosY.bin'; |
23 |
fprintf('converting %s to %s ...',fnam1,fnam2) |
24 |
vv1=rdda(fnam1,[nx1 ny1],1,prec,'b'); |
25 |
vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1; |
26 |
fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid); |
27 |
fprintf(' done\n'); |
28 |
|
29 |
fnam1='trinit1'; |
30 |
fnam2='tr1_ini.bin'; |
31 |
fprintf('converting %s to %s ...',fnam1,fnam2) |
32 |
vv1=rdda(fnam1,[nx1 ny1],1,prec,'b'); |
33 |
vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1; |
34 |
fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid); |
35 |
fprintf(' done\n'); |
36 |
|
37 |
prec='real*8'; |
38 |
fnam1='pickup'; it=200000; |
39 |
fnam2='pickup.32x32.data'; |
40 |
fprintf('converting %s to %s ...',fnam1,fnam2) |
41 |
vv1=rdmds(fnam1,it); n3d=size(vv1,3); |
42 |
vv2=zeros(nx2,ny2,n3d); vv2(2:nx2,2:ny2,:)=vv1; |
43 |
fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid); |
44 |
fprintf(' done\n'); |
45 |
|