Changeset 12303

Timestamp:

05/30/12 13:44:17 (13 years ago)

Author:

Mathieu Morlighem

Message:

Better input output

File:

• issm/trunk-jpl/src/m/utils/DataProcessing/gamv.m (modified) (3 diffs)

issm/trunk-jpl/src/m/utils/DataProcessing/gamv.m

@@ -4 +4 @@
 %   Usage:
 %      output = gamv(x,y,data,varargin)
+
+options=pairoptions(varargin{:});
+
+nlag = getfieldvalue(options,'nlag', 20);
+dlag = getfieldvalue(options,'dlag', 1000);
 
 %Write data file
@@ -26 +31 @@
 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 +46 @@
 %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')