Changeset 5035


Ignore:
Timestamp:
Feb 14, 2008, 2:17:55 PM (16 years ago)
Author:
nick
Message:

update to run_circular

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/circular_island/run_circular.py

    r5031 r5035  
    104104r2 = 2.5
    105105d=0.75
    106 res=.05
    107 #res=1.
     106#res=.05
     107res=1.
    108108L=4.11
    109109
     
    178178attribute_dic, title_index_dic = csv2dict('sqrt_table.csv')
    179179#print attribute_dic
    180 print attribute_dic['number'][10]
     180#print attribute_dic['number'][10]
    181181
    182182def test(x,y):
     
    185185#    center_y = 6.80
    186186
    187     du = (center_x-x)**2+(center_y-L-y)**2
    188 #    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)
    190190    z1=[]
    191     for d in du:
     191    for xy in list_xy:
    192192        # this finds the int that relates to the correct sqrt in the table
    193193        nn=int(round(d,2)*100)
     
    293293    #print i,poly[i]
    294294   
    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.