Ignore:
Timestamp:
08/25/20 00:32:13 (5 years ago)
Author:
jdquinn
Message:

CHG: Saving chnages so that Basile has access to potential fix to solidearthmodel class

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/basin.m

    r25065 r25455  
    117117                        %recover options
    118118                        options=pairoptions(varargin{:});
    119                         x=[]; y=[];
     119                        x=[];
     120                        y=[];
    120121
    121122                        %go through boundaries, recover edges and project them in the basin epsg reference frame:
     
    127128                                y=[y;contour.y];
    128129                        end
     130
    129131                        %close the contour:
    130132                        if x(end)~=x(1) | y(end)~=y(1),
     
    181183                                error('Basin shapefilecrop error message: epsgshapefile or projshapefile should be specified');
    182184                        end
    183 
    184185
    185186                        %create list of contours that have critical length > threshold (in lat,long):
     
    209210                        for i=1:length(contours),
    210211                                h=contours(i);
    211                                 isin=inpolygon(h.x,h.y,ba.x,ba.y); 
     212                                isin=inpolygon(h.x,h.y,ba.x,ba.y);
    212213                                if ~isempty(find(isin==0)),
    213214                                        flags(i)=1;
    214215                                end
    215216                        end
    216                         pos=find(flags);  contours(pos)=[];
     217
     218                        pos=find(flags);
     219                        contours(pos)=[];
    217220
    218221                        %Two options:
Note: See TracChangeset for help on using the changeset viewer.