Changeset 12302


Ignore:
Timestamp:
05/29/12 16:06:43 (13 years ago)
Author:
Mathieu Morlighem
Message:

Added varmap.m gamv.m

Location:
issm/trunk-jpl/src/m/utils/DataProcessing
Files:
2 added
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/utils/DataProcessing/gslib.m

    r12291 r12302  
    1 x=[];
    2 y=[];
    3 thickness=[];
    4 for i=1:5,
    5         switch(i)
    6                 case 1, filename='2002_Antarctica_P3chile';
    7                 case 2, filename='2004_Antarctica_P3chile';
    8                 case 3, filename='2009_Antarctica_DC8';
    9                 case 4, filename='2009_Antarctica_TO';
    10                 case 5, filename='2010_Antarctica_DC8';
    11                 otherwise, error('not supported');
    12         end
    13         disp(['Loading ' filename '...']);
    14         A=load(filename);
    15         [tempx tempy]=ll2xy(A.lat(:),A.lon(:),-1,0,71);
    16         x=[x;tempx];
    17         y=[y;tempy];
    18         thickness=[thickness;A.thickness];
    19 end
    20 [xl yl]=basinzoom(pairoptions('basin','Pine Island Glacier'));
    21 x=x(1:100:end);
    22 y=y(1:100:end);
    23 data=thickness(1:100:end);
     1function output = gslib(x,y,data,varargin);
     2%GSLIB - use gslib for Kriging
     3%
     4%   Usage:
     5%      output = gslib(x,y,data,varargin)
    246
    257%Output Matrix
    … …  
    4426c = (sill-nugget);
    4527a = sqrt(3)*range;
    46 
    4728
    4829%Write data file
    … …  
    10889
    10990        tic;system([issmdir() '/externalpackages/gslib/install/kb2d kb2d.par']);toc;
     91        delete('kb2d.par');
    11092
    11193        %Read output
Note: See TracChangeset for help on using the changeset viewer.