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