Changeset 3830
- Timestamp:
- Oct 19, 2006, 6:35:45 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/data_manager.py
r3826 r3830 4250 4250 4251 4251 lonlatdep = p_array.array('f') 4252 points_num = 29144252 ####points_num = 2914 4253 4253 lonlatdep.read(f, columns * points_num) 4254 4254 lonlatdep = array(lonlatdep, typecode=Float) … … 4258 4258 4259 4259 lon, lat, depth = lon_lat2grid(lonlatdep) 4260 lon_sorted = l on[:]4260 lon_sorted = list(lon) 4261 4261 lon_sorted.sort() 4262 4263 # if not lon == lon_sorted: 4262 4264 4263 if not allclose(lon, lon_sorted): 4265 4264 msg = "Longitudes in mux file are not in ascending order" 4266 4265 raise IOError, msg 4267 lat_sorted = l at[:]4268 lat_sorted.sort() 4269 # if not lat == lat_sorted: 4266 lat_sorted = list(lat) 4267 lat_sorted.sort() 4268 4270 4269 if not allclose(lat, lat_sorted): 4271 4270 msg = "Latitudes in mux file are not in ascending order" … … 4353 4352 LAT = 1 4354 4353 QUANTITY = 2 4355 4356 # print 'long_lat_dep', type(long_lat_dep), long_lat_dep.shape 4357 # print 'long_lat_dep', long_lat_dep 4354 4355 long_lat_dep = ensure_numeric(long_lat_dep, Float) 4358 4356 4359 4357 num_points = long_lat_dep.shape[0]
Note: See TracChangeset
for help on using the changeset viewer.