clear
load_constants
  LAT=[ 31.625  31.049  31.367  31.212  30.792  30.931  31.543  31.825 ...
	31.649]';
  LON=[-81.150 -80.480 -80.582 -80.163 -80.352 -80.771 -80.224 -80.304 ...
	-79.974]';
 trx0=-7.672028358888956e+06 ;lo0=-80.75;  
 try0=3.128247751710321e+06;la0=31.25;  
[Xd,Yd]=convll2m(LON,LAT,lo0,la0); 
Xd=Xd-trx0; Yd=Yd-try0; 
scll=loadgrid('../MESH/sab_clim2ll');
load results_fundy_struct.mat
outflag=0;
bfdbig=genbfd(scll,[50 50],outflag);
close 
bfdv=genbfd(scll,[25 25],outflag);
close 
for i=1:12
 RSIG{i}.d=RSIG{i}.d/1000;   % convert to kms
 RR{i}.x=RR{i}.x/1000;   % convert to kms
end


% Make surf elev plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   elev=zeta{i}';
   [c,h,hcb,DATA]=fdcontour(bfdbig,elev,'Qmin',-.75,'Qmax',-.25,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_surf_elev.epsc']);
   eval(['!convert ',cal.months{i},'_surf_elev.epsc ',cal.months{i},'_surf_elev.gif ']);
end

close all

% Make surf vel plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   plotbnd(scll); axeq;axtt;h=lcontour(scll,'z',[25 100 200]);
   vel=Vsurf{i}(:,1)+sqrt(-1)*Vsurf{i}(:,2);
   for iii=1:length(vel); 
      if (abs(vel(iii))>.45);vel(iii)=NaN;end;
      if (scll.z(iii)>300);vel(iii)=NaN;end;
   end;
   hv=fdvector(bfdv,real(vel),imag(vel),25,25,'ScaleLabel','m/s',...
'ScaleFac',.15,'Color','b','ScaleXor',-81.,'ScaleYor',33.1,'VecType','arrow');
   set(hv(2),'FontSize',11);
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_surf_vel.epsc']);
   eval(['!convert ',cal.months{i},'_surf_vel.epsc ',cal.months{i},'_surf_vel.gif ']);
end

close all

% Make depth averaged vel plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   plotbnd(scll); axeq;axtt;h=lcontour(scll,'z',[25 100 200]);
   vel=Vbar{i}(1:9606)'+sqrt(-1)*Vbar{i}(9607:end)';
   for iii=1:length(vel); 
      if (abs(vel(iii))>.5);vel(iii)=NaN;end;
      if (scll.z(iii)>250);vel(iii)=NaN;end;
   end;
   hv=fdvector(bfdv,real(vel),imag(vel),25,25,'ScaleLabel','m/s',...
'ScaleFac',.1,'Color','b','ScaleXor',-81.,'ScaleYor',33.1,'VecType','arrow');
   set(hv(2),'FontSize',11);
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_daver_vel.epsc']);
   eval(['!convert ',cal.months{i},'_daver_vel.epsc ',cal.months{i},'_daver_vel.gif ']);
end

close all

% Make (surface-bottom) density  plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   rhodiff=(RHOsurf{i}-RHObot{i})';
    [c,h,hcb,DATA]=fdcontour(bfdbig,rhodiff,'Qmin',-4,'Qmax',4,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_diff_dens.epsc']);
   eval(['!convert ',cal.months{i},'_diff_dens.epsc ',cal.months{i},'_diff_dens.gif ']);
end

close all

% Make temperature transect plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   pcolor(RSIG{i}.d,RSIG{i}.q1,RSIG{i}.q7)
   caxis([8 30]);  shading interp
   colormap(jet(16)); plotbnd(RR{i},'Color','k','LineWidth',2);
   ylim([-100 0])
   set(gca,'Ytick',[-100:20:0],'YtickLabel',[-100:20:0],'FontSize',12)
   set(gca,'Xtick',[0:20:120],'XtickLabel',[0:20:120],'FontSize',12)
   hcb=colorbar('horiz');
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.3 0.19 0.435 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.33,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_tran_temp.epsc']);
   eval(['!convert ',cal.months{i},'_tran_temp.epsc ',cal.months{i},'_tran_temp.gif ']);
end

close all

% Make salinity transect plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   pcolor(RSIG{i}.d,RSIG{i}.q1,RSIG{i}.q8)
   caxis([32 36.6]);  shading interp
   colormap(jet(16)); plotbnd(RR{i},'Color','k','LineWidth',2);
   ylim([-100 0])
   set(gca,'Ytick',[-100:20:0],'YtickLabel',[-100:20:0],'FontSize',12)
   set(gca,'Xtick',[0:20:120],'XtickLabel',[0:20:120],'FontSize',12)
   hcb=colorbar('horiz');
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.3 0.19 0.435 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.33,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_tran_sali.epsc']);
   eval(['!convert ',cal.months{i},'_tran_sali.epsc ',cal.months{i},'_tran_sali.gif ']);
end

close all

% Make density transect plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   pcolor(RSIG{i}.d,RSIG{i}.q1,RSIG{i}.q9-1000)
   caxis([22 27]);  shading interp
   colormap(jet(16)); plotbnd(RR{i},'Color','k','LineWidth',2);
   ylim([-100 0])
   set(gca,'Ytick',[-100:20:0],'YtickLabel',[-100:20:0],'FontSize',12)
   set(gca,'Xtick',[0:20:120],'XtickLabel',[0:20:120],'FontSize',12)
   hcb=colorbar('horiz');
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.3 0.19 0.435 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.33,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_tran_dens.epsc']);
   eval(['!convert ',cal.months{i},'_tran_dens.epsc ',cal.months{i},'_tran_dens.gif ']);
end

close all

% Make normal velocity transect plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   pcolor(RSIG{i}.d,RSIG{i}.q1,unnowi{i});hold on;
   caxis([-.2 .2]);  shading interp
  [trash,hcon]=contour(RSIG{i}.d,RSIG{i}.q1,unnowi{i},[0 0],'k-');hold off;
   colormap(jet(16)); plotbnd(RR{i},'Color','k','LineWidth',2);
   ylim([-100 0])
   set(gca,'Ytick',[-100:20:0],'YtickLabel',[-100:20:0],'FontSize',12)
   set(gca,'Xtick',[0:20:120],'XtickLabel',[0:20:120],'FontSize',12)
   hcb=colorbar('horiz');
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.3 0.19 0.435 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.33,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_tran_vel.epsc']);
   eval(['!convert ',cal.months{i},'_tran_vel.epsc ',cal.months{i},'_tran_vel.gif ']);
end

close all

% Make surface density  plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Rsu=(RHOsurf{i}-1000)';
    [c,h,hcb,DATA]=fdcontour(bfdbig,Rsu,'Qmin',22,'Qmax',27,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_surf_dens.epsc']);
   eval(['!convert ',cal.months{i},'_surf_dens.epsc ',cal.months{i},'_surf_dens.gif ']);
end

close all

% Make surface salinity  plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Ssu=(Ssurf{i});
    [c,h,hcb,DATA]=fdcontour(bfdbig,Ssu,'Qmin',32,'Qmax',36.6,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_surf_sali.epsc']);
   eval(['!convert ',cal.months{i},'_surf_sali.epsc ',cal.months{i},'_surf_sali.gif ']);
end

close all

% Make surface temperature plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Tsu=(Tsurf{i});
    [c,h,hcb,DATA]=fdcontour(bfdbig,Tsu,'Qmin',8,'Qmax',30,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_surf_temp.epsc']);
   eval(['!convert ',cal.months{i},'_surf_temp.epsc ',cal.months{i},'_surf_temp.gif ']);
end

close all

% Make bottom density  plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Rbo=(RHObot{i}-1000)';
    [c,h,hcb,DATA]=fdcontour(bfdbig,Rbo,'Qmin',22,'Qmax',27,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_bott_dens.epsc']);
   eval(['!convert ',cal.months{i},'_bott_dens.epsc ',cal.months{i},'_bott_dens.gif ']);
end

close all

% Make surface salinity  plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Sbo=(Sbot{i});
    [c,h,hcb,DATA]=fdcontour(bfdbig,Sbo,'Qmin',32,'Qmax',36.6,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_bott_sali.epsc']);
   eval(['!convert ',cal.months{i},'_bott_sali.epsc ',cal.months{i},'_bott_sali.gif ']);
end

close all

% Make surface temperature plot
for i=1:12
   figure
   ha=subplot(1,1,1);
   Tbo=(Tbot{i});
    [c,h,hcb,DATA]=fdcontour(bfdbig,Tbo,'Qmin',8,'Qmax',30,...
'NColMapBase',128,'NContour',32,'Colorbar','off');
   shading flat
   lndfill(scll,[1 1 1])
   plotbnd(scll);axeq;axtt;
   hc=lcontour(scll,'z',[25 100 200]);maskit300(0,'w');
   set(gca,'Ytick',[28:2:34],'YtickLabel',[28:2:34],'FontSize',12)
   set(gca,'Xtick',[-81:1:-77],'XtickLabel',[-81:1:-77],'FontSize',12)
   line(LON,LAT,'Marker','.','LineStyle','none','MarkerSize',12,'Color','k')
   hcb=colorbar('horiz');
   set(gca,'XLim',[-81.6,-77],'YLim',[27.87,34.11]);
   set(hcb,'Position',[0.3,0.09,0.435,0.025],'FontSize',12);
   set(ha,'Position',[0.13 0.19 0.775 0.815]);
   set(ha,'LineWidth',3,'TickDir','out')
   h=axes('Position',[0 0 1 1],'Visible','off');
   text(.35,.93,cal.months{i},'FontSize',14,'FontWeight','bold')
   eval(['print -depsc -painters ',cal.months{i},'_bott_temp.epsc']);
   eval(['!convert ',cal.months{i},'_bott_temp.epsc ',cal.months{i},'_bott_temp.gif ']);
end

close all