Changeset 5602 for anuga_core/source/anuga/shallow_water/data_manager.py
- Timestamp:
- Aug 4, 2008, 1:41:59 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/shallow_water/data_manager.py
r5589 r5602 4926 4926 assert nt>0,msg 4927 4927 4928 OFFSET=5 #number of site parameters p passed back with data 4929 #p=[geolat,geolon,depth,start_tstep,finish_tstep] 4930 4931 times=dt*arange(1,(data.shape[1]-OFFSET)+1) 4932 latitudes=zeros(data.shape[0],Float) 4933 longitudes=zeros(data.shape[0],Float) 4934 elevation=zeros(data.shape[0],Float) 4935 quantity=zeros((data.shape[0],data.shape[1]-OFFSET),Float) 4928 OFFSET=5 # Number of site parameters p passed back with data 4929 # p=[geolat,geolon,depth,start_tstep,finish_tstep] 4930 4931 parameters_index = data.shape[1]-OFFSET 4932 4933 times=dt*arange(parameters_index) 4934 latitudes=zeros(data.shape[0], Float) 4935 longitudes=zeros(data.shape[0], Float) 4936 elevation=zeros(data.shape[0], Float) 4937 quantity=zeros((data.shape[0], parameters_index), Float) 4938 4939 4936 4940 starttime=1e16 4937 for i in range(0,data.shape[0]): 4938 latitudes[i]=data[i][data.shape[1]-OFFSET] 4939 longitudes[i]=data[i][data.shape[1]-OFFSET+1] 4940 elevation[i]=-data[i][data.shape[1]-OFFSET+2] 4941 quantity[i]=data[i][:-OFFSET] 4942 starttime=min(dt*data[i][data.shape[1]-OFFSET+3],starttime) 4941 for i in range(0, data.shape[0]): 4942 latitudes[i]=data[i][parameters_index] 4943 longitudes[i]=data[i][parameters_index+1] 4944 elevation[i]=-data[i][parameters_index+2] 4945 quantity[i][:]=data[i][:parameters_index] # Was data[i][:-OFFSET] 4946 4947 starttime=min(dt*data[i][parameters_index+3],starttime) 4943 4948 4944 4949 return times, latitudes, longitudes, elevation, quantity, starttime 4950 4945 4951 4946 4952 def mux2sww_time(mux_times, mint, maxt):
Note: See TracChangeset
for help on using the changeset viewer.