Changes between Initial Version and Version 1 of shelf_thickness

Timestamp:

10/12/10 14:32:44 (15 years ago)

Author:

Mathieu Morlighem

Comment:

--

shelf_thickness

@@ +1 @@
+== 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);