%plot ROMS angular cross section from NJ coast to shelf break
%to compare to POM cross section

clear all; close all; %pack

%load('eastcoastmask_new.mat') %coastline

load('roms_jview_temp_greg.mat') 
romslat = lat_i(1,:);
romslon = lon_i(1,:);
%[romsxgrid,romsygrid]=meshgrid(romslon,romslat);


load bathymetry_USeastcoast
baddepth = find(depthi == 999999);
if isempty(baddepth) ~= 1
    depthi(baddepth) = NaN;
end

depthnew = interp2(loni,lati,-depthi,romslon,romslat);

close all
figure;
set(gcf,'Paperposition',[0 0 6 6]);
%pcolor(xs(125:140),zs,squeeze(Tdat(K,:,175,125:140))); shading flat; caxis([9 24]); %xslice
pcolor(lon_i,-dep_it,squeeze(tempi_dep)); shading flat; caxis([9 30]); %jslice
set(gca,'YDir','reverse')

hold on

plot(romslon,depthnew,'k','Linewidth',2)

xlim([-74.75 -72])
ylim([0 200])
%xlim([-74.75 -71])
%ylim([0 1000])
%ylim([0 4000])
%xlim([-74.75 -71.5])
colormap(jet); hc=colorbar; set(hc,'tickdir','out');
grid

set(gca,'tickdir','out','fontsize',12);

ttl{1}=['Hurricane irene09 ROMS Simulation: Initial time: 2011082700Z'];
ttl{2}=['T:     Cycle000 run'];

title(ttl,'fontsize',12);

%text(0.01,0.1,'24-hr (2011/08/26 00Z)',...
%    'fontsize',12,'fontweight','bold','color','k','units','normalized');
%text(0.01,0.05,['Xslice at lat=',num2str(ys(175))],...
%    'fontsize',12,'fontweight','bold','color','k','units','normalized');

wysiwyg
eval('print -dpng /Users/seroka/Research/data/HWRF/ROMS_T_xslice_cyc000_11082700_30C.png');
