krd=3; kgr=0;

% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/high_res_cube/matlab/v_integrate.m,v 1.1 2007/03/28 18:07:27 dimitri 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