1 |
dimitri |
1.1 |
krd=3; kgr=0; |
2 |
|
|
|
3 |
|
|
% $Header: /u/gcmpack/MITgcm/utils/matlab/cs_grid/bk_line/use_psiLine.m,v 1.1 2005/09/15 16:46:28 jmc Exp $ |
4 |
|
|
% $Name: $ |
5 |
|
|
listI=[02232 04320 06552 08712 10944 13104 15336 17568 19728 21960 24120 26352 ... |
6 |
|
|
28584 30600 32832 34992 37224 39384 41616 43848 46008 48240 50400 52632 ... |
7 |
|
|
54864 56880 59112 61272 63504 65664 67896]; |
8 |
|
|
listR=listI(end-11:end); |
9 |
|
|
deltaT=listI(end-12:end-1); deltaT=listR-deltaT; |
10 |
|
|
|
11 |
|
|
Rac='grid_files/'; |
12 |
|
|
drF=squeeze(rdmds([Rac,'DRF']))'; nr=length(drF); |
13 |
|
|
|
14 |
|
|
if rem(krd,2) == 1, |
15 |
|
|
%- load the grid dx,dy , convert to 10^6 m : |
16 |
|
|
dyg=rdmds([Rac,'DYG']); dyg=dyg*1.e-6; |
17 |
|
|
fprintf(' load dyg file\n'); |
18 |
|
|
nc=size(dyg,2); nPg=nc*nc*6; |
19 |
|
|
dyg=reshape(dyg,[nPg 1]); |
20 |
|
|
end |
21 |
|
|
|
22 |
|
|
if rem(krd,2) == 1, |
23 |
|
|
%-- load velocity : |
24 |
|
|
rac='output_files/'; |
25 |
|
|
%nit=28584; |
26 |
|
|
%uu=rdmds([rac,'UVELMASS'],nit); |
27 |
|
|
ut=zeros(nPg,1); |
28 |
|
|
for n=1:length(listR), |
29 |
|
|
nit=listR(n); |
30 |
|
|
namf=sprintf([rac,'UVELMASS.%10.10i.data'],nit); |
31 |
|
|
fprintf([' process file: ',namf,' ...']); |
32 |
|
|
for k=1:nr, |
33 |
|
|
fprintf(' %i',k); |
34 |
|
|
uk=rdda(namf,[nPg 1],k,'real*4','b'); |
35 |
|
|
uk=dyg.*uk; |
36 |
|
|
ut=ut+drF(k)*deltaT(n)*uk; |
37 |
|
|
end |
38 |
|
|
fprintf(' done\n'); |
39 |
|
|
end |
40 |
|
|
ut=ut/sum(deltaT); |
41 |
|
|
namfil='UintTrsYr.bin'; |
42 |
|
|
fprintf([' write to file: ',namfil,' ...']); |
43 |
|
|
fid=fopen(namfil,'w','b'); fwrite(fid,ut,'real*8'); fclose(fid); |
44 |
|
|
fprintf(' done\n'); |
45 |
|
|
end |
46 |
|
|
|
47 |
|
|
if krd > 1, |
48 |
|
|
clear dyg uk ut |
49 |
|
|
dxg=rdmds([Rac,'DXG']); dxg=dxg*1.e-6; |
50 |
|
|
fprintf(' load dxg file\n'); |
51 |
|
|
nc=size(dxg,2); nPg=nc*nc*6; |
52 |
|
|
dxg=reshape(dxg,[nPg 1]); |
53 |
|
|
end |
54 |
|
|
|
55 |
|
|
if krd > 1, |
56 |
|
|
%-- load velocity : |
57 |
|
|
rac='output_files/'; |
58 |
|
|
%nit=28584; |
59 |
|
|
%vv=rdmds([rac,'VVELMASS'],nit); |
60 |
|
|
vt=zeros(nPg,1); |
61 |
|
|
for n=1:length(listR), |
62 |
|
|
nit=listR(n); |
63 |
|
|
namf=sprintf([rac,'VVELMASS.%10.10i.data'],nit); |
64 |
|
|
fprintf([' process file: ',namf,' ...']); |
65 |
|
|
for k=1:nr, |
66 |
|
|
fprintf(' %i',k); |
67 |
|
|
vk=rdda(namf,[nPg 1],k,'real*4','b'); |
68 |
|
|
vk=dxg.*vk; |
69 |
|
|
vt=vt+drF(k)*deltaT(n)*vk; |
70 |
|
|
end |
71 |
|
|
fprintf(' done\n'); |
72 |
|
|
end |
73 |
|
|
vt=vt/sum(deltaT); |
74 |
|
|
namfil='VintTrsYr.bin'; |
75 |
|
|
fprintf([' write to file: ',namfil,' ...']); |
76 |
|
|
fid=fopen(namfil,'w','b'); fwrite(fid,vt,'real*8'); fclose(fid); |
77 |
|
|
fprintf(' done\n'); |
78 |
|
|
end |
79 |
|
|
|
80 |
|
|
return |