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	mlosch	1.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);