Changeset 18382


Ignore:
Timestamp:
08/13/14 13:47:29 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: more flexible call to gslib

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified issm/trunk-jpl/src/m/contrib/gslib/gslib.m

    r17677 r18382  
    1 function [B E] = gslib(x,y,data,xmin,ymin,nx,ny,deltax,deltay)
     1function [B E] = gslib(x,y,data,xmin,ymin,nx,ny,deltax,deltay,varargin)
    22%GSLIB - use gslib for Kriging
    33%
     
    55%      output = gslib(x,y,data,varargin)
    66
     7%process options
     8options = pairoptions(varargin{:});
     9
    710%Variogram
    8 nugget=10;
    9 sill  =164;
    10 range =25763;
     11nugget= getfieldvalue(options,'nugget',10);
     12sill  = getfieldvalue(options,'sill',164);
     13range = getfieldvalue(options,'range',25763);
    1114
    1215%Kriging options
    13 mindata = 1;
    14 maxdata = 50;
    15 maxsearchradius = 50000;
     16mindata = getfieldvalue(options,'mindata',1);
     17maxdata = getfieldvalue(options,'maxdata',50);
     18maxsearchradius = getfieldvalue(options,'searchrange',50000);
    1619
    1720%Some intermediaries (Convert to gslib's parameters);
Note: See TracChangeset for help on using the changeset viewer.