Index: /issm/trunk-jpl/src/m/modeldata/interpFromGeotiff.m
===================================================================
--- /issm/trunk-jpl/src/m/modeldata/interpFromGeotiff.m	(revision 27911)
+++ /issm/trunk-jpl/src/m/modeldata/interpFromGeotiff.m	(revision 27912)
@@ -1,6 +1,16 @@
-function dataout = interpFromGeotiff(geotiffname,X,Y,nanValue),
+function dataout = interpFromGeotiff(geotiffname,X,Y,nanValue,fillholes)
+%INTERPFROMGEOTIFF - interpolate field in geotiff onto list of points
+%
+%   Usage:
+%      dataout = interpFromGeotiff(geotiffname,X,Y,nanValue,fillholes)
+%      dataout = interpFromGeotiff(geotiffname,X,Y);
+
 
 if nargin < 4
 	nanValue = 10^30;
+	fillholes = false;
+end
+if nargin < 5
+	fillholes = false;
 end
 
@@ -64,5 +74,9 @@
 		data(find(data<=nanValue))=NaN;
 	end
-
+	if fillholes
+		disp('Filling holes');
+		data = inpaint_nans(data);
+		disp('done');
+	end
 end