clear
scxy=loadgrid('../MESH/sab_clim2');
sacr=loadgrid('sab_ec2000_200');
CASE=['Jan';'Feb';'Mar';'Apr';'May';'Jun';'Jul';'Aug';'Sep';'Oct';'Nov';'Dec'];
nsamp=50;
  xv=[sacr.x(2091) sacr.x(1712)];
  yv=[sacr.y(2091) sacr.y(1712)];
for i=1:12
  DATADIR=['../RESULTS/',CASE(i,1:3)];
  icqfile=[DATADIR '/F6.HOT.', CASE(i,1:3), '_nub.icq4'];
  icq4datanowi=read_icq4(icqfile); 
  nn=icq4datanowi.nn;
  nnv=icq4datanowi.nnv;

  VRT=icq42vrt(icq4datanowi);
  VRT(:,:,9)=sw_dens0(VRT(:,:,8),VRT(:,:,7));
  [RSIG{i},RR{i}]=transect(scxy,VRT,xv,yv,nsamp);
  [Tnowi{i},unnowi{i},upnowi{i}]=transport(...
                RSIG{i}.d,xv,yv,RSIG{i}.q1,RSIG{i}.q2,RSIG{i}.q3);
   zeta{i}=icq4datanowi.HMID-scxy.z';
   Vtop{i}=[(mean((icq4datanowi.UZMID(:,ceil(nnv/2):nnv))'))'... 
   	    (mean((icq4datanowi.VZMID(:,ceil(nnv/2):nnv))'))'];
   Vdow{i}=[(mean((icq4datanowi.UZMID(:,1:ceil(nnv/2)))'))'...
   	    (mean((icq4datanowi.VZMID(:,1:ceil(nnv/2)))'))'];
   Vsurf{i}=[icq4datanowi.UZMID(:,nnv) icq4datanowi.VZMID(:,nnv)];
   Vbar{i}=[icq4datanowi.UMID icq4datanowi.VMID];
   Vbot{i}=[icq4datanowi.UZMID(:,1) icq4datanowi.VZMID(:,1)];
   VR6{i}=[icq4datanowi.UZMID(8607,:) icq4datanowi.VZMID(8607,:)];
   Tsurf{i}=icq4datanowi.TMPMID(:,21);
   Ssurf{i}=icq4datanowi.SALMID(:,21);
   Tbot{i}=icq4datanowi.TMPMID(:,1);
   Sbot{i}=icq4datanowi.SALMID(:,1);
   RHOsurf{i}=VRT(nnv,:,9);
   RHObot{i}=VRT(1,:,9);
%   RHOdepth=VRT(1:nnv,:,9);
%   for kk=1:nn;
%      irhomax(kk)=find(RHOdepth(:,kk)==max(RHOdepth(:,kk)));
%      RHOmax{i}(1,kk)=RHOdepth(irhomax(kk),kk);
%      depth{i}(1,kk)=-icq4datanowi.ZMID(kk,irhomax(kk));
%   end
%   clear RHOdepth;
%   DATADIR='/afs/isis.unc.edu/depts/marine/opnml/scratch/f6_winds/';
%   icqfile=[DATADIR, CASE(i,1:3), '_coads_wind.icq4'];
%   icq4wind=read_icq4(icqfile); 
%   %SW means southwest node 2057 ||||| R6 is taken as node 8607
%   zetawind{i}=icq4wind.HMID-scxy.z';
%  Vsurfwind{i}=[icq4wind.UZMID(:,nnv) icq4wind.VZMID(:,nnv)];
%   Vbarwind{i}=[icq4wind.UMID icq4wind.VMID];
%   Vbotwind{i}=[icq4wind.UZMID(:,1) icq4wind.VZMID(:,1)];
%   VSWwind{i}=[icq4wind.UZMID(2057,:) icq4wind.VZMID(2057,:)]; 
%   VR6wind{i}=[icq4wind.UZMID(8607,:) icq4wind.VZMID(8607,:)];   
%   VRTwind=icq42vrt(icq4wind);
%   [RSIGwind{i},RRwind{i}]=transect(scxy,VRTwind,xv,yv,nsamp);
%   [Twind{i},unwind{i},upwind{i}]=transport(...
%            RSIGwind{i}.d,xv,yv,RSIGwind{i}.q1,RSIGwind{i}.q2,RSIGwind{i}.q3);
%   RSIGwind{i}.q1=-RSIGwind{i}.q1;
%   RRwind{i}.q1=-RRwind{i}.q1;
%   RRwind{i}.y=-RRwind{i}.y;
end

%clear icq4datanowi VRT CASE DATADIR i icqfile icq4wind VRTwind
clear icq4datanowi VRT CASE DATADIR i icqfile 
save results_fundy_struct.mat

%   RSIG.d=RSIG.d/1000;   % convert to kms