MITgcm_contrib/grid_bathymetry_matlab/get_etopo2.m

function [etopo2,lon,lat]=get_etopo2(lon,lat)
% Function [etopo2,lon,lat]=get_etopo2(lon,lat)
%   get etopo2 bathymetry
%   etopo2 is decimated on grid defined by input lon/lat;
%   if input lon/lat is not specified, full resolution etopo2
%   is returned. etopo2 bathymetry file obtained from
%   http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO2/ETOPO2.raw.bin
%
% INPUTS
% lon     longitude East (deg)
% lat     latitude North (deg)
%
% OUTPUTS
% lon     longitude East (deg)
% lat     latitude North (deg)
% etopo2  2-D array containing bathymetry (m)

elon=0:1/30:(360-1/30);
elat=(-90+1/30):1/30:90;

etopo2=rot90(readbin('ETOPO2.raw.bin',[10800 5400],1,'int16'))';
etopo2=etopo2([5401:10800 1:5400],:);

if nargin < 2
  lon=elon;
  lat=elat;
else
  etopo2=decim2(elon,elat,etopo2,lon,lat);
end
1	function [etopo2,lon,lat]=get_etopo2(lon,lat)
2	% Function [etopo2,lon,lat]=get_etopo2(lon,lat)
3	% get etopo2 bathymetry
4	% etopo2 is decimated on grid defined by input lon/lat;
5	% if input lon/lat is not specified, full resolution etopo2
6	% is returned. etopo2 bathymetry file obtained from
7	% http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO2/ETOPO2.raw.bin
8	%
9	% INPUTS
10	% lon longitude East (deg)
11	% lat latitude North (deg)
12	%
13	% OUTPUTS
14	% lon longitude East (deg)
15	% lat latitude North (deg)
16	% etopo2 2-D array containing bathymetry (m)
17
18	elon=0:1/30:(360-1/30);
19	elat=(-90+1/30):1/30:90;
20
21	etopo2=rot90(readbin('ETOPO2.raw.bin',[10800 5400],1,'int16'))';
22	etopo2=etopo2([5401:10800 1:5400],:);
23
24	if nargin < 2
25	lon=elon;
26	lat=elat;
27	else
28	etopo2=decim2(elon,elat,etopo2,lon,lat);
29	end