Changeset 5035
- Timestamp:
- Feb 14, 2008, 2:17:55 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_validation/circular_island/run_circular.py
r5031 r5035 104 104 r2 = 2.5 105 105 d=0.75 106 res=.05107 #res=1.106 #res=.05 107 res=1. 108 108 L=4.11 109 109 … … 178 178 attribute_dic, title_index_dic = csv2dict('sqrt_table.csv') 179 179 #print attribute_dic 180 print attribute_dic['number'][10]180 #print attribute_dic['number'][10] 181 181 182 182 def test(x,y): … … 185 185 # center_y = 6.80 186 186 187 du= (center_x-x)**2+(center_y-L-y)**2188 #print 'x',min(x),max(x)189 #print 'y',min(y),max(y)187 list_xy = (center_x-x)**2+(center_y-L-y)**2 188 print 'x',min(x),max(x) 189 print 'y',min(y),max(y) 190 190 z1=[] 191 for d in du:191 for xy in list_xy: 192 192 # this finds the int that relates to the correct sqrt in the table 193 193 nn=int(round(d,2)*100) … … 293 293 #print i,poly[i] 294 294 295 run_up, x_y= get_maximum_inundation_data(filename=model_name+'.sww',polygon=poly_segment, verbose=False)296 297 print 'maxd %.3f, %.6f, %.6f, %.5s, %.5s '%(((angle/pi)*180)+180, run_up*100, run_up-water_height*100, x_y[0],x_y[1])298 299 300 301 302 303 304 295 run_up, location = get_maximum_inundation_data(filename=model_name+'.sww',polygon=poly_segment, verbose=False) 296 297 print 'maxd %.3f, %.6f, %.6f, %.5s, %.5s '%(((angle/pi)*180)+180, run_up*100, run_up-water_height*100, location[0],location[1]) 298 299 300 301 302 303 304
Note: See TracChangeset
for help on using the changeset viewer.