Changeset 4030
- Timestamp:
- Nov 24, 2006, 9:18:19 AM (17 years ago)
- Files:
-
- 4 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2006/make_report.py
r4029 r4030 73 73 # WA DLI data 74 74 production_dirs = {#'20061008_234702': 'Mw 8-7', # 2500 res 75 #'20061004_230438': 'Mw 8-7' ,# 7500 res75 #'20061004_230438': 'Mw 8-7'} # 7500 res 76 76 '20061022_224422': 'Mw 8-7'} # refined mesh 77 77 #'20061006_062319': 'Mw 8-5'} # 2500 res … … 99 99 #reportname = 'latexoutput', 100 100 plot_quantity = ['stage', 'speed'], 101 surface = True,101 #surface = True, 102 102 time_min = None, 103 103 time_max = None, -
anuga_work/production/hobart_2006/project.py
r4029 r4030 79 79 #gauge_filename = gaugedir + 'hobart_gauges_final.csv' 80 80 #gauge_filename = gaugedir + 'hobartGauges20061019.csv' 81 gauge_filename = gaugedir + 'riskAreas231106.csv' 81 #gauge_filename = gaugedir + 'riskAreas231106.csv' 82 gauge_filename = gaugedir + 'hobart_gauges2.csv' 82 83 #gauge_filename = gaugedir + 'gauges_MRT.csv' 83 84 #gauge_filename_bom = gaugedir + 'gauges_bruny_fine.csv' -
anuga_work/production/sydney_2006/smf_volume_conservation.py
r3190 r4030 45 45 #axis([0, 5, -2, 2]) 46 46 47 figure( 1)47 figure(21) 48 48 plot(X/lam0,eta1,'b.-',X/lam0,eta2,'g.-') 49 49 title('Surface elevation of SMF: kappa = 0.83 and 1') … … 51 51 ylabel('eta(x,y=0)') 52 52 axis([0, 5, -2, 2]) 53 savefig('check_fig2') 53 54 54 55 def set_integral_zero(x,dx): … … 76 77 #print 'dx: %0.5f kappad: %f' % (dx, kappad1) 77 78 78 figure(2 )79 figure(22) 79 80 plot(dx_vec1,kappad_vec,'.-') 80 81 title('Volume conservation') … … 97 98 98 99 99 figure( 3)100 figure(23) 100 101 plot(x0_vec2,kappad_vec2,'.-') 101 102 title('Volume conservation - using xg')
Note: See TracChangeset
for help on using the changeset viewer.