clear;
prec='real*4';
ieee='b';
%% 4 degrees
% nx=90;
% ny=40;
% fid=fopen('bathymetry4.bin','r',ieee); 
%% 2 degrees
% nx=180;
% ny=80;
% fid=fopen('bathymetry2.bin','r',ieee);
%% 1 degrees
nx=360;
ny=180;
fid=fopen('bathy_1deg_pangea15z.bin','r',ieee);

F_num=fread(fid,[nx,ny],prec); 
%F_num=fread(fid,1000,prec);
fclose(fid);

%% 1 degrees
[x,y]=meshgrid(linspace(0,359,nx),linspace(-90,89,ny));
%% 2 degrees
%[x,y]=meshgrid(linspace(0,358,nx),linspace(-90,88,ny));
%%
figure
contourf(x,y,F_num')
colorbar