Changeset 24897
- Timestamp:
- 05/22/20 17:21:51 (5 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/mesh3dsurface.m
r24894 r24897 213 213 end 214 214 elseif strcmpi(geometry,'line'), 215 count er=1;215 count=1; 216 216 for i=1:self.numberofelements, 217 217 el=self.elements(i,:); 218 218 %first line: 219 contours(count er).x = [self.long(el(1)) self.long(el(2))];220 contours(count er).y = [self.lat(el(1)) self.lat(el(2))];221 contours(count er).Geometry = 'Line';219 contours(count).x = [self.long(el(1)) self.long(el(2))]; 220 contours(count).y = [self.lat(el(1)) self.lat(el(2))]; 221 contours(count).Geometry = 'Line'; 222 222 223 223 %second line: 224 contours(count er+1).x = [self.long(el(2)) self.long(el(3))];225 contours(count er+1).y = [self.lat(el(2)) self.lat(el(3))];226 contours(count er+1).Geometry = 'Line';224 contours(count+1).x = [self.long(el(2)) self.long(el(3))]; 225 contours(count+1).y = [self.lat(el(2)) self.lat(el(3))]; 226 contours(count+1).Geometry = 'Line'; 227 227 228 228 %third line: 229 contours(count er+2).x = [self.long(el(3)) self.long(el(1))];230 contours(count er+2).y = [self.lat(el(3)) self.lat(el(1))];231 contours(count er+2).Geometry = 'Line';229 contours(count+2).x = [self.long(el(3)) self.long(el(1))]; 230 contours(count+2).y = [self.lat(el(3)) self.lat(el(1))]; 231 contours(count+2).Geometry = 'Line'; 232 232 233 %increase count er:234 count er+= 3;233 %increase count: 234 count += 3; 235 235 end 236 236 elseif strcmpi(geometry,'polygon'), -
issm/trunk-jpl/src/m/classes/mesh3dsurface.py
r24896 r24897 218 218 contours[i] = contour 219 219 elif geometry == 'line': 220 count er= 0220 count = 0 221 221 for i in range(len(self.numberofelements)): 222 222 el = self.elements[i] … … 227 227 contour.y = [self.lat[el[0]], self.lat[el[1]]] 228 228 contour.Geometry = 'Line' 229 contours[count er] = contour229 contours[count] = contour 230 230 231 231 #second line: … … 234 234 contour.y = [self.lat[el[1]], self.lat[el[2]]] 235 235 contour.Geometry = 'Line' 236 contours[count er+ 1] = contour236 contours[count + 1] = contour 237 237 238 238 #second line: … … 241 241 contour.y = [self.lat[el[2]], self.lat[el[0]]] 242 242 contour.Geometry = 'Line' 243 contours[count er+ 2] = contour244 245 #increase count er:246 count er+= 3243 contours[count + 2] = contour 244 245 #increase count: 246 count += 3 247 247 elif geometry == 'polygon': 248 248 # TODO: Refactor the following to reduce repeated code, or
Note:
See TracChangeset
for help on using the changeset viewer.