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);