1 |
jmc |
1.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 |
|
|
|