mlosch/interp_llc/gen_masks.m

% This is a script that executes the commands for
% creating a "PMASK" file (notation of MITgcm Classic)
% which is a 3D field containing 1 for ocean and 0 for
% land. This mask is *NOT* necessarily the mask used
% in the model but is representative of the largest
% possible mask that the model might use (cf. hFacMin).
% We need it here in order to guarantee that there
% is valid data in the climatological input data files.
%
% Uses data files: bathymetry.bin HGRID.mat VGRID.mat
% Creates data files: pmask.bin
% Creates arrays: n/a
% Uses m-scripts: hfac
%
% Created  11/11/99 by adcroft@mit.edu
% Modified          by
% Maintained by adcroft@mit.edu, abiastoch@ucsd.edu

clear all
format compact
more off

load FMT.mat
load HGRID.mat nxc nyc
load VGRID.mat dz

load HN
%Hn=rdda('bathymetry.bin',size(hn),1,fmt,Ieee);
%if max(abs(hn(:)-Hn(:)))
% error('HN.mat and bathymetry.bin don''t match');
%end

% Generate mask and store in pmask.bin
[msk] = hfac(dz,hn,0,0);
msk(find(msk~=0))=1;

save MASKS msk
%wrda('pmask.bin',msk,1,fmt,Ieee);
1	% This is a script that executes the commands for
2	% creating a "PMASK" file (notation of MITgcm Classic)
3	% which is a 3D field containing 1 for ocean and 0 for
4	% land. This mask is NOT necessarily the mask used
5	% in the model but is representative of the largest
6	% possible mask that the model might use (cf. hFacMin).
7	% We need it here in order to guarantee that there
8	% is valid data in the climatological input data files.
9	%
10	% Uses data files: bathymetry.bin HGRID.mat VGRID.mat
11	% Creates data files: pmask.bin
12	% Creates arrays: n/a
13	% Uses m-scripts: hfac
14	%
15	% Created 11/11/99 by adcroft@mit.edu
16	% Modified by
17	% Maintained by adcroft@mit.edu, abiastoch@ucsd.edu
18
19	clear all
20	format compact
21	more off
22
23	load FMT.mat
24	load HGRID.mat nxc nyc
25	load VGRID.mat dz
26
27	load HN
28	%Hn=rdda('bathymetry.bin',size(hn),1,fmt,Ieee);
29	%if max(abs(hn(:)-Hn(:)))
30	% error('HN.mat and bathymetry.bin don''t match');
31	%end
32
33	% Generate mask and store in pmask.bin
34	[msk] = hfac(dz,hn,0,0);
35	msk(find(msk~=0))=1;
36
37	save MASKS msk
38	%wrda('pmask.bin',msk,1,fmt,Ieee);