function plot_singv_domain(outdir)

% Load shoreline data
load /nfs/see-archive-10_a26/lecimb/cathryn/gshhs/shoreline.mat
%load /nfs/see-archive-10_a26/lecimb/cathryn/gshhs/country_boundaries.mat
%cmap = cptcmap('/nfs/see-fs-02_users/earceb/matlab/cptfiles/WhBlGrYeRe.cpt');

% Read in netcdf data
ht = nc_varget('/nfs/a161/earceb/Singapore/MJO_dc/ancils/km2pt2/qrparm.orog.globe.nc','ht');
x = nc_varget('/nfs/a161/earceb/Singapore/MJO_dc/ancils/km2pt2/qrparm.orog.globe.nc','x');
y = nc_varget('/nfs/a161/earceb/Singapore/MJO_dc/ancils/km2pt2/qrparm.orog.globe.nc','y');

% Plot
contourf(x,y,ht,[0:100:2500],'LineColor','None')
caxis([0 2500])
%colormap(cmap)
colorbar
hold on
plot(shoreline(:,1),shoreline(:,2),'k-','LineWidth',1.5) 
plot([96 112],[-8 7],'k-','LineWidth',2)

plot([101.5 101.5],[-3 -1],'k-','LineWidth',2)
plot([104.5 104.5],[-3 -1],'k-','LineWidth',2)
plot([101.5 104.5],[-3 -3],'k-','LineWidth',2)
plot([101.5 104.5],[-1 -1],'k-','LineWidth',2)

plot([95 95],[-6 -1],'k-','LineWidth',2)
plot([100 100],[-6 -1],'k-','LineWidth',2)
plot([95 100],[-6 -6],'k-','LineWidth',2)
plot([95 100],[-1 -1],'k-','LineWidth',2)

set(gca,'FontSize',12,'LineWidth',1.5,'TickDir','Out')
xlim([93 112])
ylim([-8 8])
ylabel('Latitude')
xlabel('Longitude') 

printname = ['SINGV_domain_2pt2km_subdomains.png'];
eval(['print ''',outdir,'/',printname,''' -dpng -r200']);
close