== Matlab file ==

%Paths to datasets
modeldatapath='/u/astrid-r1b/larour/ModelData';
thicknesspath=[modeldatapath '/ASEThickness/ASE_thickness'];
surfacepath=[modeldatapath '/BamberDEMAntarctica1km/surface_smooth30_lowslope'];
firnpath=[modeldatapath '/BroekeFirn1km/firn.mat'];

%Read Surface and interpolate onto model
md.surface=InterpFromFile(md.x,md.y,surfacepath,10);

%Correct surface using firn layerfrom Broeke
md.firn_layer=InterpFromFile(md.x,md.y,firnpath,0);
rho_ice=917;
rho_firn=830;
firn_layer_correction=md.firn_layer*(1-rho_firn/rho_ice);
md.surface=md.surface-firn_layer_correction;

%Load thickness and interpolate onto model
load(thicknesspath);
md.thickness=InterpFromGridToMesh(x_m,y_m,thickness,md.x,md.y,10);

%Correct thickness using Bamber's dem
md=ThicknessCorrection(md,15000);