Changeset 8744

Timestamp:

06/28/11 14:33:36 (14 years ago)

Author:

jschierm

Message:

expread.m: Added closed structure field, and only delete last point if closed.

File:

• issm/trunk/src/m/utils/Exp/readwrite/expread.m (modified) (2 diffs)

issm/trunk/src/m/utils/Exp/readwrite/expread.m ¶

@@ -57 +57 @@
         %Get number of nods and density
    A=fscanf(fid,'%f %f',[1 2]);
-   if whole==1,
-       Struct(count).nods=A(1);
-   else
-       Struct(count).nods=A(1)-1;
-   end
+   Struct(count).nods=A(1);
    Struct(count).density=A(2);
 
@@ -75 +71 @@
         if(Struct(count).nods~=length(Struct(count).x))error(['Profile ' num2str(count) ' reports incorrect length']); end;
 
-        %skip last coordinate if whole=0,
-   if whole==0, fscanf(fid,'%f %f',[1 2]); end
+        %Check if closed
+        if (Struct(count).nods > 1) && ...
+           (Struct(count).x(end) == Struct(count).x(1)) && ...
+           (Struct(count).y(end) == Struct(count).y(1))
+                Struct(count).closed=true;
+                %skip last coordinate if whole=0,
+                if whole==0
+                        Struct(count).nods=Struct(count).nods-1;
+                        Struct(count).x   =Struct(count).x(1:end-1,1);
+                        Struct(count).y   =Struct(count).y(1:end-1,1);
+                end
+        else
+                Struct(count).closed=false;
+        end
 
 end