Changeset 3824
- Timestamp:
- Oct 18, 2006, 1:52:52 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/data_manager.py
r3820 r3824 4224 4224 """ 4225 4225 columns = 3 # long, lat , depth 4226 #FIXME use mux_file, not f 4226 4227 f = open(file_in, 'rb') 4227 4228 … … 4273 4274 4274 4275 for i in range(time_step_count): 4275 #Read in a time slice 4276 #Read in a time slice from mux file 4276 4277 hz_p_array = p_array.array('f') 4277 4278 hz_p_array.read(f, points_num) … … 4279 4280 hz_p = reshape(hz_p, (len(lon), len(lat))) 4280 4281 hz_p = transpose(hz_p) #mux has lat varying fastest, nc has long v.f. 4281 4282 4283 #write time slice to nc file 4282 4284 nc_file.store_timestep(hz_p) 4283 4285 #FIXME should I close the mux file here? 4286 f.close() 4284 4287 nc_file.close() 4285 4288
Note: See TracChangeset
for help on using the changeset viewer.