Changeset 4748 for anuga_validation/UQ_runup_2006/run_dam.py
- Timestamp:
- Oct 23, 2007, 2:49:17 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_validation/UQ_runup_2006/run_dam.py
r4678 r4748 31 31 32 32 slopes = [0] 33 frictions = [0.0 05, 0.01, 0.03] #can enter multiple frictions eg [0.00, 0.01]33 frictions = [0.01] #can enter multiple frictions eg [0.00, 0.01] 34 34 inital_depths = [0.2] #can enter multiple initial depths eg [0.1, 0.2] 35 35 gate_positions = [0.75] 36 36 flume_lengths = [3] #can enter multiple flume lengths eg [3.0, 5.0] 37 # NOTE flume_lengths is not used. It's hardcoded into create_mesh 37 38 38 39 for slope in slopes: … … 61 62 project.outputtimedir + 'create_mesh.py') 62 63 print'output dir', project.outputtimedir 63 64 #FIXME this isn't working65 #normal screen output is stored in66 screen_output_name = project.outputtimedir + "screen_output.txt"67 screen_error_name = project.outputtimedir + "screen_error.txt"68 64 69 65 #------------------------------------------------------------------------- … … 126 122 print 'That took %.2f seconds' %(time.time()-t0) 127 123 print 'finished' 128 129 points = [[0.45, 0.2]] 130 131 124 125 points = [[gate_position - 0.65,0.2], 126 [gate_position - 0.55,0.2], 127 [gate_position - 0.45,0.2], 128 [gate_position - 0.35,0.2], 129 [0.45,0.2], 130 [gate_position - 0.25,0.2] 131 ] 132 132 #------------------------------------------------------------------------- 133 133 # Calculate gauge info 134 134 #------------------------------------------------------------------------- 135 #interpolate_sww2csv(project.outputtimedir + project.basename \ 136 # + id+".sww", 137 # points, 138 # project.depth_filename + id + '.csv', 139 # project.velocity_x_filename + id + '.csv', 140 # project.velocity_y_filename + id + '.csv') 135 print "name", project.outputtimedir + project.basename \ 136 + id+".sww" 137 interpolate_sww2csv(project.outputtimedir + project.basename \ 138 + id+".sww", 139 points, 140 project.depth_filename + id + '.csv', 141 project.velocity_x_filename + id + '.csv', 142 project.velocity_y_filename + id + '.csv') 143 141 144 142 145 #-------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.