Index: ../trunk-jpl/src/m/contrib/gslib/gslib.m
===================================================================
--- ../trunk-jpl/src/m/contrib/gslib/gslib.m	(revision 18381)
+++ ../trunk-jpl/src/m/contrib/gslib/gslib.m	(revision 18382)
@@ -1,18 +1,21 @@
-function [B E] = gslib(x,y,data,xmin,ymin,nx,ny,deltax,deltay)
+function [B E] = gslib(x,y,data,xmin,ymin,nx,ny,deltax,deltay,varargin)
 %GSLIB - use gslib for Kriging
 %
 %   Usage:
 %      output = gslib(x,y,data,varargin)
 
+%process options
+options = pairoptions(varargin{:});
+
 %Variogram
-nugget=10;
-sill  =164;
-range =25763;
+nugget= getfieldvalue(options,'nugget',10);
+sill  = getfieldvalue(options,'sill',164);
+range = getfieldvalue(options,'range',25763);
 
 %Kriging options
-mindata = 1;
-maxdata = 50;
-maxsearchradius = 50000;
+mindata = getfieldvalue(options,'mindata',1);
+maxdata = getfieldvalue(options,'maxdata',50);
+maxsearchradius = getfieldvalue(options,'searchrange',50000);
 
 %Some intermediaries (Convert to gslib's parameters);
 c = (sill-nugget);
