Index: /issm/trunk-jpl/src/m/utils/DataProcessing/gamv.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/DataProcessing/gamv.m	(revision 12302)
+++ /issm/trunk-jpl/src/m/utils/DataProcessing/gamv.m	(revision 12303)
@@ -4,4 +4,9 @@
 %   Usage:
 %      output = gamv(x,y,data,varargin)
+
+options=pairoptions(varargin{:});
+
+nlag = getfieldvalue(options,'nlag', 20);
+dlag = getfieldvalue(options,'dlag', 1000);
 
 %Write data file
@@ -26,7 +31,7 @@
 fprintf(fid,'%-30s %s\n','-1.0e21 1.0e21'             ,'\trimming limits');
 fprintf(fid,'%-30s %s\n','gamv.out'                   ,'\file for variogram output');
-fprintf(fid,'%-30s %s\n','20'                         ,'\number of lags');
-fprintf(fid,'%-30s %s\n','5.0'                        ,'\lag separation distance');
-fprintf(fid,'%-30s %s\n','3.0'                        ,'\lag tolerance');
+fprintf(fid,'%-30s %s\n',num2str(nlag,'%i')           ,'\number of lags');
+fprintf(fid,'%-30s %s\n',num2str(dlag,'%g')           ,'\lag separation distance');
+fprintf(fid,'%-30s %s\n',num2str(dlag/2,'%g')         ,'\lag tolerance');
 fprintf(fid,'%-30s %s\n','3'                          ,'\number of directions');
 fprintf(fid,'%-30s %s\n','0.0 90.0 50.0 0.0 90.0 50.0','\azm, atol, bandh, dip, dtol, bandv');
@@ -41,2 +46,25 @@
 %Call gamv
 system([issmdir() '/externalpackages/gslib/install/gamv gamv.par']);
+delete('gamv.par');
+
+%Read output
+output   = struct('Semivariogram',[],'Covariance',[]);
+counter1 = 1;
+counter2 = 1;
+fid=fopen('gamv.out','r');
+while (~feof(fid)),
+	A=fscanf(fid,'%s',1);
+	if strcmp(A,'Covariance');
+		A=fscanf(fid,'%s',4); %Read tail:Data head:Data direction  2
+		output(counter1).Covariance=fscanf(fid,'%i %g %g %i %g %g',[6 nlag+2])';
+		counter1=counter1+1;
+	elseif strcmp(A,'Semivariogram'),
+		A=fscanf(fid,'%s',4); %Read tail:Data head:Data direction  2
+		output(counter2).Semivariogram=fscanf(fid,'%i %g %g %i %g %g',[6 nlag+2])';
+		counter2=counter2+1;
+	else
+		%do nothing
+	end
+end
+fclose(fid);
+delete('gamv.out')