% input files
mesh_file = 'my_mesh.mat'
tpxoh = 'h_**_tpxo9_atlas_30_v5.nc';
tpxou = 'u_**_tpxo9_atlas_30_v5.nc';

% load the mesh
load(mesh_file)

% add the sponge layer (width 10% of M2) 
sponge_opt = 1; % for tidal boundary conditions
m = Calc_Sponge(m,[12.42*3600, 0.1],sponge_opt);

%% Make the fort.15 struct
TS = start_time; %[date-string]
TE = end_time; %[date-string]
DT = time_step; %[s]
CONST = "major8"; % tidal constituent information
m = Make_f15(m,TS,TE,DT,'const',CONST,'tidal_database',tpxoh);

% add the sponge info for this tidal
m = Make_f5354(m, tpxoh, tpxou);
% write out mesh files..
write(m,'sponge_mesh')