@@ -13,7 +13,7 @@ function segmentstobasin(basin,varargin)
1313% See also EXPREAD
1414
1515 % some checks
16- if exist(basin ),
16+ if exist(basin )
1717 % choice=input(['A file ' basin ' already exists, do you want to modify it? (y/n)'],'s');
1818 % if ~strcmpi(choice,'y'),
1919 % disp('no modification done ... exiting');
@@ -22,15 +22,16 @@ function segmentstobasin(basin,varargin)
2222 end
2323
2424 % go through the list of basins
25- if mod(length(varargin ),2 )~=0 ,
25+ if mod(length(varargin ),2 )~=0
2626 error(' an even number of arguments should be provided after the basin name' );
2727 end
2828
2929 domain.x= []; domain.y= []; domain.nods= 1 ;
30- for i= 1 : nargin / 2 ,
30+ for i= 1 : nargin / 2
3131 expfile= varargin{(i - 1 )*2 + 1 };
3232 invert= varargin{(i - 1 )*2 + 2 };
33- if isexp(expfile ),
33+ [path ,name ,ext ]=fileparts(expfile );
34+ if strcmpi(ext ,' .exp' ),
3435 expstruct= expread(expfile ,' invert' ,invert );
3536 else
3637 expstruct= shpread(expfile ,' invert' ,invert );
@@ -45,7 +46,8 @@ function segmentstobasin(basin,varargin)
4546 domain.y= [domain .y ;domain .y(1 )];
4647 domain.Geometry= ' Polygon' ;
4748
48- if isexp(basin ),
49+ [path ,name ,ext ]=fileparts(basin );
50+ if strcmpi(ext ,' .exp' )
4951 expwrite(domain ,basin );
5052 else
5153 shpwrite(domain ,basin );
0 commit comments