Changeset 5590
- Timestamp:
- Aug 1, 2008, 10:05:08 AM (15 years ago)
- Location:
- anuga_work/development/Hinwood_2008
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/Hinwood_2008/calc_norm.py
r5577 r5590 14 14 from os import sep, path 15 15 from csv import writer 16 from time import localtime, strftime 17 18 # Related major packages 16 19 from Numeric import arange, array, zeros, Float, where, greater, less, \ 17 20 compress, argmin, choose, searchsorted, sqrt, sum 18 21 19 # Related major packages20 22 import project 21 23 from os import sep … … 142 144 if is_interactive: 143 145 ion() 144 146 time_date = strftime('plot date: %d/%m/%Y Time: %H:%M:%S', 147 localtime()) 145 148 subplot(212) 146 plot_title = id + " Root Mean Square Deviation comparison" 149 plot_title = id + " Root Mean Square Deviation comparison" + '\n' \ 150 + time_date 151 147 152 title(plot_title) 148 153 y_label = "RMSD" 149 154 ylabel(y_label) 155 xlabel("x location, m") 150 156 grid(True) 151 157 … … 164 170 rmsd_list = [float(x) for x in simulation['rmsd']] 165 171 lines.append(plot(locations, rmsd_list)) 172 173 174 for break_x in run_data['break_xs']: 175 axvspan(break_x-0.001,break_x+0.001, facecolor='g') 166 176 167 177 figlegend(lines, outputdir_tags,'upper left') … … 273 283 #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!! 274 284 #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!! 275 for outputdir_tag in outputdir_tags:276 auto_rrms(outputdir_tag, scenarios, "stage")285 #for outputdir_tag in outputdir_tags: 286 # auto_rrms(outputdir_tag, scenarios, "stage") 277 287 278 288 #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!! 279 #auto_plot_rrms(outputdir_tags, scenarios, "stage")280 plot_settings(outputdir_tags, scenarios, "stage")289 auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage") 290 #plot_settings(outputdir_tags, scenarios, "stage") -
anuga_work/development/Hinwood_2008/run_dam.py
r5577 r5590 112 112 finaltime = None 113 113 maximum_triangle_area=0.0001 114 #outputdir_name += '_good' 115 116 elif run_type == 7: 117 # this is not a test 118 # Output will go to a file 119 # The sww file will be interpolated 120 yieldstep = 0.01 121 finaltime = None 122 maximum_triangle_area=0.00001 123 #outputdir_name += '_good' 124 125 elif run_type == 8: 126 # this is not a test 127 # Output will go to a file 128 # The sww file will be interpolated 129 yieldstep = 0.05 130 finaltime = None 131 maximum_triangle_area=0.00001 114 132 #outputdir_name += '_good' 115 133 … … 232 250 Bd = Dirichlet_boundary([0.3,0,0]) 233 251 Bts = Time_boundary(domain, function) 252 #Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function) 234 253 domain.set_boundary( {'wall': Br, 'wave': Bts} ) 235 254 #domain.set_boundary( {'wall': Br, 'wave': Bd} ) … … 287 306 # 5 is 0.001 288 307 # 6 is 0.0001 308 # 7 is 0.00001 289 309 290 310 #run_type = 1 291 run_type = 5311 run_type = 8 292 312 #for run_data in [scenarios[5]]: 293 313 #scenarios = scenarios[2:] … … 295 315 width = 1.0 296 316 width = 0.1 317 #width = 0.01 297 318 for run_data in scenarios: 298 319 pro_instance = main( run_data['scenario_id'] + '_boundary.tsm' , … … 302 323 outputdir_name=run_data['scenario_id'], 303 324 use_limits=False, 304 friction=0.012, 325 friction=0.012, 305 326 end_tag='_G') 306 327 #gauges_for_slope(pro_instance.outputdir,[run_data]) -
anuga_work/development/Hinwood_2008/slope.py
r5577 r5590 674 674 #scenarios = [scenarios[0]] 675 675 outputdir_tag = "_good_tri_area_0.01_limiterD" 676 outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.0 1_F"676 outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.0001_F" 677 677 slope_tag = "" 678 678 #outputdir_tag = "_test_limiterC" 679 #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!679 scenarios = [scenarios[4]] # !!!!!!!!!!!!!!!!!!!!!! 680 680 #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!! 681 681 682 gauges_for_slope(slope_tag, outputdir_tag, scenarios)682 #gauges_for_slope(slope_tag, outputdir_tag, scenarios) 683 683 #auto_graph_slopes(outputdir_tag, scenarios) #, is_interactive=True) 684 auto_find_min_slopes(slope_tag, outputdir_tag, scenarios)684 #auto_find_min_slopes(slope_tag, outputdir_tag, scenarios) 685 685 #auto_graph_froudes(outputdir_tag, scenarios) 686 #auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios)686 auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios) 687 687 #g = Get_file_name(scenarios[0], outputdir_tag, slope_tag) 688 688 #for wave_file, save_as, wave_number in Get_file_name(
Note: See TracChangeset
for help on using the changeset viewer.