1 |
|
2 |
nx=20; ny=16; nit=62; |
3 |
|
4 |
nrc=nx*ny*nit; |
5 |
%- check size: |
6 |
!ls -l kf_fu kf_fv |
7 |
fprintf('expected size to read : %i\n',nrc*4) |
8 |
|
9 |
%- read wind stress |
10 |
namF='kf_fu'; fprintf([' read file: ',namF,' ...']); |
11 |
fid=fopen(namF,'r','b'); tx=fread(fid,nrc,'real*4'); fclose(fid); |
12 |
% tx=rdda(namF,[nx ny nit],1,'real*4','b'); |
13 |
fprintf(' done\n'); |
14 |
|
15 |
namF='kf_fv'; fprintf([' read file: ',namF,' ...']); |
16 |
fid=fopen(namF,'r','b'); ty=fread(fid,nrc,'real*4'); fclose(fid); |
17 |
%ty=rdda(namF,[nx ny nit],1,'real*4','b'); |
18 |
fprintf(' done\n'); |
19 |
|
20 |
%- from "usual" old convention to MITgcm sign convention: flux > 0 is upward |
21 |
% reverse sign of tx,ty |
22 |
tx=-tx; |
23 |
ty=-ty; |
24 |
|
25 |
namF='kf_tx'; fprintf([' write tx to file: ',namF,' ...']); |
26 |
fid=fopen(namF,'w','b'); fwrite(fid,tx,'real*4'); fclose(fid); |
27 |
fprintf(' done\n'); |
28 |
|
29 |
namF='kf_ty'; fprintf([' write ty to file: ',namF,' ...']); |
30 |
fid=fopen(namF,'w','b'); fwrite(fid,ty,'real*4'); fclose(fid); |
31 |
fprintf(' done\n'); |
32 |
|