Changeset 5592


Ignore:
Timestamp:
Aug 1, 2008, 1:33:26 PM (16 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/perth
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/export_results.py

    r5569 r5592  
    1010#time_dir = '20080714_095738_run_final_0.6_exmouth_kvanputt'
    1111#time_dir = '20080714_134935_run_final_0_exmouth_kvanputt'
    12 time_dir = '20080718_112228_run_final_0_exmouth_kvanputt'
     12#time_dir = '20080718_112228_run_final_0_exmouth_kvanputt'
     13#time_dir = '20080725_173911_run_final_0.6_polyline_alpha0.1_kvanputt'
     14#time_dir = '20080728_112519_run_final_0.6_polyline_alpha0.2_kvanputt'
     15time_dir = '20080728_154123_run_final_0.6_polyline_alpha0.15_kvanputt'
    1316
    1417cellsize = 25
  • anuga_work/production/perth/get_gauges.py

    r5581 r5592  
    3636    for i, name in enumerate(quantity_names):
    3737        quantities[name] = fid.variables[name][:]
    38         if inside_points is not None:
    39             quantities[name] = take(quantities[name],indices,1)
    40         else:
    41             msg = 'No gauges found inside polygon'
    42             raise msg
    43            
     38##        if inside_points is not None:
     39##            quantities[name] = take(quantities[name],indices,1)
     40##        else:
     41##            msg = 'No gauges found inside polygon'
     42##            raise msg
     43##           
    4444    for j in range(len(x)-1):
    4545        locx=int(x[j])
     
    6363polygon=project.poly_all
    6464points,quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir,project.scenario_name),polygon,verbose=False)
    65 
     65print quantities['stage'][0,:]
     66print len(points), len(elevation), len(quantities['stage'][0,:])
    6667assert len(points)==len(elevation)==len(quantities['stage'][0,:])
  • anuga_work/production/perth/get_timeseries.py

    r5576 r5592  
    1717
    1818#timestamp='20080724_121200_run_trial_0.6_polyline_alpha0.1_kvanputt'
    19 timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
     19#timestamp='20080724_161830_run_final_0.6_polyline_alpha0.1_kvanputt'
     20timestamp='20080725_173911_run_final_0.6_polyline_alpha0.1_kvanputt'
    2021
    2122filename=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
     
    2324
    2425sww2csv_gauges(filename,
    25                 project.gauges_dir_name,
     26                #project.gauges_dir_name,
     27                project.gauges_dir_name2,
    2628                quantities = ['stage','speed','depth','elevation'],
    2729                verbose=True)
  • anuga_work/production/perth/project.py

    r5581 r5592  
    3333#Making assumptions about the location of scenario data
    3434state = 'western_australia'
     35#scenario_name = 'perth_44unitsources'
    3536scenario_name = 'perth'
    3637scenario = 'perth_tsunami_scenario'
     
    3839tide = 0.0 #0.6
    3940
    40 alpha = 0.1
     41alpha = 0.15
    4142friction=0.01
    4243starttime=0
     
    152153#gauges
    153154gauge_name = 'perth.csv'
     155gauge_name2 = 'thinned_MGA50+50.csv'
    154156
    155157gauges_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'gauges'+sep
    156158beach_gauges = gauges_dir + 'beach_gauges.csv'
    157159gauges_dir_name = gauges_dir + gauge_name
     160gauges_dir_name2 = gauges_dir + gauge_name2
    158161
    159162buildings_filename = gauges_dir + 'Perth_resA.csv'
Note: See TracChangeset for help on using the changeset viewer.