
% Read input file
inp = ncget('REF_20060116_18');

% Plot continental boundaries (rotated latitude/longitude)
figure;
inp.COAST_RLON.data ( inp.COAST_RLON.data< -900 ) = NaN;
inp.COAST_RLAT.data ( inp.COAST_RLAT.data< -900 ) = NaN;
plot(inp.COAST_RLON.data,inp.COAST_RLAT.data,'k');

% Plot continental boundaries (x/y)
figure;
inp.COAST_X.data ( inp.COAST_X.data< -5000 ) = NaN;
inp.COAST_Y.data ( inp.COAST_Y.data< -5000 ) = NaN;
plot(inp.COAST_X.data,inp.COAST_Y.data,'k');
