Changeset 5147 for anuga_validation/circular_island/run_circular.py
 Timestamp:
 Mar 10, 2008, 3:43:22 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_validation/circular_island/run_circular.py
r5131 r5147 30 30 time = strftime('%Y%m%d_%H%M%S',localtime()) 31 31 32 res=0.0 02532 res=0.05 33 33 34 34 dir_comment = time+'_res_'+str(res)+'_'+str(user) … … 113 113 center_y = 13.80 114 114 115 #d=0.75116 115 #L=3.6 and therefore the gauges 1,2,3 and 4 are located 1.8 from the toe of the island 117 116 # 13.8  3.6(radii of island)  1.8 = 8.4 … … 126 125 angle = ((angle180)/180.0)*pi 127 126 128 p = get_xy(center_x,center_y,4. 0,angle)127 p = get_xy(center_x,center_y,4.5,angle) 129 128 poly1.append([p[0],p[1]]) 130 129 … … 136 135 #convert Degs to Rads 137 136 angle = ((angle180)/180.0)*pi 138 p = get_xy(center_x,center_y, 2.6,angle)137 p = get_xy(center_x,center_y,3.2,angle) 139 138 poly2.append([p[0],p[1]]) 140 139 … … 179 178 180 179 #attribute_dic, title_index_dic = csv2dict('sqrt_table.csv') 181 #print attribute_dic182 #print attribute_dic['number'][10]183 180 184 181 def circular_island_elevation(x,y): 185 182 water_depth = 0.32 186 # list_xy_square = (center_xx)**2+(center_yLy)**2187 183 list_xy = sqrt((center_xx)**2+(center_yLy)**2) 188 184 print 'x',min(x),max(x) … … 192 188 for distance in list_xy: 193 189 194 # z = distance*.25+0.9195 # z = distance*.25+0.578190 # The cone has a 1/4 slope and a 7.2 diameter this make the height 0.9m 191 # to get the stage at 0 elevation + cone height and  water depth 196 192 z = distance*.25+0.9water_depth 197 # z = zz*.25+0.8975198 193 199 194 if z <0water_depth: … … 201 196 if z > .625water_depth: 202 197 z=0.625water_depth 203 # print 'hello',z204 198 list_z.append(z) 205 199 return list_z … … 286 280 # print p1,p2,p3,p4,angle 287 281 poly_segment =[[p1[0],p1[1]],[p2[0],p2[1]],[p3[0],p3[1]],[p4[0],p4[1]]] 288 # print i, poly_segment, angle289 #print i,poly[i]290 282 291 283 run_up, location = get_maximum_inundation_data(filename=sww_file,polygon=poly_segment, verbose=False) … … 297 289 sww2csv_gauges(sww_file=sww_file, 298 290 gauge_file=anuga_dir+'gauges'+sep+'gauges.csv', 299 gauge_file='gauges.csv',291 # gauge_file='gauges.csv', 300 292 quantities = ['stage','depth', 'elevation', 'xmomentum', 'ymomentum'], 301 293 verbose=True,
Note: See TracChangeset
for help on using the changeset viewer.