Index: /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMartos2017.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMartos2017.m	(revision 25269)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMartos2017.m	(revision 25269)
@@ -0,0 +1,18 @@
+function out = interpMartos2017(X,Y)
+%INTERPMARTOS2017 - interpolate geothermal heat flux
+%
+%   Usage:
+%      out = interpMartos2017(X,Y)
+
+switch oshostname(),
+	case {'ronne'}
+		gtfpath='/home/ModelData/Antarctica/GeothermalMartos/Antarctic_GHF.xyz';
+	otherwise
+		error('machine not supported yet');
+end
+
+%Load data
+data = load(gtfpath);
+
+%Interpolate using nearest neighbor (dataset stops at ocean boundary!)
+out = Kriging(data(:,1),data(:,2),data(:,3),X,Y,'output','nearestneighbor')/1e3; %from mW/m2 to W/m2
Index: /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMouginotAnt2019.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMouginotAnt2019.m	(revision 25268)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpMouginotAnt2019.m	(revision 25269)
@@ -24,5 +24,5 @@
 id2y=min(numel(ydata),posy(end)+offset);
 
-disp(['   -- Mouginot 2017: loading velocities']);
+disp(['   -- Mouginot 2019: loading velocities']);
 vxdata = double(ncread(nc,'VX',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
 vydata = double(ncread(nc,'VY',[id1x id1y],[id2x-id1x+1 id2y-id1y+1],[1 1]))';
@@ -30,5 +30,5 @@
 ydata=ydata(id1y:id2y);
 
-disp(['   -- Mouginot 2017: interpolating ']);
+disp(['   -- Mouginot 2019: interpolating ']);
 vxout = InterpFromGrid(xdata,ydata,vxdata,double(X),double(Y));
 vyout = InterpFromGrid(xdata,ydata,vydata,double(X),double(Y));