high_res_cube/matlab/v_integrate.m

krd=3; kgr=0;

% $Header: /u/gcmpack/MITgcm/utils/matlab/cs_grid/bk_line/use_psiLine.m,v 1.1 2005/09/15 16:46:28 jmc Exp $
% $Name:  $
listI=[02232 04320 06552 08712 10944 13104 15336 17568 19728 21960 24120 26352 ...
       28584 30600 32832 34992 37224 39384 41616 43848 46008 48240 50400 52632 ...
       54864 56880 59112 61272 63504 65664 67896];
listR=listI(end-11:end);
deltaT=listI(end-12:end-1); deltaT=listR-deltaT;

Rac='grid_files/';
drF=squeeze(rdmds([Rac,'DRF']))'; nr=length(drF);

if rem(krd,2) == 1,
%- load the grid dx,dy , convert to 10^6 m :
 dyg=rdmds([Rac,'DYG']); dyg=dyg*1.e-6;
 fprintf(' load dyg file\n');
 nc=size(dyg,2); nPg=nc*nc*6;
 dyg=reshape(dyg,[nPg 1]);
end

if rem(krd,2) == 1,
%-- load velocity :
 rac='output_files/';
%nit=28584;
%uu=rdmds([rac,'UVELMASS'],nit);
 ut=zeros(nPg,1);
 for n=1:length(listR),
  nit=listR(n);
  namf=sprintf([rac,'UVELMASS.%10.10i.data'],nit);
  fprintf([' process file: ',namf,' ...']);
  for k=1:nr,
   fprintf(' %i',k);
   uk=rdda(namf,[nPg 1],k,'real*4','b');
   uk=dyg.*uk;
   ut=ut+drF(k)*deltaT(n)*uk;
  end
  fprintf(' done\n');
 end
 ut=ut/sum(deltaT);
 namfil='UintTrsYr.bin';
 fprintf([' write to file: ',namfil,' ...']);
 fid=fopen(namfil,'w','b'); fwrite(fid,ut,'real*8'); fclose(fid);
 fprintf(' done\n');
end

if krd > 1,
 clear dyg uk ut 
 dxg=rdmds([Rac,'DXG']); dxg=dxg*1.e-6;
 fprintf(' load dxg file\n');
 nc=size(dxg,2); nPg=nc*nc*6;
 dxg=reshape(dxg,[nPg 1]);
end

if krd > 1,
%-- load velocity :
 rac='output_files/';
%nit=28584;
%vv=rdmds([rac,'VVELMASS'],nit);
 vt=zeros(nPg,1);
 for n=1:length(listR),
  nit=listR(n);
  namf=sprintf([rac,'VVELMASS.%10.10i.data'],nit);
  fprintf([' process file: ',namf,' ...']);
  for k=1:nr,
   fprintf(' %i',k);
   vk=rdda(namf,[nPg 1],k,'real*4','b');
   vk=dxg.*vk;
   vt=vt+drF(k)*deltaT(n)*vk;
  end
  fprintf(' done\n');
 end
 vt=vt/sum(deltaT);
 namfil='VintTrsYr.bin';
 fprintf([' write to file: ',namfil,' ...']);
 fid=fopen(namfil,'w','b'); fwrite(fid,vt,'real*8'); fclose(fid);
 fprintf(' done\n');
end

return
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=ncnc6;
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=ncnc6;
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