Changeset 4949
- Timestamp:
- Jan 17, 2008, 11:49:00 AM (17 years ago)
- Files:
-
- 6 added
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/documentation/user_manual/anuga_user_manual.tex
r4948 r4949 1328 1328 \file{GetTimeseries.py} which generates the csv files for each point location. The csv files 1329 1329 can then be used in \code{csv2timeseries_graphs} to create the timeseries plot for each desired 1330 quantity. \code{csv2timeseries_graphs} relies on \code{pylab} to be install which is not part1330 quantity. \code{csv2timeseries_graphs} relies on \code{pylab} to be installed which is not part 1331 1331 of the standard \code{anuga} release, however it can be downloaded and installed from \code{http://matplotlib.sourceforge.net/} 1332 1332 … … 1340 1340 Figure \ref{fig:reef} shows the time series for the quantity stage for the 1341 1341 Elford Reef location for each scenario (the elevation at this location is negative, 1342 therefore stage is the mo st appropriate to use). Note the large negative stage value when the slide was1342 therefore stage is the more appropriate quantity to plot). Note the large negative stage value when the slide was 1343 1343 introduced. This is due to the double gaussian form of the initial surface 1344 1344 displacement of the slide. By contrast, the time series for depth is shown for the onshore location of the Cairns -
anuga_core/documentation/user_manual/demos/cairns/GetTimeseries.py
r4948 r4949 43 43 verbose=True) 44 44 except ImportError: 45 #ANUGA don't need pylab to work so the system doesn't 46 #rely on pylab being installed 47 print 'must have pylab install to generate plots' 45 #ANUGA does not rely on pylab to work 46 print 'must have pylab installed to generate plots' 48 47 49 48 -
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r4936 r4949 1235 1235 units = 'm / sec' 1236 1236 if which_quantity == 'bearing': 1237 #due_east = 90.0*ones(shape(model_time[0:n[j]-1,k,j],Float))1238 #due_west = 270.0*ones(shape(model_time[0:n[j]-1,k,j],Float))1239 print 'hello', bearings1240 print 'east', due_east1241 1237 plot(model_time[0:n[j]-1,k,j], bearings[0:n[j]-1,k,j], '-', 1242 1238 model_time[0:n[j]-1,k,j], due_west[0:n[j]-1], '-.', -
anuga_work/production/pt_hedland_2006/convert2latlon.py
r4927 r4949 31 31 32 32 timedir='20070613_044204_run_final_3.6_exmouth_nbartzis' 33 filename = project.outputdir + timedir + sep + 'stage0_interp.csv' 34 #filename = project.outputdir + timedir + sep + 'stage5_interp.csv' 33 #timedir='20070613_061053_run_final_3.6_onslow_nbartzis' 34 35 #filename = project.outputdir + timedir + sep + 'stage0_interp.csv' 36 filename = project.outputdir + timedir + sep + 'stage5_interp.csv' 35 37 #filename = project.outputdir+ timedir + sep + 'stage20_interp.csv' 36 #filename = project.outputdir+ timedir + sep + 'stage 80_interp.csv'37 filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'38 #filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'38 #filename = project.outputdir+ timedir + sep + 'stage50_interp.csv' 39 #filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv' 40 filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv' 39 41 #filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv' 40 #filename_out = project.outputdir + timedir + sep + 'stage 80_convert.csv'42 #filename_out = project.outputdir + timedir + sep + 'stage50_convert.csv' 41 43 alter_file(filename,filename_out) -
anuga_work/production/pt_hedland_2006/plot_contour_data.py
r4927 r4949 8 8 directory = project.outputdir 9 9 timedir='20070613_044204_run_final_3.6_exmouth_nbartzis' 10 #timedir='20070613_061053_run_final_3.6_onslow_nbartzis' 10 11 11 12 def get_data(filename): … … 19 20 return max(stage), min(stage) 20 21 21 files = ['stage0_convert.csv','stage5_convert.csv','stage20_convert.csv','stage80_convert.csv'] 22 depth = [0,5,20,80] 22 files = ['stage0_convert.csv','stage5_convert.csv','stage20_convert.csv','stage50_convert.csv'] 23 #files = ['stage0_convert.csv','stage20_convert.csv','stage50_convert.csv'] 24 depth = [0,5,20,50] 25 #depth = [0,20,50] 23 26 stage_max = [] 24 27 stage_min = [] … … 33 36 from anuga.abstract_2d_finite_volumes.util import greens_law 34 37 from Numeric import arange 35 d1 = 80.38 d1 = 50. 36 39 d2 = arange(d1,0.15,-0.1) 37 40 h1 = max_stage … … 53 56 #Magnitude 9.0 event (Sumba) at Pt Hedland') 54 57 legend(['ANUGA','Green\'s law']) 55 axis([-5, 85,min(stage_max)*0.9,max(stage_max)*1.1])58 axis([-5,55,min(stage_max)*0.9,max(stage_max)*1.1]) 56 59 grid(True) 57 60 #savefig('stratification_onslow_gun')
Note: See TracChangeset
for help on using the changeset viewer.