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	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