Index: /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpShepherd2019.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpShepherd2019.m	(revision 24874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpShepherd2019.m	(revision 24874)
@@ -0,0 +1,53 @@
+function dhdt=interpShepherd2019(X,Y,string,varargin)
+%INTERPSHEPHERD2019 - interpolate Shepherd2019 data
+%
+%   Available data:
+%      1.  dhdt_1992_1996
+%      2.  dhdt_1997_2001
+%      3.  dhdt_2002_2006
+%      4.  dhdt_2007_2011
+%      5.  dhdt_2012_2016
+%      6.  dhdt_1992_2017
+%      7.  uncert_1992_1996
+%      8.  uncert_1997_2001
+%      9.  uncert_2002_2006
+%      10.  uncert_2007_2011
+%      11.  uncert_2012_2016
+%      12.  uncert_1992_2017
+%
+%   Usage:
+%      [dataout] = interpShepherd2019(X,Y,'dhdt_1992_2017')
+
+options={'dhdt_1992_1996','dhdt_1997_2001','dhdt_2002_2006','dhdt_2007_2011','dhdt_2012_2016','dhdt_1992_2017',...
+			'uncert_1992_1996','uncert_1997_2001','uncert_2002_2006','uncert_2007_2011','uncert_2012_2016','uncert_1992_2017'};
+tf=strcmp(string,options);
+
+if ~any(tf)
+	disp('String not available!');
+   disp('The options are:');
+   disp(options);
+   error('String not available. See message above.');
+end
+
+switch oshostname(),
+   case {'ronne'}
+		nc='/home/ModelData/Antarctica/DHDTShepherd/antarctic_dhdt_5km_grid_1992_2017.nc';
+	case {'recruta'}
+		nc='/home/santos/ModelData/CPOM_dhdt_shepherd_2019/antarctic_dhdt_5km_grid_1992_2017.nc';
+	otherwise
+      error('machine not supported yet');
+end
+
+if nargin==3,
+   method='linear';% default
+else
+   method=varargin{1};
+end
+
+xdata=double(ncread(nc,'x'));
+ydata=double(ncread(nc,'y'));
+data=double(ncread(nc,string))';
+
+dhdt=InterpFromGrid(xdata,ydata,data,X,Y,method);
+
+end
Index: /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpSmith2020.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpSmith2020.m	(revision 24874)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpSmith2020.m	(revision 24874)
@@ -0,0 +1,45 @@
+function dhdt=interpSmith2020(X,Y,string,varargin)
+%INTERPSMITH2020 - interpolate Smith2020 data
+%
+%	Available data:
+%		Filtered mass-change maps, for display only, units of m(ice-equivalent)/yr:
+%		ais_floating_filt.tif
+%		ais_grounded_filt.tif
+%		gris_filt.tif
+
+%		Raw mass-change maps, suitable for generation of basin-by-basin mass-change estimates, units of m(ice-equivalent)/yr:
+%		ais_floating.tif
+%		ais_grounded.tif
+%		gris.tif
+%
+%   Usage:
+%      [dataout] = interpSmith2020(X,Y,'ais_floating_filt')
+
+options={'ais_floating_filt','ais_grounded_filt','gris_filt','ais_floating','ais_grounded','gris'};
+tf=strcmp(string,options);
+
+if ~any(tf)
+	disp('String not available!');
+	disp('The options are:');
+	disp(options);
+	error('String not available. See message above.');
+end
+	
+switch oshostname(),
+	case {'ronne'}
+		if strcmp(string,'gris_filt') | strcmp(string,'gris')
+			path='/home/ModelData/Greeland/DHDTSmith/'
+		else
+			path='/home/ModelData/Antarctica/DHDTSmith/'
+		end
+	case {'recruta'}
+		path='/home/santos/ModelData/ICESat1_ICESat2_mass_change/';
+	otherwise
+		error('machine not supported yet');
+end
+
+file=strcat(path,string,'.tif');
+
+dhdt=interpFromGeotiff(file,X,Y);
+
+end