-
Notifications
You must be signed in to change notification settings - Fork 61
/
Copy pathm_gshhs_h.m
41 lines (32 loc) · 1.22 KB
/
m_gshhs_h.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
function m_gshhs_h(varargin);
% M_GSHHS_H Add a coastline to a given map using the 'high' resolution of
% the Global Self-consistant Hierarchical High-resolution
% Shorelines.
%
% M_GSHHS_H( (standard line option,...,...) ) draws the coastline
% as a simple line.
% M_GSHHS_H('patch' ( ,standard patch options,...,...) ) draws the
% coastline as a number of patches.
%
% M_GSHHS_H('save',FILENAME) saves the extracted coastline data
% for the current projection in a file FILENAME. This allows
% speedier replotting using M_USERCOAST(FILENAME).
%
% See also M_PROJ, M_GRID, M_COAST, M_GSHHS_L, M_GSHHS_I, M_GSHHS_C
% M_USERCOAST
% Rich Pawlowicz (rich@ocgy.ubc.ca) 15/June/98
%
%
% This software is provided "as is" without warranty of any kind. But
% it's mine, so you can't sell it.
FILNAME='private/gshhs_h.b';
% Set current projection to geographic
Currentmap=m_coord('set');
m_coord('geographic');
if length(varargin)>1 & strcmp(varargin{1},'save'),
[ncst,Area,k]=mu_coast('h',FILNAME);
eval(['save ' varargin{2} ' ncst k Area']);
else
mu_coast('h',FILNAME,varargin{:},'tag','m_gshhs_h');
end;
m_coord(Currentmap.name);