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